A broad deglacial δ13C minimum event in planktonic foraminiferal records in the Okinawa Trough

The equatorial Pacific upwelling zone has been suspected of playing an important role in the global atmospheric CO2 changes associated with glacial-interglacial cycles.In order to assess the influencing scope of the surface water deglacial δ13C minimum in the tropical Iow-latitude Pacific,the core DGKS9603, collected from the middle Okinawa Trough, was examined for δ13C records of planktonic foraminifera N. dutertrei and G. ruber. The planktonic foraminiferal δ13C records show a clear decreasing event from 20 to 6 cal. kaBP., which is characterized by long duration of about 14 ka and amplitude shift of 0.4 × 10-3. Its minimum value occurred at 15.7 cai kaBP. The event shows fairly synchrony with the surface water deglacial δ13C minimum identiffed in the tropical Pacific and its marginal seas. Because there is no evidence in planktonic foraminiferal fauna and δ18O records for upwelling and river runoff enhancement,the broad deglacial δ13C minimum event in planktonic foraminiferal records revealed in core DGKS9603 might have been the direct influencing result of the deglacial surface water of the tropical Pacific. The identification for the event in the Okinawa Trough provides new evidence that the water evolution in the tropical low-latitude Pacific plays a key role in large regional, even global carbon cycle.     

Early Pleistocene formation of the asymmetric east-west pattern of upper water structure in the equatorial Pacific Ocean

Surface-and subsurface-dwelling planktonic foraminifera from the upper 43 m of Hole A at the Ocean Drilling Program (ODP) Site 807,which was recovered from the western Pacific warm pool during ODP Leg 130,were analyzed for stable oxygen and carbon isotopes.By comparing these results with data from ODP Site 851 in the eastern equatorial Pacific,this study has reconstructed the paleoceanographic changes in upper ocean waters in the equatorial Pacific since 2.5 Ma.During the period from 1.6-1.4 Ma,the oxygen isotopes of surface and subsurface waters were found to markedly change in the western and eastern equatorial Pacific,further confirming the final formation of the well-defined asymmetric east-west (E-W) pattern at that time.This feature was similar to the zonal temperature gradient (sea surface temperature is higher in the west and lower in the east) and the asymmetric upper water structure (thermocline depth is deeper in the west and shallower in the east) in the modern equatorial Pacific.The zonal gradient change of subsurface water δ18O was greater than that of surface water δ18O,indicating that the formation of the asymmetric E-W pattern in the equatorial Pacific should be much more related to the shoaled thermocline and markedly decreased subsurface water temperature in the eastern equatorial Pacific.Moreover,since ～1.6 Ma,the carbon isotopic differences between surface and subsurface waters clearly decreased in the equatorial Pacific,and their long-term eccentricity periods changed from 400 ka to ～500 ka,reflecting the reorganization of the ocean carbon reservoir.This probably resulted from the deep water reorganization in the Southern Ocean at that time and its enhanced influence on the tropical Pacific (especially subsurface water).Our study demonstrates that the tropical ocean plays an important role in global climate change.     

C－isotope composition and significance of the Sinian on the Tarim plate

102 samples of limestones,mudrocks shales and tillites from the Sinian in Kuluketage have been analyzed ,Four formations in the Sinian belong to glacial deposits in which the differences in δ^13C are distinct (-14.7‰- 4.2‰),Distinctly negative δ^13C excursions for three times support the viwe of three glacial deposits in this area suggested by previous authors,except for the durations of them ,δ^13C values are generally positive in the Beiyixi Formation of lower part of the Lower Sinian ,representing probably the rock records of a continental rift related to break-up of the Rodinian supercontinent,There are three times of distinct excursions of δ^13C curve in all strata,The first one is from positive values at the bottom to about -5‰ at the top of the Beiyixi Formation ,The second one,from 0‰-3‰ in the Arletonggou Formation of upper part of the Lower Sinian to -3.4‰---14.4‰ in the Teruiaiken Formation.The third one,from positive values in the Zhamoketi Formation of lower part of the Upper Sinian to negative values in glacial varves at the top of the Hankeerqiaoke Formation.     

Different climatic controls of soil δ~(13)C_(org) in three mid-latitude regions of the Northern Hemisphere since the Last Glacial period

Paleoecological records of soil δ13Corg from three regions in the middle latitudes of the Northern Hemisphere, including the Chinese Loess Plateau (CLP), the Great Plains and adjacent areas of North America and northwestern Europe, showed different variations since the Last Glacial period. An at-tempt was made to evaluate the causes for the difference in δ13Corg on the basis of the modern climatic data collected in these regions and of the modern C3 and C4 plant distributions. The analysis indicates that temperature, especially the growing season temperature, has a dominant control on the growth of C4 plants. When the mean annual or growing season temperatures are below the "threshold value" , the growth of C4 plants is limited. When the temperature is above the "threshold value" , C4 plants can grow under a wide range of precipitation. However, when the precipitation is high enough to favor the growth of trees, the proportions of C4 plants in local biomass will decline. The implicit control factor recovered by sedimentary records is consistent with the control factor on modern C3/C4 distribution. Pure C3 plants have been dominating the local biomass since the Last Glacial period in European loess region, mainly owing to the low local temperature. The increases in C4 plants from the late Pleistocene to the Holocene in the Chinese Loess Plateau, the Great Plains and adjacent areas, mainly reflect the influ-ence of increasing temperature.     

Different climatic controls of soil δ 13Corg in three mid-latitude regions of the Northern Hemisphere since the Last Glacial period

Paleoecological records of soil δ^13Corg from three regions in the middle latitudes of the Northern Hemisphere, including the Chinese Loess Plateau （CLP）, the Great Plains and adjacent areas of North America and northwestern Europe, showed different variations since the Last Glacial period. An attempt was made to evaluate the causes for the difference in δ^13Corg on the basis of the modern climatic data collected in these regions and of the modern C3 and C4 plant distributions. The analysis indicates that temperature, especially the growing season temperature, has a dominant control on the growth of C4 plants. When the mean annual or growing season temperatures are below the ＂threshold value＂, the growth of C4 plants is limited. When the temperature is above the ＂threshold value＂, C4 plants can grow under a wide range of precipitation. However, when the precipitation is high enough to favor the growth of trees, the proportions of C4 plants in local biomass will decline. The implicit control factor recovered by sedimentary records is consistent with the control factor on modern C3/C4 distribution. Pure C3 plants have been dominating the local biomass since the Last Glacial period in European loess region, mainly owing to the low local temperature. The increases in C4 plants from the late Pleistocene to the Holocene in the Chinese Loess Plateau, the Great Plains and adjacent areas, mainly reflect the influence of increasing temperature.     

Relationship between the tropical cyclone genesis over the Northwest Pacific and the sea surface temperature anomalies

Using the tropical cyclone (TC) data derived from the Joint Typhoon Warning Center (JTWC) and the sea surface temperature data derived from the Joint Environmental Data Analysis Center (JEDAC) at the Scripps Institute of Oceanography from January 1955 to December 2000, we analyzed the relationship between the TC genesis over the Northwest Pacific (NWP) and the sea surface temperature anomalies (SSTA) over the Pacific basin. A long-term trend indicated that the highest frequency of monthly TC genesis appeared earlier and the annual genesis sum increased gradually during the last half century with some oscillations. No significant synchronous correlation was found between the NWP TC events and the SSTA over the Pacific basin, while the annual sum of TC genesis was closely related with the SSTA averaged from the first three months (January, February and March) of the year in the equatorial western and eastern Pacific and over mid-high latitudes of the North Pacific. The results implied that there are an interannual El Nio SSTA mode in the equatorial western and eastern Pacific and an interdecadal SSTA mode in the northern Pacific, which affected the TC genesis. A regression analysis between the first three-month SSTA and the annual TC sum based on two time scales was conducted. The correlation coefficient between simulated and observed TC sums reached a high value of 0.77.     

Signal propagations and linkages of subsurface temperature anomalies in the tropical Pacific and Indian Ocean

The propagation characteristics of signals along different zonal-time profiles are analyzed using surface and subsurface temperature anomalies over the tropical Pacific and Indian oceans. Analyses show that there are intrinsic relationships between El Nio events in the eastern equatorial Pacific and dipole events in the equatorial Indian Ocean. In the region of tropical North Pacific between the equator and 16°N, there is a circle of propagation of subsurface temperature anomalies. El Nio events only happen when the warm subsurface signals reach the eastern equatorial Pacific. Dipole events are characterized when a warm subsurface signal travels along off-equatorial Indian Ocean to the western boundary. From these analyses, we believe that subsurface temperature anomalies can be considered to be the oceanographic early signal to forecast El Nio events in Pacific Ocean and dipole events in Indian Ocean, respectively.     

Pulleniatina Minimum Event during the last deglaciation in the southern South China Sea

The planktonic foraminiferal faunal census of core MD 05-2894 （7°2.25′N, 111°33.11′E, water depth 1982 m）, retrieved from the southern South China Sea （SCS） during the ＂Marco Polo＂ cruise in 2005, was performed to investigate the abundance changes of a subsurface dweller, Pulleniatina obliquiloculata. The results display that the abundance of P. obliquiloculata nearly declines to zero during 16.0--14.9 ka, corresponding to the Heinrich 1 （H1） cold interval. The unexpected decrease of P. obliquiloculata occurs in the adjacent cores, roughly between 17 and 14.8 ka based on the previous studies. Accordingly, the Pulleniatina Minimum Event in the last deglaciation can serve as a good stratigraphical indicator, at least in the southern SCS. To further explore the changes of sea surface temperature （SST） and subsurface seawater temperature （SSST）, we made parallel Mg/Ca measurements on surface dweller Globigerinoides tuber and subsurface dweller P. obliquiloculata tests. Since the last deglaciation, the SSTs show a continuous increasing trend towards the late Holocene, while the warming of the subsurface water is punctuated by a 2℃-cooling interval across the deglacial Pulleniatina Minimum Event. Both increased 5180 differences between G. ruber and P. obliquiloculata, and increased temperature differences between surface and subsurface water suggest a shoaling of the mixed layer during the deglacial Pulleniatina Minimum Event. Therefore, we consider that the significant changes in the upper ocean structure are responsible for the Pulleniatina Minimum Event during the last deglaciation in the southern SCS.     

Relationship between the tropical cyclone genesis over the Northwest Pacific and the sea surface temperature anomalies

Using the tropical cyclone (TC) data derived from the Joint Typhoon Warning Center (JTWC) and the sea surface temperature data derived from the Joint Environmental Data Analysis Center (JEDAC) at the Scripps Institute of Oceanography from January 1955 to December 2000, we analyzed the relationship between the TC genesis over the Northwest Pacific (NWP) and the sea surface temperature anomalies (SSTA) over the Pacific basin. A long-term trend indicated that the highest frequency of monthly TC genesis appeared earlier and the annual genesis sum increased gradually during the last half century with some oscillations. No significant synchronous correlation was found between the NWP TC events and the SSTA over the Pacific basin, while the annual sum of TC genesis was closely related with the SSTA averaged from the first three months (January, February and March) of the year in the equatorial western and eastern Pacific and over mid-high latitudes of the North Pacific. The results implied that there are an interannual El Nio SSTA mode in the equatorial western and eastern Pacific and an interdecadal SSTA mode in the northern Pacific, which affected the TC genesis. A regression analysis between the first three-month SSTA and the annual TC sum based on two time scales was conducted. The correlation coefficient between simulated and observed TC sums reached a high value of 0.77.     

Carbon isotopic record from Upper Devonian carbonates at Dongcun in Guilin, southern China, supporting the worldwide pattern of carbon isotope excursions during Frasnian-Famennian transition

Two positive δ^13C excursions are presented in records from the Frasnlan-Famenninn (F-F) marine carbonate sediments in Europe, America, Africa, and Australia,having been considered as a worldwide pattern, and attributed to enhanced organic carbon burial during the F-F biological mass extinction, However, this worldwide pattern has not been revealed from the well-deposited Late Devonian sequences in southern China. In this paper, a detailed investigation has been made on the Late Devonian section at Dongcun, Guilin, southern China to constrain perturbationsin δ^13C of carbonates in the F-F deposited sequence. The result from this section also indicates two positive δ^13C excursions during the F-F transition. The first excursion with an amplitude of 1.5‰ occurred at the bottom of linguiformis Zone, later than the early excursion existing in the Late rhenana Zone of the Late Devonian profiles in other continents,especially, in central Europe. This difference has been expected to be a result as conodont Palmatolepis linguiformis occurred earlier in southern China than other sites. The second excursion with an amplitude of 2.1‰ is located at the F-F boundary, same as the records from other continents.This result strongly supports the view that two carbon isotope positive excursions during the F-F transition are common in carbonate sediments, resulting from worldwide increases of organic carbon burial intensity.     

Sea surface temperature east of Australia： A predictor of tropical cyclone frequency over the western North Pacific？

The relationship between sea surface temperature (SST) east of Australia and tropical cyclone frequency over the western North Pacific (WNPTCF) is analyzed by use of observation data.The WNPTCF from June to October is correlated negatively to spring SST east of Australia.When the spring SST is in the positive phase,a cyclonic circulation anomaly in the upper troposphere and an anticyclonic circulation anomaly in the lower troposphere prevail over the western North Pacific from June to October,concurrent with an anomalous atmospheric subsidence and an enlarged vertical zonal wind shear.These conditions are unfavorable for tropical cyclone genesis,and thus WNPTCF decreases.The negative phase of the spring SST east of Australia leads to more tropical cyclones over the western North Pacific.The spring SST east of Australia may give rise to simultaneous change in tropical atmospheric circulation via the teleconnection wave train,and then subsequently affect atmospheric circulation variation over the western North Pacific.     

Reconstruction of surface ocean water <Emphasis Type="Italic">p</Emphasis>CO<Subscript>2</Subscript>(aq) in Nansha area,the South China Sea during the last 30 ka

The reconstruction of pCO2 in the tropic ocean is one of the most important issues to understand global climatic changes.In this study,the high-resolution stratigraphic analysis of core 17962 was conducted,which is Iocated in the southern South China Sea(SCS).The contents of sedimentary organic matter,the stable carbon isotopic composition of sedimentary organic matter,and the δ^13C values of black carbon and terrigenons n-alkanes were determined.And the δ^13Cwc value of carbon derived from aquatic was calculated.On the basis of δ^13Corg-pCO2equation proposed by Popp et al.(1989),we estimated the pCO2 in the Nansha area,SCS,since the last glaciation using δ^13Cwc instead of δ^13Corg.The results show that the average pCO2 was estimated at 240 ppm V during the last glaciation,and at 320ppm V in the Holocene.A comparison of surface sea pCO2 with the atmosphere CO2 recorded in the Vostok ice core,indicates that CO2 in surface water of the southern SCS could emit into atmosphere during the last 30ka.     

Sources and transportation of suspended matter and sediment in the southern Yellow Sea： Evidence from stable carbon isotopes

The concentrations of total suspended matter (TSM) and the compositions of organic stable carbon isotopes of TSM and bottom sediments were analyzed to study the sources of TSM and sediments and the transportation processes. For this study, 284 TSM samples and 64 sediment ones taken from 67 stations along 7 transects and in 5 layers were collected in the southern Yellow Sea on the cruise in May, 1998. The main sediment transportation pattern in the southern Yellow Sea was obtained by analyzing the distribution characteristics of TSM concentration and particulate organic carbon δ^13C values. It was confirmed from the pattern that the bottom layer plays a more important role than the surface one in the transportation processes of terrigenous material to the central deep-water area of the southern Yellow Sea. The Yellow Sea circulation is an important control factor in determining the sediment transportation pattern in the southern Yellow Sea. The carbon isotope signals of sedimentary organic matter confirmed that the main material in sediments with high sedimentation rate in the Shandong subaqueous delta originated from the modern Yellow River.The terrigenous sediments in the deep-water area of the southern Yellow Sea are mainly from the abandoned Yellow River and the modern Yellow River, and a small portion of them are from the modern Yangtze material. The amount of terrigenous material from Korean Peninsula and its influence range are relatively limited. The conclusions derived from TSM and stable carbon isotopes were further confirmed by another independent material source tracer--PAHs.     

An improvement of the too cold tongue in the tropical Pacific with the development of an ocean-wave-atmosphere coupled numerical model

A common problem in the application of the coupled ocean-atmosphere general circulation models (CGCMs) without flux correction is that the simulated equatorial cold tongue in general tends to be too strong, narrow, and extending too far west. The causes are not well understood yet. One possible reason may be the simulated mixed layer depth (MLD) is too shallow in the tropical Pacific due to insufficient vertical mixing in the OGCM. It is believed that the wave-induced vertical mixing can greatly improve the simulation of the MLD and thermocline structure. In this study, the coupled ocean-atmosphere general circulation model (FGCM-0) incorporated with wave-induced mixing has been employed to simulate the tropical Pacific sea surface temperature (SST). Generally, the wave-induced mixing lowers the SST in the OGCM because the strengthened vertical mixing can bring more cold water upward. However, in the coupled model, the non-uniformity of the space distribution in SST drop generates a horizontal gradient of the sea surface air pressure, and thus yields surface wind field anomaly. The wind anomaly leads to both ocean surface circulation anomaly and downwelling anomaly, which can restrain the overly-westward extension of the cold tongue in the tropical Pacific. Compared with the model results from the original model (FGCM-0), the modeled SST is higher by more than 0.8℃, with a maximum of 1.2℃in the western Pacific (160—180°E, 0—3°N) . The eastern boundary of the isotherm of 26. 0℃also moves from 165°E to 180°E. The overly westward extension of the simulated equatorial cold tongue is suppressed with the incorporation of the wave-induced mixing in the coupled model. The simulated SST shows in general improved results with a maximum improvement of more than 1.0℃. The simulated SST improvement in the north tropical Pacific is much better than that of the south tropical Pacific.     

An improvement of the too cold tongue in the tropical Pacific with the development of an ocean-wave-atmosphere coupled numerical model

A common problem in the application of the coupled ocean-atmosphere general circulation models (CGCMs) without flux correction is that the simulated equatorial cold tongue in general tends to be too strong, narrow, and extending too far west. The causes are not well understood yet. One possible reason may be the simulated mixed layer depth (MLD) is too shallow in the tropical Pacific due to insufficient vertical mixing in the OGCM. It is believed that the wave-induced vertical mixing can greatly improve the simulation of the MLD and thermocline structure. In this study, the coupled ocean-atmosphere general circulation model (FGCM-0) incorporated with wave-induced mixing has been employed to simulate the tropical Pacific sea surface temperature (SST). Generally, the wave-induced mixing lowers the SST in the OGCM because the strengthened vertical mixing can bring more cold water upward. However, in the coupled model, the non-uniformity of the space distribution in SST drop generates a horizontal gradient of the sea surface air pressure, and thus yields surface wind field anomaly. The wind anomaly leads to both ocean surface circulation anomaly and downwelling anomaly, which can restrain the overly-westward extension of the cold tongue in the tropical Pacific. Compared with the model results from the original model (FGCM-0), the modeled SST is higher by more than 0.8℃, with a maximum of 1.2℃ in the western Pacific (160—180°E, 0—3°N). The eastern boundary of the isotherm of 26.0℃ also moves from 165°E to 180°E. The overly westward extension of the simulated equatorial cold tongue is suppressed with the incorporation of the wave-induced mixing in the coupled model. The simulated SST shows in general improved results with a maximum improvement of more than 1.0℃. The simulated SST improvement in the north tropical Pacific is much better than that of the south tropical Pacific.     

Exploring cyclic changes of the ocean carbon reservoir

A 5-Ma record from ODP Site 1143 has re-vealed the long-term cycles of 400—500 ka in the carbon isotope variations. The periodicity is correlatable all over the global ocean and hence indicative of low-frequency changes in the ocean carbon reservoir. As the same periodicity is also found in carbonate and eolian dust records in the tropical ocean, it may have been caused by such low-latitude proc-esses like monsoon. According to the Quaternary records from Site 1143 and elsewhere, major ice-sheet expansion and major transition in glacial cyclicity (such as the Mid-Brunhes Event and the Mid-Pleistocene Revolution ) were all pre-ceded by reorganization in the ocean carbon reservoir ex-pressed as an episode of carbon isotope maximum (d 13Cmax), implying the role of carbon cycling in modulating the glacial periodicity. The Quaternary glacial cycles, therefore, should no more be ascribed to the physical response to insolation changes at the Northern Hemisphere high latitudes alone; rather, they have been driven by the 揹ouble forcing? a combination of processes at both high and low latitudes, and of processes in both physical (ice-sheet) and biogeochemical (carbon cycling) realms. As the Earth is now passing through a new carbon isotope maximum, it is of vital impor-tance to understand the cyclic variations in the ocean carbon reservoir and its climate impact. The Pre-Quaternary varia-tions in carbon and oxygen isotopes are characterized by their co-variations at the 400-ka eccentricity band, but the response of d 13C and d 18O to orbital forcing in the Quater-nary became diverged with the growth of the Arctic ice-sheet. The present paper is the second summary report of ODP Leg 184 to the South China Sea.     

^226Ra evidence for the ecosystem shift over the past 40 years in the North Pacific Subtropical Gyre

Surface seawater was collected for ~(226)Ra measurement in the North Pacific Subtropical Gyre from July to October, 1999 and October to December, 2003. Combined with the historical data reported for this sea area, a declined trend of surface ~(226)Ra concentrations was observed since 1960s, indicating the ecosystem shift in response to global warming. On one side, the enhanced stratification of the upper water column resulting from global warming reduced the ~(226)Ra input from the depth, on the other, the temporal increase of biological production resulting from the climate-related ecosystem structure change strengthened the ~(226)Ra removal from the surface ocean. Both the physical and biological processes resulted in the decrease of surface ~(226)Ra concentrations in the North Pacific Subtropical Gyre. The temporal trend of surface ~(226)Ra concentrations was consistent with the trends of chlorophyll a, silicate, phosphate and primary production previously reported. This study provided ~(226)Ra evidence for the ecosystem shift under global change.     

Analysis of carbon isotopes in airborne carbonate and implications for aeolian sources

Methods were developed to determine the mass ratios of carbon isotopes in trace amounts of aerosol carbonate. A Finnigan MAT 252 mass spectrometer fitted with an on-line Kiel device was to determine the ^13C/^12C ratio in CO2 produced from the carbonate. A study using these methods was conducted to characterize the carbonate carbon isotopes in aerosol samples collected in Xi‘an on dusty and normal days during March and April 2002. Results of the study demonstrate that insights into the origin of the dust can be deduced from its isotopic composition. That is, the δ^13C ofcarbonate for dust storm samples ranged from -1.4‰ to -4.2‰, and this is consistent with sandy materials in dust source regions upwind. In contrast, for non-dusty days δ^13C ranged from -7.5% to -9.3‰, which is more similar to fine particles emitted from local surface soils. Comparisons of dust storm aerosols with surface soils from source regions and with aerosol samples collected downwind indicate that the δ^13C values did not change appreciably during longrange transport. Therefore, carbon isotopes have the potential for distinguishing among source materials, and this approach provides a powerful new tool for identifying dust provenance.     

Solar forced transient evolution of Pacific upper water thermal structure during the Holocene in an earth system model of intermediate complexity

Forced by transient solar activities since 7 ka,the thermal structures of the Pacific upper water at boreal winter are featured by an enhanced response of 3-dimensional Western Pacific Warm Pool(WPWP) in an Earth system model of Intermediate Complexity at centennial scales.During solar maximum periods,the magnitude of surface ocean temperature variations is 30% larger in the western tropical Pacific than in the Ni o3 region,while at subsurface,it is 40% larger in the subtropical North Pacific than in the western Equatorial Pacific.They compromise stronger zonal and meridional thermal gradients in surface and subsurface Pacific respectively which are both linearly responded to solar forcing at centennial periods.The surface gradient is most sensitive at 208-year period while the subsurface gradient shows more significance at periods longer than 208-year.Also noteworthy are two differences:(1) the phase lags at these periods of surface gradient are slightly smaller than that of subsurface;(2) the 148-year and 102-year periods in surface gradient are lost in subsurface gradient.These modeled features preliminary confirm the centennial fluctuations of WPWP in paleo-proxies and a potential solar forcing during the Holocene.     

Relative humidity reconstruction for northwestern China's Altay Mountains using tree-ring δ~(18)O

Relative humidity is an important factor in water and water vapor feedback cycles.In this study,we established a 222-year annual tree-ring δ~(18)O chronology for Siberian larch(Larix sibirica Ldb.)from the Altay Mountains in northwestern China.Climate response analyses revealed that the relative humidity was the primary factor limiting tree-ring δ~(18)O fractionation.Based on our analysis,tree-ring δ~(18)O can be used to reconstruct the July–August relative humidity based on both a reasonable mechanism of tree-ring δ~(18)O fractionation and a statistically significant regression model.We used this model to reconstruct variations in the July–August relative humidity,and the model explained 47.4% of the total variation in the measured relative humidity data from 1961 to 2011.The relative humidity in the study area increased from 1900 to the 1990s and decreased thereafter.Two regime-shift dry periods were detected during the study period(one from 1817 to 1830 and the other from 2004 to 2011).