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: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 125(6), (2020): e2019JB019239, doi:10.1029/2019JB019239.

Description: P‐to‐S‐converted waves observed in controlled‐source multicomponent ocean bottom seismometer (OBS) records were used to derive the Vp/Vs structure of Cascadia Basin sediments. We used P‐to‐S waves converted at the basement to derive an empirical function describing the average Vp/Vs of Cascadia sediments as a function of sediment thickness. We derived one‐dimensional interval Vp/Vs functions from semblance velocity analysis of S‐converted intrasediment and basement reflections, which we used to define an empirical Vp/Vs versus burial depth compaction trend. We find that seaward from the Cascadia deformation front, Vp/Vs structure offshore northern Oregon and Washington shows little variability along strike, while the structure of incoming sediments offshore central Oregon is more heterogeneous and includes intermediate‐to‐deep sediment layers of anomalously elevated Vp/Vs. These zones with elevated Vp/Vs are likely due to elevated pore fluid pressures, although layers of high sand content intercalated within a more clayey sedimentary sequence, and/or a higher content of coarser‐grained clay minerals relative to finer‐grained smectite could be contributing factors. We find that the proto‐décollement offshore central Oregon develops within the incoming sediments at a low‐permeability boundary that traps fluids in a stratigraphic level where fluid overpressure exceeds 50% of the differential pressure between the hydrostatic pressure and the lithostatic pressure. Incoming sediments with the highest estimated fluid overpressures occur offshore central Oregon where deformation of the accretionary prism is seaward vergent. Conversely, landward vergence offshore northern Oregon and Washington correlates with more moderate pore pressures and laterally homogeneous Vp/Vs functions of Cascadia Basin sediments.

Description: This research was funded by National Science Foundation (NSF) Grant OCE‐1657237 to J. P. C, OCE‐1657839 to A. F. A. and S. H., and OCE‐1657737 to S. M. C. Data used in this study were acquired with funding from NSF Grants OCE‐1029305 and OCE‐1249353. Data used in this research were provided by instruments from the Ocean Bottom Seismic Instrument Center (http://obsic.whoi.edu, formerly OBSIP), which is funded by the NSF. OBSIC/OBSIP data are archived at the IRIS Data Management Center (http://www.iris.edu) under network code X6 (https://doi.org/10.7914/SN/X6_2012). Data processing was conducted with Emerson‐Paradigm Software package Echos licensed to Woods Hole Oceanographic Institution under Paradigm Academic Software Program and MATLAB package SeismicLab of the University of Alberta, Canada (http://seismic-lab.physics.ualberta.ca), under GNU General Public License (MATLAB® is a registered trademark of MathWorks).

Description: 2020-11-28

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(8), (2020): e2020JC016068, doi:10.1029/2020JC016068.

Description: Labrador Sea Water (LSW) is a major component of the deep limb of the Atlantic Meridional Overturning Circulation, yet LSW transport pathways and their variability lack a complete description. A portion of the LSW exported from the subpolar gyre is advected eastward along the North Atlantic Current and must contend with the Mid‐Atlantic Ridge before reaching the eastern basins of the North Atlantic. Here, we analyze observations from a mooring array and satellite altimetry, together with outputs from a hindcast ocean model simulation, to estimate the mean transport of LSW across the Charlie‐Gibbs Fracture Zone (CGFZ), a primary gateway for the eastward transport of the water mass. The LSW transport estimated from the 25‐year altimetry record is 5.3 ± 2.9 Sv, where the error represents the combination of observational variability and the uncertainty in the projection of the surface velocities to the LSW layer. Current velocities modulate the interannual to higher‐frequency variability of the LSW transport at the CGFZ, while the LSW thickness becomes important on longer time scales. The modeled mean LSW transport for 1993–2012 is higher than the estimate from altimetry, at 8.2 ± 4.1 Sv. The modeled LSW thickness decreases substantially at the CGFZ between 1996 and 2009, consistent with an observed decline in LSW volume in the Labrador Sea after 1994. We suggest that satellite altimetry and continuous hydrographic measurements in the central Labrador Sea, supplemented by profiles from Argo floats, could be sufficient to quantify the LSW transport at the CGFZ.

Description: A. G. N. appreciates conversations with Kathy Donohue, Tom Rossby and Lisa Beal, which helped to interpret the results. J. B. P. acknowledges support from NSF through Grant OCE‐1947829. The authors thank all colleagues and ship crew involved in the R/V Meteor cruise M‐82/2 and Maria S. Merian cruise MSM‐21/2. The mooring data presented in this paper were funded by NSF through Grant OCE‐0926656.

Description: 2021-01-03

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geophysical Research Letters 47(3), (2020): e2019GL086703, doi:10.1029/2019GL086703.

Description: Salt marsh assessments focus on vertical metrics such as accretion or lateral metrics such as open‐water conversion, without exploration of how the dimensions are related. We exploited a novel geospatial data set to explore how elevation is related to the unvegetated‐vegetated marsh ratio (UVVR), a lateral metric, across individual marsh “units” within four estuarine‐marsh systems. We find that elevation scales consistently with the UVVR across systems, with lower elevation units demonstrating more open‐water conversion and higher UVVRs. A normalized elevation‐UVVR relationship converges across systems near the system‐mean elevation and a UVVR of 0.1, a critical threshold identified by prior studies. This indicates that open‐water conversion becomes a dominant lateral instability process at a relatively conservative elevation threshold. We then integrate the UVVR and elevation to yield lifespan estimates, which demonstrate that higher elevation marshes are more resilient to internal deterioration, with an order‐of‐magnitude longer lifespan than predicted for lower elevation marshes.

Description: This study was supported by the USGS through the Coastal Marine Hazards/Resources Program, the National Park Service through the Natural Resource Preservation Program, and the U.S. Fish and Wildlife Service through the Science Support Partnership. Erika Lentz, Elizabeth Pendleton, Meagan Gonneea, Joel Carr, and two anonymous reviewers provided constructive advice on the study. S.F. was partly supported by US National Science Foundation award 1637630 (PIE LTER), 1832221 (VCR LTER). The geospatial data used in this study are published in the Coastal Wetlands Synthesis Products catalog on ScienceBase (https://www.sciencebase.gov/catalog/item/5b73325ee4b0f5d5787c5ff3).

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research- Biogeosciences 125(4), (2020): e2019JG005158, doi:10.1029/2019JG005158.

Description: Long‐term soil warming can decrease soil organic matter (SOM), resulting in self‐reinforcing feedback to the global climate system. We investigated additional consequences of SOM reduction for soil water holding capacity (WHC) and soil thermal and hydrological buffering. At a long‐term soil warming experiment in a temperate forest in the northeastern United States, we suspended the warming treatment for 104 days during the summer of 2017. The formerly heated plot remained warmer (+0.39 °C) and drier (−0.024 cm3 H2O cm−3 soil) than the control plot throughout the suspension. We measured decreased SOM content (−0.184 g SOM g−1 for O horizon soil, −0.010 g SOM g−1 for A horizon soil) and WHC (−0.82 g H2O g−1 for O horizon soil, −0.18 g H2O g−1 for A horizon soil) in the formerly heated plot relative to the control plot. Reduced SOM content accounted for 62% of the WHC reduction in the O horizon and 22% in the A horizon. We investigated differences in SOM composition as a possible explanation for the remaining reductions with Fourier transform infrared (FTIR) spectra. We found FTIR spectra that correlated more strongly with WHC than SOM, but those particular spectra did not differ between the heated and control plots, suggesting that SOM composition affects WHC but does not explain treatment differences in this study. We conclude that SOM reductions due to soil warming can reduce WHC and hydrological and thermal buffering, further warming soil and decreasing SOM. This feedback may operate in parallel, and perhaps synergistically, with carbon cycle feedbacks to climate change.

Description: We would like to acknowledge Jeffery Blanchard, Priya Chowdhury, Kristen DeAngelis, Luiz Dominguez‐Horta, Kevin Geyer, Rachelle Lacroix, Xaiojun Liu, William Rodriguez, and Alexander Truchonand and for assistance with field sampling. We would like to acknowledge Michael Bernard for assistance with field sampling and lab work. We would like to acknowledge Aaron Ellison for statistical consultation. This research was financially supported by the U.S. National Science Foundation's Long Term Ecological Research Program (NSF‐DEB‐0620443 and NSF‐DEB‐1237491), the Long Term Research in Environmental Biology Program (NSF DEB‐1456528) , and the U.S. Department of Energy (DOE‐DE‐SC0005421 and DOE‐DE‐SC0010740). Data used in this study are available from the Harvard Forest Data Archive (Datasets HF018‐03, HF018‐04, and HF018‐13), accessible at https://harvardforest.fas.harvard.edu/harvard‐forest‐data‐archive.

Description: 2020-10-04

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(8), (2020): e2020JC016445, doi:10.1029/2020JC016445.

Description: The Mid‐Atlantic Bight (MAB) Cold Pool is a bottom‐trapped, cold (temperature below 10°C) and fresh (practical salinity below 34) water mass that is isolated from the surface by the seasonal thermocline and is located over the midshelf and outer shelf of the MAB. The interannual variability of the Cold Pool with regard to its persistence time, volume, temperature, and seasonal along‐shelf propagation is investigated based on a long‐term (1958–2007) high‐resolution regional model of the northwest Atlantic Ocean. A Cold Pool Index is defined and computed in order to quantify the strength of the Cold Pool on the interannual timescale. Anomalous strong, weak, and normal years are categorized and compared based on the Cold Pool Index. A detailed quantitative study of the volume‐averaged heat budget of the Cold Pool region (CPR) has been examined on the interannual timescale. Results suggest that the initial temperature and abnormal warming/cooling due to advection are the primary drivers in the interannual variability of the near‐bottom CPR temperature anomaly during stratified seasons. The long persistence of temperature anomalies from winter to summer in the CPR also suggests a potential for seasonal predictability.

Description: This work was funded by the National Oceanic and Atmospheric Administration through Awards NOAA‐NA‐15OAR4310133 and NOAA‐NA‐13OAR4830233 and the National Science Foundation Awards OCE‐1049088, OCE‐1419584, and OCE‐0961545.

Description: 2021-02-03

Description: Author Posting. © American Geophysical Union, 2020. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 125(8), (2020): e2020JC016197, doi:10.1029/2020JC016197.

Description: Synoptic shipboard measurements, together with historical hydrographic data and satellite data, are used to elucidate the detailed structure of the Atlantic Water (AW) boundary current system in the southern Canada Basin and its connection to the upstream source of AW in the Chukchi Borderland. Nine high‐resolution occupations of a transect extending from the Beaufort shelf to the deep basin near 152°W, taken between 2003 and 2018, reveal that there are two branches of the AW boundary current that flow beneath and counter to the Beaufort Gyre. Each branch corresponds to a warm temperature core and transports comparable amounts of Fram Strait Branch Water between roughly 200–700 m depth, although they are characterized by a different temperature/salinity (T/S) structure. The mean volume flux of the combined branches is 0.87 ± 0.13 Sv. Using the historical hydrographic data, the two branches are tracked upstream by their temperature cores and T/S signatures. This sheds new light on how the AW negotiates the Chukchi Borderland and why two branches emerge from this region. Lastly, the propagation of warm temperature anomalies through the region is quantified and shown to be consistent with the deduced circulation scheme.

Description: This work was funded by the following sources: National Science Foundation Grants PLR‐1504333, OPP‐1733564, and OPP‐1504394; National Oceanic and Atmospheric Administration Grant NA14OAR4320158; and National Aeronautics and Space Administration Grant NNX10AF42G.

Description: 2021-01-27

Description: Vegetation biomass is a globally important climate-relevant terrestrial carbon pool. Landsat, Sentinel-2 and Sentinel-1 satellite missions provide a landscape-level opportunity to upscale tundra vegetation communities and biomass in high latitude terrestrial environments. We assessed the applicability of landscape-level remote sensing for the low Arctic Lena Delta region in Northern Yakutia, Siberia, Russia. The Lena Delta is the largest delta in the Arctic and is located North of the treeline and the 10 °C July isotherm at 72° Northern Latitude in the Laptev Sea region. During the LENA2018 expedition, we set up plots for plant projective cover and Above Ground Biomass (AGB) and sampled shrubs for shrub-ring analyses. AGB is providing the magnitude of the carbon flux, whereas stand age is irreplaceable to provide the cycle rate. AGB data and shrub age data clearly show a separation between i) low disturbance landscape types with dominant AGB moss contribution, but always low vascular plant AGB (〈0.5 kg m-2) characterised by old shrubs of several decades of stand age versus ii) a much higher vascular plant AGB contribution (〉 0.5 kg m-2) with only young shrubs in high disturbance regimes. The low disturbance regimes are represented on the Holocene and Pleistocene delta terraces in form of azonal polygonal tundra complexes and softly dissected valleys with zonal tussock tundra. In contrast, the high disturbance regimes are sites of thermo-erosion such as along thermo-erosional valleys and on floodplains. We upscaled AGB and above ground carbon pool ages using a Sentinel-2 satellite acquisition from early August 2018. We classified via classification training using Elementary Sampling Units that are the 30 m x 30 m vegetation field plots. We then used the land cover classes and grouped them according to their settings either in high disturbance or low disturbance regimes with each associated AGB value ranges and shrub age regimes. We also evaluated circum-Arctic harmonized ESA GlobPermafrost land cover and vegetation height remote sensing products covering subarctic to Arctic land cover types for the central Lena Delta. The products are freely available and published in the PANGAEA data repository under https://doi.org/10.1594/PANGAEA.897916 and https://doi.org/10.1594/PANGAEA.897045. ESA GlobPermafrost land cover and vegetation height remote sensing products and our Sentinel-2 derived AGB product for the central Lena Delta shows realistic spatial patterns of landcover classes and biomass distribution at landscape level. However, in all products, the high biomass patches of high shrubs in the tundra landscape could not spatially be resolved as they are confined to patchy and linear distribution, not representing large enough areas suitable for upscaling. We found that high disturbance regimes with linked high and rapid AGB fluxes are distributed mainly on the floodplains and as patches along thermoerosioal features, e.g. valleys. Whereas the low disturbance landscapes on Yedoma upland tundra and Holocene terraces occur with larger area coverage representing decades slower and in magnitude smaller AGB fluxes.