Predicted present-day evolution patterns of ice thickness and bedrock elevation over Greenland and Antarctica

This paper discusses predicted evolution patterns of present-day changes of ice thickness, surface elevation, and bedrock elevation over the Greenland and Antarctic continents. These were obtained from calculations with dynamic 3-D ice sheet models which were coupled to a visco-elastic solid Earth model. The experiments were initialized over the last two glacial cycles and subsequently averaged over the last 200 years to obtain the current evolution. The calculations indicate that the Antarctic Ice Sheet is still adjusting to the last glacial-interglacial transition yielding a decreasing ice volume and a rising bedrock elevation of the order of several centimetres per year. The Greenland Ice Sheet was found to be close to a stationary state with a mean thickness change of only a few millimetres per year, but the calculations revealed large spatial differences. Predicted patterns over Greenland are characterized by a small thickening over the ice sheet interior and a general thinning of the ablation area. In Antarctica, almost all of the predicted changes are concentrated in the West Antarctic Ice Sheet, which is still retreating at both the Weddell and Ross Sea margins. Over most of both ice sheets, the model indicates that the surface elevation trend is dominated by ice thickness changes rather than by bedrock elevation changes.     

Chemistry of rainwater in the Massif Central (France): a strontium isotope and major element study

Atmospheric aerosols (sea salt, crustal dust, and biogenic aerosols) are the primary source of dissolved species in rainwater as well as one of the sources of dissolved species in river water. Chemical weathering studies require quantification of this atmospheric input. The crustal component of atmospheric input can have various origins, both distant and local. The proportions of the various inputs (marine, distant or local) are determined in this study.Strontium isotope ratios and Ca, Na, K, Mg, Al, Cl, SO₄, NO₃ and Sr concentrations were measured in rainwater samples collected in the Massif Central (France) over a period of one year. Each sample, collected automatically, represents a monthly series of rain events. Chemical composition of the rainwater samples varied considerably and the ⁸⁷Sr/⁸⁶Sr ratios ranged between 0.709198 and 0.713143.Using Na as an indicator of marine origin, and Al for the crustal input in rain samples, the proportion of marine and crustal elements was estimated from elemental ratios. A marine origin of 4 to 100% of Cl, of 0.6 to 20% of the SO₄, of <1 to 10% of Ca, <1 to 40% of K, 4 to 100% of Mg and 1 to 44% of Sr was determined.Strontium isotopes were used to characterize the crustal sources. The ⁸⁷Sr/⁸⁶Sr ratios of the crustal sources varied considerably from 0.7092 to 0.71625 and indicate the occurrence of multiple sources for the crustal component in the analysed rainwaters.     

Robust dual‐sensor data sea‐surface ghost attenuation for quality control purposes

Three‐dimensional receiver ghost attenuation (deghosting) of dual‐sensor towed‐streamer data is straightforward, in principle. In its simplest form, it requires applying a three‐dimensional frequency–wavenumber filter to the vertical component of the particle motion data to correct for the amplitude reduction on the vertical component of non‐normal incidence plane waves before combining with the pressure data. More elaborate techniques use three‐dimensional filters to both components before summation, for example, for ghost wavelet dephasing and mitigation of noise of different strengths on the individual components in optimum deghosting. The problem with all these techniques is, of course, that it is usually impossible to transform the data into the crossline wavenumber domain because of aliasing. Hence, usually, a two‐dimensional version of deghosting is applied to the data in the frequency–inline wavenumber domain. We investigate going down the “dimensionality ladder” one more step to a one‐dimensional weighted summation of the records of the collocated sensors to create an approximate deghosting procedure. We specifically consider amplitude‐balancing weights computed via a standard automatic gain control before summation, reminiscent of a diversity stack of the dual‐sensor recordings. This technique is independent of the actual streamer depth and insensitive to variations in the sea‐surface reflection coefficient. The automatic gain control weights serve two purposes: (i) to approximately correct for the geometric amplitude loss of the Z data and (ii) to mitigate noise strength variations on the two components. Here, Z denotes the vertical component of the velocity of particle motion scaled by the seismic impedance of the near‐sensor water volume. The weights are time‐varying and can also be made frequency‐band dependent, adapting better to frequency variations of the noise. The investigated process is a very robust, almost fully hands‐off, approximate three‐dimensional deghosting step for dual‐sensor data, requiring no spatial filtering and no explicit estimates of noise power. We argue that this technique performs well in terms of ghost attenuation (albeit, not exact ghost removal) and balancing the signal‐to‐noise ratio in the output data. For instances where full three‐dimensional receiver deghosting is the final product, the proposed technique is appropriate for efficient quality control of the data acquired and in aiding the parameterisation of the subsequent deghosting processing.     

Simulating rainfall time-series: how to account for statistical variability at multiple scales?

Daily rainfall is a complex signal exhibiting alternation of dry and wet states, seasonal fluctuations and an irregular behavior at multiple scales that cannot be preserved by stationary stochastic simulation models. In this paper, we try to investigate some of the strategies devoted to preserve these features by comparing two recent algorithms for stochastic rainfall simulation: the first one is the modified Markov model, belonging to the family of Markov-chain based techniques, which introduces non-stationarity in the chain parameters to preserve the long-term behavior of rainfall. The second technique is direct sampling, based on multiple-point statistics, which aims at simulating a complex statistical structure by reproducing the same data patterns found in a training data set. The two techniques are compared by first simulating a synthetic daily rainfall time-series showing a highly irregular alternation of two regimes and then a real rainfall data set. This comparison allows analyzing the efficiency of different elements characterizing the two techniques, such as the application of a variable time dependence, the adaptive kernel smoothing or the use of low-frequency rainfall covariates. The results suggest, under different data availability scenarios, which of these elements are more appropriate to represent the rainfall amount probability distribution at different scales, the annual seasonality, the dry-wet temporal pattern, and the persistence of the rainfall events.     

First Multi‐Isotopic (Pb‐Nd‐Sr‐Zn‐Cu‐Fe) Characterisation of Dust Reference Materials (ATD and BCR‐723): A Multi‐Column Chromatographic Method Optimised to Trace Mineral and Anthropogenic Dust Sources

Atmospheric dust is an integral component of the Earth system with major implications for the climate, biosphere and public health. In this context, identifying and quantifying the provenance and the processes generating the various types of dust found in the atmosphere is paramount. Isotopic signatures of Pb, Nd, Sr, Zn, Cu and Fe are commonly used as sensitive geochemical tracers. However, their combined use is limited by the lack of (a) a dedicated chromatographic protocol to separate the six elements of interest for low‐mass samples and (b) specific reference materials for dust. Indeed, our work shows that USGS rock reference materials BHVO‐2, AGV‐2 and G‐2 are not applicable as substitute reference materials for dust. We characterised the isotopic signatures of these six elements in dust reference materials ATD and BCR‐723, representatives of natural and urban environments, respectively. To achieve this, we developed a specific procedure for dust, applicable in the 4–25 mg mass range, to separate the six elements using a multi‐column ion‐exchange chromatographic method and MC‐ICP‐MS measurements.     

Towards Improved Remedial Outcomes in Categorical Aquifers with an Iterative Ensemble Smoother

Categorical parameter distributions consisting of geologic facies with distinct properties, for example, high-permeability channels embedded in a low-permeability matrix, are common at contaminated sites. At these sites, low-permeability facies store solute mass, acting as secondary sources to higher-permeability facies, sustaining concentrations for decades while increasing risk and cleanup costs. Parameter estimation is difficult in such systems because the discontinuities in the parameter space hinder the inverse problem. This paper presents a novel approach based on Traveling Pilot Points (TRIPS) and an iterative ensemble smoother (IES) to solve the categorical inverse problem. Groundwater flow and solute transport in a hypothetical aquifer with a categorical parameter distribution are simulated using MODFLOW 6. Heads and concentrations are recorded at multiple monitoring locations. IES is used to generate posterior ensembles assuming a TRIPS prior and an approximate multi-Gaussian prior. The ensembles are used to predict solute concentrations and mass into the future. The evaluation also includes an assessment of how the number of measurements and the choice of the geological prior determine the characteristics of the posterior ensemble and the resulting predictions. The results indicate that IES was able to efficiently sample the posterior distribution and showed that even with an approximate geological prior, a high degree of parameterization and history matching could lead to parameter ensembles that can be useful for making certain types of predictions (heads, concentrations). However, the approximate geological prior was insufficient for predicting mass. The analysis demonstrates how decision-makers can quantify uncertainty and make informed decisions with an ensemble-based approach.     

Sediment flux-driven channel geometry adjustment of bedrock and mixed gravel–bedrock rivers

Sediment supply (Q_s) is often overlooked in modelling studies of landscape evolution, despite sediment playing a key role in the physical processes that drive erosion and sedimentation in river channels. Here, we show the direct impact of the supply of coarse-grained, hard sediment on the geometry of bedrock channels from the Rangitikei River, New Zealand. Channels receiving a coarse bedload sediment supply are systematically (up to an order of magnitude) wider than channels with no bedload sediment input for a given discharge. We also present physical model experiments of a bedrock river channel with a fixed water discharge (1.5 l min⁻¹) under different Q_s (between 0 and 20 g l⁻¹) that allow the quantification of the role of sediment in setting the width and slope of channels and the distribution of shear stress within channels. The addition of bedload sediment increases the width, slope and width-to-depth ratio of the channels, and increasing sediment loads promote emerging complexity in channel morphology and shear stress distributions. Channels with low Q_s are characterized by simple in-channel morphologies with a uniform distribution of shear stress within the channel while channels with high Q_s are characterized by dynamic channels with multiple active threads and a non-uniform distribution of shear stress. We compare bedrock channel geometries from the Rangitikei and the experiments to alluvial channels and demonstrate that the behaviour is similar, with a transition from single-thread and uniform channels to multiple threads occurring when bedload sediment is present. In the experimental bedrock channels, this threshold Q_s is when the input sediment supply exceeds the transport capacity of the channel. Caution is required when using the channel geometry to reconstruct past environmental conditions or to invert for tectonic uplift rates, because multiple configurations of channel geometry can exist for a given discharge, solely due to input Q_s. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Multiple-point statistics using multi-resolution images

Stochastic Environmental Research and Risk Assessment - Multiple-point statistics (MPS) is a simulation technique allowing to generate images that reproduce the spatial features present in a...     

Coleopteran evidence for a mosaic of environments at high altitude in the eastern Pyrénées,France, during the climatic transition between the Allerød and Younger Dryas

Late‐glacial environmental and climatic implications are inferred from an insect fauna from organic sediments infilling a palaeochannel on the banks of the River Têt, eastern Pyrénées, France. A pine cone in association with the insect fauna has been radiocarbon dated to 10 920 ± 60 yr BP, namely close to the Allerød – Younger Dryas boundary. Two distinct insect associations appear to be recognisable here. One is an assemblage typical of the high altitude forest and a second is characteristic of an alpine grassland. The close coexistence of these two assemblages is attributed to the climatic cooling towards the start of the Younger Dryas Stadial, when the forest cover broke up into remnant patches interspersed by alpine grassland. It is suggested that in a region of such high relief a mosaic of habitats may have been caused by patchy differences in insolation aspect, especially during a period of climatic deterioration. Copyright © 1999 John Wiley & Sons, Ltd.     

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.