Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

Past and future evolution of Abies alba forests in Europe - comparison of a dynamic vegetation model with palaeo data and observations

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling–Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe.

Modeling the marine aragonite cycle: changes under rising carbon dioxide and its role in shallow water CaCO3 dissolution

The marine aragonite cycle has been included in the global biogeochemical model PISCES to study the role of aragonite in shallow water CaCO3 dissolution. Aragonite production is parameterized as a function of mesozooplankton biomass and aragonite saturation state of ambient waters. Observation-based estimates of marine carbonate production and dissolution are well reproduced by the model and about 60% of the combined CaCO3 water column dissolution from aragonite and calcite is simulated above 2000 m. In contrast, a calcite-only version yields a much smaller fraction. This suggests that the aragonite cycle should be included in models for a realistic representation of CaCO3 dissolution and alkalinity. For the SRES A2 CO2 scenario, production rates of aragonite are projected to notably decrease after 2050. By the end of this century, global aragonite production is reduced by 29% and total CaCO3 production by 19% relative to pre-industrial. Geographically, the effect from increasing atmospheric CO2, and the subsequent reduction in saturation state, is largest in the subpolar and polar areas where the modeled aragonite production is projected to decrease by 65% until 2100.

Modeling the stratospheric warming following the Mt. Pinatubo eruption: uncertainties in aerosol extinctions

In terms of atmospheric impact, the volcanic eruption of Mt. Pinatubo (1991) is the best characterized large eruption on record. We investigate here the model-derived stratospheric warming following the Pinatubo eruption as derived from SAGE II extinction data including recent improvements in the processing algorithm. This method, termed SAGE_4λ, makes use of the four wavelengths (385, 452, 525 and 1024 nm) of the SAGE II data when available, and uses a data-filling procedure in the opacity-induced "gap" regions. Using SAGE_4λ, we derived aerosol size distributions that properly reproduce extinction coefficients also at much longer wavelengths. This provides a good basis for calculating the absorption of terrestrial infrared radiation and the resulting stratospheric heating. However, we also show that the use of this data set in a global chemistry–climate model (CCM) still leads to stronger aerosol-induced stratospheric heating than observed, with temperatures partly even higher than the already too high values found by many models in recent general circulation model (GCM) and CCM intercomparisons. This suggests that the overestimation of the stratospheric warming after the Pinatubo eruption may not be ascribed to an insufficient observational database but instead to using outdated data sets, to deficiencies in the implementation of the forcing data, or to radiative or dynamical model artifacts. Conversely, the SAGE_4λ approach reduces the infrared absorption in the tropical tropopause region, resulting in a significantly better agreement with the post-volcanic temperature record at these altitudes.

Estimated strength of the Atlantic overturning circulation during the last deglaciation

The Atlantic meridional overturning circulation affects the latitudinal distribution of heat, and is a key component of the climate system. Proxy reconstructions, based on sedimentary Pa-231/Th-230 ratios and the difference between surface-and deep-water radiocarbon ages, indicate that during the last glacial period, the overturning circulation was reduced during millennial-scale periods of cooling(1-5). However, much debate exists over the robustness of these proxies(6-8). Here e combine proxy reconstructions of sea surface and air temperatures and a global climate model to quantitatively estimate changes in the strength of the Atlantic meridional overturning circulation during the last glacial period. We find that, relative to the Last Glacial Maximum, the overturning circulation was reduced by approximately 14 Sv during the cold Heinrich event 1. During the Younger Dryas cold event, the overturning circulation was reduced by approximately 12 Sv, relative to the preceding warm interval. These changes are consistent with qualitative estimates of the overturning circulation from sedimentary Pa-231/Th-230 ratios. In addition, we find that the strength of the overturning circulation during the Last Glacial Maximum and the Holocene epoch are indistinguishable within the uncertainty of the reconstruction.

Impact of an abrupt cooling event on interglacial methane emissions in northern peatlands

Rapid changes in atmospheric methane (CH4), temperature and precipitation are documented by Greenland ice core data both for glacial times (the so called Dansgaard-Oeschger (D-O) events) as well as for a cooling event in the early Holocene (the 8.2 kyr event). The onsets of D-O warm events are paralleled by abrupt increases in CH4 by up to 250 ppb in a few decades. Vice versa, the 8.2 kyr event is accompanied by an intermittent decrease in CH4 of about 80 ppb over 150 yr. The abrupt CH4 changes are thought to mainly originate from source emission variations in tropical and boreal wet ecosystems, but complex process oriented bottom-up model estimates of the changes in these ecosystems during rapid climate changes are still missing. Here we present simulations of CH4 emissions from northern peatlands with the LPJ-Bern dynamic global vegetation model. The model represents CH4 production and oxidation in soils and transport by ebullition, through plant aerenchyma, and by diffusion. Parameters are tuned to represent site emission data as well as inversion-based estimates of northern wetland emissions. The model is forced with climate input data from freshwater hosing experiments using the NCAR CSM1.4 climate model to simulate an abrupt cooling event. A concentration reduction of ~10 ppb is simulated per degree K change of mean northern hemispheric surface temperature in peatlands. Peatland emissions are equally sensitive to both changes in temperature and in precipitation. If simulated changes are taken as an analogy to the 8.2 kyr event, boreal peatland emissions alone could only explain 23 of the 80 ppb decline in atmospheric methane concentration. This points to a significant contribution to source changes from low latitude and tropical wetlands to this event.

Is there a common pattern of future gas-phase air pollution in Europe under diverse climate change scenarios?

Climate change alone influences future levels of tropospheric ozone and their precursors through modifications of gas-phase chemistry, transport, removal, and natural emissions. The goal of this study is to determine at what extent the modes of variability of gas-phase pollutants respond to different climate change scenarios over Europe. The methodology includes the use of the regional modeling system MM5 (regional climate model version)-CHIMERE for a target domain covering Europe. Two full-transient simulations covering from 1991–2050 under the SRES A2 and B2 scenarios driven by ECHO-G global circulation model have been compared. The results indicate that the spatial patterns of variability for tropospheric ozone are similar for both scenarios, but the magnitude of the change signal significantly differs for A2 and B2. The 1991–2050 simulations share common characteristics for their chemical behavior. As observed from the NO2 and α-pinene modes of variability, our simulations suggest that the enhanced ozone chemical activity is driven by a number of parameters, such as the warming-induced increase in biogenic emissions and, to a lesser extent, by the variation in nitrogen dioxide levels. For gas-phase pollutants, the general increasing trend for ozone found under A2 and B2 forcing is due to a multiplicity of climate factors, such as increased temperature, decreased wet removal associated with an overall decrease of precipitation in southern Europe, increased photolysis of primary and secondary pollutants as a consequence of lower cloudiness and increased biogenic emissions fueled by higher temperatures.

Land use change and nitrogen feedbacks constrain the trajectory of the land carbon sink

Our understanding of Earth's carbon climate system depends critically upon interactions between rising atmospheric CO2, changing land use, and nitrogen limitation on vegetation growth. Using a global land model, we show how these factors interact locally to generate the global land carbon sink over the past 200 years. Nitrogen constraints were alleviated by N2 fixation in the tropics and by atmospheric nitrogen deposition in extratropical regions. Nonlinear interactions between land use change and land carbon and nitrogen cycling originated from three major mechanisms: (i) a sink foregone that would have occurred without land use conversion; (ii) an accelerated response of secondary vegetation to CO2 and nitrogen, and (iii) a compounded clearance loss from deforestation. Over time, these nonlinear effects have become increasingly important and reduce the present-day net carbon sink by ~40% or 0.4 PgC yr−1.

Recent trends in Inner Asian forest dynamics to temperature and precipitation indicate high sensitivity to climate change

Semi-arid ecosystems play an important role in regulating global climate with the fate of these ecosystems in the Anthropocene depending upon interactions among temperature, precipitation, and CO2. However, in cool-arid environments, precipitation is not the only limitation to forest productivity. Interactions between changes in precipitation and air temperature may enhance soil moisture stress while simultaneously extending growing season length, with unclear consequences for net carbon uptake. This study evaluates recent trends in productivity and phenology of Inner Asian forests (in Mongolia and Northern China) using satellite remote sensing, dendrochronology, and dynamic global vegetation model (DGVM) simulations to quantify the sensitivity of forest dynamics to decadal climate variability and trends. Trends in photosynthetically active radiation fraction (FPAR) between 1982 and 2010 show a greening of about 7% of the region in spring (March, April, May), and 3% of the area ‘browning’ during summertime (June, July, August). These satellite observations of FPAR are corroborated by trends in NPP simulated by the LPJ DGVM. Spring greening trends in FPAR are mainly explained by long-term trends in precipitation whereas summer browning trends are correlated with decreasing precipitation. Tree ring data from 25 sites confirm annual growth increments are mainly limited by summer precipitation (June, July, August) in Mongolia, and spring precipitation in northern China (March, April, May), with relatively weak prior-year lag effects. An ensemble of climate projections from the IPCC CMIP3 models indicates that warming temperatures (spring, summer) are expected to be associated with higher summer precipitation, which combined with CO2 causes large increases in NPP and possibly even greater forest cover in the Mongolian steppe. In the absence of a strong direct CO2 fertilization effect on plant growth (e.g., due to nutrient limitation), water stress or decreased carbon gain from higher autotrophic respiration results in decreased productivity and loss of forest cover. The fate of these semi-arid ecosystems thus appears to hinge upon the magnitude and subtleties of CO2 fertilization effects, for which experimental observations in arid systems are needed to test and refine vegetation models.

Past and future carbon fluxes from land use change, shifting cultivation and wood harvest

Carbon emissions from anthropogenic land use (LU) and land use change (LUC) are quantified with a Dynamic Global Vegetation Model for the past and the 21st century following Representative Concentration Pathways (RCPs). Wood harvesting and parallel abandonment and expansion of agricultural land in areas of shifting cultivation are explicitly simulated (gross LUC) based on the Land Use Harmonization (LUH) dataset and a proposed alternative method that relies on minimum input data and generically accounts for gross LUC. Cumulative global LUC emissions are 72 GtC by 1850 and 243 GtC by 2004 and 27–151 GtC for the next 95 yr following the different RCP scenarios. The alternative method reproduces results based on LUH data with full transition information within <0.1 GtC/yr over the last decades and bears potential for applications in combination with other LU scenarios. In the last decade, shifting cultivation and wood harvest within remaining forests including slash each contributed 19% to the mean annual emissions of 1.2 GtC/yr. These factors, in combination with amplification effects under elevated CO2, contribute substantially to future emissions from LUC in all RCPs.

Imminent ocean acidification in the Arctic projected with the NCAR global coupled carbon cycle-climate model

Ocean acidification from the uptake of anthropogenic carbon is simulated for the industrial period and IPCC SRES emission scenarios A2 and B1 with a global coupled carbon cycle-climate model. Earlier studies identified seawater saturation state with respect to aragonite, a mineral phase of calcium carbonate, as a key variable governing impacts on corals and other shell-forming organisms. Globally in the A2 scenario, water saturated by more than 300%, considered suitable for coral growth, vanishes by 2070 AD (CO2≈630 ppm), and the ocean volume fraction occupied by saturated water decreases from 42% to 25% over this century. The largest simulated pH changes worldwide occur in Arctic surface waters, where hydrogen ion concentration increases by up to 185% (ΔpH=−0.45). Projected climate change amplifies the decrease in Arctic surface mean saturation and pH by more than 20%, mainly due to freshening and increased carbon uptake in response to sea ice retreat. Modeled saturation compares well with observation-based estimates along an Arctic transect and simulated changes have been corrected for remaining model-data differences in this region. Aragonite undersaturation in Arctic surface waters is projected to occur locally within a decade and to become more widespread as atmospheric CO2 continues to grow. The results imply that surface waters in the Arctic Ocean will become corrosive to aragonite, with potentially large implications for the marine ecosystem, if anthropogenic carbon emissions are not reduced and atmospheric CO2 not kept below 450 ppm.

Oxygen and indicators of stress for marine life in multi-model global warming projections

Decadal-to-century scale trends for a range of marine environmental variables in the upper mesopelagic layer (UML, 100–600 m) are investigated using results from seven Earth System Models forced by a high greenhouse gas emission scenario. The models as a class represent the observation-based distribution of oxygen (O2) and carbon dioxide (CO2), albeit major mismatches between observation-based and simulated values remain for individual models. By year 2100 all models project an increase in SST between 2 °C and 3 °C, and a decrease in the pH and in the saturation state of water with respect to calcium carbonate minerals in the UML. A decrease in the total ocean inventory of dissolved oxygen by 2% to 4% is projected by the range of models. Projected O2 changes in the UML show a complex pattern with both increasing and decreasing trends reflecting the subtle balance of different competing factors such as circulation, production, remineralization, and temperature changes. Projected changes in the total volume of hypoxic and suboxic waters remain relatively small in all models. A widespread increase of CO2 in the UML is projected. The median of the CO2 distribution between 100 and 600m shifts from 0.1–0.2 mol m−3 in year 1990 to 0.2–0.4 mol m−3 in year 2100, primarily as a result of the invasion of anthropogenic carbon from the atmosphere. The co-occurrence of changes in a range of environmental variables indicates the need to further investigate their synergistic impacts on marine ecosystems and Earth System feedbacks.