Géosciences Environnement Toulouse

Key Points GRACE detects the groundwater depletion in North China. Official report underestimates the groundwater depletion rate in North China. High groundwater depletion exists in Beijing and Hebei province.

We review recent progress in our understanding of the global cycling of mercury (Hg), including best estimates of Hg concentrations and pool sizes in major environmental compartments and exchange processes within and between these reservoirs. Recent advances include the availability of new global datasets covering areas of the world where environmental Hg data were previously lacking; integration of these data into global and regional models is continually improving estimates of global Hg cycling. New analytical techniques, such as Hg stable isotope characterization, provide novel constraints of sources and transformation processes. The major global Hg reservoirs that are, and continue to be, affected by anthropogenic activities include the atmosphere (4.4–5.3 Gt), terrestrial environments (particularly soils: 250–1000 Gg), and aquatic ecosystems (e.g., oceans: 270–450 Gg). Declines in anthropogenic Hg emissions between 1990 and 2010 have led to declines in atmospheric Hg0 concentrations and HgII wet deposition in Europe and the US (− 1.5 to − 2.2% per year). Smaller atmospheric Hg0 declines (− 0.2% per year) have been reported in high northern latitudes, but not in the southern hemisphere, while increasing atmospheric Hg loads are still reported in East Asia. New observations and updated models now suggest high concentrations of oxidized HgII in the tropical and subtropical free troposphere where deep convection can scavenge these HgII reservoirs. As a result, up to 50% of total global wet HgII deposition has been predicted to occur to tropical oceans. Ocean Hg0 evasion is a large source of present-day atmospheric Hg (approximately 2900 Mg/year; range 1900–4200 Mg/year). Enhanced seawater Hg0 levels suggest enhanced Hg0 ocean evasion in the intertropical convergence zone, which may be linked to high HgII deposition. Estimates of gaseous Hg0 emissions to the atmosphere over land, long considered a critical Hg source, have been revised downward, and most terrestrial environments now are considered net sinks of atmospheric Hg due to substantial Hg uptake by plants. Litterfall deposition by plants is now estimated at 1020–1230 Mg/year globally. Stable isotope analysis and direct flux measurements provide evidence that in many ecosystems Hg0 deposition via plant inputs dominates, accounting for 57–94% of Hg in soils. Of global aquatic Hg releases, around 50% are estimated to occur in China and India, where Hg drains into the West Pacific and North Indian Oceans. A first inventory of global freshwater Hg suggests that inland freshwater Hg releases may be dominated by artisanal and small-scale gold mining (ASGM; approximately 880 Mg/year), industrial and wastewater releases (220 Mg/year), and terrestrial mobilization (170–300 Mg/year). For pelagic ocean regions, the dominant source of Hg is atmospheric deposition; an exception is the Arctic Ocean, where riverine and coastal erosion is likely the dominant source. Ocean water Hg concentrations in the North Atlantic appear to have declined during the last several decades but have increased since the mid-1980s in the Pacific due to enhanced atmospheric deposition from the Asian continent. Finally, we provide examples of ongoing and anticipated changes in Hg cycling due to emission, climate, and land use changes. It is anticipated that future emissions changes will be strongly dependent on ASGM, as well as energy use scenarios and technology requirements implemented under the Minamata Convention. We predict that land use and climate change impacts on Hg cycling will be large and inherently linked to changes in ecosystem function and global atmospheric and ocean circulations. Our ability to predict multiple and simultaneous changes in future Hg global cycling and human exposure is rapidly developing but requires further enhancement.

Ecologists have long sought to understand the factors controlling the structure of savanna vegetation. Using data from 2154 sites in savannas across Africa, Australia, and South America, we found that increasing moisture availability drives increases in fire and tree basal area, whereas fire reduces tree basal area. However, among continents, the magnitude of these effects varied substantially, so that a single model cannot adequately represent savanna woody biomass across these regions. Historical and environmental differences drive the regional variation in the functional relationships between woody vegetation, fire, and climate. These same differences will determine the regional responses of vegetation to future climates, with implications for global carbon stocks.

As the world warms, there is a profound need to improve projections of climate change. Although the latest Earth system models offer an unprecedented number of features, fundamental uncertainties continue to cloud our view of the future. Past climates provide the only opportunity to observe how the Earth system responds to high carbon dioxide, underlining a fundamental role for paleoclimatology in constraining future climate change. Here, we review the relevancy of paleoclimate information for climate prediction and discuss the prospects for emerging methodologies to further insights gained from past climates. Advances in proxy methods and interpretations pave the way for the use of past climates for model evaluation-a practice that we argue should be widely adopted.

Abstract Combined analyses of deep tow magnetic anomalies and International Ocean Discovery Program Expedition 349 cores show that initial seafloor spreading started around 33 Ma in the northeastern South China Sea (SCS), but varied slightly by 1–2 Myr along the northern continent‐ocean boundary (COB). A southward ridge jump of ∼20 km occurred around 23.6 Ma in the East Subbasin; this timing also slightly varied along the ridge and was coeval to the onset of seafloor spreading in the Southwest Subbasin, which propagated for about 400 km southwestward from ∼23.6 to ∼21.5 Ma. The terminal age of seafloor spreading is ∼15 Ma in the East Subbasin and ∼16 Ma in the Southwest Subbasin. The full spreading rate in the East Subbasin varied largely from ∼20 to ∼80 km/Myr, but mostly decreased with time except for the period between ∼26.0 Ma and the ridge jump (∼23.6 Ma), within which the rate was the fastest at ∼70 km/Myr on average. The spreading rates are not correlated, in most cases, to magnetic anomaly amplitudes that reflect basement magnetization contrasts. Shipboard magnetic measurements reveal at least one magnetic reversal in the top 100 m of basaltic layers, in addition to large vertical intensity variations. These complexities are caused by late‐stage lava flows that are magnetized in a different polarity from the primary basaltic layer emplaced during the main phase of crustal accretion. Deep tow magnetic modeling also reveals this smearing in basement magnetizations by incorporating a contamination coefficient of 0.5, which partly alleviates the problem of assuming a magnetic blocking model of constant thickness and uniform magnetization. The primary contribution to magnetic anomalies of the SCS is not in the top 100 m of the igneous basement.

Abiotic release of nitrous acid (HONO) in equilibrium with soil nitrite (NO2(-)) was suggested as an important contributor to the missing source of atmospheric HONO and hydroxyl radicals (OH). The role of total soil-derived HONO in the biogeochemical and atmospheric nitrogen cycles, however, has remained unknown. In laboratory experiments, we found that for nonacidic soils from arid and arable areas, reactive nitrogen emitted as HONO is comparable with emissions of nitric oxide (NO). We show that ammonia-oxidizing bacteria can directly release HONO in quantities larger than expected from the acid-base and Henry's law equilibria of the aqueous phase in soil. This component of the nitrogen cycle constitutes an additional loss term for fixed nitrogen in soils and a source for reactive nitrogen in the atmosphere.

Global plastic litter pollution has been increasing alongside demand since plastic products gained commercial popularity in the 1930's. Current plastic pollutant research has generally assumed that once plastics enter the ocean they are there to stay, retained permanently within the ocean currents, biota or sediment until eventual deposition on the sea floor or become washed up onto the beach. In contrast to this, we suggest it appears that some plastic particles could be leaving the sea and entering the atmosphere along with sea salt, bacteria, virus' and algae. This occurs via the process of bubble burst ejection and wave action, for example from strong wind or sea state turbulence. In this manuscript we review evidence from the existing literature which is relevant to this theory and follow this with a pilot study which analyses microplastics (MP) in sea spray. Here we show first evidence of MP particles, analysed by μRaman, in marine boundary layer air samples on the French Atlantic coast during both onshore (average of 2.9MP/m3) and offshore (average of 9.6MP/m3) winds. Notably, during sampling, the convergence of sea breeze meant our samples were dominated by sea spray, increasing our capacity to sample MPs if they were released from the sea. Our results indicate a potential for MPs to be released from the marine environment into the atmosphere by sea-spray giving a globally extrapolated figure of 136000 ton/yr blowing on shore.

Abstract. Natural and human-ignited fires affect all major biomes, altering ecosystem structure, biogeochemical cycles and atmospheric composition. Satellite observations provide global data on spatiotemporal patterns of biomass burning and evidence for the rapid changes in global fire activity in response to land management and climate. Satellite imagery also provides detailed information on the daily or sub-daily position of fires that can be used to understand the dynamics of individual fires. The Global Fire Atlas is a new global dataset that tracks the dynamics of individual fires to determine the timing and location of ignitions, fire size and duration, and daily expansion, fire line length, speed, and direction of spread. Here, we present the underlying methodology and Global Fire Atlas results for 2003–2016 derived from daily moderate-resolution (500 m) Collection 6 MCD64A1 burned-area data. The algorithm identified 13.3 million individual fires over the study period, and estimated fire perimeters were in good agreement with independent data for the continental United States. A small number of large fires dominated sparsely populated arid and boreal ecosystems, while burned area in agricultural and other human-dominated landscapes was driven by high ignition densities that resulted in numerous smaller fires. Long-duration fires in boreal regions and natural landscapes in the humid tropics suggest that fire season length exerts a strong control on fire size and total burned area in these areas. In arid ecosystems with low fuel densities, high fire spread rates resulted in large, short-duration fires that quickly consumed available fuels. Importantly, multiday fires contributed the majority of burned area in all biomass burning regions. A first analysis of the largest, longest and fastest fires that occurred around the world revealed coherent regional patterns of extreme fires driven by large-scale climate forcing. Global Fire Atlas data are publicly available through http://www.globalfiredata.org (last access: 9 August 2018) and https://doi.org/10.3334/ORNLDAAC/1642, and individual fire information and summary data products provide new information for benchmarking fire models within ecosystem and Earth system models, understanding vegetation–fire feedbacks, improving global emissions estimates, and characterizing the changing role of fire in the Earth system.

Key Points Remote sensing validation of a hydrologic‐hydrodynamic model of the Amazon Uncertainty of precipitation and river‐floodplain parameters cause model errors Importance of floodplains and backwater effects on flood waves traveling

Environmental context Natural organic matter exerts a powerful control on chemical conditions in waters and soils, affecting pH and influencing the biological availability, transport and retention of metals. To quantify the reactions, we collated a wealth of laboratory data covering 40 metals and acid–base reactions, and used them to parameterise the latest in a series of Humic Ion-Binding Models. Model VII is now available to interpret field data, and contribute to the prediction of environmental chemistry. Abstract Humic Ion-Binding Model VII aims to predict the competitive reactions of protons and metals with natural organic matter in soils and waters, based on laboratory results with isolated humic and fulvic acids (HA and FA). Model VII is simpler in its postulated multidentate metal binding sites than the previous Model VI. Three model parameters were eliminated by using a formal relationship between monodentate binding to strong- and weak-acid oxygen-containing ligands, and removing factors that provide ranges of ligand binding strengths. Thus Model VII uses a single adjustable parameter, the equilibrium constant for monodentate binding to strong-acid (carboxylate) groups (KMA), for each metallic cation. Proton-binding parameters, and mean values of log KMA were derived by fitting 248 published datasets (28 for protons, 220 for cationic metals). Default values of log KMA for FA were obtained by combining the fitted values for FA, results for HA, and the relationship for different metals between log KMA and equilibrium constants for simple oxygen-containing ligands. The equivalent approach was used for HA. The parameterised model improves on Model VI by incorporating more metals (40), providing better descriptions of metal binding at higher pH, and through more internally consistent parameter values.

H2O, CO2, SO2, O2, H2, H2S, HCl, chlorinated hydrocarbons, NO, and other trace gases were evolved during pyrolysis of two mudstone samples acquired by the Curiosity rover at Yellowknife Bay within Gale crater, Mars. H2O/OH-bearing phases included 2:1 phyllosilicate(s), bassanite, akaganeite, and amorphous materials. Thermal decomposition of carbonates and combustion of organic materials are candidate sources for the CO2. Concurrent evolution of O2 and chlorinated hydrocarbons suggests the presence of oxychlorine phase(s). Sulfides are likely sources for sulfur-bearing species. Higher abundances of chlorinated hydrocarbons in the mudstone compared with Rocknest windblown materials previously analyzed by Curiosity suggest that indigenous martian or meteoritic organic carbon sources may be preserved in the mudstone; however, the carbon source for the chlorinated hydrocarbons is not definitively of martian origin.

Abstract Gravimetry is a well-established technique for the determination of sub-surface mass distribution needed in several fields of geoscience, and various types of gravimeters have been developed over the last 50 years. Among them, quantum gravimeters based on atom interferometry have shown top-level performance in terms of sensitivity, long-term stability and accuracy. Nevertheless, they have remained confined to laboratories due to their complex operation and high sensitivity to the external environment. Here we report on a novel, transportable, quantum gravimeter that can be operated under real world conditions by non-specialists, and measure the absolute gravitational acceleration continuously with a long-term stability below 10 nm.s −2 (1 μ Gal). It features several technological innovations that allow for high-precision gravity measurements, while keeping the instrument light and small enough for field measurements. The instrument was characterized in detail and its stability was evaluated during a month-long measurement campaign.

Abstract. Polar amplification – the phenomenon where external radiative forcing produces a larger change in surface temperature at high latitudes than the global average – is a key aspect of anthropogenic climate change, but its causes and consequences are not fully understood. The Polar Amplification Model Intercomparison Project (PAMIP) contribution to the sixth Coupled Model Intercomparison Project (CMIP6; Eyring et al., 2016) seeks to improve our understanding of this phenomenon through a coordinated set of numerical model experiments documented here. In particular, PAMIP will address the following primary questions: (1) what are the relative roles of local sea ice and remote sea surface temperature changes in driving polar amplification? (2) How does the global climate system respond to changes in Arctic and Antarctic sea ice? These issues will be addressed with multi-model simulations that are forced with different combinations of sea ice and/or sea surface temperatures representing present-day, pre-industrial and future conditions. The use of three time periods allows the signals of interest to be diagnosed in multiple ways. Lower-priority tier experiments are proposed to investigate additional aspects and provide further understanding of the physical processes. These experiments will address the following specific questions: what role does ocean–atmosphere coupling play in the response to sea ice? How and why does the atmospheric response to Arctic sea ice depend on the pattern of sea ice forcing? How and why does the atmospheric response to Arctic sea ice depend on the model background state? What have been the roles of local sea ice and remote sea surface temperature in polar amplification, and the response to sea ice, over the recent period since 1979? How does the response to sea ice evolve on decadal and longer timescales? A key goal of PAMIP is to determine the real-world situation using imperfect climate models. Although the experiments proposed here form a coordinated set, we anticipate a large spread across models. However, this spread will be exploited by seeking “emergent constraints” in which model uncertainty may be reduced by using an observable quantity that physically explains the intermodel spread. In summary, PAMIP will improve our understanding of the physical processes that drive polar amplification and its global climate impacts, thereby reducing the uncertainties in future projections and predictions of climate change and variability.

We present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Frchet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 'Sesame' of our widely used open source spectral-element package SPECFEM3D.

Atmospheric correction over inland and coastal waters is one of the major remaining challenges in aquatic remote sensing, often hindering the quantitative retrieval of biogeochemical variables and analysis of their spatial and temporal variability within aquatic environments. The Atmospheric Correction Intercomparison Exercise (ACIX-Aqua), a joint NASA – ESA activity, was initiated to enable a thorough evaluation of eight state-of-the-art atmospheric correction (AC) processors available for Landsat-8 and Sentinel-2 data processing. Over 1000 radiometric matchups from both freshwaters (rivers, lakes, reservoirs) and coastal waters were utilized to examine the quality of derived aquatic reflectances (ρ̂w). This dataset originated from two sources: Data gathered from the international scientific community (henceforth called Community Validation Database, CVD), which captured predominantly inland water observations, and the Ocean Color component of AERONET measurements (AERONET-OC), representing primarily coastal ocean environments. This volume of data permitted the evaluation of the AC processors individually (using all the matchups) and comparatively (across seven different Optical Water Types, OWTs) using common matchups. We found that the performance of the AC processors differed for CVD and AERONET-OC matchups, likely reflecting inherent variability in aquatic and atmospheric properties between the two datasets. For the former, the median errors in ρ̂w560 and ρ̂w664 were found to range from 20 to 30% for best-performing processors. Using the AERONET-OC matchups, our performance assessments showed that median errors within the 15–30% range in these spectral bands may be achieved. The largest uncertainties were associated with the blue bands (25 to 60%) for best-performing processors considering both CVD and AERONET-OC assessments. We further assessed uncertainty propagation to the downstream products such as near-surface concentration of chlorophyll-a (Chla) and Total Suspended Solids (TSS). Using satellite matchups from the CVD along with in situ Chla and TSS, we found that 20–30% uncertainties in ρ̂w490≤λ≤743nm yielded 25–70% uncertainties in derived Chla and TSS products for top-performing AC processors. We summarize our results using performance matrices guiding the satellite user community through the OWT-specific relative performance of AC processors. Our analysis stresses the need for better representation of aerosols, particularly absorbing ones, and improvements in corrections for sky- (or sun-) glint and adjacency effects, in order to achieve higher quality downstream products in freshwater and coastal ecosystems.

Abstract A multimodel, multiresolution set of simulations over the period 1950–2014 using a common forcing protocol from CMIP6 HighResMIP have been completed by six modeling groups. Analysis of tropical cyclone performance using two different tracking algorithms suggests that enhanced resolution toward 25 km typically leads to more frequent and stronger tropical cyclones, together with improvements in spatial distribution and storm structure. Both of these factors reduce typical GCM biases seen at lower resolution. Using single ensemble members of each model, there is little evidence of systematic improvement in interannual variability in either storm frequency or accumulated cyclone energy as compared with observations when resolution is increased. Changes in the relationships between large-scale drivers of climate variability and tropical cyclone variability in the Atlantic Ocean are also not robust to model resolution. However, using a larger ensemble of simulations (of up to 14 members) with one model at different resolutions does show evidence of increased skill at higher resolution. The ensemble mean correlation of Atlantic interannual tropical cyclone variability increases from ~0.5 to ~0.65 when resolution increases from 250 to 100 km. In the northwestern Pacific Ocean the skill keeps increasing with 50-km resolution to 0.7. These calculations also suggest that more than six members are required to adequately distinguish the impact of resolution within the forced signal from the weather noise.

The Gravity Recovery and Climate Experiment (GRACE) satellite mission, which was in operation from March 2002 to June 2017, was the first remote sensing mission to provide temporal variations of Terrestrial Water Storage (TWS), which is the sum of the water masses that were contained in the soil column (i.e., snow, surface water, soil moisture, and groundwater), at a spatial resolution of a few hundred kilometers. As in situ level measurements are generally not sufficiently available for monitoring groundwater changes at the regional-scale, this unique dataset, combined with external information, is widely used to quantify the interannual variations of groundwater storage in the world’s major aquifers. GRACE-based groundwater changes revealed significant aquifer depletion over large regions, such as the Middle East, the northwest India aquifer, the North China Plain aquifer, the Murray-Darling Basin in Australia, the High Plains, and the California Central Valley aquifers in the United States of America (USA), but were also used to estimate groundwater-related parameters such as the specific yield, which relates groundwater level to storage, or to define the indices of groundwater depletion and stress. In this review, the approaches used for estimating groundwater storage variations are presented along with the main applications of GRACE data for groundwater monitoring. Issues that were related to the use of GRACE-based TWS are also addressed.

The long-term isolation of South America during most of the Cenozoic produced a highly peculiar terrestrial vertebrate biota, with a wide array of mammal groups, among which caviomorph rodents and platyrrhine primates are Mid-Cenozoic immigrants. In the absence of indisputable pre-Oligocene South American rodents or primates, the mode, timing and biogeography of these extraordinary dispersals remained debated. Here, we describe South America's oldest known rodents, based on a new diverse caviomorph assemblage from the late Middle Eocene (approx. 41 Ma) of Peru, including five small rodents with three stem caviomorphs. Instead of being tied to the Eocene/Oligocene global cooling and drying episode (approx. 34 Ma), as previously considered, the arrival of caviomorphs and their initial radiation in South America probably occurred under much warmer and wetter conditions, around the Mid-Eocene Climatic Optimum. Our phylogenetic results reaffirm the African origin of South American rodents and support a trans-Atlantic dispersal of these mammals during Middle Eocene times. This discovery further extends the gap (approx. 15 Myr) between first appearances of rodents and primates in South America.

Field projects are beginning to demonstrate the potential for carbon storage in basaltic rocks.

Gaseous elemental mercury (GEM) is the dominant form of mercury in the atmosphere. Its conversion into oxidized gaseous and particulate forms is thought to drive atmospheric mercury wet deposition to terrestrial and aquatic ecosystems, where it can be subsequently transformed into toxic methylmercury. The contribution of mercury dry deposition is however largely unconstrained. Here we examine mercury mass balance and mercury stable isotope composition in a peat bog ecosystem. We find that isotope signatures of living sphagnum moss (Δ(199)Hg = -0.11 ± 0.09‰, Δ(200)Hg = 0.03 ± 0.02‰, 1σ) and recently accumulated peat (Δ(199)Hg = -0.22 ± 0.06‰, Δ(200)Hg = 0.00 ± 0.04‰, 1σ) are characteristic of GEM (Δ(199)Hg = -0.17 ± 0.07‰, Δ(200)Hg = -0.05 ± 0.02‰, 1σ), and differs from wet deposition (Δ(199)Hg = 0.73 ± 0.15‰, Δ(200)Hg = 0.21 ± 0.04‰, 1σ). Sphagnum covered during three years by transparent and opaque surfaces, which eliminate wet deposition, continue to accumulate Hg. Sphagnum Hg isotope signatures indicate accumulation to take place by GEM dry deposition, and indicate little photochemical re-emission. We estimate that atmospheric mercury deposition to the peat bog surface is dominated by GEM dry deposition (79%) rather than wet deposition (21%). Consequently, peat deposits are potential records of past atmospheric GEM concentrations and isotopic composition.