ALBERT

Description: We present a workflow to estimate geostatistical aquifer parameters from pumping test data using the Python package welltestpy. The procedure of pumping test analysis is exemplified for two data sets from the Horkheimer Insel site and from the Lauswiesen site, Germany. The analysis is based on a semi‐analytical drawdown solution from the upscaling approach Radial Coarse Graining, which enables to infer log‐transmissivity variance and horizontal correlation length, beside mean transmissivity, and storativity, from pumping test data. We estimate these parameters of aquifer heterogeneity from type‐curve analysis and determine their sensitivity. This procedure, implemented in welltestpy, is a template for analyzing any pumping test. It goes beyond the possibilities of standard methods, for example, based on Theis' equation, which are limited to mean transmissivity and storativity. A sensitivity study showed the impact of observation well positions on the parameter estimation quality. The insights of this study help to optimize future test setups for geostatistical aquifer analysis and provides guidance for investigating pumping tests with regard to aquifer statistics using the open‐source software package welltestpy.

Description: Article impact statement: We present a workflow to infer parameters of subsurface heterogeneity from pumping test data exemplified at two sites using welltestpy.

Description: German Federal Environmental Foundation (DBU) http://dx.doi.org/10.13039/100007636

Description: In designed experiments, different sources of variability and an adequate scale of measurement need to be considered, but not all approaches in common usage are equally valid. In order to elucidate the importance of sources of variability and choice of scale, we conducted an experiment where the effects of biochar and slurry applications on soil properties related to soil fertility were studied for different designs: (a) for a field‐scale sampling design with either a model soil (without natural variability) as an internal control or with composited soils, (b) for a design with a focus on amendment variabilities, and (c) for three individual field‐scale designs with true field replication and a combined analysis representative of the population of loess‐derived soils. Three silty loam sites in Germany were sampled and the soil macroaggregates were crushed. For each design, six treatments (0, 0.15 and 0.30 g slurry‐N kg−1 with and without 30 g biochar kg−1) were applied before incubating the units under constant soil moisture conditions for 78 days. CO2 fluxes were monitored and soils were analysed for macroaggregate yields and associated organic carbon (C). Mixed‐effects models were used to describe the effects. For all soil properties, results for the loess sites differed with respect to significant contributions of fixed effects for at least one site, suggesting the need for a general inclusion of different sites. Analysis using a multilevel model allowed generalizations for loess soils to be made and showed that site:slurry:biochar and site:slurry interactions were not negligible for macroaggregate yields. The use of a model soil as an internal control enabled observation of variabilities other than those related to soils or amendments. Experiments incorporating natural variability in soils or amendments resulted in partially different outcomes, indicating the need to include all important sources of variability. Highlights Effects of biochar and slurry applications were studied for different designs and mixed‐effects models were used to describe the effects. Including an internal control allowed observation of, e.g., methodological and analytical variabilities. The results suggested the need for a general inclusion of different sites. Analysis using a multilevel model allowed generalizations for loess soils. The results indicated the need to include all important sources of variability.

Description: Temperate forest soils are often considered as an important sink for atmospheric carbon (C), thereby buffering anthropogenic CO2 emissions. However, the effect of tree species composition on the magnitude of this sink is unclear. We resampled a tree species common garden experiment (six sites) a decade after initial sampling to evaluate whether forest floor (FF) and topsoil organic carbon (Corg) and total nitrogen (Nt) stocks changed in dependence of tree species (Norway spruce—Picea abies L., European beech—Fagus sylvatica L., pedunculate oak—Quercus robur L., sycamore maple—Acer pseudoplatanus L., European ash—Fraxinus excelsior L. and small‐leaved lime—Tilia cordata L.). Two groups of species were identified in terms of Corg and Nt distribution: (1) Spruce with high Corg and Nt stocks in the FF developed as a mor humus layer which tended to have smaller Corg and Nt stocks and a wider Corg:Nt ratio in the mineral topsoil, and (2) the broadleaved species, of which ash and maple distinguished most clearly from spruce by very low Corg and Nt stocks in the FF developed as mull humus layer, had greater Corg and Nt stocks, and narrow Corg:Nt ratios in the mineral topsoil. Over 11 years, FF Corg and Nt stocks increased most under spruce, while small decreases in bulk mineral soil (esp. in 0–15 cm and 0–30 cm depth) Corg and Nt stocks dominated irrespective of species. Observed decadal changes were associated with site‐related and tree species‐mediated soil properties in a way that hinted towards short‐term accumulation and mineralisation dynamics of easily available organic substances. We found no indication for Corg stabilisation. However, results indicated increasing Nt stabilisation with increasing biomass of burrowing earthworms, which were highest under ash, lime and maple and lowest under spruce. Highlights We studied if tree species differences in topsoil Corg and Nt stocks substantiate after a decade. The study is unique in its repeated soil sampling in a multisite common garden experiment. Forest floors increased under spruce, but topsoil stocks decreased irrespective of species. Changes were of short‐term nature. Nitrogen was most stable under arbuscular mycorrhizal species.

Description: Deutsche Forschungsgemeinschaff (DFG)

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Sadai, S., Condron, A., DeConto, R., & Pollard, D. Future climate response to Antarctic Ice Sheet melt caused by anthropogenic warming. Science Advances, 6(39), (2020): eaaz1169, doi:10.1126/sciadv.aaz1169.

Description: Meltwater and ice discharge from a retreating Antarctic Ice Sheet could have important impacts on future global climate. Here, we report on multi-century (present–2250) climate simulations performed using a coupled numerical model integrated under future greenhouse-gas emission scenarios IPCC RCP4.5 and RCP8.5, with meltwater and ice discharge provided by a dynamic-thermodynamic ice sheet model. Accounting for Antarctic discharge raises subsurface ocean temperatures by 〉1°C at the ice margin relative to simulations ignoring discharge. In contrast, expanded sea ice and 2° to 10°C cooler surface air and surface ocean temperatures in the Southern Ocean delay the increase of projected global mean anthropogenic warming through 2250. In addition, the projected loss of Arctic winter sea ice and weakening of the Atlantic Meridional Overturning Circulation are delayed by several decades. Our results demonstrate a need to accurately account for meltwater input from ice sheets in order to make confident climate predictions.

Description: This research was supported by the NSF Office of Polar Programs through NSF grant 1443347, the Biological and Environmental Research (BER) division of the U.S. Department of Energy through grant DE-SC0019263, the NSF through ICER 1664013, and by a grant to the NASA Sea Level Science Team 80NSSC17K0698.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Foukal, N. P., Gelderloos, R., & Pickart, R. S. A continuous pathway for fresh water along the East Greenland shelf. Science Advances, 6(43), (2020): eabc4254, doi:10.1126/sciadv.abc4254.

Description: Export from the Arctic and meltwater from the Greenland Ice Sheet together form a southward-flowing coastal current along the East Greenland shelf. This current transports enough fresh water to substantially alter the large-scale circulation of the North Atlantic, yet the coastal current’s origin and fate are poorly known due to our lack of knowledge concerning its north-south connectivity. Here, we demonstrate how the current negotiates the complex topography of Denmark Strait using in situ data and output from an ocean circulation model. We determine that the coastal current north of the strait supplies half of the transport to the coastal current south of the strait, while the other half is sourced from offshore via the shelfbreak jet, with little input from the Greenland Ice Sheet. These results indicate that there is a continuous pathway for Arctic-sourced fresh water along the entire East Greenland shelf from Fram Strait to Cape Farewell.

Description: Funding for this work comes from the NSF under grant numbers OCE-1756361 and OCE-1558742 (N.P.F. and R.S.P.) and grant numbers OCE-1756863 and OAC-1835640 (R.G.).

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Liu, C.-Z., Dick, H. J. B., Mitchell, R. N., Wei, W., Zhang, Z.-Y., Hofmann, A. W., Yang, J.-F., & Li, Y. Archean cratonic mantle recycled at a mid-ocean ridge. Science Advances, 8(22), (2022): eabn6749, https://doi.org/10.1126/sciadv.abn6749.

Description: Basalts and mantle peridotites of mid-ocean ridges are thought to sample Earth’s upper mantle. Osmium isotopes of abyssal peridotites uniquely preserve melt extraction events throughout Earth history, but existing records only indicate ages up to ~2 billion years (Ga) ago. Thus, the memory of the suspected large volumes of mantle lithosphere that existed in Archean time (〉2.5 Ga) has apparently been lost somehow. We report abyssal peridotites with melt-depletion ages up to 2.8 Ga, documented by extremely unradiogenic 187Os/188Os ratios (to as low as 0.1095) and refractory major elements that compositionally resemble the deep keels of Archean cratons. These oceanic rocks were thus derived from the once-extensive Archean continental keels that have been dislodged and recycled back into the mantle, the feasibility of which we confirm with numerical modeling. This unexpected connection between young oceanic and ancient continental lithosphere indicates an underappreciated degree of compositional recycling over time.

Description: This study was financially supported by the National Science Fund for Distinguished Young Scholars 42025201 (to C.-Z.L.), the National Key Research and Development Project of China 2020YFA0714801 (to C.-Z.L.), the Strategic Priority Research Program of the Chinese Academy of Sciences XDA13010106 (to C.-Z.L.), the Strategic Priority Research Program of the Chinese Academy of Sciences XDB42020301 (to C.-Z.L.), and NSF grants 2114652 and 1657983 (to H.J.B.D.).

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Gomaa, F., Utter, D. R., Powers, C., Beaudoin, D. J., Edgcomb, V. P., Filipsson, H. L., Hansel, C. M., Wankel, S. D., Zhang, Y., & Bernhard, J. M. Multiple integrated metabolic strategies allow foraminiferan protists to thrive in anoxic marine sediments. Science Advances, 7(22), (2021): eabf1586, https://doi.org/10.1126/sciadv.abf1586.

Description: Oceanic deoxygenation is increasingly affecting marine ecosystems; many taxa will be severely challenged, yet certain nominally aerobic foraminifera (rhizarian protists) thrive in oxygen-depleted to anoxic, sometimes sulfidic, sediments uninhabitable to most eukaryotes. Gene expression analyses of foraminifera common to severely hypoxic or anoxic sediments identified metabolic strategies used by this abundant taxon. In field-collected and laboratory-incubated samples, foraminifera expressed denitrification genes regardless of oxygen regime with a putative nitric oxide dismutase, a characteristic enzyme of oxygenic denitrification. A pyruvate:ferredoxin oxidoreductase was highly expressed, indicating the capability for anaerobic energy generation during exposure to hypoxia and anoxia. Near-complete expression of a diatom’s plastid genome in one foraminiferal species suggests kleptoplasty or sequestration of functional plastids, conferring a metabolic advantage despite the host living far below the euphotic zone. Through a unique integration of functions largely unrecognized among “typical” eukaryotes, benthic foraminifera represent winning microeukaryotes in the face of ongoing oceanic deoxygenation.

Description: his project was funded by the U.S. NSF IOS 1557430 and 1557566. H.L.F. acknowledges support from the Swedish Research Council VR (grant number 2017-04190).

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Le Roux, V., Urann, B. M., Brunelli, D., Bonatti, E., Cipriani, A., Demouchy, S., & Monteleone, B. D. Postmelting hydrogen enrichment in the oceanic lithosphere. Science Advances, 7(24), (2021): eabf6071, https://doi.org/10.1126/sciadv.abf6071.

Description: The large range of H2O contents recorded in minerals from exhumed mantle rocks has been challenging to interpret, as it often records a combination of melting, metasomatism, and diffusional processes in spatially isolated samples. Here, we determine the temporal variations of H2O contents in pyroxenes from a 24-Ma time series of abyssal peridotites exposed along the Vema fracture zone (Atlantic Ocean). The H2O contents of pyroxenes correlate with both crustal ages and pyroxene chemistry and increase toward younger and more refractory peridotites. These variations are inconsistent with residual values after melting and opposite to trends often observed in mantle xenoliths. Postmelting hydrogen enrichment occurred by ionic diffusion during cryptic metasomatism of peridotite residues by low-degree, volatile-rich melts and was particularly effective in the most depleted peridotites. The presence of hydrous melts under ridges leads to widespread hydrogen incorporation in the oceanic lithosphere, likely lowering mantle viscosity compared to dry models.

Description: Funding for this study was supported by NSF EAR-P&G 1524311 and 1839128 to V.L.R. and the Andrew W. Mellon Foundation Award for Innovative Research to V.L.R. A.C. and D.B. were funded by the Italian Programma di Rilevante Interesse Nazionale PRIN 20178LPCPW and PRIN2017KY5ZX8, respectively. Revisions were performed within the duration of a “Visiting Scholar at SCIENCE 2020” award to V.L.R. (University of Copenhagen, Denmark), with support from the Department of Geosciences and Natural Resource Management, Section for Geology.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Clemens, S. C., Yamamoto, M., Thirumalai, K., Giosan, L., Richey, J. N., Nilsson-Kerr, K., Rosenthal, Y., Anand, P., & McGrath, S. M. Remote and local drivers of Pleistocene South Asian summer monsoon precipitation: a test for future predictions. Science Advances, 7(23), (2021): eabg3848, https://doi.org/10.1126/sciadv.abg3848.

Description: South Asian precipitation amount and extreme variability are predicted to increase due to thermodynamic effects of increased 21st-century greenhouse gases, accompanied by an increased supply of moisture from the southern hemisphere Indian Ocean. We reconstructed South Asian summer monsoon precipitation and runoff into the Bay of Bengal to assess the extent to which these factors also operated in the Pleistocene, a time of large-scale natural changes in carbon dioxide and ice volume. South Asian precipitation and runoff are strongly coherent with, and lag, atmospheric carbon dioxide changes at Earth’s orbital eccentricity, obliquity, and precession bands and are closely tied to cross-equatorial wind strength at the precession band. We find that the projected monsoon response to ongoing, rapid high-latitude ice melt and rising carbon dioxide levels is fully consistent with dynamics of the past 0.9 million years.

Description: S.C.C. and S.M.M. were supported by U.S. NSF OCE1634774. M.Y. was funded by JSPS grants JPMXS05R2900001 and 19H05595 and JAMSTEC Exp. 353 postcruise study. K.N.-K. and P.A. were supported by UK-IODP, Open University, and NERC (NE/L002493/1), K.T. was supported by the Technology and Research Initiative Fund, Arizona Board of Regents.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Seltzer, A. M., Bekaert, D. V., Barry, P. H., Durkin, K. E., Mace, E. K., Aalseth, C. E., Zappala, J. C., Mueller, P., Jurgens, B., & Kulongoski, J. T. Groundwater residence time estimates obscured by anthropogenic carbonate. Science Advances, 7(17), (2021): eabf3503, https://doi.org/10.1126/sciadv.abf3503.

Description: Groundwater is an important source of drinking and irrigation water. Dating groundwater informs its vulnerability to contamination and aids in calibrating flow models. Here, we report measurements of multiple age tracers (14C, 3H, 39Ar, and 85Kr) and parameters relevant to dissolved inorganic carbon (DIC) from 17 wells in California’s San Joaquin Valley (SJV), an agricultural region that is heavily reliant on groundwater. We find evidence for a major mid-20th century shift in groundwater DIC input from mostly closed- to mostly open-system carbonate dissolution, which we suggest is driven by input of anthropogenic carbonate soil amendments. Crucially, enhanced open-system dissolution, in which DIC equilibrates with soil CO2, fundamentally affects the initial 14C activity of recently recharged groundwater. Conventional 14C dating of deeper SJV groundwater, assuming an open system, substantially overestimates residence time and thereby underestimates susceptibility to modern contamination. Because carbonate soil amendments are ubiquitous, other groundwater-reliant agricultural regions may be similarly affected.

Description: his work was conducted as a part of the USGS National Water Quality Assessment Program (NAWQA) Enhanced Trends Project (https://water.usgs.gov/nawqa/studies/gwtrends/). Measurements at Argonne National Laboratory were supported by Department of Energy, Office of Science under contract DE-AC02-06CH11357. Measurements at Pacific Northwest National Laboratory were part of the Ultra-Sensitive Nuclear Measurements Initiative conducted under the Laboratory Directed Research and Development Program. PNNL is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. This work was also partially supported by NSF award OCE-1923915 (to A.M.S. and P.H.B. at WHOI).

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Zhang, Y., Gazel, E., Gaetani, G. A., & Klein, F. Serpentinite-derived slab fluids control the oxidation state of the subarc mantle. Science Advances, 7(48), (2021): eabj2515, https://doi.org/10.1126/sciadv.abj2515.

Description: Recent geochemical evidence confirms the oxidized nature of arc magmas, but the underlying processes that regulate the redox state of the subarc mantle remain yet to be determined. We established a link between deep subduction-related fluids derived from dehydration of serpentinite ± altered oceanic crust (AOC) using B isotopes and B/Nb as fluid proxies, and the oxidized nature of arc magmas as indicated by Cu enrichment during magma evolution and V/Yb. Our results suggest that arc magmas derived from source regions influenced by a greater serpentinite (±AOC) fluid component record higher oxygen fugacity. The incorporation of this component into the subarc mantle is controlled by the subduction system’s thermodynamic conditions and geometry. Our results suggest that the redox state of the subarc mantle is not homogeneous globally: Primitive arc magmas associated with flat, warm subduction are less oxidized overall than those generated in steep, cold subduction zones.

Description: Y.Z. acknowledges funding from the National Science Foundation of China (91958213), the Chinese Academy of Sciences (XDB42020402), and the Shandong Provincial Natural Science Foundation, China (ZR2020QD068). This study was supported in part by the U.S. National Science Foundation NSF EAR 1826673 to E.G. and G.A.G. and OCE 1756349 to E.G.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Trembath-Reichert, E., Shah Walter, S. R., Ortiz, M. A. F., Carter, P. D., Girguis, P. R., & Huber, J. A. Multiple carbon incorporation strategies support microbial survival in cold subseafloor crustal fluids. Science Advances, 7(18), (2021): eabg0153, https://doi.org/10.1126/sciadv.abg0153.

Description: Biogeochemical processes occurring in fluids that permeate oceanic crust make measurable contributions to the marine carbon cycle, but quantitative assessments of microbial impacts on this vast, subsurface carbon pool are lacking. We provide bulk and single-cell estimates of microbial biomass production from carbon and nitrogen substrates in cool, oxic basement fluids from the western flank of the Mid-Atlantic Ridge. The wide range in carbon and nitrogen incorporation rates indicates a microbial community well poised for dynamic conditions, potentially anabolizing carbon and nitrogen at rates ranging from those observed in subsurface sediments to those found in on-axis hydrothermal vent environments. Bicarbonate incorporation rates were highest where fluids are most isolated from recharging bottom seawater, suggesting that anabolism of inorganic carbon may be a potential strategy for supplementing the ancient and recalcitrant dissolved organic carbon that is prevalent in the globally distributed subseafloor crustal environment.

Description: The Gordon and Betty Moore Foundation sponsored most of the observatory components at North Pond through grant GBMF1609. This work was supported by the National Science Foundation through grants NSF OCE-1745589, OCE-1635208, and OCE-1062006 to J.A.H. and NSF OCE-1635365 to P.R.G. and S.R.S.W.; NASA Postdoctoral Fellowship with the NASA Astrobiology Institute to E.T.-R.; L’Oréal USA For Women in Science Fellowship to E.T.-R.; and Woods Hole Partnership Education Program, sponsored by the Woods Hole Diversity Initiative to M.A.F.O. The Center for Dark Energy Biosphere Investigations (C-DEBI OCE-0939564) also supported the participation of J.A.H. and P.D.C. This is C-DEBI contribution number 564.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Johnson, A. C., Ostrander, C. M., Romaniello, S. J., Reinhard, C. T., Greaney, A. T., Lyons, T. W., & Anbar, A. D. Reconciling evidence of oxidative weathering and atmospheric anoxia on Archean Earth. Science Advances, 7(40), (2021): eabj0108, https://doi.org/10.1126/sciadv.abj0108.

Description: Evidence continues to emerge for the production and low-level accumulation of molecular oxygen (O2) at Earth’s surface before the Great Oxidation Event. Quantifying this early O2 has proven difficult. Here, we use the distribution and isotopic composition of molybdenum in the ancient sedimentary record to quantify Archean Mo cycling, which allows us to calculate lower limits for atmospheric O2 partial pressures (PO2) and O2 production fluxes during the Archean. We consider two end-member scenarios. First, if O2 was evenly distributed throughout the atmosphere, then PO2 〉 10–6.9 present atmospheric level was required for large periods of time during the Archean eon. Alternatively, if O2 accumulation was instead spatially restricted (e.g., occurring only near the sites of O2 production), then O2 production fluxes 〉0.01 Tmol O2/year were required. Archean O2 levels were vanishingly low according to our calculations but substantially above those predicted for an abiotic Earth system.

Description: We would like to thank our funding sources, including FESD “Dynamics of Earth System Oxygenation” (NSF EAR 1338810 to A.D.A.), NASA Earth and Space Science Fellowship awarded to A.C.J. (80NSSC17K0498), NSF EAR PF to A.C.J. (1952809), and WHOI Postdoctoral Fellowship to C.M.O. C.T.R. acknowledges support from the NASA Astrobiology Institute. We also acknowledge support from the Metal Utilization and Selection across Eons (MUSE) Interdisciplinary Consortium for Astrobiology Research, sponsored by the National Aeronautics and Space Administration Science Mission Directorate (19-ICAR19_2-0007).

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bekaert, D. V., Auro, M., Shollenberger, Q. R., Liu, M.-C., Marschall, H., Burton, K. W., Jacobsen, B., Brennecka, G. A., McPherson, G. J., von Mutius, R., Sarafian, A., & Nielsen, S. G. Fossil records of early solar irradiation and cosmolocation of the CAI factory: a reappraisal. Science Advances, 7(40), (2021): eabg8329, https://doi.org/10.1126/sciadv.abg8329.

Description: Calcium-aluminum–rich inclusions (CAIs) in meteorites carry crucial information about the environmental conditions of the nascent Solar System prior to planet formation. Based on models of 50V–10Be co-production by in-situ irradiation, CAIs are considered to have formed within ~0.1 AU from the proto-Sun. Here, we present vanadium (V) and strontium (Sr) isotopic co-variations in fine- and coarse-grained CAIs and demonstrate that kinetic isotope effects during partial condensation and evaporation best explain V isotope anomalies previously attributed to solar particle irradiation. We also report initial excesses of 10Be and argue that CV CAIs possess essentially a homogeneous level of 10Be, inherited during their formation. Based on numerical modeling of 50V–10Be co-production by irradiation, we show that CAI formation during protoplanetary disk build-up likely occurred at greater heliocentric distances than previously considered, up to planet-forming regions (~1AU), where solar particle fluxes were sufficiently low to avoid substantial in-situ irradiation of CAIs.

Description: This study was funded by NASA Emerging Worlds grant NNX16AD36G to S.G.N. and prepared by LLNL under contract DE-AC52-07NA27344 with release number LLNL-JRNL-819045. M.C.L acknowledges the support by the NASA grant 80NSSC20K0759. The UCLA ion microprobe facility is partially supported by a grant from the NSF Instrumentation and Facilities program.

Description: Little research attention has been given to validating clusters obtained from the groundwater geochemistry of the waterworks' capture zone with a prevailing lake‐groundwater exchange. To address this knowledge gap, we proposed a new scheme whereby Gaussian finite mixture modeling (GFMM) and Spike‐and‐Slab Bayesian (SSB) algorithms were utilized to cluster the groundwater geochemistry while quantifying the probability of the resulting cluster membership against each other. We applied GFMM and SSB to 13 geochemical parameters collected during different sampling periods at 13 observation points across the Barnim Highlands plateau located in the northeast of Berlin, Germany; this included 10 observation wells, two lakes, and a gallery of drinking production wells. The cluster analysis of GFMM yielded nine clusters, either with a probability ≥0.8, while the SSB produced three hierarchical clusters with a probability of cluster membership varying from 〈0.2 to 〉0.8. The findings demonstrated that the clustering results of GFMM were in good agreement with the classification as per the principal component analysis and Piper diagram. By superimposing the parameter clustering onto the observation clustering, we could identify discrepancies that exist among the parameters of a certain cluster. This enables the identification of different factors that may control the geochemistry of a certain cluster, although parameters of that cluster share a strong similarity. The GFMM results have shown that from 2002, there has been active groundwater inflow from the lakes towards the capture zone. This means that it is necessary to adopt appropriate measures to reverse the inflow towards the lakes.

Description: Article impact statement: The probability of cluster membership quantified using an algorithm should be validated against another probabilistic‐based classifier.

Description: Federal Ministry of Education and Research http://dx.doi.org/10.13039/501100002347

Description: Copper (Cu) is an essential element for plants and microorganisms and at larger concentrations a toxic pollutant. A number of factors controlling Cu dynamics have been reported, but information on quantitative relationships is scarce. We aimed to (i) quantitatively describe and predict soil Cu concentrations (CuAR) in aqua regia considering site‐specific effects and effects of pH, soil organic carbon (SOC) and cation exchange capacity (CEC), and (ii) study the suitability of mixed‐effects modelling and rule‐based models for the analysis of long‐term soil monitoring data. Thirteen uncontaminated long‐term monitoring soil profiles in southern Germany were analysed. Since there was no measurable trend of increasing CuAR concentrations with time in the respective depth ranges of the sites, data from different sampling dates were combined and horizon‐specific regression analyses including model simplifications were carried out for 10 horizons. Fixed‐ and mixed‐effects models with the site as a random effect were useful for the different horizons and significant contributions (either of main effects or interactions) of SOC, CEC and pH were present for 9, 8 and 7 horizons, respectively. Horizon‐specific rule‐based cubist models described the CuAR data similarly well. Validations of cubist models and mixed‐effects models for the CuAR concentrations in A horizons were successful for the given population after random splitting into calibration and validation samples, but not after independent validations with random splitting according to sites. Overall, site, CEC, SOC and pH provide important information for a description of CuAR concentrations using the different regression approaches. Highlights: Information on quantitative relationships for factors controlling Cu dynamics is scarce. Site, CEC, SOC and pH provide important information for a description of Cu concentrations. Validations of cubist models and mixed‐effects models for A horizons were successful for a closed population of sites.

Description: Bavarian State Ministry of the Environment and Consumer Protection http://dx.doi.org/10.13039/501100010219

Description: Ministry of Agriculture and Environment Mecklenburg‐Western Pomerania

Description: Geological archives record multiple reversals of Earth's magnetic poles, but the global impacts of these events, if any, remain unclear. Uncertain radiocarbon calibration has limited investigation of the potential effects of the last major magnetic inversion, known as the Laschamps Excursion 41 to 42 thousand years ago (ka). We use ancient New Zealand kauri trees (Agathis australis) to develop a detailed record of atmospheric radiocarbon levels across the Laschamps Excursion. We precisely characterize the geomagnetic reversal and perform global chemistry-climate modeling and detailed radiocarbon dating of paleoenvironmental records to investigate impacts. We find that geomagnetic field minima ~42 ka, in combination with Grand Solar Minima, caused substantial changes in atmospheric ozone concentration and circulation, driving synchronous global climate shifts that caused major environmental changes, extinction events, and transformations in the archaeological record.

Description: In recent years, German cities were heavily impacted by pluvial flooding and related damage is projected to increase due to climate change and urbanisation. It is important to ask how to improve urban pluvial flood risk management. To understand the current state of property level adaptation, a survey was conducted in four municipalities that had recently been impacted by pluvial flooding. A hybrid framework based on the Protection Motivation Theory (PMT) and the Protection Action Decision Model (PADM) was used to investigate drivers of adaptive behaviour through both descriptive and regression analyses. Descriptive statistics revealed that participants tended to instal more low‐ and medium‐cost measures than high‐cost measures. Regression analyses showed that coping appraisal increased protection motivation, but that the adaptive behaviour also depends on framing factors, particularly homeownership. We further found that, while threat appraisal solely affects protection motivation and responsibility appraisal affects solely maladaptive thinking, coping appraisal affects both. Our results indicate that PMT is a solid starting point to study adaptive behaviours in the context of pluvial flooding, but we need to go beyond that by, for instance, considering factors of the PADM, such as responsibility, ownership, or respondent age, to fully understand this complex decision‐making process.

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Soil aeration is a critical factor for oxygen‐limited subsoil processes, as transport by diffusion and advection is restricted by the long distance to the free atmosphere. Oxygen transport into the soil matrix is highly dependent on its connectivity to larger pore channels like earthworm and root colonised biopores. Here we hypothesize that the soil matrix around biopores represents different connectivity depending on biopore genesis and actual coloniser. We analysed the soil pore system of undisturbed soil core samples around biopores generated or colonised by roots and earthworms and compared them with the pore system of soil, not in the immediacy of a biopore. Oxygen partial pressure profiles and gas relative diffusion was measured in the rhizosphere and drilosphere from the biopore wall into the bulk soil with microelectrodes. The measurements were linked with structural features such as porosity and connectivity obtained from X‐ray tomography and image analysis. Aeration was enhanced in the soil matrix surrounding biopores in comparison to the bulk soil, shown by higher oxygen concentrations and higher relative diffusion coefficients. Biopores colonised by roots presented more connected lateral pores than earthworm colonised ones, which resulted in enhanced aeration of the rhizosphere compared to the drilosphere. This has influenced biotic processes (microbial turnover/mineralization or root respiration) at biopore interfaces and highlights the importance of microstructural features for soil processes and their dependency on the biopore's coloniser.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Freeman, D. H., & Ward, C. P. Sunlight-driven dissolution is a major fate of oil at sea. Science Advances, 8(7), (2022): eabl7605, https://doi.org/10.1126/sciadv.abl7605.

Description: Oxygenation reactions initiated by sunlight can transform insoluble components of crude oil at sea into water-soluble products, a process called photo-dissolution. First reported a half century ago, photo-dissolution has never been included in spill models because key parameters required for rate modeling were unknown, including the wavelength and photon dose dependence. Here, we experimentally quantified photo-dissolution as a function of wavelength and photon dose, making possible a sensitivity analysis of environmental variables in hypothetical spill scenarios and a mass balance assessment for the 2010 Deepwater Horizon (DwH) spill. The sensitivity analysis revealed that rates were most sensitive to oil slick thickness, season/latitude, and wavelength and less sensitive to photon dose. We estimate that 3 to 17% (best estimate 8%) of DwH surface oil was subject to photo-dissolution, comparable in magnitude to other widely recognized fate processes. Our findings invite a critical reevaluation of surface oil budgets for both DwH and future spills at sea.

Description: This work was supported by the Fisheries and Oceans Canada Multi-Partner Research Initiative award to C.P.W. (project #1.06), the NSF Graduate Research Fellowship awarded to D.H.F. (award #174530), and NSF-OCE grant #1841092 to C.P.W.

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Peng, Q., Xie, S.-P., Wang, D., Huang, R. X., Chen, G., Shu, Y., Shi, J.-R., & Liu, W. Surface warming-induced global acceleration of upper ocean currents. Science Advances, 8(16), (2022): eabj8394, https://doi.org/10.1126/sciadv.abj8394.

Description: How the ocean circulation changes in a warming climate is an important but poorly understood problem. Using a global ocean model, we decompose the problem into distinct responses to changes in sea surface temperature, salinity, and wind. Our results show that the surface warming effect, a robust feature of anthropogenic climate change, dominates and accelerates the upper ocean currents in 77% of the global ocean. Specifically, the increased vertical stratification intensifies the upper subtropical gyres and equatorial currents by shoaling these systems, while the differential warming between the Southern Ocean upwelling zone and the region to the north accelerates surface zonal currents in the Southern Ocean. In comparison, the wind stress and surface salinity changes affect regional current systems. Our study points a way forward for investigating ocean circulation change and evaluating the uncertainty.

Description: Q.P. is supported by the National Natural Science Foundation of China (42005035), the Science and Technology Planning Project of Guangzhou (202102020935), and the Independent Research Project Program of State Key Laboratory of Tropical Oceanography (LTOZZ2102). D.W. is supported by the National Natural Science Foundation of China (92158204), and the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (311020004). S.-P.X. is supported by the National Science Foundation (AGS-1934392). Y.S. is supported by the National Key Research and Development Program of China (2016YFC1401702). G.C. is supported by National Natural Science Foundation of China (41822602). The numerical simulation is supported by the High-Performance Computing Division and HPC managers of W. Zhou and D. Sui in the South China Sea Institute of Oceanology.

Description: Application of farmyard manure (FYM) is common practice to improve physical and chemical properties of arable soil and crop yields. However, studies on effects of FYM application mainly focussed on topsoils, whereas subsoils have rarely been addressed so far. We, therefore, investigated the effects of 36‐year FYM application with different rates of annual organic carbon (OC) addition (0, 469, 938 and 1875 g C m−2 a−1) on OC contents of a Chernozem in 0–30 cm (topsoil) and 35–45 cm (subsoil) depth. We also investigated its effects on soil structure and hydraulic properties in subsoil. X‐ray computed tomography was used to analyse the response of the subsoil macropore system (≥19 μm) and the distribution of particulate organic matter (POM) to different FYM applications, which were related to contents in total OC (TOC) and water‐extractable OC (WEOC). We show that FYM‐C application of 469 g C m−2 a−1 caused increases in TOC and WEOC contents only in the topsoil, whereas rates of ≥938 g C m−2 a−1 were necessary for TOC enrichment also in the subsoil. At this depth, the subdivision of TOC into different OC sources shows that most of the increase was due to fresh POM, likely by the stimulation of root growth and bioturbation. The increase in subsoil TOC went along with increases in macroporosity and macropore connectivity. We neither observed increases in plant‐available water capacity nor in unsaturated hydraulic conductivity. In conclusion, only very high application of FYM over long periods can increase OC content of subsoil at our study site, but this increase is largely based on fresh, easily degradable POM and likely accompanied by high C losses when considering the discrepancy between OC addition rate by FYM and TOC response in soil. Highlights A new image processing procedure to distinguish fresh and decomposed POM. The increase of subsoil C stock based to a large extend on fresh, labile POM. Potential of arable subsoils for long‐term C storage by large FYM application rates is limited. The increase in TOC has no effect on hydraulic properties of the subsoil.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in DiBenedetto, M., Qin, Z., & Suckale, J. Crystal aggregates record the pre-eruptive flow field in the volcanic conduit at Kilauea, Hawaii. Science Advances, 6(49), (2020): eabd4850, doi:10.1126/sciadv.abd4850.

Description: Developing reliable, quantitative conduit models that capture the physical processes governing eruptions is hindered by our inability to observe conduit flow directly. The closest we get to direct evidence is testimony imprinted on individual crystals or bubbles in the conduit and preserved by quenching during the eruption. For example, small crystal aggregates in products of the 1959 eruption of Kīlauea Iki, Hawaii contain overgrown olivines separated by large, hydrodynamically unfavorable angles. The common occurrence of these aggregates calls for a flow mechanism that creates this crystal misorientation. Here, we show that the observed aggregates are the result of exposure to a steady wave field in the conduit through a customized, process-based model at the scale of individual crystals. We use this model to infer quantitative attributes of the flow at the time of aggregate formation; notably, the formation of misoriented aggregates is only reproduced in bidirectional, not unidirectional, conduit flow.

Description: M.D. acknowledges support the Stanford Gerald J. Lieberman Fellowship and the Postdoctoral Scholarship from Woods Hole Oceanographic Institution.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Praetorius, S. K., Condron, A., Mix, A. C., Walczak, M. H., McKay, J. L., & Du, J. The role of northeast pacific meltwater events in deglacial climate change. Science Advances, 6(9), (2020): eaay2915, doi:10.1126/sciadv.aay2915.

Description: Columbia River megafloods occurred repeatedly during the last deglaciation, but the impacts of this fresh water on Pacific hydrography are largely unknown. To reconstruct changes in ocean circulation during this period, we used a numerical model to simulate the flow trajectory of Columbia River megafloods and compiled records of sea surface temperature, paleo-salinity, and deep-water radiocarbon from marine sediment cores in the Northeast Pacific. The North Pacific sea surface cooled and freshened during the early deglacial (19.0-16.5 ka) and Younger Dryas (12.9-11.7 ka) intervals, coincident with the appearance of subsurface water masses depleted in radiocarbon relative to the sea surface. We infer that Pacific meltwater fluxes contributed to net Northern Hemisphere cooling prior to North Atlantic Heinrich Events, and again during the Younger Dryas stadial. Abrupt warming in the Northeast Pacific similarly contributed to hemispheric warming during the Bølling and Holocene transitions. These findings underscore the importance of changes in North Pacific freshwater fluxes and circulation in deglacial climate events.

Description: The research was partly supported by the NSF through grants ARC-257 1204045 and PLR-1417667. The numerical model simulations used resources from the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under contract no. DE-AC02-05CH11231.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Torres, J. P., Lin, Z., Watkins, M., Salcedo, P. F., Baskin, R. P., Elhabian, S., Safavi-Hemami, H., Taylor, D., Tun, J., Concepcion, G. P., Saguil, N., Yanagihara, A. A., Fang, Y., McArthur, J. R., Tae, H. S., Finol-Urdaneta, R. K., Özpolat, B. D., Olivera, B. M., & Schmidt, E. W. Small-molecule mimicry hunting strategy in the imperial cone snail, Conus imperialis. Science Advances, 7(11), (2021): eabf2704, https://doi.org/10.1126/sciadv.abf2704.

Description: Venomous animals hunt using bioactive peptides, but relatively little is known about venom small molecules and the resulting complex hunting behaviors. Here, we explored the specialized metabolites from the venom of the worm-hunting cone snail, Conus imperialis. Using the model polychaete worm Platynereis dumerilii, we demonstrate that C. imperialis venom contains small molecules that mimic natural polychaete mating pheromones, evoking the mating phenotype in worms. The specialized metabolites from different cone snails are species-specific and structurally diverse, suggesting that the cones may adopt many different prey-hunting strategies enabled by small molecules. Predators sometimes attract prey using the prey’s own pheromones, in a strategy known as aggressive mimicry. Instead, C. imperialis uses metabolically stable mimics of those pheromones, indicating that, in biological mimicry, even the molecules themselves may be disguised, providing a twist on fake news in chemical ecology.

Description: Research reported in this publication was supported by NIH R35GM12252, with contributions to biological work from NIH Fogarty International Center U19TW008163, NIH P01GM48677, and DOD CDMRP W81XWH-17-1-0413. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Guillermic, M., Cameron, L. P., De Corte, I., Misra, S., Bijma, J., de Beer, D., Reymond, C. E., Westphal, H., Ries, J. B., & Eagle, R. A. Thermal stress reduces pocilloporid coral resilience to ocean acidification by impairing control over calcifying fluid chemistry. Science Advances, 7(2), (2021): eaba9958, https://doi.org/10.1126/sciadv.aba9958.

Description: The combination of thermal stress and ocean acidification (OA) can more negatively affect coral calcification than an individual stressors, but the mechanism behind this interaction is unknown. We used two independent methods (microelectrode and boron geochemistry) to measure calcifying fluid pH (pHcf) and carbonate chemistry of the corals Pocillopora damicornis and Stylophora pistillata grown under various temperature and pCO2 conditions. Although these approaches demonstrate that they record pHcf over different time scales, they reveal that both species can cope with OA under optimal temperatures (28°C) by elevating pHcf and aragonite saturation state (Ωcf) in support of calcification. At 31°C, neither species elevated these parameters as they did at 28°C and, likewise, could not maintain substantially positive calcification rates under any pH treatment. These results reveal a previously uncharacterized influence of temperature on coral pHcf regulation—the apparent mechanism behind the negative interaction between thermal stress and OA on coral calcification.

Description: R.A.E. and J.B.R. acknowledge support from National Science Foundation grants OCE-1437166 and OCE-1437371. The work was also supported by the “Laboratoire d’Excellence” LabexMER (ANR-10-LABX-19), cofunded by a grant from the French government under the program “Investissements d’Avenir,” and an IAGC student grant 2017. R.A.E. acknowledges financial and logistical support from the Pritzker Endowment to UCLA IoES, and J.B.R. acknowledges support from the ZMT and the Hanse-Wissenschaftskolleg Fellowship Program and the NSF OCE award #1437371.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Law, K. L., Starr, N., Siegler, T. R., Jambeck, J. R., Mallos, N. J., & Leonard, G. H. The United States' contribution of plastic waste to land and ocean. Science Advances, 6(44), (2020): eabd0288, doi:10.1126/sciadv.abd0288.

Description: Plastic waste affects environmental quality and ecosystem health. In 2010, an estimated 5 to 13 million metric tons (Mt) of plastic waste entered the ocean from both developing countries with insufficient solid waste infrastructure and high-income countries with very high waste generation. We demonstrate that, in 2016, the United States generated the largest amount of plastic waste of any country in the world (42.0 Mt). Between 0.14 and 0.41 Mt of this waste was illegally dumped in the United States, and 0.15 to 0.99 Mt was inadequately managed in countries that imported materials collected in the United States for recycling. Accounting for these contributions, the amount of plastic waste generated in the United States estimated to enter the coastal environment in 2016 was up to five times larger than that estimated for 2010, rendering the United States’ contribution among the highest in the world.

Description: This work was funded by Ocean Conservancy through support from the Arthur Vining Davis Foundations.

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Skinner, C., Mill, A. C., Fox, M. D., Newman, S. P., Zhu, Y., Kuhl, A., & Polunin, N. V. C. Offshore pelagic subsidies dominate carbon inputs to coral reef predators. Science Advances, 7(8), (2021): eabf3792, https://doi.org/10.1126/sciadv.abf3792.

Description: Coral reefs were traditionally perceived as productive hot spots in oligotrophic waters. While modern evidence indicates that many coral reef food webs are heavily subsidized by planktonic production, the pathways through which this occurs remain unresolved. We used the analytical power of carbon isotope analysis of essential amino acids to distinguish between alternative carbon pathways supporting four key reef predators across an oceanic atoll. This technique separates benthic versus planktonic inputs, further identifying two distinct planktonic pathways (nearshore reef-associated plankton and offshore pelagic plankton), and revealing that these reef predators are overwhelmingly sustained by offshore pelagic sources rather than by reef sources (including reef-associated plankton). Notably, pelagic reliance did not vary between species or reef habitats, emphasizing that allochthonous energetic subsidies may have system-wide importance. These results help explain how coral reefs maintain exceptional productivity in apparently nutrient-poor tropical settings, but also emphasize their susceptibility to future ocean productivity fluctuations.

Description: Sample analysis funding was provided by NERC LSMSF grant BRIS/102/0717 and BRIS/125/1418. C.S. was supported by a Newcastle University SAgE DTA studentship and a cooperative agreement with Banyan Tree.

Description: © The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Richter, M., Nebel, O., Maas, R., Mather, B., Nebel-Jacobsen, Y., Capitanio, F. A., Dick, H. J. B., & Cawood, P. A. An early cretaceous subduction-modified mantle underneath the ultraslow spreading Gakkel Ridge, Arctic Ocean. Science Advances, 6(44), (2020): eabb4340, doi:10.1126/sciadv.abb4340.

Description: Earth’s upper mantle, as sampled by mid-ocean ridge basalts (MORBs) at oceanic spreading centers, has developed chemical and isotopic heterogeneity over billions of years through focused melt extraction and re-enrichment by recycled crustal components. Chemical and isotopic heterogeneity of MORB is dwarfed by the large compositional spectrum of lavas at convergent margins, identifying subduction zones as the major site for crustal recycling into and modification of the mantle. The fate of subduction-modified mantle and if this heterogeneity transmits into MORB chemistry remains elusive. Here, we investigate the origin of upper mantle chemical heterogeneity underneath the Western Gakkel Ridge region in the Arctic Ocean through MORB geochemistry and tectonic plate reconstruction. We find that seafloor lavas from the Western Gakkel Ridge region mirror geochemical signatures of an Early Cretaceous, paleo-subduction zone, and conclude that the upper mantle can preserve a long-lived, stationary geochemical memory of past geodynamic processes.

Description: O.N. was supported by the Australian Research Council (grant FT140101062). P.A.C. was supported by the Australian Research Council (grant FL160100168). H.J.B.D. was supported by the NSF (grants PLR 9912162, PLR 0327591, OCE 0930487, and OCE 1434452). M.R. was supported by a graduate scholarship of Monash University and the SEAE.