Comparison of climate field reconstruction techniques: Application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.

Finding Homogeneity in Heterogeneity - A New Approach to Quantifying Landscape Mosaics Developed for the Lao PDR

A key challenge for land change science in general and research on swidden agriculture in particular, is linking land cover information to human–environment interactions over larger spatial areas. In Lao PDR, a country facing rapid and multi-level land change processes, this hinders informed policy- and decision-making. Crucial information on land use types and people involved is still lacking. This article proposes an alternative approach for the description of landscape mosaics. Instead of analyzing local land use combinations, we studied land cover mosaics at a meso-level of spatial scale and interpreted these in terms of human–environmental interactions. These landscape mosaics were then overlaid with population census data. Results showed that swidden agricultural landscapes, involving 17% of the population, dominate 29% of the country, while permanent agricultural landscapes involve 74% of the population in 29% of the territory. Forests still form an important component of these landscape mosaics.

Arctic Winds in the “Twentieth Century Reanalysis”

Climate in the European part of the Arctic underwent a rapid warming between the 1910s and the 1930s. Previous studies have addressed the role of atmospheric circulation in this period based on geopotential height fields because observations of upper-level winds in the Arctic are rare. Here we analyse winds over the Arctic and specifically over Spitsbergen in the “Twentieth Century Reanalyses” (20CR). We compare in situ upper-air wind measurements performed in 1912 and 1913 in Spitsbergen with six-hourly 20CR data. Furthermore, we compare monthly-to-seasonal 20CR winds at 700 hPa over the European Arctic with statistically reconstructed winds at 3 km altitude. Finally, we analyse long-term trends in Arctic winds in 20CR. The general agreement between observed upper-air winds and 20CR on the day-to-day scale is rather poor, which is not surprising given the paucity of observations in the Arctic at that time that constrain 20CR. In contrast, the seasonally averaged winds (which represent a larger spatial scale) in 20CR compare well with statistically reconstructed winds. The analysis of long term near-surface wind time series in 20CR shows arguably artificial trends from 1871 to around the 1950s over sparsely observed regions, particularly oceanic regions. Densely observed regions such as Europe or the USA show no such trends. This analysis shows that great care needs to be taken when working with 20CR in the Arctic and other sparsely observed regions.

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake-stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecologi- cal speciation. Ecological speciation is well studied in parts of the native range of the three-spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving differ- ent stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajec- tories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known eco- logical relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.

Transport of 234U in the Opalinus Clay on centimetre to decimetre scales

The Opalinus Clay formation in North Switzerland is a potential host rock for a deep underground radioactive waste repository. The distribution of U-238, U-234 and Th-230 was studied in rock samples of the Opalinus Clay from an exploratory borehole at Benken (Canton of Zurich) using MC-ICP-MS. The aim of U-234 was to assess the in situ, long-term migration behaviour in this rock. Very low hydraulic conductivities of the Opalinus Clay, reducing potential of the pore water and its chemical equilibrium with the host rock are expected to render both U-238 and Th-230 immobile. If U is heterogeneously distributed in the Opalinus Clay, gradients in the supply of U-234 from the rock matrix to the pore water by the decay of U-238 will be established. Diffusive redistribution separates U-234 from its immobile parent U-238 resulting in bulk rock U-234/U-238 activity disequilibria. These may provide a means of estimating the mobility of U-234 in the rock if the diffusion rate of U-234 is significant compared to its decay rate. Sampling was carried out on two scales. Drilling of cm-spaced samples from the drill-core was done to study mobility over short distances and elucidate possible small-scale lithological control. Homogenized 25-cm-long portions of a 2-m-long drill-core section were prepared to provide information on transport over a longer distance. Variations in U and/or Th content on the cm-scale between clays and carbonate-sandy layers are revealed by beta-scanning, which shows that the (dominant) clay is richer in both elements. Samples were digested using aqua regia followed by total HF dissolution, yielding two fractions. in all studied samples U was found to be concentrated in the HF digestion fraction. It has a high U/Th ratio and a study by SEM-EDS points to sub-mu m up to several mu m in size zircon grains as the main U-rich phase. This fraction consistently has U-234/U-238 activity ratios below unity. The minute zircon grains constitute the major reservoir of U in the rock and act as constant rate suppliers of U-234 into the rock matrix and the pore water. The aqua regia leach fraction was found to be enriched in Th, and complementary to the HF fraction, having U-234/U-238 activity ratios above unity. It is believed that these U activity ratios reflect the surplus of having U-234 delivered from the zircon grains. Some cm-spaced samples show bulk rock U-234/U-238 activity ratios that are markedly out of equilibrium. In most of them a striking negative correlation between the total U content and the bulk rock U-234/U-238 activity ratios is observed. This is interpreted to indicate net U-234 transfer from regions of higher supply of U-234 towards those of lower supply which is, in most cases, equivalent to transfer from clayey towards carbonate/sandy portions of the rock. In contrast, the 25 cm averaged samples all have uniform bulk rock U-234/U-238 activity ratios in equilibrium, indicating U immobility in the last 1-1.5 Ma on this spatial scale. It is concluded that the small-scale lithological variations which govern U spatial distribution in the Opalinus Clay are the major factor determining U-234 in situ supply rates, regulating its diffusive fluxes and controlling the observed bulk rock U-234/U-238 activity ratios. A simple box-model is presented to simulate the measured bulk rock U-234/U-238 activity ratios and to give an additional insight into the studied system. (C) 2008 Elsevier Ltd. All rights reserved.

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.

Pollen-based quantitative reconstructions of Holocene regional vegetation cover (plant-functional types and land-cover types) in Europe suitable for climate modelling

We present quantitative reconstructions of regional vegetation cover in north-western Europe, western Europe north of the Alps, and eastern Europe for five time windows in the Holocene around 6k, 3k, 0.5k, 0.2k, and 0.05k calendar years before present (bp)] at a 1 degrees x1 degrees spatial scale with the objective of producing vegetation descriptions suitable for climate modelling. The REVEALS model was applied on 636 pollen records from lakes and bogs to reconstruct the past cover of 25 plant taxa grouped into 10 plant-functional types and three land-cover types evergreen trees, summer-green (deciduous) trees, and open land]. The model corrects for some of the biases in pollen percentages by using pollen productivity estimates and fall speeds of pollen, and by applying simple but robust models of pollen dispersal and deposition. The emerging patterns of tree migration and deforestation between 6k bp and modern time in the REVEALS estimates agree with our general understanding of the vegetation history of Europe based on pollen percentages. However, the degree of anthropogenic deforestation (i.e. cover of cultivated and grazing land) at 3k, 0.5k, and 0.2k bp is significantly higher than deduced from pollen percentages. This is also the case at 6k in some parts of Europe, in particular Britain and Ireland. Furthermore, the relationship between summer-green and evergreen trees, and between individual tree taxa, differs significantly when expressed as pollen percentages or as REVEALS estimates of tree cover. For instance, when Pinus is dominant over Picea as pollen percentages, Picea is dominant over Pinus as REVEALS estimates. These differences play a major role in the reconstruction of European landscapes and for the study of land cover-climate interactions, biodiversity and human resources.

Gene flow and the maintenance of species boundaries

Hybrid zones are regions where individuals from genetically differentiated populations meet and mate, resulting in at least some offspring of mixed ancestry. Patterns of gene flow (introgression) in hybrid zones vary across the genome, allowing assessment of the role of individual genes or genome regions in reproductive isolation. Here, we document patterns of introgression between two recently diverged species of field crickets. We sampled at a very fine spatial scale and genotyped crickets for 110 highly differentiated single nucleotide polymorphisms (SNPs) identified through transcriptome scans. Using both genomic and geographic cline analysis, we document remarkably abrupt transitions (<100 m) in allele frequencies for 50 loci, despite high levels of gene flow at other loci. These are among the steepest clines documented for any hybridizing taxa. Furthermore, the cricket hybrid zone provides one of the clearest examples of the semi-permeability of species boundaries. Comparisons between data from the fine-scale transect and data (for the same set of markers) from sampling a much larger area in a different region of the cricket hybrid zone reveal consistent patterns of introgression for individual loci. The consistency in patterns of introgression between these two distant and distinct regions of the hybrid zone suggests that strong selection is acting to maintain abrupt discontinuities within the hybrid zone and that genomic regions with restricted introgression likely include genes that contribute to nonecological prezygotic barriers.

Pits formation from volatile outgassing on 67P/Churyumov-Gerasimenko

We investigate the thermal evolution of comet 67P/Churyumov-Gerasimenko's subsurface in the Seth_01 region, where active pits have been observed by the ESA/Rosetta mission. Our simulations show that clathrate destabilization and amorphous ice crystallization can occur at depths corresponding to those of the observed pits in a timescale shorter than 67P/Churyumov-Gerasimenko's lifetime in the comet's activity zone in the inner solar system. Sublimation of crystalline ice down to such depths is possible only in the absence of a dust mantle, which requires the presence of dust grains in the matrix small enough to be dragged out by gas from the pores. Our results are consistent with both pits formation via sinkholes or subsequent to outbursts, the dominant process depending on the status of the subsurface porosity. A sealed dust mantle would favor episodic and disruptive outgassing as a result of increasing gas pressure in the pores, while high porosity should allow the formation of large voids in the subsurface due to the continuous escape of volatiles. We finally conclude that the subsurface of 67P/Churyumov-Gerasimenko is not uniform at a spatial scale of similar to 100-200 m.

Extracting more out of relocation data: building movement models as mixtures of random walks

We present a framework for fitting multiple random walks to animal movement paths consisting of ordered sets of step lengths and turning angles. Each step and turn is assigned to one of a number of random walks, each characteristic of a different behavioral state. Behavioral state assignments may be inferred purely from movement data or may include the habitat type in which the animals are located. Switching between different behavioral states may be modeled explicitly using a state transition matrix estimated directly from data, or switching probabilities may take into account the proximity of animals to landscape features. Model fitting is undertaken within a Bayesian framework using the WinBUGS software. These methods allow for identification of different movement states using several properties of observed paths and lead naturally to the formulation of movement models. Analysis of relocation data from elk released in east-central Ontario, Canada, suggests a biphasic movement behavior: elk are either in an "encamped" state in which step lengths are small and turning angles are high, or in an "exploratory" state, in which daily step lengths are several kilometers and turning angles are small. Animals encamp in open habitat (agricultural fields and opened forest), but the exploratory state is not associated with any particular habitat type.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2004)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2004 just prior to mowing (during peak standing biomass in late May and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Stable isotope ratios and chemistry on corals from the Caribbean Sea

Instrumental climate data are limited in length and only available with low spatial coverage before the middle of the 20th century. This is too short to reliably determine and interpret decadal and longer scale climate variability and to understand the underlying mechanisms with sufficient accuracy. A proper knowledge of past variability of the climate system is needed to assess the anthropogenic impact on climate and ecosystems, and also important with regard to long-range climate forecasting. Highly-resolved records of past climate variations that extend beyond pre-industrial times can significantly help to understand long-term climate changes and trends. Indirect information on past environmental and climatic conditions can be deduced from climate-sensitive proxies. Large colonies of massive growing tropical reef corals have been proven to sensitively monitor changes in ambient seawater. Rapid skeletal growth, typically ranging between several millimeters to centimeters per year, allows the development of proxy records at sub-seasonal resolution. Stable oxygen isotopic composition and trace elemental ratios incorporated in the aragonitic coral skeleton can reveal a detailed history of past environmental conditions, e.g., sea surface temperature (SST). In general, coral-based reconstructions from the tropical Atlantic region have lagged behind the extensive work published using coral records from the Indian and Pacific Oceans. Difficulties in the analysis of previously utilized coral archives from the Atlantic, typically corals of the genera Montastrea and Siderastrea, have so far exacerbated the production of long-term high-resolution proxy records. The objective of this study is the evaluation of massive fast-growing corals of the species Diploria strigosa as a new marine archive for climate reconstructions from the tropical Atlantic region. For this purpose, coral records from two study sites in the eastern Caribbean Sea (Guadeloupe, Lesser Antilles; and Archipelago Los Roques, Venezuela) were examined. At Guadeloupe, a century-long monthly resolved multi-proxy coral record was generated. Results present the first d18O (Sr/Ca)-SST calibration equations for the Atlantic braincoral Diploria strigosa, that are robust and consistent with previously published values using other coral species from different regions. Both proxies reflect local variability of SST on a sub-seasonal scale, which is a precondition for studying seasonally phase-locked climate variations, as well as track variability on a larger spatial scale (i.e., in the Caribbean and tropical North Atlantic). Coral Sr/Ca reliably records local annual to interannual temperature variations and is higher correlated to in-situ air temperature than to grid-SST. The warming calculated from coral Sr/Ca is concurrent with the strong surface temperature increase at the study site during the past decades. Proxy data show a close relationship to major climate signals from the tropical Pacific and North Atlantic (the El Niño Southern Oscillation (ENSO) and the North Atlantic Oscillation (NAO)) affecting the seasonal cycle of SST in the North Tropical Atlantic (NTA). Coral oxygen isotopes are also influenced by seawater d18O (d18Osw) which is linked to the hydrological cycle, and capture large-scale climate variability in the NTA region better than Sr/Ca. Results from a quantitative comparison between extreme events in the two most prominent modes of external forcing, namely the ENSO and NAO, and respective events recorded in seasonal coral d18O imply that SST variability at the study site is highly linked to Pacific and North Atlantic variability, by this means supporting the assumptions of observational- and model-based studies which suggest a strong impact of ENSO and NAO forcings onto the NTA region through a modulation of trade wind strength in winter. Results from different spectral analysis tools suggest that interannual climate variability recorded by the coral proxies is II largely dictated by Pacific ENSO forcing, whereas at decadal and longer timescales the influence of the NAO is dominan. tThe Archipelago Los Roques is situated in the southeastern Caribbean Sea, north of the Venezuelan coast. Year-to-year variations in monthly resolved coral d18O of a nearcentury- long Diploria strigosa record are significantly correlated with SST and show pronounced multidecadal variations. About half of the variance in coral d18O can be explained by variations in seawater d18O, which can be estimated by calculating the d18Oresidual via subtracting the SST component from measured coral d18O. The d18Oresidual and a regional precipitation index are highly correlated at low frequencies, suggesting that d18Osw variations are primarily atmospheric-driven. Warmer SSTs at Los Roques broadly coincide with higher precipitation in the southeastern Caribbean at multidecadal time scales, effectively strengthening the climate signal in the coral d18O record. The Los Roques coral d18O record displays a strong and statistically significant relationship to different indices of hurricane activity during the peak of the Atlantic hurricane season in boreal summer and is a particularly good indicator of decadal-multidecadal swings in the latter indices. In general, the detection of long-term changes and trends in Atlantic hurricane activity is hampered due to the limited length of the reliable instrumental record and the known inhomogeneity in the observational databases which result from changes in observing practice and technology over the years. The results suggest that coral-derived proxy data from Los Roques can be used to infer changes in past hurricane activity on timescales that extend well beyond the reliable record. In addition, the coral record exhibits a clear negative trend superimposed on the decadal to multidecadal cycles, indicating a significant warming and freshening of surface waters in the genesis region of tropical cyclones during the past decades. The presented coral d18O time series provides the first and, so far, longest continuous coral-based record of hurricane activity. It appears that the combination of both signals (SST and d18Osw) in coral d18O leads to an amplification of large-scale climate signals in the record, and makes coral d18O even a better proxy for hurricane activity than SST alone. Atlantic hurricane activity naturally exhibits strong multidecadal variations that are associated with the Atlantic Multidecadal Oscillation (AMO), the major mode of lowfrequency variability in the North Atlantic Ocean. However, the mechanisms underlying this multidecadal variability remain controversial, primarily because of the limited instrumental record. The Los Roques coral d18O displays strong multidecadal variability with a period of approximately 60 years that is closely related to the AMO, making the Archipelago Los Roques a very sensitive location for studying low-frequency climate variability in the Atlantic Ocean. In summary, the coral records presented in this thesis capture different key climate variables in the north tropical Atlantic region very well, indicating that fast-growing Diploria strigosa corals represent a promising marine archive for further proxy-based reconstructions of past climate variability on a range of time scales.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2007)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2007 just prior to mowing (during peak standing biomass in early June and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four (May) or three (August) rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Aboveground community and species-specific plant biomass from the Jena Experiment (Main Experiment, year 2006)

This data set contains aboveground community biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice in 2006 just prior to mowing (during peak standing biomass in early June and in late August) on all experimental plots of the main experiment. This was done by clipping the vegetation at 3 cm above ground in four rectangles of 0.2 x 0.5 m per large plot. The location of these rectangles was assigned prior to each harvest by random selection of coordinates within the core area of the plots (i.e. the central 10 x 15 m). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material (i.e., dead plant material in the data file), and remaining plant material that could not be assigned to any category (i.e., unidentified plant material in the data file). All biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The data for individual samples and the mean over samples for the biomass measures on the community level are given. Overall, analyses of the community biomass data have identified species richness as well as functional group composition as important drivers of a positive biodiversity-productivity relationship.

Reefs shift from net accretion to net erosion along a natural environmental gradient

Coral reefs persist in an accretion-erosion balance and ocean acidification resulting from anthropogenic CO2 emissions threatens to shift this balance in favor of net reef erosion. Corals and calcifying algae, largely responsible for reef accretion, are vulnerable to environmental changes associated with ocean acidification, but the direct effects of lower pH on reef erosion has received less attention, particularly in the context of known drivers of bioerosion and natural variability. This study examines the balance between reef accretion and erosion along a well-characterized natural environmental gradient in Kane'ohe Bay, Hawai'i using experimental blocks of coral skeleton. Comparing before and after micro-computed tomography (µCT) scans to quantify net accretion and erosion, we show that, at the small spatial scale of this study (tens of meters), pH was a better predictor of the accretion-erosion balance than environmental drivers suggested by prior studies, including resource availability, temperature, distance from shore, or depth. In addition, this study highlights the fine-scale variation of pH in coastal systems and the importance of microhabitat variation for reef accretion and erosion processes. We demonstrate significant changes in both the mean and variance of pH on the order of meters, providing a local perspective on global increases in pCO2. Our findings suggest that increases in reef erosion, combined with expected decreases in calcification, will accelerate the shift of coral reefs to an erosion-dominated system in a high-CO2 world. This shift will make reefs increasingly susceptible to storm damage and sea-level rise, threatening the maintenance of the ecosystem services that coral reefs provide.