Institut des Géosciences de l'Environnement

Agricultural food production and agriculturally-related change in land use substantially contribute to greenhouse-gas
\nemissions worldwide. Four-fifths of agricultural emissions arise from the livestock sector. Although livestock products
\nare a source of some essential nutrients, they provide large amounts of saturated fat, which is a known risk factor for
\ncardiovascular disease. We considered potential strategies for the agricultural sector to meet the target recommended by the UK Committee on Climate Change to reduce UK emissions from the concentrations recorded in 1990 by 80% by 2050, which would require a 50% reduction by 2030. With use of the UK as a case study, we identified that a
\ncombination of agricultural technological improvements and a 30% reduction in livestock production would be needed to meet this target; in the absence of good emissions data from Brazil, we assumed for illustrative purposes that the required reductions would be the same for our second case study in São Paulo city. We then used these data to model the potential benefits of reduced consumption of livestock products on the burden of ischaemic heart disease: disease burden would decrease by about 15% in the UK (equivalent to 2850 disability-adjusted life-years [DALYs] per million
\npopulation in 1 year) and 16% in São Paulo city (equivalent to 2180 DALYs per million population in 1 year). Although likely to yield benefits to health, such a strategy will probably encounter cultural, political, and commercial resistance, and face technical challenges. Coordinated intersectoral action is needed across agricultural, nutritional, public health,
\nand climate change communities worldwide to provide affordable, healthy, low-emission diets for all societies.

We present the global general circulation model IPSL-CM5 developed to study the long-term response of the climate system to natural and anthropogenic forcings as part of the 5th Phase of the Coupled Model Intercomparison Project (CMIP5). This model includes an interactive carbon cycle, a representation of tropospheric and stratospheric chemistry, and a comprehensive representation of aerosols. As it represents the principal dynamical, physical, and bio-geochemical processes relevant to the climate system, it may be referred to as an Earth System Model. However, the IPSL-CM5 model may be used in a multitude of configurations associated with different boundary conditions and with a range of complexities in terms of processes and interactions. This paper presents an overview of the different model components and explains how they were coupled and used to simulate historical climate changes over the past 150 years and different scenarios of future climate change. A single version of the IPSL-CM5 model (IPSL-CM5A-LR) was used to provide climate projections associated with different socio-economic scenarios, including the different Representative Concentration Pathways considered by CMIP5 and several scenarios from the Special Report on Emission Scenarios considered by CMIP3. Results suggest that the magnitude of global warming projections primarily depends on the socio-economic scenario considered, that there is potential for an aggressive mitigation policy to limit global warming to about two degrees, and that the behavior of some components of the climate system such as the Arctic sea ice and the Atlantic Meridional Overturning Circulation may change drastically by the end of the twenty-first century in the case of a no climate policy scenario. Although the magnitude of regional temperature and precipitation changes depends fairly linearly on the magnitude of the projected global warming (and thus on the scenario considered), the geographical pattern of these changes is strikingly similar for the different scenarios. The representation of atmospheric physical processes in the model is shown to strongly influence the simulated climate variability and both the magnitude and pattern of the projected climate changes.

, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

Abstract This study presents the global climate model IPSL‐CM6A‐LR developed at Institut Pierre‐Simon Laplace (IPSL) to study natural climate variability and climate response to natural and anthropogenic forcings as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This article describes the different model components, their coupling, and the simulated climate in comparison to previous model versions. We focus here on the representation of the physical climate along with the main characteristics of the global carbon cycle. The model's climatology, as assessed from a range of metrics (related in particular to radiation, temperature, precipitation, and wind), is strongly improved in comparison to previous model versions. Although they are reduced, a number of known biases and shortcomings (e.g., double Intertropical Convergence Zone [ITCZ], frequency of midlatitude wintertime blockings, and El Niño–Southern Oscillation [ENSO] dynamics) persist. The equilibrium climate sensitivity and transient climate response have both increased from the previous climate model IPSL‐CM5A‐LR used in CMIP5. A large ensemble of more than 30 members for the historical period (1850–2018) and a smaller ensemble for a range of emissions scenarios (until 2100 and 2300) are also presented and discussed.

We use updated drainage inventory, ice thickness, and ice velocity data to calculate the grounding line ice discharge of 176 basins draining the Antarctic Ice Sheet from 1979 to 2017. We compare the results with a surface mass balance model to deduce the ice sheet mass balance. The total mass loss increased from 40 ± 9 Gt/y in 1979-1990 to 50 ± 14 Gt/y in 1989-2000, 166 ± 18 Gt/y in 1999-2009, and 252 ± 26 Gt/y in 2009-2017. In 2009-2017, the mass loss was dominated by the Amundsen/Bellingshausen Sea sectors, in West Antarctica (159 ± 8 Gt/y), Wilkes Land, in East Antarctica (51 ± 13 Gt/y), and West and Northeast Peninsula (42 ± 5 Gt/y). The contribution to sea-level rise from Antarctica averaged 3.6 ± 0.5 mm per decade with a cumulative 14.0 ± 2.0 mm since 1979, including 6.9 ± 0.6 mm from West Antarctica, 4.4 ± 0.9 mm from East Antarctica, and 2.5 ± 0.4 mm from the Peninsula (i.e., East Antarctica is a major participant in the mass loss). During the entire period, the mass loss concentrated in areas closest to warm, salty, subsurface, circumpolar deep water (CDW), that is, consistent with enhanced polar westerlies pushing CDW toward Antarctica to melt its floating ice shelves, destabilize the glaciers, and raise sea level.

We reconstruct the mass balance of the Greenland Ice Sheet using a comprehensive survey of thickness, surface elevation, velocity, and surface mass balance (SMB) of 260 glaciers from 1972 to 2018. We calculate mass discharge, D, into the ocean directly for 107 glaciers (85% of D) and indirectly for 110 glaciers (15%) using velocity-scaled reference fluxes. The decadal mass balance switched from a mass gain of +47 ± 21 Gt/y in 1972-1980 to a loss of 51 ± 17 Gt/y in 1980-1990. The mass loss increased from 41 ± 17 Gt/y in 1990-2000, to 187 ± 17 Gt/y in 2000-2010, to 286 ± 20 Gt/y in 2010-2018, or sixfold since the 1980s, or 80 ± 6 Gt/y per decade, on average. The acceleration in mass loss switched from positive in 2000-2010 to negative in 2010-2018 due to a series of cold summers, which illustrates the difficulty of extrapolating short records into longer-term trends. Cumulated since 1972, the largest contributions to global sea level rise are from northwest (4.4 ± 0.2 mm), southeast (3.0 ± 0.3 mm), and central west (2.0 ± 0.2 mm) Greenland, with a total 13.7 ± 1.1 mm for the ice sheet. The mass loss is controlled at 66 ± 8% by glacier dynamics (9.1 mm) and 34 ± 8% by SMB (4.6 mm). Even in years of high SMB, enhanced glacier discharge has remained sufficiently high above equilibrium to maintain an annual mass loss every year since 1998.

Increased concentrations of atmospheric greenhouse gases have led to a global mean surface temperature 1.0°C higher than during the pre-industrial period. We expand on the recent IPCC Special Report on global warming of 1.5°C and review the additional risks associated with higher levels of warming, each having major implications for multiple geographies, climates, and ecosystems. Limiting warming to 1.5°C rather than 2.0°C would be required to maintain substantial proportions of ecosystems and would have clear benefits for human health and economies. These conclusions are relevant for people everywhere, particularly in low- and middle-income countries, where the escalation of climate-related risks may prevent the achievement of the United Nations Sustainable Development Goals.

Abstract The European Project for Ice Coring in Antarctica Dome ice core from Dome C (EDC) has allowed for the reconstruction of atmospheric CO 2 concentrations for the last 800,000 years. Here we revisit the oldest part of the EDC CO 2 record using different air extraction methods and sections of the core. For our established cracker system, we found an analytical artifact, which increases over the deepest 200 m and reaches 10.1 ± 2.4 ppm in the oldest/deepest part. The governing mechanism is not yet fully understood, but it is related to insufficient gas extraction in combination with ice relaxation during storage and ice structure. The corrected record presented here resolves partly ‐ but not completely ‐ the issue with a different correlation between CO 2 and Antarctic temperatures found in this oldest part of the records. In addition, we provide here an update of 800,000 years atmospheric CO 2 history including recent studies covering the last glacial cycle.

Abstract. The recent evolution of Pamir-Karakoram-Himalaya (PKH) glaciers, widely acknowledged as valuable high-altitude as well as mid-latitude climatic indicators, remains poorly known. To estimate the region-wide glacier mass balance for 9 study sites spread from the Pamir to the Hengduan Shan (eastern Himalaya), we compared the 2000 Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) to recent (2008–2011) DEMs derived from SPOT5 stereo imagery. During the last decade, the region-wide glacier mass balances were contrasted with moderate mass losses in the eastern and central Himalaya (−0.22 ± 0.12 m w.e. yr−1 to −0.33 ± 0.14 m w.e. yr−1) and larger losses in the western Himalaya (−0.45 ± 0.13 m w.e. yr−1). Recently reported slight mass gain or balanced mass budget of glaciers in the central Karakoram is confirmed for a larger area (+0.10 ± 0.16 m w.e. yr−1) and also observed for glaciers in the western Pamir (+0.14 ± 0.13 m w.e. yr−1). Thus, the "Karakoram anomaly" should be renamed the "Pamir-Karakoram anomaly", at least for the last decade. The overall mass balance of PKH glaciers, −0.14 ± 0.08 m w.e. yr−1, is two to three times less negative than the global average for glaciers distinct from the Greenland and Antarctic ice sheets. Together with recent studies using ICESat and GRACE data, DEM differencing confirms a contrasted pattern of glacier mass change in the PKH during the first decade of the 21st century.

. Here we use field observations and air-quality modelling to quantify the major primary and secondary sources of particulate matter and of oxidative potential in Europe. We find that secondary inorganic components, crustal material and secondary biogenic organic aerosols control the mass concentration of particulate matter. By contrast, oxidative potential concentration is associated mostly with anthropogenic sources, in particular with fine-mode secondary organic aerosols largely from residential biomass burning and coarse-mode metals from vehicular non-exhaust emissions. Our results suggest that mitigation strategies aimed at reducing the mass concentrations of particulate matter alone may not reduce the oxidative potential concentration. If the oxidative potential can be linked to major health impacts, it may be more effective to control specific sources of particulate matter rather than overall particulate mass.

Abstract. The mountain cryosphere of mainland Europe is recognized to have important impacts on a range of environmental processes. In this paper, we provide an overview on the current knowledge on snow, glacier, and permafrost processes, as well as their past, current, and future evolution. We additionally provide an assessment of current cryosphere research in Europe and point to the different domains requiring further research. Emphasis is given to our understanding of climate–cryosphere interactions, cryosphere controls on physical and biological mountain systems, and related impacts. By the end of the century, Europe's mountain cryosphere will have changed to an extent that will impact the landscape, the hydrological regimes, the water resources, and the infrastructure. The impacts will not remain confined to the mountain area but also affect the downstream lowlands, entailing a wide range of socioeconomical consequences. European mountains will have a completely different visual appearance, in which low- and mid-range-altitude glaciers will have disappeared and even large valley glaciers will have experienced significant retreat and mass loss. Due to increased air temperatures and related shifts from solid to liquid precipitation, seasonal snow lines will be found at much higher altitudes, and the snow season will be much shorter than today. These changes in snow and ice melt will cause a shift in the timing of discharge maxima, as well as a transition of runoff regimes from glacial to nival and from nival to pluvial. This will entail significant impacts on the seasonality of high-altitude water availability, with consequences for water storage and management in reservoirs for drinking water, irrigation, and hydropower production. Whereas an upward shift of the tree line and expansion of vegetation can be expected into current periglacial areas, the disappearance of permafrost at lower altitudes and its warming at higher elevations will likely result in mass movements and process chains beyond historical experience. Future cryospheric research has the responsibility not only to foster awareness of these expected changes and to develop targeted strategies to precisely quantify their magnitude and rate of occurrence but also to help in the development of approaches to adapt to these changes and to mitigate their consequences. Major joint efforts are required in the domain of cryospheric monitoring, which will require coordination in terms of data availability and quality. In particular, we recognize the quantification of high-altitude precipitation as a key source of uncertainty in projections of future changes. Improvements in numerical modeling and a better understanding of process chains affecting high-altitude mass movements are the two further fields that – in our view – future cryospheric research should focus on.

Abstract. The Earth system model EC-Earth3 for contributions to CMIP6 is documented here, with its flexible coupling framework, major model configurations, a methodology for ensuring the simulations are comparable across different high-performance computing (HPC) systems, and with the physical performance of base configurations over the historical period. The variety of possible configurations and sub-models reflects the broad interests in the EC-Earth community. EC-Earth3 key performance metrics demonstrate physical behavior and biases well within the frame known from recent CMIP models. With improved physical and dynamic features, new Earth system model (ESM) components, community tools, and largely improved physical performance compared to the CMIP5 version, EC-Earth3 represents a clear step forward for the only European community ESM. We demonstrate here that EC-Earth3 is suited for a range of tasks in CMIP6 and beyond.

Abstract. The aim of this paper is to provide the community with a comprehensive overview of the studies of glaciers in the tropical Andes conducted in recent decades leading to the current status of the glaciers in the context of climate change. In terms of changes in surface area and length, we show that the glacier retreat in the tropical Andes over the last three decades is unprecedented since the maximum extension of the Little Ice Age (LIA, mid-17th–early 18th century). In terms of changes in mass balance, although there have been some sporadic gains on several glaciers, we show that the trend has been quite negative over the past 50 yr, with a mean mass balance deficit for glaciers in the tropical Andes that is slightly more negative than the one computed on a global scale. A break point in the trend appeared in the late 1970s with mean annual mass balance per year decreasing from −0.2 m w.e. in the period 1964–1975 to −0.76 m w.e. in the period 1976–2010. In addition, even if glaciers are currently retreating everywhere in the tropical Andes, it should be noted that this is much more pronounced on small glaciers at low altitudes that do not have a permanent accumulation zone, and which could disappear in the coming years/decades. Monthly mass balance measurements performed in Bolivia, Ecuador and Colombia show that variability of the surface temperature of the Pacific Ocean is the main factor governing variability of the mass balance at the decadal timescale. Precipitation did not display a significant trend in the tropical Andes in the 20th century, and consequently cannot explain the glacier recession. On the other hand, temperature increased at a significant rate of 0.10 °C decade−1 in the last 70 yr. The higher frequency of El Niño events and changes in its spatial and temporal occurrence since the late 1970s together with a warming troposphere over the tropical Andes may thus explain much of the recent dramatic shrinkage of glaciers in this part of the world.

Globally, thermodynamics explains an increase in atmospheric water vapor with warming of around 7%/°C near to the surface. In contrast, global precipitation and evaporation are constrained by the Earth's energy balance to increase at ∼2-3%/°C. However, this rate of increase is suppressed by rapid atmospheric adjustments in response to greenhouse gases and absorbing aerosols that directly alter the atmospheric energy budget. Rapid adjustments to forcings, cooling effects from scattering aerosol, and observational uncertainty can explain why observed global precipitation responses are currently difficult to detect but are expected to emerge and accelerate as warming increases and aerosol forcing diminishes. Precipitation increases with warming are expected to be smaller over land than ocean due to limitations on moisture convergence, exacerbated by feedbacks and affected by rapid adjustments. Thermodynamic increases in atmospheric moisture fluxes amplify wet and dry events, driving an intensification of precipitation extremes. The rate of intensification can deviate from a simple thermodynamic response due to in-storm and larger-scale feedback processes, while changes in large-scale dynamics and catchment characteristics further modulate the frequency of flooding in response to precipitation increases. Changes in atmospheric circulation in response to radiative forcing and evolving surface temperature patterns are capable of dominating water cycle changes in some regions. Moreover, the direct impact of human activities on the water cycle through water abstraction, irrigation, and land use change is already a significant component of regional water cycle change and is expected to further increase in importance as water demand grows with global population.

The future international Surface Water and Ocean Topography (SWOT) Mission, planned for launch in 2021, will make high-resolution 2D observations of sea-surface height using SAR radar interferometric techniques. SWOT will map the global and coastal oceans up to 77.6° latitude every 21 days over a swath of 120 km (20 km nadir gap). Today’s 2D mapped altimeter data can resolve ocean scales of 150 km wavelength whereas the SWOT measurement will extend our 2D observations down to 15-30 km, depending on sea state. SWOT will offer new opportunities to observe the oceanic dynamic processes at these scales, that are important in the generation and dissipation of kinetic energy in the ocean, and act as one of the main gateways connecting the interior of the ocean to the upper layer. The active vertical exchanges linked to these scales have impacts on the local and global budgets of heat and carbon, and on nutrients for biogeochemical cycles. This review paper highlights the issues being addressed by the SWOT science community to understand SWOT’s very precise SSH / surface pressure observations, and it explores how SWOT data will be combined with other satellite and in-situ data and models to better understand the upper ocean 4D circulation (x,y,z,t) over the next decade. SWOT’s new SAR-interferometry technology aims to observe ocean SSH scales down to 15-30 km in wavelength. At these scales, SSH includes “balanced” geostrophic eddy motions and high-frequency internal tides and internal waves. This presents both a challenge in reconstructing the 4D upper ocean circulation, or in the assimilation of SSH in models, but also an opportunity to have global observations of the 2D structure of these phenomena, and to learn more about their interactions. At these small scales, the ocean dynamics evolve rapidly, and combining SWOT 2D SSH data with other satellite or in-situ data with different space-time coverage is also a challenge. SWOT’s new technology will be a forerunner for the future altimetric observing system, and so advancing on these issues today will pave the way for our future.

-g carbon) loss to 70-Pg C gain. For the RCP8.5 projection, losses in soil carbon varied between 74 and 652 Pg C (mean loss, 341 Pg C). For the RCP4.5 projection, gains in vegetation carbon were largely responsible for the overall projected net gains in ecosystem carbon by 2299 (8- to 244-Pg C gains). In contrast, for the RCP8.5 projection, gains in vegetation carbon were not great enough to compensate for the losses of carbon projected by four of the five models; changes in ecosystem carbon ranged from a 641-Pg C loss to a 167-Pg C gain (mean, 208-Pg C loss). The models indicate that substantial net losses of ecosystem carbon would not occur until after 2100. This assessment suggests that effective mitigation efforts during the remainder of this century could attenuate the negative consequences of the permafrost carbon-climate feedback.

Abstract. The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120–800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on δ15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores.

Abstract. A long term data acquisition effort of profile soil moisture is under way in southwestern France at 13 automated weather stations. This ground network was developed in order to validate remote sensing and model soil moisture estimates. In this paper, both those in situ observations and a synthetic data set covering continental France are used to test a simple method to retrieve root zone soil moisture from a time series of surface soil moisture information. A recursive exponential filter equation using a time constant, T, is used to compute a soil water index. The Nash and Sutcliff coefficient is used as a criterion to optimise the T parameter for each ground station and for each model pixel of the synthetic data set. In general, the soil water indices derived from the surface soil moisture observations and simulations agree well with the reference root-zone soil moisture. Overall, the results show the potential of the exponential filter equation and of its recursive formulation to derive a soil water index from surface soil moisture estimates. This paper further investigates the correlation of the time scale parameter T with soil properties and climate conditions. While no significant relationship could be determined between T and the main soil properties (clay and sand fractions, bulk density and organic matter content), the modelled spatial variability and the observed inter-annual variability of T suggest that a weak climate effect may exist.

Abstract. Simulations of climate over the Last Millennium (850–1850 CE) have been incorporated into the third phase of the Paleoclimate Modelling Intercomparison Project (PMIP3). The drivers of climate over this period are chiefly orbital, solar, volcanic, changes in land use/land cover and some variation in greenhouse gas levels. While some of these effects can be easily defined, the reconstructions of solar, volcanic and land use-related forcing are more uncertain. We describe here the approach taken in defining the scenarios used in PMIP3, document the forcing reconstructions and discuss likely implications.

The European Project for Ice Coring in Antarctica Dome C ice core enables us to extend existing records of atmospheric methane (CH4) and nitrous oxide (N2O) back to 650,000 years before the present. A combined record of CH4 measured along the Dome C and the Vostok ice cores demonstrates, within the resolution of our measurements, that preindustrial concentrations over Antarctica have not exceeded 773 +/- 15 ppbv (parts per billion by volume) during the past 650,000 years. Before 420,000 years ago, when interglacials were cooler, maximum CH4 concentrations were only about 600 ppbv, similar to lower Holocene values. In contrast, the N2O record shows maximum concentrations of 278 +/- 7 ppbv, slightly higher than early Holocene values.