Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP)

Global wetlands are believed to be climate sensitive, and are the largest natural emitters of methane (CH4). Increased wetland CH4 emissions could act as a positive feedback to future warming. The Wetland and Wetland CH4 Inter-comparison of Models Project (WETCHIMP) investigated our present ability to simulate large-scale wetland characteristics and corresponding CH4 emissions. To ensure inter-comparability, we used a common experimental protocol driving all models with the same climate and carbon dioxide (CO2) forcing datasets. The WETCHIMP experiments were conducted for model equilibrium states as well as transient simulations covering the last century. Sensitivity experiments investigated model response to changes in selected forcing inputs (precipitation, temperature, and atmospheric CO2 concentration). Ten models participated, covering the spectrum from simple to relatively complex, including models tailored either for regional or global simulations. The models also varied in methods to calculate wetland size and location, with some models simulating wetland area prognostically, while other models relied on remotely sensed inundation datasets, or an approach intermediate between the two. Four major conclusions emerged from the project. First, the suite of models demonstrate extensive disagreement in their simulations of wetland areal extent and CH4 emissions, in both space and time. Simple metrics of wetland area, such as the latitudinal gradient, show large variability, principally between models that use inundation dataset information and those that independently determine wetland area. Agreement between the models improves for zonally summed CH4 emissions, but large variation between the models remains. For annual global CH4 emissions, the models vary by ±40% of the all-model mean (190 Tg CH4 yr−1). Second, all models show a strong positive response to increased atmospheric CO2 concentrations (857 ppm) in both CH4 emissions and wetland area. In response to increasing global temperatures (+3.4 °C globally spatially uniform), on average, the models decreased wetland area and CH4 fluxes, primarily in the tropics, but the magnitude and sign of the response varied greatly. Models were least sensitive to increased global precipitation (+3.9 % globally spatially uniform) with a consistent small positive response in CH4 fluxes and wetland area. Results from the 20th century transient simulation show that interactions between climate forcings could have strong non-linear effects. Third, we presently do not have sufficient wetland methane observation datasets adequate to evaluate model fluxes at a spatial scale comparable to model grid cells (commonly 0.5°). This limitation severely restricts our ability to model global wetland CH4 emissions with confidence. Our simulated wetland extents are also difficult to evaluate due to extensive disagreements between wetland mapping and remotely sensed inundation datasets. Fourth, the large range in predicted CH4 emission rates leads to the conclusion that there is both substantial parameter and structural uncertainty in large-scale CH4 emission models, even after uncertainties in wetland areas are accounted for.

Sea-surface temperature records of Termination 1 in the Gulf of California: Challenges for seasonal and interannual analogues of tropical Pacific climate change

Centennial-scale records of sea-surface temperature and opal composition spanning the Last Glacial Maximum and Termination 1 (circa 25–6 ka) are presented here from Guaymas Basin in the Gulf of California. Through the application of two organic geochemistry proxies, the U37K′ index and the TEX86H index, we present evidence for rapid, stepped changes in temperatures during deglaciation. These occur in both temperature proxies at 13 ka (∼3°C increase in 270 years), 10.0 ka (∼2°C decrease over ∼250 years) and at 8.2 ka (3°C increase in <200 years). An additional rapid warming step is also observed in TEX86H at 11.5 ka. In comparing the two temperature proxies and opal content, we consider the potential for upwelling intensity to be recorded and link this millennial-scale variability to shifting Intertropical Convergence Zone position and variations in the strength of the Subtropical High. The onset of the deglacial warming from 17 to 18 ka is comparable to a “southern hemisphere” signal, although the opal record mimics the ice-rafting events of the north Atlantic (Heinrich events). Neither the modern seasonal cycle nor El Niño/Southern Oscillation patterns provide valid analogues for the trends we observe in comparison with other regional records. Fully coupled climate model simulations confirm this result, and in combination we question whether the seasonal or interannual climate variations of the modern climate are valid analogues for the glacial and deglacial tropical Pacific.

Upland streamwater nitrate dynamics across decadal to sub-daily timescales: a case study of Plynlimon, Wales

Streamwater nitrate dynamics in the River Hafren, Plynlimon, mid-Wales were investigated over decadal to sub-daily timescales using a range of statistical techniques. Long-term data were derived from weekly grab samples (1984–2010) and high-frequency data from 7-hourly samples (2007–2009) both measured at two sites: a headwater stream draining moorland and a downstream site below plantation forest. This study is one of the first to analyse upland streamwater nitrate dynamics across such a wide range of timescales and report on the principal mechanisms identified. The data analysis provided no clear evidence that the long-term decline in streamwater nitrate concentrations was related to a decline in atmospheric deposition alone, because nitrogen deposition first increased and then decreased during the study period. Increased streamwater temperature and denitrification may also have contributed to the decline in stream nitrate concentrations, the former through increased N uptake rates and the latter resultant from increased dissolved organic carbon concentrations. Strong seasonal cycles, with concentration minimums in the summer, were driven by seasonal flow minimums and seasonal biological activity enhancing nitrate uptake. Complex diurnal dynamics were observed, with seasonal changes in phase and amplitude of the cycling, and the diurnal dynamics were variable along the river. At the moorland site, a regular daily cycle, with minimum concentrations in the early afternoon, corresponding with peak air temperatures, indicated the importance of instream biological processing. At the downstream site, the diurnal dynamics were a composite signal, resultant from advection, dispersion and nitrate processing in the soils of the lower catchment. The diurnal streamwater nitrate dynamics were also affected by drought conditions. Enhanced diurnal cycling in Spring 2007 was attributed to increased nitrate availability in the post-drought period as well as low flow rates and high temperatures over this period. The combination of high-frequency short-term measurements and long-term monitoring provides a powerful tool for increasing understanding of the controls of element fluxes and concentrations in surface waters.

The influence of vegetation, fire spread and fire behaviour on biomass burning and trace gas emissions: results from a process-based model

A process-based fire regime model (SPITFIRE) has been developed, coupled with ecosystem dynamics in the LPJ Dynamic Global Vegetation Model, and used to explore fire regimes and the current impact of fire on the terrestrial carbon cycle and associated emissions of trace atmospheric constituents. The model estimates an average release of 2.24 Pg C yr−1 as CO2 from biomass burning during the 1980s and 1990s. Comparison with observed active fire counts shows that the model reproduces where fire occurs and can mimic broad geographic patterns in the peak fire season, although the predicted peak is 1–2 months late in some regions. Modelled fire season length is generally overestimated by about one month, but shows a realistic pattern of differences among biomes. Comparisons with remotely sensed burnt-area products indicate that the model reproduces broad geographic patterns of annual fractional burnt area over most regions, including the boreal forest, although interannual variability in the boreal zone is underestimated.

Global patterns of vegetation response to millennial-scale variability and rapid climate change during the last glacial period

Ninety-four sites worldwide have sufficient resolution and dating to document the impact of millennial-scale climate variability on vegetation and fire regimes during the last glacial period. Although Dansgaard–Oeschger (D–O) cycles all show a basically similar gross structure, they vary in the magnitude and the length of the warm and cool intervals. We illustrate the geographic patterns in the climate-induced changes in vegetation by comparing D–O 6, D–O 8 and D–O 19. There is a strong response to both D–O warming events and subsequent cooling, most marked in the northern extratropics. Pollen records from marine cores from the northern extratropics confirm that there is no lag between the change in climate and the vegetation response, within the limits of the dating resolution (50–100 years). However, the magnitude of the change in vegetation is regionally specific and is not a simple function of either the magnitude or the duration of the change in climate as registered in Greenland ice cores. Fire regimes also show an initial immediate response to climate changes, but during cooling intervals there is a slow recovery of biomass burning after the initial reduction, suggesting a secondary control through the recovery of vegetation productivity. In the extratropics, vegetation changes are largely determined by winter temperatures while in the tropics they are largely determined by changes in plant-available water. Tropical vegetation records show changes corresponding to Heinrich Stadials but the response to D–O warming events is less marked than in the northern extratropics. There are very few high-resolution records from the Southern Hemisphere extratropics, but these records also show both a vegetation and fire response to millennial-scale climate variability. It is not yet possible to determine unequivocally whether terrestrial records reflect the asynchroneity apparent in the ice-core records.

Paralog re-emergence: a novel, historically contingent mechanism in the evolution of antimicrobial resistance

Evolution of resistance to drugs and pesticides poses a serious threat to human health and agricultural production. CYP51 encodes the target site of azole fungicides, widely used clinically and in agriculture. Azole resistance can evolve due to point mutations or overexpression of CYP51, and previous studies have shown that fungicide-resistant alleles have arisen by de novo mutation. Paralogs CYP51A and CYP51B are found in filamentous ascomycetes, but CYP51A has been lost from multiple lineages. Here, we show that in the barley pathogen Rhynchosporium commune, re-emergence of CYP51A constitutes a novel mechanism for the evolution of resistance to azoles. Pyrosequencing analysis of historical barley leaf samples from a unique long-term experiment from 1892 to 2008 indicates that the majority of the R. commune population lacked CYP51A until 1985, after which the frequency of CYP51A rapidly increased. Functional analysis demonstrates that CYP51A retains the same substrate as CYP51B, but with different transcriptional regulation. Phylogenetic analyses show that the origin of CYP51A far predates azole use, and newly sequenced Rhynchosporium genomes show CYP51A persisting in the R. commune lineage rather than being regained by horizontal gene transfer; therefore, CYP51A re-emergence provides an example of adaptation to novel compounds by selection from standing genetic variation.

Coupled versus uncoupled hindcast simulations of the Madden-Julian oscillation in the Year of Tropical Convection

This study investigates the impact of a full interactive ocean on daily initialised 15 day hindcasts of the Madden-Julian oscillation (MJO), measured against a Met Office Unified Model (MetUM) atmosphere control simulation (AGCM) during a 3 month period of the Year of Tropical Convection (YOTC). Results indicated that the coupled configuration (CGCM) extends MJO predictability over that of the AGCM, by up to 3-5 days. Propagation is improved in the CGCM, which we partly attribute to a more realistic phase relationship between sea surface temperature (SST) and convection. In addition, the CGCM demonstrates skill in representing downwelling oceanic Kelvin and Rossby waves which warm SSTs along their trajectory, with the potential to feed back on the atmosphere. These results imply that an ocean model capable of simulating internal ocean waves may be required to capture the full effect of air-sea coupling for the MJO.

Association of vitamin D status with arterial blood pressure and hypertension risk: a mendelian randomisation study

BACKGROUND: Low plasma 25-hydroxyvitamin D (25[OH]D) concentration is associated with high arterial blood pressure and hypertension risk, but whether this association is causal is unknown. We used a mendelian randomisation approach to test whether 25(OH)D concentration is causally associated with blood pressure and hypertension risk. METHODS: In this mendelian randomisation study, we generated an allele score (25[OH]D synthesis score) based on variants of genes that affect 25(OH)D synthesis or substrate availability (CYP2R1 and DHCR7), which we used as a proxy for 25(OH)D concentration. We meta-analysed data for up to 108 173 individuals from 35 studies in the D-CarDia collaboration to investigate associations between the allele score and blood pressure measurements. We complemented these analyses with previously published summary statistics from the International Consortium on Blood Pressure (ICBP), the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and the Global Blood Pressure Genetics (Global BPGen) consortium. FINDINGS: In phenotypic analyses (up to n=49 363), increased 25(OH)D concentration was associated with decreased systolic blood pressure (β per 10% increase, -0·12 mm Hg, 95% CI -0·20 to -0·04; p=0·003) and reduced odds of hypertension (odds ratio [OR] 0·98, 95% CI 0·97-0·99; p=0·0003), but not with decreased diastolic blood pressure (β per 10% increase, -0·02 mm Hg, -0·08 to 0·03; p=0·37). In meta-analyses in which we combined data from D-CarDia and the ICBP (n=146 581, after exclusion of overlapping studies), each 25(OH)D-increasing allele of the synthesis score was associated with a change of -0·10 mm Hg in systolic blood pressure (-0·21 to -0·0001; p=0·0498) and a change of -0·08 mm Hg in diastolic blood pressure (-0·15 to -0·02; p=0·01). When D-CarDia and consortia data for hypertension were meta-analysed together (n=142 255), the synthesis score was associated with a reduced odds of hypertension (OR per allele, 0·98, 0·96-0·99; p=0·001). In instrumental variable analysis, each 10% increase in genetically instrumented 25(OH)D concentration was associated with a change of -0·29 mm Hg in diastolic blood pressure (-0·52 to -0·07; p=0·01), a change of -0·37 mm Hg in systolic blood pressure (-0·73 to 0·003; p=0·052), and an 8·1% decreased odds of hypertension (OR 0·92, 0·87-0·97; p=0·002). INTERPRETATION: Increased plasma concentrations of 25(OH)D might reduce the risk of hypertension. This finding warrants further investigation in an independent, similarly powered study.

Global analysis of radiative forcing from fire-induced shortwave albedo change

Land surface albedo, a key parameter to derive Earth's surface energy balance, is used in the parameterization of numerical weather prediction, climate monitoring and climate change impact assessments. Changes in albedo due to fire have not been fully investigated on a continental and global scale. The main goal of this study, therefore, is to quantify the changes in instantaneous shortwave albedo produced by biomass burning activities and their associated radiative forcing. The study relies on the MODerate-resolution Imaging Spectroradiometer (MODIS) MCD64A1 burned-area product to create an annual composite of areas affected by fire and the MCD43C2 bidirectional reflectance distribution function (BRDF) albedo snow-free product to compute a bihemispherical reflectance time series. The approximate day of burning is used to calculate the instantaneous change in shortwave albedo. Using the corresponding National Centers for Environmental Prediction (NCEP) monthly mean downward solar radiation flux at the surface, the global radiative forcing associated with fire was computed. The analysis reveals a mean decrease in shortwave albedo of −0.014 (1σ = 0.017), causing a mean positive radiative forcing of 3.99 Wm−2 (1σ = 4.89) over the 2002–20012 time period in areas affected by fire. The greatest drop in mean shortwave albedo change occurs in 2002, which corresponds to the highest total area burned (378 Mha) observed in the same year and produces the highest mean radiative forcing (4.5 Wm−2). Africa is the main contributor in terms of burned area, but forests globally give the highest radiative forcing per unit area and thus give detectable changes in shortwave albedo. The global mean radiative forcing for the whole period studied (~0.0275 Wm−2) shows that the contribution of fires to the Earth system is not insignificant.

Simulation of tree-ring widths with a model for primary production, carbon allocation, and growth

We present a simple, generic model of annual tree growth, called "T". This model accepts input from a first-principles light-use efficiency model (the "P" model). The P model provides values for gross primary production (GPP) per unit of absorbed photosynthetically active radiation (PAR). Absorbed PAR is estimated from the current leaf area. GPP is allocated to foliage, transport tissue, and fine-root production and respiration in such a way as to satisfy well-understood dimensional and functional relationships. Our approach thereby integrates two modelling approaches separately developed in the global carbon-cycle and forest-science literature. The T model can represent both ontogenetic effects (the impact of ageing) and the effects of environmental variations and trends (climate and CO2) on growth. Driven by local climate records, the model was applied to simulate ring widths during the period 1958–2006 for multiple trees of Pinus koraiensis from the Changbai Mountains in northeastern China. Each tree was initialised at its actual diameter at the time when local climate records started. The model produces realistic simulations of the interannual variability in ring width for different age cohorts (young, mature, and old). Both the simulations and observations show a significant positive response of tree-ring width to growing-season total photosynthetically active radiation (PAR0) and the ratio of actual to potential evapotranspiration (α), and a significant negative response to mean annual temperature (MAT). The slopes of the simulated and observed relationships with PAR0 and α are similar; the negative response to MAT is underestimated by the model. Comparison of simulations with fixed and changing atmospheric CO2 concentration shows that CO2 fertilisation over the past 50 years is too small to be distinguished in the ring-width data, given ontogenetic trends and interannual variability in climate.

Causal relationships vs. emergent patterns in the global controls of fire frequency

Global controls on month-by-month fractional burnt area (2000–2005) were investigated by fitting a generalised linear model (GLM) to Global Fire Emissions Database (GFED) data, with 11 predictor variables representing vegetation, climate, land use and potential ignition sources. Burnt area is shown to increase with annual net primary production (NPP), number of dry days, maximum temperature, grazing-land area, grass/shrub cover and diurnal temperature range, and to decrease with soil moisture, cropland area and population density. Lightning showed an apparent (weak) negative influence, but this disappeared when pure seasonal-cycle effects were taken into account. The model predicts observed geographic and seasonal patterns, as well as the emergent relationships seen when burnt area is plotted against each variable separately. Unimodal relationships with mean annual temperature and precipitation, population density and gross domestic product (GDP) are reproduced too, and are thus shown to be secondary consequences of correlations between different controls (e.g. high NPP with high precipitation; low NPP with low population density and GDP). These findings have major implications for the design of global fire models, as several assumptions in current models – most notably, the widely assumed dependence of fire frequency on ignition rates – are evidently incorrect.

Christian Gottlob Heyne and the changing fortunes of the commentary in the age of Altertumswissenschaft

This chapter discusses the pedagogic and scholarly priorities that informed Heyne’s commentaries on Tibullus (1755), Virgil (1767-75) and Homer (1802), as well as their initial critical reception. Like those of his teachers, Gesner and Ernesti, Heyne’s works eschew detailed textual scholarship in favour of aesthetic and historicizing appreciation of literary works as wholes. Their formal innovations – most notably the relegation of advanced philological discussions to endnotes and the inclusion of excursuses on significant historical and cultural questions – are an attempt to tailor a traditional format to the demands of an Enlightened age and the cultural-historical interests of the new Altertumswissenschaft. The chapter discusses their contrasting critical receptions in order to raise questions about the viability of Heyne’s endeavours to make a traditional medium fit new concerns.

Vertical structure and physical processes of the Madden-Julian oscillation: Linking hindcast fidelity to simulated diabatic heating and moistening

Many theories for the Madden-Julian oscillation (MJO) focus on diabatic processes, particularly the evolution of vertical heating and moistening. Poor MJO performance in weather and climate models is often blamed on biases in these processes and their interactions with the large-scale circulation. We introduce one of three components of a model-evaluation project, which aims to connect MJO fidelity in models to their representations of several physical processes, focusing on diabatic heating and moistening. This component consists of 20-day hindcasts, initialised daily during two MJO events in winter 2009-10. The 13 models exhibit a range of skill: several have accurate forecasts to 20 days' lead, while others perform similarly to statistical models (8-11 days). Models that maintain the observed MJO amplitude accurately predict propagation, but not vice versa. We find no link between hindcast fidelity and the precipitation-moisture relationship, in contrast to other recent studies. There is also no relationship between models' performance and the evolution of their diabatic-heating profiles with rain rate. A more robust association emerges between models' fidelity and net moistening: the highest-skill models show a clear transition from low-level moistening for light rainfall to mid-level moistening at moderate rainfall and upper-level moistening for heavy rainfall. The mid-level moistening, arising from both dynamics and physics, may be most important. Accurately representing many processes may be necessary, but not sufficient for capturing the MJO, which suggests that models fail to predict the MJO for a broad range of reasons and limits the possibility of finding a panacea.

Vertical structure and physical processes of the Madden-Julian oscillation: synthesis and summary

The "Vertical structure and physical processes of the Madden-Julian oscillation (MJO)" project comprises three experiments, designed to evaluate comprehensively the heating, moistening and momentum associated with tropical convection in general circulation models (GCMs). We consider here only those GCMs that performed all experiments. Some models display relatively higher or lower MJO fidelity in both initialized hindcasts and climate simulations, while others show considerable variations in fidelity between experiments. Fidelity in hindcasts and climate simulations are not meaningfully correlated. The analysis of each experiment led to the development of process-oriented diagnostics, some of which distinguished between GCMs with higher or lower fidelity in that experiment. We select the most discriminating diagnostics and apply them to data from all experiments, where possible, to determine if correlations with MJO fidelity hold across scales and GCM states. While normalized gross moist stability had a small but statistically significant correlation with MJO fidelity in climate simulations, we find no link with fidelity in medium-range hindcasts. Similarly, there is no association between timestep-to-timestep rainfall variability, identified from short hindcasts, and fidelity in medium-range hindcasts or climate simulations. Two metrics that relate precipitation to free-tropospheric moisture--the relative humidity for extreme daily precipitation, and variations in the height and amplitude of moistening with rain rate--successfully distinguish between higher- and lower-fidelity GCMs in hindcasts and climate simulations. To improve the MJO, developers should focus on relationships between convection and both total moisture and its rate of change. We conclude by offering recommendations for further experiments.

Vertical structure and physical processes of the Madden-Julian oscillation: exploring key model physics in climate simulations

Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one-fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models, but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high and low rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.