Mean δ18O Research Articles

Challenging the Atacama desert: Agronomic and water conditions for pre-Hispanic maize agriculture in hyper arid environments inferred by δ18O isotopes

During the Formative Period (ca.2400–950 years BP), pre-Hispanic farmers in Pampa del Tamarugal developed a complex agricultural system in the hyper-arid Atacama Desert in Tarapacá, northern Chile. This system involved numerous agricultural fields congregated near the Tarapacá Valley’s perennial stream and the Guatacondo ravine’s ephemeral stream. Well-established villages such as Caserones, Pircas, Ramaditas, and Guatacondo accompanied these developments. However, the importance of understanding the water sources has been neglected despite the efforts to understand the relationships between pre-Hispanic agriculture and their living environment. Thus, we presumed the use of local water sources without a clear understanding of their exploitation methods or the associated technological implications. Furthermore, there is limited research on groundwater use in Tarapacá. This research aims to understand water sources used in the cultivation of maize (Zea mays) using δ18O isotope values obtained from pre-Hispanic maize kernels in Tarapacá. We compared these values with published δ18O values of water sources and applied a generalized linear model (GLM) with a Gaussian distribution, performing a Tukey’s post hoc test for multiple comparisons of means with heteroscedasticity-consistent covariance estimation. The best-fit model was identified using a stepwise model selection procedure based on the Bayesian information criterion (BIC). Our results indicate that mean δ18O values of organic matter from maize kernels range from 24.73 ‰ to 31.65 ‰. The best performing model on δ18O only included Group (BIC=298.2) as the explanatory variable as Period, Site, and Weight had no effects. These findings point towards a significant statistical relationship between the δ18O values of organic matter derived from maize kernels and the specific geographic regions they originate from. These values also show an enrichment of δ18O isotope in Tarapacá samples, except for Pica 8. These results indicate diverse agricultural strategies that utilized different water sources including perennial flow in the Tarapacá River, ephemeral runoff at Guatacondo, and groundwater in the Pampa del Tamarugal. Also, due to the enrichment of δ18O, we suggest that the circular structures found among the fields were used as water reservoirs leading to the observed enrichment of δ18O isotope values. We discuss the role of ancient agriculture technology in water management, the role of maize, cultural strategies, and the relationships with their environment. In conclusion, pre-Hispanic farmers managed limited water sources successfully despite intermittent drought for hundreds of years until the 1970 s, when agriculture was abandoned due to multiple factors.

Increasing effect of biocrusts on evaporation is evidenced by simulating evaporation and diffusion experiments and water stable isotope analysis

Evaporation, one of the main components of the water and energy balance between soil and atmosphere, plays a key role in hydrological processes particularly in semi-arid and arid climate regions. As critical living organism communities inhabiting the soil-atmosphere boundary, biocrusts are capable of modifying near-surface soil properties and structure in drylands. However, the roles of biocrusts in mediating soil evaporation still remain greatly controversial, and the literature yet lacks in a systematic assessment. In this study, evaporation experiments were conducted under simulated isothermal and non-isothermal conditions for bare soil and three types of biocrusts (cyanobacterial, moss, and their mixture) sampled from a semi-arid climate region of the Chinese Loess Plateau. During the evaporation, soil temperature and moisture were measured at depths 2, 5, 10, and 15 cm for those biocrusts and bare soil. Their water vapor diffusion process and stable isotopic composition were also analyzed in the laboratory. We found that the measured evaporation yielded 4.9%, 17.5%, and 31.6% higher average evaporation rate for cyanobacterial, cyano-moss, and moss crusts, respectively, when compared to the bare soil. In comparison to bare soil, all types of biocrusts markedly reduced the rate of temperature rise, prolonged the duration of wetting time, and delayed the formation of dry surface layer. The largest regulating effects on soil temperature and moisture were found in the moss crusts instead of cyanobacterial or cyano-moss crusts. The increasing evaporation in biocrusts was attributed to their higher fine particle content, lower bulk density, and greater water holding capacity. Moreover, the water vapor diffusion rate and diffusivity of biocrusts were 18.2% and 16.4% higher than those of bare soil, respectively, which is attributable to the increased total porosity and better soil structure. Additionally, the isotopic composition of soil water (0–50 cm) in bare soil and moss crusts was distributed below the local meteoric water line. The mean δ18O and δ2H in soil water gradually decreased with increasing soil depth, and the moss crusts were more isotopically enriched than bare soil. Based on the enhanced soil evaporation rate, improved vapor diffusivity, and enriched isotopic composition of biocrusts in contrast to bare soil, we conclude that biocrusts exacerbate vapor transport and evaporative loss, and consequently they are crucial surface cover in rebalancing surface ecohydrological and biochemical processes in global drylands. These findings improve our understanding of the complexity of dryland hydrology, and a full consideration should be given in future biocrust studies.

Altitude effect of precipitation isotopes in an arid mountain-basin system: Observation and modelling around the world’s second-largest shifting desert

The stable water isotopes of precipitation provide important information about the hydrological circulation. In the arid mountain-basin system in central Asia, the altitude effect of precipitation isotopes has been a controversial topic in recent years, but the sample availability in extreme environments constrains the accurate understanding of the relationship between altitude and stable isotopes in precipitation. Based on the observation of precipitation isotopes around the Tarim Basin covered by the world’s second-largest shifting desert, we examined the relationship between altitude and isotope composition. There is an altitude effect of precipitation isotopes between the basin and the surrounding mountains, with the modelled gradient for annual mean δ18O being approximately 1.96 ‰ per 1000 m, which is weaker than the observed gradient focusing on the oases (2.75 ‰ per 1000 m). The largest modelled difference in δ18O between 1000–2000 m and 2000–3000 m above sea level occurs in August and September. The periods with a larger gradient of altitude effect usually have higher temperature and more precipitation. Across the westerlies-dominated central Asia, the below-cloud evaporation enhances the altitude gradient of precipitation isotopes for most areas. The findings are useful to understand the local and remote drivers of precipitation isotopes and the paleoaltimetry of stable isotopes in climate proxies.

Summer climate information recorded in tree-ring oxygen isotope chronologies from seven locations in the Republic of Korea

The Republic of Korea is characterized by its north-to-south stretch and high mountain ranges along the eastern coast, resulting in terrain with higher elevation in the east and lower in the west. These geographical features typically lead to regional climate differences, either based on latitude or from east to west. In the present study, for effectiveness, the entire Korean peninsula was divided into four regions based on the geographical features: The Northeast Coast (NEC), Central Inland (MI), Southeast Coast (SEC), and South Coast (SC). Two test sites were chosen from each region, except for the SC. The linear relationship between the altitude of sites and the mean oxygen isotope ratio (δ18O) revealed a negative correlation; the highest (1,447 m a.s.l.) and the lowest altitude (86 m a.s.l.) sites had a mean δ18O of 27.03‰ and 29.67‰, respectively. The sites selected from the same region exhibited stronger correlation coefficients (0.75–0.79) and Glk (Gleichläufigkeit) (74–83%) between the tree-ring oxygen isotope chronologies (δ18OTR chronologies) than those from different regions (0.60–0.69/70–79%). However, subtle variations in pattern were observed in the comparison period during a few selected intervals (approximately 10 years). All the regional δ18OTR chronologies exhibited positive correlations with either June or July temperatures over Korea, whereas negative correlations with regional summer precipitation and SPEI-3. Moreover, the chronologies showed notable negative correlations with the water condition of western Japan. The findings of this study can be used as a scientific reference for the study of variations of rainfall in East Asia using δ18OTR chronology.

Transmission characteristics and the factors influencing stable oxygen isotopes in precipitation, soil water, and drip water in Remi Cave, Western Hunan, China

Understanding the modern transport processes of δ18O between different waters in cave systems can provide a basis for the paleoenvironmental interpretation of records of stalagmite δ18O. We conducted monthly in-situ monitoring of local precipitation, soil water at three different soil depths, five drip water sites and environmental parameters in Remi cave, Western Hunan Province, Central China, from December 2020 to January 2024. During the monitoring period, due to differences in geological structure and ventilation effect at different locations within the cave, cave air temperature and relative humidity gradually remained constant from the cave entrance to the deep parts of the cave. On intra-annual timescales, precipitation δ18Op were depleted in the wet season and enriched in the dry season, which were mainly controlled by changes in moisture sources, upstream convective activity, and rainout effects. With the increasing thickness of the soil layer above the cave, soil water δ18Osoil gradually became negative and leveled off, possibly because shallow soil water responded sensitively to local year-round precipitation, accompanied by surface evaporation. While, deep soil water with more negative δ18O, may be mainly replenished by heavy precipitation during the wet season, but the effect of water retained in the soil from preceding years couldn’t be excluded and would be needed to verify by conducting more monitoring studies on the variation of δ18Osoil covering the long period. The mean δ18O values of drip water at sites D1 and D2 near the cave entrance were heavier than that of shallow soil water and precipitation, possibly due to evaporation from the thin overlying soil layer, together with ventilation effect. There were seasonal variations of δ18Od values at site D5 in the deep parts of the cave, which may be associated with the seasonal variation of drip rate at this site. Due to the smoothing effect of soil layer and epikarst zone, there were no seasonality in δ18Od values at sites D3 and D4, in the middle and deep parts of the cave. The mean δ18Od at these two sites were lighter than amount-weighted annual mean δ18Op, and close to the deep δ18Osoil and the wet-season amount-weighted mean δ18Op values. This suggests that the thick soil layer overlying the cave raises the precipitation threshold, causing these sites to be mainly recharged by heavy rainfall in the wet season. Therefore, the δ18O of stalagmites growing under these drip sites may mainly inherit the wet-season δ18Op signal. The short monitoring period of this study (three years) and the lack of modern calcite δ18O data, mean that caution is need in interpreting stalagmite δ18O on a timescale longer than that of this study.

Diagenetic evolution and origin of cements of the Cretaceous carbonates in the Akal Block, southeast Gulf of Mexico

The Cretaceous succession of deep marine facies (580 m-thick) of the Akal Block, southeast of the Gulf of Mexico, is dominated by lime mudstones that have been extensively affected by dolomitization. Petrographic examination, including cathodoluminescence (CL), and microthermometry reveal three phases of dolomitization that occurred as replacement and pore-/fracture-filling cements. These are (1) an early replacement dolomicrite (D1, <50 μm) with dull orange CL, (2) eu-to subhedral dolomite (D2, 60 μm–150 μm) with distinctive cloudy cores (dull CL) and clean edges (non-luminescent to dull CL), and intercrystalline porosity (6 %), and (3) anhedral vug-/fracture-filling dolomite (D3, 200 μm–500 μm) exhibiting CL zonation. The mean δ18O values of those dolomites show a trend of slight decrease from D1 to D3 (D1 = 0.0 ± 1.3‰ VPDB, D2 = −0.2 ± 1.1‰ VPDB, D3 = −1.2 ± 1.6‰ VPDB), which reflects the influence of progressive burial. The low Sr contents (101 ± 29 ppm) of D1, along with its near-micritic grain size and estimated δ18OD1fluid of parent dolomitizing fluid (−3.6‰–1.9‰ SMOW), support early-stage dolomitization at low temperature of near-surface conditions from solutions of mixed meteoric and seawaters. The higher estimated δ18Ofluid values of D2 and D3 (3.0–13.9‰ SMOW) are expected for fluids in burial settings of high temperatures. The homogenization temperatures of the primary two-phase fluid inclusions of D2 (Th = 99.1 ± 7.9 °C) and D3 (Th = 89.5 ± 11.2 °C) suggest precipitation at deeper burial settings, which is also consistent with their lower Sr contents (80 ± 14 ppm and 98 ± 33 ppm, respectively). The lower mean Th value of D3 than that of D2 likely reflects the influence of tectonic uplift that affected the Akal Block. The shale-normalized (REESN) patterns of the dolomites mimic that of modern seawater and those of D1 and D2, in particular, coincide with that of the lime mudstones (C1), thus suggesting that C1 was their precursor and also D3 originated almost from the same parent fluid that its geochemical composition evolved with progressive burial. The latest fracture-filling calcite cement (C2) postdated D3. The investigation of the products of burial diagenesis, such as fracturing, dissolution, and dolomitization events and their relationship with the tectonic uplift of the Akal Block allows a better understanding of the factors controlling the quality of Cretaceous reservoir rocks.

Spatiotemporal Absorption Features of Yellow Willow Water Usage on the Southern Edge of the Semi-Arid Hunshandak Sandland in China

The improvement of water usage efficiency and productivity, as well as the development of effective water management plans, necessitates a comprehensive understanding of how water utilization patterns in different soil layers within arid and semi-arid climates impact the capacity of plant roots to absorb water. However, there is currently no knowledge regarding the water use strategies employed by artificial yellow willow. So, we conducted a study on the hydrogen and oxygen isotopic composition of rainfall in yellow willow (Salix gordejevii) from the semi-arid region located at the southern edge of the Hunshandak Sandland in China. This study utilized measured data on xylem water, groundwater, soil moisture, and rainfall. By employing a combination of the direct comparison method and the MixSIAR model, we investigated the water uptake strategies employed by yellow willow throughout its growing season. The findings revealed that the mean δ D was highest in precipitation and lowest in groundwater, whereas the mean δ18O was highest in stem water and lowest in groundwater. The δ D and δ18O fluctuated significantly in precipitation but were steady in groundwater. Because precipitation was significantly less than evaporation, the slope and intercept were lower for the local than global atmospheric precipitation line. Water availability steadily declined with increasing depth. Lower δ18O values were caused by precipitation diluting the soil water. The MixSIAR results indicated that the primary source in May, September, and October was utilized at 19%, 18%, and 18%, respectively. In contrast, the utilization rate of each source varied considerably in June, July, and August (the primary source was utilized at 19%, 18%, and 18%, respectively). Comparatively high rates of water absorption and utilization were observed in June (19% of the total water source), July (18%), and August (23%). Therefore, the vertical distribution of the root system and variations in the soil water content regulate water usage for the yellow willow. To prevent excessive water usage and promote ecosystem restoration with artificial yellow willow plantations in water-limited desert settings, policy makers should consider the patterns of plant water use and soil water availability. By selecting drought-adapted plant species and optimizing irrigation management, it is possible to reduce water wastage and ensure that water is used efficiently for revegetation and ecosystem restoration, avoiding overuse of water and maintaining the sustainability of revegetation in water-stressed desert areas.

Shallow carbonate geochemistry in the Bahamas since the last interglacial period

The carbon isotope composition (δ13C) of ancient shallow-water carbonates frequently is used to reconstruct changes to Earth's global carbon cycle and to perform chemostratigraphic correlation. However, previous work demonstrates that local banktop processes also exert an important control on shallow carbonate δ13C as well as other isotope systems like δ18O. To effectively interpret ancient δ13C records, we must understand how both global carbon cycle perturbations and changes to local conditions are translated to the stratigraphic record. Modern environments, while imperfect analogues, can serve as a guide for interpreting physical and geochemical records of more ancient environments. Shallow carbonate strata from the most recent Pleistocene glacial cycles, which drove high-amplitude perturbations to sea level, temperature, and pCO2 without significantly altering the δ13C of global-mean seawater DIC, present an opportunity to begin untangling signals of global and local processes. However, the geochemistry of Pleistocene platform carbonates largely was overprinted by dissolution and meteoric diagenesis during glacial sea level lowstands. To understand how shallow carbonate geochemistry has changed during the Pleistocene, we instead look to the periplatformal slope and proximal basins. These deep environments serve as a refuge for carbonate produced on the shelf and exported to the slope, and contain a record of shallow carbonate geochemistry that persists across glacial cycles. We study 21 short piston cores from around Bahamian platforms to quantify differences in banktop production and geochemistry between the Holocene interglacial, the last glacial period, and the last interglacial (LIG) period. We show that mud production persists on the periplatformal slopes during the last glacial period, but differences in geochemistry between glacial and interglacial carbonates are a complex function of surface conditions and diagenesis. In contrast, Holocene and LIG carbonates show no evidence of post-depositional alteration, and offer the chance to study differences in δ13C and δ18O between interglacials. We find that while the δ13C of aragonite mud is the same during the Holocene and LIG, the LIG carbonate factory may have delivered more aragonite mud to the periplatform. In addition, the mean δ18O of this mud is elevated compared to the Holocene. We posit that these differences are caused by changes to regional climate and LIG surface conditions.

The Chemical Properties and Formation Mechanisms of Shallow Groundwater in the Guohe River Basin, China

An investigation of the distribution and control factors of groundwater is significant for the rational exploitation and utilization of groundwater resources. This study analyzes the hydrogeochemical processes and control factors of shallow groundwater in the Guohe River Basin. A hydrogeological survey was conducted, and hydrochemical and hydrogen–oxygen isotopic data of 125 samples of surface water and groundwater were analyzed. The results showed that the total dissolved solid (TDS) content in shallow groundwater was 138–2967 mg/L, with an average of 831 mg/L. A decline in the TDS was observed from the upper to the lower reaches. The contents of the anions and cations in the shallow groundwater were in the order HCO3− &gt; Cl− &gt; SO42− and Na+ and K+ &gt; Mg2+ &gt; Ca2+, respectively. The cation exchange increased the aqueous concentrations of Na+ and K+, and the TDS content was highly correlated with the contents of Na+, Cl−, and SO42− ions. The δD and δ18O values in shallow groundwater increased from the upper to the lower reaches, with the mean δD values being −59.72‰, −53.58‰, and −47.17‰ and the mean δ18O values being −8.33‰, −7.37‰, and −6.43‰. The contribution rates of the recharge source, evaporation, and water–rock interaction to the groundwater TDS concentration were 20.4%, 29.5%, and 50.1%, respectively. The water–rock interaction dominated the formation of shallow groundwater in the Guohe River Basin. The dissolution of salt rock and gypsum contributed to ion formation in shallow groundwater. The research findings can be used to improve the groundwater quality in the Guohe River Basin.

Hydrochemical and isotopic studies providing a new functional model for the coastal aquifers in Douala Coastal Sedimentary Basin (DCSB)/Cameroon

A new conceptual model of the hydrogeological systems in Cameroon's Douala Coastal Sedimentary Basin (DCSB) was constructed. The model is based upon the basin's known geology, plus data from recent field campaigns that allowed the collection of rainwater and groundwater samples for analyses of stable isotopes (δ2H, δ18O, δ13C), radiogenic isotopes (3H, 14C), and water chemistry. Aquifer characteristics that were thereby deciphered include recharge, isotopic distributions, residence times, and mixing processes. Rainfall samples (mean δ18O = −2.0 ‰; mean δ2H = −6.80 ‰; weighted mean = −2.4 ‰ δ18O, −9.85 ‰ δ2H) scatter along two distinct lines, thus indicating that local rainfall events undergo processes during convective events, variability in humidity, amount effects, and seasonal variations. Stable isotope values of river water samples are close to the weighted mean of local precipitation, with some downstream enrichment. The Quaternary/Mio-Pliocene superficial aquifer system (depth < 70 m) and the intermediate Oligocene/Upper Eocene aquifer system (depth: 70 to 200 m) exhibit evidence of similar fractionation processes through an enrichment gradient of δ-values. The enrichment is more pronounced at the top of the superficial aquifer, which is very exposed to direct rainfall water infiltration, evaporation, and amount effects. The depth profiles of δ-values coupled to water chemistry and tritium contents, evidence leakage between (i) the superficial system's Quaternary alluvium sands and Mio-Pliocene sands; and (ii) the superficial and intermediate systems. Thus, the aquifers that contain modern, post nuclear groundwater are characterized by flow exchanges and direct recharge from rainfall events. In contrast, the Upper Eocene system has depleted δ-values and lower bicarbonate contents, suggesting not only that this system was recharged by rapid infiltration (with limited effect of evaporation), but that this recharge occurred during a cooler time in the past. The residence times (computed from 14C dates) indicate uncorrected ages ranging from hundreds to thousands of years.

O- and H-isotope record of Cape Town rainfall from 1996 to 2022: the effect of increasing temperature, and the ‘water crisis’ of 2015 to 2018

Abstract The O- and H-isotope composition of rainfall collected monthly at the University of Cape Town (UCT) between 2009 and 2022 has been added to a previously published data set from 1996 to 2008 to make a continuous 27 year record. Monthly rainfall over the 27 year period has a range in δD and δ18O values from -57 to +21‰ and -8.1 to +3.5‰, respectively, and shows limited but discernible temperature and amount effects. The 27 year rainfall record defines a Local Meteoric Water Line (LMWL) whose equation is δD = 5.96*δ18O + 7.00 (r = 0.88), a slight change from the LMWL of 1996 to 2009 (6.41*δ18O + 8.89). Annual rainfall at UCT has varied from ~992 mm (2017) to 1996 mm (2013), with no systematic change in annual rainfall amount (r = −0.16). However, from 2015 to 2022, the average annual rainfall of 1 145 mm has been below the 27 year average of 1 313 mm. Mean monthly temperature has increased from 1996 to 2022 (r = 0.53), and the weighted mean annual δD and δ18O values have increased by ~4‰ (r = 0.53) and ~0.5‰ (r = 0.64), respectively, over the 27 years. The UCT data and the data for 1961 to 1974 from Cape Town International Airport plot around the same LMWL, with an average deuterium excess (d-excess) of 13.78 and 12.95, respectively. Natural springs in the area plot close to the LMWL with an average d-excess of 14.15, whereas local well-point and borehole water samples generally plot below the LMWL with an average d-excess of 10.65. These differences can be explained by relatively rapid recharge of springs and slower recharge of groundwater, with the latter containing an additional component, that could either be rainwater that fell during a period of hotter drier climate or, more probably, municipal mains water. Long-term monitoring of groundwater and spring water as well as rainwater would be of great help in assessing the sustainability of groundwater use, among other important questions.

Altitude Correction of GCM-Simulated Precipitation Isotopes in a Valley Topography of the Chinese Loess Plateau

Altitude is one of the important factors influencing the spatial distribution of precipitation, especially in a complex topography, and simulations of isotope-enabled climate models can be improved by altitude correlation. Here we compiled isotope observations at 12 sites in Lanzhou, and examined the relationship between isotope error and altitude in this valley in the Chinese Loess Plateau using isoGSM2 isotope simulations. Before altitude correction, the performance using the nearest four grid boxes to the target site is better than that using the nearest box; the root mean square error in δ18O using the nearest four grid boxes averagely decreases by 0.37‰ compared to that using the nearest grid boxes, and correlation coefficient increases by 0.05. The influences of altitude on precipitation isotope errors were examined, and the linear relationship between altitude error and isotope simulations was calculated. The strongest altitude isotopic gradient between δ18O mean bias error and altitude error is in summer, and the weakest is in winter. The regression relationships were used to correct the simulated isotope composition. After altitude correction, the root mean square error decreases by 1.21‰ or 0.86‰ using the nearest one or four grid boxes, respectively, and the correlation coefficient increases by 0.13 or 0.08, respectively. The differences between methods using the nearest one or four grids are also weakened, and the differences are 0.02‰ for root mean square error and −0.01 for the correlation coefficient. The altitude correction of precipitation isotopes should be considered to downscale the simulations of climate models, especially in complex topography.

Moisture sources and climatic effects controlling precipitation stable isotope composition in a western Mediterranean island (Pianosa, Italy)

The Mediterranean basin is indicated as a hot spot of climate change, which is an area whose climate is especially responsive to variations. The insular environment is one of the most threatened by the current climate change, especially in terms of drought events, with serious consequences for water scarcity and water stress. This issue is even enhanced in small islands, whose ecosystems are among more sensitive to climatic changes and water availability. The stable isotope composition of hydrogen (δ2H) and oxygen (δ18O) in precipitation is globally recognized as a powerful natural tracer in the water cycle and represents the starting point to investigate hydrological processes. The understanding of the prevailing factors that drive the isotopic variability of precipitation in the Mediterranean is therefore essential to unravel the hydrological processes and to ensure proper and sustainable management of potentially vulnerable resources to climate change. Here, we discuss the results of multi-year isotopic monitoring in the period 2014–2021 of monthly precipitation collected on Pianosa Island (Italy), a small island located in the northern Tyrrhenian (western Mediterranean). The lower slope and intercept of the Local Meteoric Water Line of the island compared to the Global Meteoric Water Line indicated warmer and drier climatic conditions, suggesting the existence of sub-cloud evaporation processes of raindrops during precipitation, especially in summer. The mean δ18O of precipitation was lower with respect to other sites placed at higher elevation in this Mediterranean region, due to the lack of summer precipitation which were generally enriched in heavy isotopes. Temperature and amount effects may explain part of the δ18O variability observed at the monthly and seasonal scale. An HYSPLIT-based moisture uptake analysis indicated the area between the western Mediterranean basin, Italy, and the Adriatic Sea as the region that supplied most of the humidity associated with monthly precipitation samples on Pianosa Island. Less moisture was picked from the northwestern areas of Europe, the North Atlantic Ocean, the proximal Atlantic Ocean, the Iberian Peninsula and North Africa. Consistently with the rainout effect, the higher the moisture fraction picked from the more proximal regions, the more positive the δ18O of precipitation occurring on Pianosa Island; conversely, the higher the percentage of moisture sourced from more distal regions, the more negative the δ18O. A multiple linear model was proposed to predict the δ18O of monthly precipitation from temperature, precipitation amount and moisture origin data, which explained 45% of the δ18O variability. The deuterium excess variability on the island was partly controlled by the local climatic variables, whose effect potentially modifies the original d-excess signature imprinted at the moisture source. No relationship was found between the precipitation deuterium excess and moisture sources, suggesting that more attention should be paid when using the deuterium excess as a tracer of moisture origin, especially in the Mediterranean.

Apatite U-Pb geochronological and geochemical constraints on the Mazhala Au-Sb deposit in South Tibet, China

Due to lack of suitable datable minerals, timing of the Mazhala Au-Sb deposit in South Tibet remains unresolved. Also, previous geochemical studies suggest an ambiguous conclusion on the source of ore-forming fluid, with diverse viewpoints of magmatic water, mixed magmatic–meteoric water or mixed metamorphic-meteoric water, respectively. In this study, hydrothermal apatite U-Pb dating constrains the Mazhala Au-Sb deposits to have formed at ca. 35 Ma, consistent with the timing of Eocene leucogranite intrusions in South Tibet. The Mazhala Au-Sb mineralization occurred during the main collision stage to late collision stage in Eocene to Oligocene, and predates the Cuonadong Sn mineralization in South Tibet which occurred during the crustal extension stage in Miocene. An EPMA analysis shows the apatite belongs to fluorapatite that have F contents in the dominant range of 2–4%. The Mazhala apatite displays an overall MREE enrichment and LREE- and HREE depletion distribution pattern with positive Eu anomalies, and have δ18O values dominantly in the range of 16-19‰ with a mean δ18O‰ value of 17.74 ± 0.15‰ (2σ), which suggest a relatively homogenous chemical composition for the ore forming fluids with a relatively stable temperature regime. In combination with previous temperature estimates, the ore forming fluids are inferred to have a dominant δ18Ofluid value of 13-16‰, which overlap with metamorphic water.The combined age and geochemical data suggest Eocene leucogranite intrusion induced metamorphic fluid circulation within Paleozoic-Mesozoic metamorphic sedimentary rocks around the study area, which prompted precipitation of fluorapatite together with Au-Sb bearing pyrite and arsenopyrite in the Mazhala Au-Sb deposit. This study puts an emphasis on the contribution of metamorphic water as the dominant ore-forming fluid origin in the Mazhala Au-Sb deposit. In contrast, the close spatial and temporal association of the Cuonadong Sn polymetallic mineralization with the Miocene leucogranite intrusion-induced hydrothermal activity in the Cuonadong dome demonstrates that their genetic links, as was also evidenced by C-H-O-S-Pb isotope data that point to a dominant contritution of Miocene leucogranite to ore-forming fluid and materials for the Cuonadong Sn polymetallic ore deposits. As such, from a perspective of ore-forming model in South Tibet it be seen that during the main collision stage to late collision stage ore-forming fluid and materials were mainly derived from metamorphic water, whereas during the crustal extension stage, ore-forming fluid and materials were dominantly sourced from magmatic water.

Uncorrected soil water isotopes through cryogenic vacuum distillation may lead to a false estimation on plant water sources

Abstract Successful use of stable isotopes (δ2H and δ18O) in ecohydrological studies relies on the accurate extraction of unfractionated water from different types of soil samples. Cryogenic vacuum distillation (CVD) is a common laboratory‐based technique used for soil water extraction; however, the reliability of this technique in reflecting soil water δ2H and δ18O is still of concern. This study examines the reliability of a newly developed automatic cryogenic vacuum distillation (ACVD) system. We further assessed the impacts of extraction parameters (i.e. extraction time, temperature and vacuum) and soil properties on the recovery of soil water δ2H and δ18O for the ACVD and traditional cryogenic vacuum distillation (TCVD) systems. Finally, we investigated the potential influence of CVD (ACVD and TCVD) technique on the prediction of plant water uptake through a sensitivity analysis. Both ACVD and TCVD similarly extracted water from the rewetted soils, but none of the CVD systems successfully recovered the isotopic signatures of doped water from soil materials. Mean δ2H offsets of extracted soil water were −2.6 ± 1.3‰ and −2.4 ± 1.7‰ for ACVD and TCVD, respectively, while mean δ18O offsets were −0.16 ± 0.14‰ and −0.39 ± 0.37‰. The isotopic offsets of CVD systems were positively correlated with soil clay content, and negatively correlated with soil water content. Using corrected soil data (with CVD offsets) could improve the prediction of plant water uptake based on its high correlation with the environmental factors. This study identifies the isotopic offsets of CVD systems (i.e. ACVD and TCVD) and provides possible solutions for better predicting plant water sources. Even though, the wide use of CVD techniques probably induce noticeable uncertainties in the prediction of plants water uptake depths. The dataset of soil water extraction in this study will have implications for the technological development of CVD techniques.

Western South Atlantic seasonal variability recorded in a mid-deglacial bivalve from the outer Uruguayan continental shelf

The oceanographic dynamics on the continental shelf off southeastern South America are primarily controlled by the southward-flowing warm Brazil Current, converging with the northward-directed cold Malvinas (Falkland) Current, and interacting with the continental discharge of the Plata River. The seasonally reversing regional wind field together with the seasonal cycle of riverine discharge, determines which of these three components provides the dominant forcing. The Uruguayan shelf is thus located in a transitional zone that extends from the region influenced by the Brazil-Malvinas Confluence (BMC) in the open ocean to the Subtropical Shelf Front (STSF) on the continental shelf. Understanding how the resulting oceanic seasonal variability responded to different climatic boundary conditions may shed light on its future behavior. This study presents the first reconstruction of mid-deglacial seasonal hydrographic variability on the continental shelf off southeastern South America in seasonal resolution based on stable oxygen and carbon isotopes (δ18O, δ13C) from a thick-walled shell of a long-living bivalve. The mid-deglacial (14.3 cal ka BP; Bølling-Allerød interstadial) Retrotapes exalbidus bivalve shows a mean δ18O of 3.27 ± 0.42‰ (2.50 ± 0.42‰ when corrected for changes in global ice volume) and a seasonal δ18O amplitude of 1.69‰ for raw isotopic excursions. Moreover, the δ13C exhibits abrupt negative peaks coincident with more negative δ18O values that indicate seasons of elevated freshwater discharge. Finally, the growth rate of the bivalve suggests that the specimen was closer to the metabolically optimum than modern individual of this species from southern South America. Combining biogeographic and ecologic information with these isotopic data, the results point to colder waters and a slightly lower mid-deglacial seasonal amplitude in temperature compared with modern conditions at this shelf site. Because of the northward-displaced Plata River mouth during deglacial times, negative δ13C peaks are expected to reflect an influence of non-point freshwater sources in the form of small fluvial distributaries along the paleo-coast. Most of this signal may, however, be driven by seasonal metabolic effects associated with low ambient water temperatures related to a shallow-water environment located closer to the respective paleo-coastline due to the low sea level at those times.

Syn-collisional extension and Ni-Cu sulfide-bearing mafic magma emplacement along the Irtysh Shear Zone in the Central Asian Orogenic Belt

The Central Asian Orogenic Belt was characterized by a long-lived accretionary history from Neoproterozoic to Paleozoic, followed by a collisional phase of orogeny in the latest Paleozoic with the closure of the Paleo-Asian Ocean. However, Ni-Cu sulfide mineralization associated with mantle-derived mafic magmas in the Central Asian Orogenic Belt was relatively short-lived. Thus, the coherence between the protracted accretionary and collisional history of the Central Asian Orogenic Belt and the transient Ni-Cu sulfide-bearing mafic magmatism is difficult to understand. We investigated both Ni-Cu sulfide-bearing and sulfide-barren mafic intrusions along a major suture zone of the Irtysh Shear Zone within the Central Asian Orogenic Belt to understand the origin of short-lived Ni-Cu sulfide mineralization in the context of the collision of two arc systems of the Chinese Altai and the East Junggar regions. Zircon U-Pb dating results for the Ni-Cu sulfide-bearing mafic intrusions show that they were only emplaced briefly from 290 Ma to 280 Ma after subduction terminated (ca. 312 Ma), substantially earlier than Ni-Cu sulfide-barren mafic intrusions (ca. 274−271 Ma). Zircons from all Ni-Cu sulfide-bearing/barren mafic intrusions yield arithmetic mean δ18O values ranging from 6.0 to 7.4‰, higher than the typical mantle value (∼5.3‰), which indicates their derivation from a metasomatized mantle with inherited subduction signatures. The δ11B values of the rocks from these intrusions range from −17.2‰ to −2.3‰, significantly lower than those of typical volcanic arcs (up to +18‰), which suggests that the metasomatic components in the mantle source were dominated by sediment melts rather than fluids. The pressure calculated using amphibole thermobarometry shows that the parental magmas of the Ni-Cu sulfide-bearing mafic intrusions were emplaced shallower than those of Ni-Cu sulfide-barren mafic intrusions. This is supported by the fact that the country rocks of the Ni-Cu sulfide-bearing mafic intrusions are shallow crustal lithologies of sedimentary rocks, whereas the country rocks of the barren intrusions are relatively deeper crustal lithologies of high-grade metamorphic rocks and granites. Olivine from the Ni-Cu sulfide-barren intrusions has forsterite (Fo) contents ranging from 51.1 mol% to 76.2 mol% and Ni from 196 ppm to 1312 ppm. These values are much lower than those of the Ni-Cu sulfide-bearing intrusions (Fo = 72.6−79.4 mol%; Ni = 1022−1925 ppm), which indicates that the parental magmas of the barren intrusions may have experienced extensive olivine crystallization. Notably, the emplacement of the Ni-Cu sulfide-bearing mafic intrusions (290−280 Ma) coincidently overlapped with a transient period of orogen-parallel extension (ca. 295−280 Ma) in response to the collision of the Chinese Altai and East Junggar. Syn-collisional extension may lead to the rapid ascent and emplacement of mantle-derived mafic magmas at a shallower crustal level, preventing significant loss of Ni from olivine crystallization at depths, and thus providing metal-rich mafic magmas for potential Ni-Cu sulfide mineralization. Our results highlight that syn-collisional extension is an essential geodynamic mechanism that controls the emplacement of Ni-Cu sulfide-bearing mafic intrusions in orogenic belts.

Tracking Magma‐Crust‐Fluid Interactions at High Temporal Resolution: Oxygen Isotopes in Young Silicic Magmas of the Taupō Volcanic Zone

AbstractOxygen isotopes are useful for tracing interactions between magmas, crustal rocks and surface‐derived waters. We use them here to consider links between voluminous silicic magmatism and large‐scale hydrothermal circulation in New Zealand's central Taupō Volcanic Zone (TVZ). We present &gt;350 measurements of plagioclase, quartz, hornblende and groundmass glass δ18O values from 40 eruptions from three discrete magmatic systems: Ōkataina and Taupō calderas, and the smaller Northeast Dome system. For each mineral, mean δ18O values vary by ∼1‰ (δ18Oplag = +6.7–7.8‰, δ18Oqtz = +7.7–+8.7‰, δ18Ohbl = +5.4–+6.4‰, δ18Oglass = +7.1–+8.0‰), and inter‐mineral fractionations mostly reflect high‐temperature equilibria. Outliers (e.g., ∼+6‰ or &gt;+10‰ plagioclase) represent contaminants incorporated on short‐enough timescales to preserve disequilibrium (∼102 yrs for plagioclase). Melt δ18O values calculated from phenocrysts are ∼+7.3–+8.0‰. Where multiple magmas were involved in the same eruption their δ18Omelt values are indistinguishable, implying that their parental mushes were isotopically well‐mixed. However, small (≤0.5‰) but consistent δ18Omelt value gradients occur over millennial timescales at Ōkataina and Taupō, with short‐term ∼0.4–0.5‰ decreases in δ18Omelt values over successive post‐caldera eruptions correlating with increases in 87Sr/86Sr. These changes reflect tens of percent assimilation of a mixture of hydrothermally altered silicic plutonic material and higher‐87Sr/86Sr greywacke. These examples represent the first evidence for assimilation of altered crust into TVZ magmas. The subtle and short‐lived isotopic signals of these interactions are only recognized through the high temporal resolution of the TVZ eruptive record and complementary radiogenic isotope data. Similar interactions may have been obscured in other nominally high‐ or normal‐δ18O magmatic systems.

Late Neoarchean crust–mantle interaction and tectonic implications in western Shandong Province, North China Craton: Evidence from a granodiorite pluton and associated magmatic enclaves

This study presents original data on the microstructure, mineralogy, whole–rock geochemistry, and isotopes of a Late Neoarchean granodiorite pluton and associated mafic microgranular enclaves (MMEs) in western Shandong Province of the eastern part of the North China Craton (NCC). The MMEs hosted in the granitic pluton are more mafic, and amphibole and biotite are the predominant mafic minerals. The sharp contacts and fine–grained textures of the MMEs indicated a rapid quenching process during a magma mingling event. Zircon and titanite UPb dating indicated that the host granodiorite pluton formed at circa 2546 Ma, and the MMEs formed at approximately 2538 Ma. The dating results indicated a crust–mantle interaction. Geochemically, the host pluton has high silica contents (SiO2 = 67.68–68.51 wt%), Mg# values (51–53) and low Aluminum Saturation Index ratios (A/CNK = 0.92–0.95), which are indicative of I–type granitic affinity. Moreover, they show arc–type trace element features characterized by enrichment of light rare earth elements (LREEs) and large ion lithophile elements (LILEs) (e.g., Cs, Rb, K, and Pb), but depletion in high field strength elements (HFSEs) (e.g., Nb, Ta, P, and Ti). In contrast, the MMEs have lower silica contents (SiO2 = 52.07–56.03 wt%), but higher Mg# values (56–68) and Cr, Ni and Cu contents of 454–633, 159–225 and 98.3–207 ppm, respectively, corresponding to a mantle component. Similar to the host pluton, the MMEs show enrichment of LREEs and LILEs, but depletion of HFSEs, reflecting an arc–type magma origin. Isotopically, the host granodiorite and MMEs both have relatively depleted NdHf isotopic compositions, with εNd(t) values of −0.67–1.22 and − 0.31–1.89, and εHf(t) values of 5.63–2.03 and − 0.18–4.48, respectively. Zircon oxygen isotopes reveal that the host granodiorite and MMEs have slightly higher oxygen isotopes than the normal mantle zircon (5.3‰), with mean δ18O values of 6.35‰ and 6.81‰, respectively, indicating assimilation of minor supracrustal materials in the magma source. However, a combination of Nd–Hf–O isotopes suggests that the MMEs haveslightly enriched isotopic compositions than the host granodiorite pluton, suggesting a “reversed isotope” feature. This important finding also enlightens that “reversed isotope” phenomenon represented magma mingling between isotopically relatively enriched mafic and relatively depleted felsic magmas. Furthermore, the Ce4+/Ce3+ ratios in zircon from the host granodiorite and MMEs indicated a low oxygen fugacity, with mean values of 24.25 and 39.69, respectively. Although the host granodiorite and MMEs have the similar oxygen fugacity ranges, the MMEs have high Cu concentrations of 98.3–207 ppm (mean value = 151 ppm), indicating a Cu–rich magma source. Integrating the petrogenesis of the Late Neoarchean granodiorite pluton and associated MMEs, it is inferred that the host granodiorite pluton was derived from partial melting of mafic lower crust with minor insignificant ancient materials, and the MMEs were derived from partial melting of depleted mantle metasomatized by sulfur–rich sediment melts/fluids (chalcophile element). Tectonically, our samples reveal an affinity with an arc in a syn–collisional setting. Together with regional rock assemblages, we suggest that the onset of plate tectonics (subduction) was to occur by the end of Neoarchean in the NCC.

Water sources and recharge mechanisms of the Yarlung Zangbo River in the Tibetan Plateau: Constraints from hydrogen and oxygen stable isotopes

Investigations of the hydrological processes of the Yarlung Zangbo River (YZR), a sizeable highland river on the southern margin of the Tibetan Plateau (TP), face challenging constraints owing to the complexities in the water sources and cycle pathway. A total of 156 stream water samples from the YZR Basin were collected and analyzed to elucidate our understanding of regional hydroclimate. The isotopic values were more negative in the uppermost reaches of the YZR and progressively increased along the river, with mean δ18O and δD values at − 15.16 ‰ and − 111.58 ‰, respectively. A clear negative gradient of − 0.19 ‰/100 m was obtained from the datasets, which is consistent with that of the southeastern TP. The linear regression for all samples was δD = 8.64 × δ18O + 19.80 (R2 = 0.93, n = 192), and the slope and intercept were highly similar to those of the meteoric water line in the TP, suggesting that precipitation plays a dominant role in stream recharge sources. Combining the results from the air mass trajectory model and the d-excess values, we found that moisture in the YZR originates mainly from the Bay of Bengal, with westerly mixing in the upper reaches. In addition to precipitation, glacial meltwater and geothermal waters were considered significant water sources for the YZR based on geographic and geological conditions. The proportion of glacial meltwater in the uppermost reaches was calculated to range from 68.2 % to 70.4 % using the end-member mixing analysis method. The contribution ratios of precipitation, glacial meltwater, and geothermal water were almost equal in the middle of the basin, while glacial meltwater contribution dominated in the downstream regions. These findings will significantly impact understanding of the hydrological cycle in the study area and provide vital baseline knowledge about the contribution proportion of different recharge sources in the YZR basin.