Early Cretaceous (late Berriasian to early Aptian) palaeoceanographic change along the northwestern Tethyan margin (Vocontian Trough,southeastern France): δ13C, δ18O and Sr-isotope belemnite and whole-rock records

Stable carbon, oxygen, and strontium isotope records were obtained from uppermost Hauterivian to lowermost Aptian belemnite rostra, which were collected in well-dated sections from the Vocontian Trough (southeastern France). This data set complements previously published belemnite-isotope records from the uppermost Berriasian-Hauterivian interval from the same basin. The belemnite carbon and oxygen isotope record is compared to the carbonate bulk-rock isotope record from the same sections, and from additional Italian sections. With regards to their long-term trends, both belemnite and whole-rock δ¹⁸O records are well correlated, except for the uppermost Hauterivian-lower Barremian interval, within which they deviate. This discrepancy is interpreted to be linked to the latest Hauterivian Faraoni oceanic anoxic event and its early Barremian aftermath. The Faraoni level is characterized by enhanced sea-water stratification, probably induced by the onset of a warmer and more humid climate along the northern Tethyan margin. The early Barremian was characterized by stronger vertical sea-water mixing reflected by a decrease in density contrast between sea-surface and deeper waters. The belemnite oxygen isotope record shows a more stable evolution with smaller fluctuations than its bulk-rock counterpart, which indicates that deeper water masses were not as much subjected to density fluctuations as sea-surface water. The comparison of belemnite and bulk-rock carbon isotope records allows observing the impact of regional influence exerted by platform carbonate ooze shedding on the carbon cycle. Discrepancies in the two records are observed during time of photozoan carbonate platform growth. The strontium isotopic record shows a gradual increase from the uppermost Berriasian to the uppermost lower Barremian followed by a rapid decrease until the uppermost Barremian and a renewed small increase within the lowermost Aptian. The major inflection point in the uppermost lower Barremian appears to predate the onset in the formation of the Ontong-Java volcanic plateau.     

Assessing Greenland ice mass loss by means of point-mass modeling: a viable methodology

Greenland ice mass loss is one of the most serious phenomena of present-day global climate change. In this context, both the quantification of overall deglaciation rates and its spatial localization are highly significant. We have thoroughly investigated the technique of point-mass modeling in order to derive mass-balance patterns from GRACE (Gravity Recovery And Climate Experiment) gravimetry. The method infers mass variations on the Earth’s surface from gravitational signals at satellite altitude. In order to solve for point-mass changes, we applied least-squares adjustment. Due to downward continuation, numerical stabilization of the inversion process gains particular significance. We stabilized the ill-posed problem by Tikhonov regularization. Our simulation and real data experiments show that point-mass modeling provides both rational deglaciation rates and high-resolution spatial mass variation patterns.     

Contrasting time trends of organic contaminants in Antarctic pelagic and benthic food webs

We demonstrate that pelagic Antarctic seabirds show significant decreases in concentrations of some persistent organic pollutants. Trends in Adélie penguins and Southern fulmars fit in a general pattern revealed by a broad literature review. Downward trends are also visible in pelagic fish, contrasting sharply with steady or increasing concentrations in Antarctic benthic organisms. Transfer of contaminants between Antarctic pelagic and benthic food webs is associated with seasonal sea-ice dynamics which may influence the balance between the final receptors of contaminants under different climatic conditions. This complicates the predictability of future trends of emerging compounds in the Antarctic ecosystem, such as of the brominated compounds that we detected in Antarctic petrels. The discrepancy in trends between pelagic and benthic organisms shows that Antarctic biota are still final receptors of globally released organic contaminants and it remains questionable whether the total environmental burden of contaminants in the Antarctic ecosystem is declining.     

Thermodynamics, structure, dynamics, and freezing of Mg2SiO4 liquid at high pressure

We perform first principles molecular dynamics simulations of Mg₂SiO₄ liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).     

Seasonal and short-term controls of riparian oxygen dynamics and the implications for redox processes

Riparian zones are highly-dynamic transition zones between surface water (SW) and groundwater (GW) and function as key biogeochemical-reactors for solutes transitioning between both compartments. Infiltration of SW rich in dissolved oxygen (DO) into the riparian aquifer can supress removal processes of redox sensitive compounds like NO₃⁻, a nutrient harmful for the aquatic ecosystem at high concentrations. Seasonal and short-term variations of temperature and hydrologic conditions can influence biogeochemical reaction rates and thus the prevailing redox conditions in the riparian zone. We combined GW tracer-tests and a 1-year high-frequency dataset of DO with data-driven simulations of DO consumption to assess the effects of seasonal and event-scale variations in temperature and transit-times on the reactive transport of DO. Damköhler numbers for DO consumption (DA_DO) were used to characterize the system in terms of DO turnover potential. Our results suggest that seasonal and short-term variations in temperature are major controls for DO turnover and the resulting concentrations at our field site, while transit-times are of minor importance. Seasonal variations of temperature in GW lead to shifts from transport-limited (DA_DO > 1) to reaction-limited conditions (DA_DO < 1), while short-term events were found to have minor impacts on the state of the system, only resulting in slightly less transport-limited conditions due to decreasing temperature and transit-times. The data-driven analyses show that assuming constant water temperature along a flowpath can lead to an over- or underestimation of reaction rates by a factor of 2–3 due to different infiltrating water temperature at the SW–GW interface, whereas the assumption of constant transit-times results in incorrect estimates of NO₃⁻ removal potential based on DA_DO approach (40%–50% difference).     

Inventory and evolution of glacial lakes since the Little Ice Age: Lessons from the case of Switzerland

Retreating glaciers give way to new landscapes with lakes as an important element. In this study, we combined available data on lake outlines with historical orthoimagery and glacier outlines for six time periods since the end of the Little Ice Age (LIA; ~1850). We generated a glacial lake inventory for modern times (2016) and traced the evolution of glacial lakes that formed in the deglaciated area since the LIA. In this deglaciated area, a total of 1192 lakes formed over the period of almost 170 years, 987 of them still in existence in 2016. Their total water surface in 2016 was 6.22 ± 0.25 km². The largest lakes are > 0.4 km² (40 ha) in size, while the majority (> 90%) are smaller than 0.01 km². Annual increase rates in area and number peaked in 1946–1973, decreased towards the end of the 20th century, and reached a new high in the latest period 2006–2016. For a period of 43 years (1973–2016), we compared modelled overdeepenings from previous studies to actual lake genesis. For a better prioritization of formation probability, we included glacier-morphological criteria such as glacier width and visible crevassing. About 40% of the modelled overdeepened area actually got covered by lakes. The inclusion of morphological aspects clearly aided in defining a lake formation probability to be linked to each modelled overdeepening. Additional morphological variables, namely dam material and type, surface runoff, and freeboard, were compiled for a subset of larger and ice-contact lakes in 2016, constituting a basis for future hazard assessment.     

An improved sampling rule for mapping geopotential functions of a planet from a near polar orbit

One of the limiting factors in the determination of gravity field solutions is the spatial sampling. Especially during phases, when the satellite repeats its own track after a short time, the spatial resolution will be limited. The Nyquist rule-of-thumb for mapping geopotential functions of a planet, also referred to as the Colombo–Nyquist rule-of-thumb, provides a limit for the maximum achievable degree of a spherical harmonic development for repeat orbits. We show in this paper that this rule is too conservative, and solutions with better spatial resolutions are possible. A new rule is introduced which limits the maximum achievable order (not degree!) to be smaller than the number of revolutions if the difference between the number of revolutions and the number of nodal days is of odd parity and to be smaller than half the number of revolutions if the difference is of even parity. The dependence on the parity is reflected in the eigenvalue spectrum of the normal matrix and becomes especially important in the presence of noise. The rule is based on applying the Nyquist sampling theorem separately in North–South and East–West direction. This is only possible for satellites in highly inclined orbits like champ and grace. Tables for these two satellite missions are also provided which indicate the passed and (in case of grace) expected repeat cycles and possible degradations in the quality of the gravity field solutions.     

Modelling potential movement in constrained travel environments using rough space–time prisms

The widespread adoption of location-aware technologies (LATs) has afforded analysts new opportunities for efficiently collecting trajectory data of moving individuals. These technologies enable measuring trajectories as a finite sample set of time-stamped locations. The uncertainty related to both finite sampling and measurement errors makes it often difficult to reconstruct and represent a trajectory followed by an individual in space–time. Time geography offers an interesting framework to deal with the potential path of an individual in between two sample locations. Although this potential path may be easily delineated for travels along networks, this will be less straightforward for more nonnetwork-constrained environments. Current models, however, have mostly concentrated on network environments on the one hand and do not account for the spatiotemporal uncertainties of input data on the other hand. This article simultaneously addresses both issues by developing a novel methodology to capture potential movement between uncertain space–time points in obstacle-constrained travel environments.