Climate-induced interannual variability of marine primary and export production in three global coupled climate carbon cycle models

Fully coupled climate carbon cycle models are sophisticated tools that are used to predict future climate change and its impact on the land and ocean carbon cycles. These models should be able to adequately represent natural variability, requiring model validation by observations. The present study focuses on the ocean carbon cycle component, in particular the spatial and temporal variability in net primary productivity (PP) and export production (EP) of particulate organic carbon (POC). Results from three coupled climate carbon cycle models (IPSL, MPIM, NCAR) are compared with observation-based estimates derived from satellite measurements of ocean colour and results from inverse modelling (data assimilation). Satellite observations of ocean colour have shown that temporal variability of PP on the global scale is largely dominated by the permanently stratified, low-latitude ocean (Behrenfeld et al., 2006) with stronger stratification (higher sea surface temperature; SST) being associated with negative PP anomalies. Results from all three coupled models confirm the role of the low-latitude, permanently stratified ocean for anomalies in globally integrated PP, but only one model (IPSL) also reproduces the inverse relationship between stratification (SST) and PP. An adequate representation of iron and macronutrient co-limitation of phytoplankton growth in the tropical ocean has shown to be the crucial mechanism determining the capability of the models to reproduce observed interactions between climate and PP.

Summer temperature in the eastern part of southern South America: its variability in the twentieth century and a teleconnection with Oceania

The 1907–2001 summer-to-summer surface air temperature variability in the eastern part of southern South America (SSA, partly including Patagonia) is analysed. Based on records from instruments located next to the Atlantic Ocean (36°S–55°S), we define indices for the interannual and interdecadal timescales. The main interdecadal mode reflects the late-1970s cold-to-warm climate shift in the region and a warm-to-cold transition during early 1930s. Although it has been in phase with the Pacific Decadal Oscillation (PDO) index since the 1960s, they diverged in the preceding decades. The main interannual variability index exhibits high spectral power at ~3.4 years and is representative of temperature variability in a broad area in the southern half of the continent. Eleven-years running correlation coefficients between this index and December-to-February (DJF) Niño3.4 show significant decadal fluctuations, out-of-phase with the running correlation with a DJF index of the Southern Annular Mode. The main interannual variability index is associated with a barotropic wavetrain-like pattern extending over the South Pacific from Oceania to SSA. During warm (cold) summers in SSA, significant anticyclonic (cyclonic) anomalies tend to predominate over eastern Australia, to the north of the Ross Sea, and to the east of SSA, whereas anomalous cyclonic (anticyclonic) circulation is observed over New Zealand and west of SSA. This teleconnection links warm (cold) SSA anomalies with dry (wet) summers in eastern Australia. The covariability seems to be influenced by the characteristics of tropical forcing; indeed, a disruption has been observed since late 1970s, presumably due to the PDO warm phase.

The dominant role of semi-arid ecosystems in the trend and variability of the land CO2 sink

The growth rate of atmospheric carbondioxide(CO2) concentrations since industrialization is characterized by large interannual variability, mostly resulting from variability in CO 2 uptake by terrestrial ecosystems (typically termed carbon sink). However, the contributions of regional ecosystems to that variability are not well known. Using an ensemble of ecosystem and land-surface models and an empirical observation-based product of global gross primary production, we show that the mean sink, trend, and interannual variability in CO2 uptake by terrestrial ecosystems are dominated by distinct biogeographic regions. Whereas the mean sink is dominated by highly productive lands (mainly tropical forests), the trend and interannual variability of the sink are dominated by semi-arid ecosystems whose carbon balance is strongly associated with circulation-driven variations in both precipitation and temperature.

Projected 21st century decrease in marine productivity: a multi-model analysis

Changes in marine net primary productivity (PP) and export of particulate organic carbon (EP) are projected over the 21st century with four global coupled carbon cycle-climate models. These include representations of marine ecosystems and the carbon cycle of different structure and complexity. All four models show a decrease in global mean PP and EP between 2 and 20% by 2100 relative to preindustrial conditions, for the SRES A2 emission scenario. Two different regimes for productivity changes are consistently identified in all models. The first chain of mechanisms is dominant in the low- and mid-latitude ocean and in the North Atlantic: reduced input of macro-nutrients into the euphotic zone related to enhanced stratification, reduced mixed layer depth, and slowed circulation causes a decrease in macro-nutrient concentrations and in PP and EP. The second regime is projected for parts of the Southern Ocean: an alleviation of light and/or temperature limitation leads to an increase in PP and EP as productivity is fueled by a sustained nutrient input. A region of disagreement among the models is the Arctic, where three models project an increase in PP while one model projects a decrease. Projected changes in seasonal and interannual variability are modest in most regions. Regional model skill metrics are proposed to generate multi-model mean fields that show an improved skill in representing observation-based estimates compared to a simple multi-model average. Model results are compared to recent productivity projections with three different algorithms, usually applied to infer net primary production from satellite observations.

Predictability of biomass burning in response to climate changes

Climate is an important control on biomass burning, but the sensitivity of fire to changes in temperature and moisture balance has not been quantified. We analyze sedimentary charcoal records to show that the changes in fire regime over the past 21,000 yrs are predictable from changes in regional climates. Analyses of paleo- fire data show that fire increases monotonically with changes in temperature and peaks at intermediate moisture levels, and that temperature is quantitatively the most important driver of changes in biomass burning over the past 21,000 yrs. Given that a similar relationship between climate drivers and fire emerges from analyses of the interannual variability in biomass burning shown by remote-sensing observations of month-by-month burnt area between 1996 and 2008, our results signal a serious cause for concern in the face of continuing global warming.

Effects of oligotrophication on primary production in peri-alpine lakes

During the second half of the 20th century untreated sewage released from housing and industry into natural waters led to a degradation of many freshwater lakes and reservoirs worldwide. In order to mitigate eutrophication, wastewater treatment plants, including Fe-induced phosphorus precipitation, were implemented throughout the industrialized world, leading to reoligotrophication in many freshwater lakes. To understand and assess the effects of reoligotrophication on primary productivity, we analyzed 28 years of 14C assimilation rates, as well as other biotic and abiotic parameters, such as global radiation, nutrient concentrations and plankton densities in peri-alpine Lake Lucerne, Switzerland. Using a simple productivity-light relationship, we estimated continuous primary production and discussed the relation between productivity and observed limnological parameters. Furthermore, we assessed the uncertainty of our modeling approach based on monthly 14C assimilation measurements using Monte Carlo simulations. Results confirm that monthly sampling of productivity is sufficient for identifying long-term trends in productivity and that conservation management has successfully improved water quality during the past three decades via reducing nutrients and primary production in the lake. However, even though nutrient concentrations have remained constant in recent years, annual primary production varies significantly from year to year. Despite the fact that nutrient concentrations have decreased by more than an order of magnitude, primary production has decreased only slightly. These results suggest that primary production correlates well to nutrients availability but meteorological conditions lead to interannual variability regardless of the trophic status of the lake. Accordingly, in oligotrophic freshwaters meteorological forcing may reduce productivity impacting on the entire food chain of the ecosystem.

IMILAST: A community effort to intercompare extratropical cyclone detection and tracking algorithms: assessing method-related uncertainties

The variability of results from different automated methods of detection and tracking of extratropical cyclones is assessed in order to identify uncertainties related to the choice of method. Fifteen international teams applied their own algorithms to the same dataset - the period 1989-2009 of interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERAInterim) data. This experiment is part of the community project Intercomparison of Mid Latitude Storm Diagnostics (IMILAST; see www.proclim.ch/imilast/index.html). The spread of results for cyclone frequency, intensity, life cycle, and track location is presented to illustrate the impact of using different methods. Globally, methods agree well for geographical distribution in large oceanic regions, interannual variability of cyclone numbers, geographical patterns of strong trends, and distribution shape for many life cycle characteristics. In contrast, the largest disparities exist for the total numbers of cyclones, the detection of weak cyclones, and distribution in some densely populated regions. Consistency between methods is better for strong cyclones than for shallow ones. Two case studies of relatively large, intense cyclones reveal that the identification of the most intense part of the life cycle of these events is robust between methods, but considerable differences exist during the development and the dissolution phases.

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.

Land-atmosphere coupling over North America in CRCM5

Land-atmosphere coupling and its impact on extreme precipitation and temperature events over North America are studied using the fifth generation of the Canadian Regional Climate Model (CRCM5). To this effect, two 30 year long simulations, spanning the 1981–2010 period, with and without land-atmosphere coupling, have been performed with CRCM5, driven by the European Centre for Medium-Range Weather Forecasts reanalysis at the boundaries. In the coupled simulation, the soil moisture interacts freely with the atmosphere at each time step, while in the uncoupled simulation, soil moisture is replaced with its climatological value computed from the coupled simulation, thus suppressing the soil moisture-atmosphere interactions. Analyses of the coupled and uncoupled simulations, for the summer period, show strong soil moisture-temperature coupling over the Great Plains, consistent with previous studies. The maxima of soil moisture-precipitation coupling is more spread out and covers the semiarid regions of the western U.S. and parts of the Great Plains. However, the strength of soil moisture-precipitation coupling is found to be generally weaker than that of soil moisture-temperature coupling. The study clearly indicates that land-atmosphere coupling increases the interannual variability of the seasonal mean daily maximum temperature in the Great Plains. Land-atmosphere coupling is found to significantly modulate selected temperature extremes such as the number of hot days, frequency, and maximum duration of hot spells over the Great Plains. Results also suggest additional hot spots, where soil moisture modulates extreme events. These hot spots are located in the southeast U.S. for the hot days/hot spells and in the semiarid regions of the western U.S. for extreme wet spells. This study thus demonstrates that climatologically wet/dry regions can become hot spots of land-atmosphere coupling when the soil moisture decreases/increases to an intermediate transitional level where evapotranspiration becomes moisture sensitive and large enough to affect the climate.

Trends in mean and extreme precipitation in the Mount Kenya region from observations and reanalyses

Land and water management in semi-arid regions requires detailed information on precipitation distribution, including extremes, and changes therein. Such information is often lacking. This paper describes statistics of mean and extreme precipitation in a unique data set from the Mount Kenya region, encompassing around 50 stations with at least 30 years of data. We describe the data set, including quality control procedures and statistical break detection. Trends in mean precipitation and extreme indices calculated from these data for individual rainy seasons are compared with corresponding trends in reanalysis products. From 1979 to 2011, mean precipitation decreased at 75% of the stations during the ‘long rains’ (March to May) and increased at 70% of the stations during the ‘short rains’ (October to December). Corresponding trends are found in the number of heavy precipitation days, and maximum of consecutive 5-day precipitation. Conversely, an increase in consecutive dry days within both main rainy seasons is found. However, trends are only statistically significant in very few cases. Reanalysis data sets agree with observations with respect to interannual variability, while correlations are considerably lower for monthly deviations (ratios) from the mean annual cycle. While some products well reproduce the rainfall climatology and some the spatial trend pattern, no product reproduces both.

Contrasting interannual and multidecadal NAO variability

Decadal and longer timescale variability in the winter North Atlantic Oscillation (NAO) has considerable impact on regional climate, yet it remains unclear what fraction of this variability is potentially predictable. This study takes a new approach to this question by demonstrating clear physical differences between NAO variability on interannual-decadal (<30 year) and multidecadal (>30 year) timescales. It is shown that on the shorter timescale the NAO is dominated by variations in the latitude of the North Atlantic jet and storm track, whereas on the longer timescale it represents changes in their strength instead. NAO variability on the two timescales is associated with different dynamical behaviour in terms of eddy-mean flow interaction, Rossby wave breaking and blocking. The two timescales also exhibit different regional impacts on temperature and precipitation and different relationships to sea surface temperatures. These results are derived from linear regression analysis of the Twentieth Century and NCEP-NCAR reanalyses and of a high-resolution HiGEM General Circulation Model control simulation, with additional analysis of a long sea level pressure reconstruction. Evidence is presented for an influence of the ocean circulation on the longer timescale variability of the NAO, which is particularly clear in the model data. As well as providing new evidence of potential predictability, these findings are shown to have implications for the reconstruction and interpretation of long climate records.

Decadal-scale variability in hazardous winds in northern Switzerland since end of the 19th century

The decadal-scale variability in winter hazardous winds in northern Switzerland from 1871 to present is investigated in the Twentieth Century Reanalysis (20CR). Independent wind speed measurements taken at Zurich climate station show that the interannual and decadal variability in hazardous winds in northern Switzerland is realistically represented in the 20CR. Both time series exhibit pronounced decadal-scale variability with periods between approximately 36 and 47 years. At these periodicities, the hazardous wind variability in northern Switzerland is positively correlated with the variability in the North Atlantic Oscillation, however the strength and statistical significance of their co-variability varies over time.