Robust dynamics of Amazon dieback to climate change with perturbed ecosystem model parameters

Climate change science is increasingly concerned with methods for managing and integrating sources of uncertainty from emission storylines, climate model projections, and ecosystem model parameterizations. In tropical ecosystems, regional climate projections and modeled ecosystem responses vary greatly, leading to a significant source of uncertainty in global biogeochemical accounting and possible future climate feedbacks. Here, we combine an ensemble of IPCC-AR4 climate change projections for the Amazon Basin (eight general circulation models) with alternative ecosystem parameter sets for the dynamic global vegetation model, LPJmL. We evaluate LPJmL simulations of carbon stocks and fluxes against flux tower and aboveground biomass datasets for individual sites and the entire basin. Variability in LPJmL model sensitivity to future climate change is primarily related to light and water limitations through biochemical and water-balance-related parameters. Temperature-dependent parameters related to plant respiration and photosynthesis appear to be less important than vegetation dynamics (and their parameters) for determining the magnitude of ecosystem response to climate change. Variance partitioning approaches reveal that relationships between uncertainty from ecosystem dynamics and climate projections are dependent on geographic location and the targeted ecosystem process. Parameter uncertainty from the LPJmL model does not affect the trajectory of ecosystem response for a given climate change scenario and the primary source of uncertainty for Amazon 'dieback' results from the uncertainty among climate projections. Our approach for describing uncertainty is applicable for informing and prioritizing policy options related to mitigation and adaptation where long-term investments are required.

Quantifying carbon-cycle feedbacks

Perturbations to the carbon cycle could constitute large feedbacks on future changes in atmospheric CO2 concentration and climate. This paper demonstrates how carbon cycle feedback can be expressed in formally similar ways to climate feedback, and thus compares their magnitudes. The carbon cycle gives rise to two climate feedback terms: the concentration–carbon feedback, resulting from the uptake of carbon by land and ocean as a biogeochemical response to the atmospheric CO2 concentration, and the climate–carbon feedback, resulting from the effect of climate change on carbon fluxes. In the earth system models of the Coupled Climate–Carbon Cycle Model Intercomparison Project (C4MIP), climate–carbon feedback on warming is positive and of a similar size to the cloud feedback. The concentration–carbon feedback is negative; it has generally received less attention in the literature, but in magnitude it is 4 times larger than the climate–carbon feedback and more uncertain. The concentration–carbon feedback is the dominant uncertainty in the allowable CO2 emissions that are consistent with a given CO2 concentration scenario. In modeling the climate response to a scenario of CO2 emissions, the net carbon cycle feedback is of comparable size and uncertainty to the noncarbon–climate response. To quantify simulated carbon cycle feedbacks satisfactorily, a radiatively coupled experiment is needed, in addition to the fully coupled and biogeochemically coupled experiments, which are referred to as coupled and uncoupled in C4MIP. The concentration–carbon and climate–carbon feedbacks do not combine linearly, and the concentration–carbon feedback is dependent on scenario and time.

Modelling monsoons: understanding and predicting current and future behaviour

The global monsoon system is so varied and complex that understanding and predicting its diverse behaviour remains a challenge that will occupy modellers for many years to come. Despite the difficult task ahead, an improved monsoon modelling capability has been realized through the inclusion of more detailed physics of the climate system and higher resolution in our numerical models. Perhaps the most crucial improvement to date has been the development of coupled ocean-atmosphere models. From subseasonal to interdecadal time scales, only through the inclusion of air-sea interaction can the proper phasing and teleconnections of convection be attained with respect to sea surface temperature variations. Even then, the response to slow variations in remote forcings (e.g., El Niño—Southern Oscillation) does not result in a robust solution, as there are a host of competing modes of variability that must be represented, including those that appear to be chaotic. Understanding the links between monsoons and land surface processes is not as mature as that explored regarding air-sea interactions. A land surface forcing signal appears to dominate the onset of wet season rainfall over the North American monsoon region, though the relative role of ocean versus land forcing remains a topic of investigation in all the monsoon systems. Also, improved forecasts have been made during periods in which additional sounding observations are available for data assimilation. Thus, there is untapped predictability that can only be attained through the development of a more comprehensive observing system for all monsoon regions. Additionally, improved parameterizations - for example, of convection, cloud, radiation, and boundary layer schemes as well as land surface processes - are essential to realize the full potential of monsoon predictability. A more comprehensive assessment is needed of the impact of black carbon aerosols, which may modulate that of other anthropogenic greenhouse gases. Dynamical considerations require ever increased horizontal resolution (probably to 0.5 degree or higher) in order to resolve many monsoon features including, but not limited to, the Mei-Yu/Baiu sudden onset and withdrawal, low-level jet orientation and variability, and orographic forced rainfall. Under anthropogenic climate change many competing factors complicate making robust projections of monsoon changes. Absent aerosol effects, increased land-sea temperature contrast suggests strengthened monsoon circulation due to climate change. However, increased aerosol emissions will reflect more solar radiation back to space, which may temper or even reduce the strength of monsoon circulations compared to the present day. Precipitation may behave independently from the circulation under warming conditions in which an increased atmospheric moisture loading, based purely on thermodynamic considerations, could result in increased monsoon rainfall under climate change. The challenge to improve model parameterizations and include more complex processes and feedbacks pushes computing resources to their limit, thus requiring continuous upgrades of computational infrastructure to ensure progress in understanding and predicting current and future behaviour of monsoons.

Using idealized snow forcing to test teleconnections with the Indian summer monsoon in the Hadley Centre GCM

Anomalous heavy snow during winter or spring has long been regarded as a possible precursor of deficient Indian monsoon rainfall during the subsequent summer. However previous work in this field is inconclusive, in terms of the mechanism that communicates snow anomalies to the monsoon summer, and even the region from which snow has the most impact. In this study we explore these issues in coupled and atmosphere-only versions of the Hadley Centre model. A 1050-year control integration of the HadCM3 coupled model, which well represents the seasonal cycle of snow cover over the Eurasian continent, is analysed and shows evidence for weakened monsoons being preceded by strong snow forcing (in the absence of ENSO) over either the Himalaya/Tibetan Plateau or north/west Eurasia regions. However, empirical orthogonal function (EOF) analysis of springtime interannual variability in snow depth shows the leading mode to have opposite signs between these two regions, suggesting that competing mechanisms may be possible. To determine the dominant region, ensemble integrations are carried out using HadAM3, the atmospheric component of HadCM3, and a variety of anomalous snow forcing initial conditions obtained from the control integration of the coupled model. Forcings are applied during spring in separate experiments over the Himalaya/Tibetan Plateau and north/west Eurasia regions, in conjunction with climatological SSTs in order to avoid the direct effects of ENSO. With the aid of idealized forcing conditions in sensitivity tests, we demonstrate that forcing from the Himalaya region is dominant in this model via a Blanford-type mechanism involving reduced surface sensible heat and longwave fluxes, reduced heating of the troposphere over the Tibetan Plateau and consequently a reduced meridional tropospheric temperature gradient which weakens the monsoon during early summer. Snow albedo is shown to be key to the mechanism, explaining around 50% of the perturbation in sensible heating over the Tibetan Plateau, and accounting for the majority of cooling through the troposphere.

The sea level conundrum: insights from paleo studies

Empirical Constraints on Future Sea Level Rise; Bern, Switzerland, 25–29 August 2008; Eustatic sea level (ESL) rise during the 21st century is perhaps the greatest threat from climate change, but its magnitude is contested. Geological records identify examples of nonlinear ice sheet response to climate forcing, suggesting a strategy for refining estimates of 21st-century sea level change. In August 2008, Past Global Changes (PAGES), International Marine Past Global Change Study (IMAGES), and the University of Bern cosponsored a workshop to address this possibility. The workshop highlighted several ways that paleoceanography studies can place limits on future sea level rise, and these are enlarged upon here.

Thresholds for irreversible decline of the Greenland ice sheet

The Greenland ice sheet will decline in volume in a warmer climate. If a sufficiently warm climate is maintained for a few thousand years, the ice sheet will be completely melted. This raises the question of whether the decline would be reversible: would the ice sheet regrow if the climate cooled down? To address this question, we conduct a number of experiments using a climate model and a high-resolution ice-sheet model. The experiments are initialised with ice sheet states obtained from various points during its decline as simulated in a high-CO2 scenario, and they are then forced with a climate simulated for pre-industrial greenhouse gas concentrations, to determine the possible trajectories of subsequent ice sheet evolution. These trajectories are not the reverse of the trajectory during decline. They converge on three different steady states. The original ice-sheet volume can be regained only if the volume has not fallen below a threshold of irreversibility, which lies between 80 and 90% of the original value. Depending on the degree of warming and the sensitivity of the climate and the ice-sheet, this point of no return could be reached within a few hundred years, sooner than CO2 and global climate could revert to a pre-industrial state, and in that case global sea level rise of at least 1.3 m would be irreversible. An even larger irreversible change to sea level rise of 5 m may occur if ice sheet volume drops below half of its current size. The set of steady states depends on the CO2 concentration. Since we expect the results to be quantitatively affected by resolution and other aspects of model formulation, we would encourage similar investigations with other models.

Lava production at Soufriere Hills Volcano, Montserrat: 1995-2009

[1] We estimate that about 1 km3 of andesitic lava has been produced at Soufrière Hills Volcano, Montserrat from 1995 to 2009. There were three major episodes of extrusion, each lasting about 2 to 3.5 years and producing about 280 to 340 M m3 of lava, and one minor episode. Our estimates account for the dense rock equivalent volumetric contributions from the core and talus components of the lava dome, pyroclastic flow deposits and air-fall deposits. By 2005 at least two thirds of the erupted mass has already entered the sea. The average lava flux across the major extrusion episodes has been 3–5 m3s−1, with short-period (10–15 days) pulses up to 10–20 m3s−1. The first and third episodes of extrusion show similar flux histories suggesting similar behaviour of the system ten years apart. Waning flux towards the end of each episode may be caused by declining overpressure in the magma reservoir.

Modelling the inorganic nitrogen behaviour in a small Mediterranean forested catchment, Fuirosos (Catalonia)

The aim of this work was to couple a nitrogen (N) sub-model to already existent hydrological lumped (LU4-N) and semi-distributed (LU4-R-N and SD4-R-N) conceptual models, to improve our understanding of the factors and processes controlling nitrogen cycling and losses in Mediterranean catchments. The N model adopted provides a simplified conceptualization of the soil nitrogen cycle considering mineralization, nitrification, immobilization, denitrification, plant uptake, and ammonium adsorption/desorption. It also includes nitrification and denitrification in the shallow perched aquifer. We included a soil moisture threshold for all the considered soil biological processes. The results suggested that all the nitrogen processes were highly influenced by the rain episodes and that soil microbial processes occurred in pulses stimulated by soil moisture increasing after rain. Our simulation highlighted the riparian zone as a possible source of nitrate, especially after the summer drought period, but it can also act as an important sink of nitrate due to denitrification, in particular during the wettest period of the year. The riparian zone was a key element to simulate the catchment nitrate behaviour. The lumped LU4-N model (which does not include the riparian zone) could not be validated, while both the semi-distributed LU4-R-N and SD4-R-N model (which include the riparian zone) gave satisfactory results for the calibration process and acceptable results for the temporal validation process.

Tracer correlations in the northern high latitude lowermost stratosphere: Influence of cross‐tropopause mass exchange

We present an analysis of trace gas correlations in the lowermost stratosphere. In‐situ aircraft measurements of CO, N2O, NOy and O3, obtained during the STREAM 1997 winter campaign, have been used to investigate the role of cross‐tropopause mass exchange on tracer‐tracer relations. At altitudes several kilometers above the local tropopause, undisturbed stratospheric air was found with NOy/NOy * ratios close to unity, NOy/O3 about 0.003–0.006 and CO mixing ratios as low as 20 ppbv (NOy * is a proxy for total reactive nitrogen derived from NOy–N2O relations measured in the stratosphere). Mixing of tropospheric air into the lowermost stratosphere has been identified by enhanced ratios of NOy/NOy * and NOy/O3, and from scatter plots of CO versus O3. The enhanced NOy/O3 ratio in the lowermost stratospheric mixing zone points to a reduced efficiency of O3 formation from aircraft NOx emissions.

The three-dimensional vortical nature of atmospheric and oceanic turbulent flows

Using a novel numerical method at unprecedented resolution, we demonstrate that structures of small to intermediate scale in rotating, stratified flows are intrinsically three-dimensional. Such flows are characterized by vortices (spinning volumes of fluid), regions of large vorticity gradients, and filamentary structures at all scales. It is found that such structures have predominantly three-dimensional dynamics below a horizontal scale LLR, where LR is the so-called Rossby radius of deformation, equal to the characteristic vertical scale of the fluid H divided by the ratio of the rotational and buoyancy frequencies f/N. The breakdown of two-dimensional dynamics at these scales is attributed to the so-called "tall-column instability" [D. G. Dritschel and M. de la Torre Juárez, J. Fluid. Mech. 328, 129 (1996)], which is active on columnar vortices that are tall after scaling by f/N, or, equivalently, that are narrow compared with LR. Moreover, this instability eventually leads to a simple relationship between typical vertical and horizontal scales: for each vertical wave number (apart from the vertically averaged, barotropic component of the flow) the average horizontal wave number is equal to f/N times the vertical wave number. The practical implication is that three-dimensional modeling is essential to capture the behavior of rotating, stratified fluids. Two-dimensional models are not valid for scales below LR. ©1999 American Institute of Physics.

Glacier shrinkage and climatic change in the Russian Altai from the mid-20th century: An assessment using remote sensing and PRECIS regional climate model

This paper examines changes in the surface area of glaciers in the North and South Chuya Ridges, Altai Mountains in 1952-2004 and their links with regional climatic variations. The glacier surface areas for 2004 were derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. Data from the World Glacier Inventory (WGI)dating to 1952 and aerial photographs from 1952 were used to estimate the changes. 256 glaciers with a combined area of 253±5.1 km2 have been identified in the region in 2004. Estimation of changes in extent of 126 glaciers with the individual areas not less than 0.5 km2 in 1952 revealed a 19.7±5.8% reduction. The observed glacier retreat is primarily driven by an increase in summer temperatures since the 1980s when air temperatures were increasing at a rate of 0.10 - 0.13oC a-1 at the glacier tongue elevation. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in the Altai in 2071-2100 by 6-7oC and 3-5oC respectively in comparison with 1961-1990 while annual precipitation will increase by 15% and 5%. The length of the ablation season will extend from June-August to the late April – early October. The projected increases in precipitation will not compensate for the projected warming and glaciers will continue to retreat in the 21st century under both B2 and A2 scenarios.

The role of eddies in driving the tropospheric response to stratospheric heating perturbations

A simplified general circulation model has been used to investigate the chain of causality whereby changes in tropospheric circulation and temperature are produced in response to stratospheric heating perturbations. Spinup ensemble experiments have been performed to examine the evolution of the tropospheric circulation in response to such perturbations. The primary aim of these experiments is to investigate the possible mechanisms whereby a tropospheric response to changing solar activity over the 11-yr solar cycle could be produced in response to heating of the equatorial lower stratosphere. This study therefore focuses on a stratospheric heating perturbation in which the heating is largest in the tropics. For comparison, experiments are also performed in which the stratosphere is heated uniformly at all latitudes and in which it is heated preferentially in the polar region. Thus, the mechanisms discussed have a wider relevance for the impact of stratospheric perturbations on the troposphere. The results demonstrate the importance of changing eddy momentum fluxes in driving the tropospheric response. This is confirmed by the lack of a similar response in a zonally symmetric model with fixed eddy forcing. Furthermore, it is apparent that feedback between the tropospheric eddy fluxes and tropospheric circulation changes is required to produce the full model response. The quasigeostrophic index of refraction is used to diagnose the cause of the changes in eddy behavior. It is demonstrated that the latitudinal extent of stratospheric heating is important in determining the direction of displacement of the tropospheric jet and storm track.