A comprehensive benchmarking system for evaluating global vegetation models.

We present a benchmark system for global vegetation models. This system provides a quantitative evaluation of multiple simulated vegetation properties, including primary production; seasonal net ecosystem production; vegetation cover, composition and 5 height; fire regime; and runoff. The benchmarks are derived from remotely sensed gridded datasets and site-based observations. The datasets allow comparisons of annual average conditions and seasonal and inter-annual variability, and they allow the impact of spatial and temporal biases in means and variability to be assessed separately. Specifically designed metrics quantify model performance for each process, 10 and are compared to scores based on the temporal or spatial mean value of the observations and a “random” model produced by bootstrap resampling of the observations. The benchmark system is applied to three models: a simple light-use efficiency and water-balance model (the Simple Diagnostic Biosphere Model: SDBM), and the Lund-Potsdam-Jena (LPJ) and Land Processes and eXchanges (LPX) dynamic global 15 vegetation models (DGVMs). SDBM reproduces observed CO2 seasonal cycles, but its simulation of independent measurements of net primary production (NPP) is too high. The two DGVMs show little difference for most benchmarks (including the interannual variability in the growth rate and seasonal cycle of atmospheric CO2), but LPX represents burnt fraction demonstrably more accurately. Benchmarking also identified 20 several weaknesses common to both DGVMs. The benchmarking system provides a quantitative approach for evaluating how adequately processes are represented in a model, identifying errors and biases, tracking improvements in performance through model development, and discriminating among models. Adoption of such a system would do much to improve confidence in terrestrial model predictions of climate change 25 impacts and feedbacks.

A comprehensive benchmarking system for evaluating global vegetation models

We present a benchmark system for global vegetation models. This system provides a quantitative evaluation of multiple simulated vegetation properties, including primary production; seasonal net ecosystem production; vegetation cover; composition and height; fire regime; and runoff. The benchmarks are derived from remotely sensed gridded datasets and site-based observations. The datasets allow comparisons of annual average conditions and seasonal and inter-annual variability, and they allow the impact of spatial and temporal biases in means and variability to be assessed separately. Specifically designed metrics quantify model performance for each process, and are compared to scores based on the temporal or spatial mean value of the observations and a "random" model produced by bootstrap resampling of the observations. The benchmark system is applied to three models: a simple light-use efficiency and water-balance model (the Simple Diagnostic Biosphere Model: SDBM), the Lund-Potsdam-Jena (LPJ) and Land Processes and eXchanges (LPX) dynamic global vegetation models (DGVMs). In general, the SDBM performs better than either of the DGVMs. It reproduces independent measurements of net primary production (NPP) but underestimates the amplitude of the observed CO2 seasonal cycle. The two DGVMs show little difference for most benchmarks (including the inter-annual variability in the growth rate and seasonal cycle of atmospheric CO2), but LPX represents burnt fraction demonstrably more accurately. Benchmarking also identified several weaknesses common to both DGVMs. The benchmarking system provides a quantitative approach for evaluating how adequately processes are represented in a model, identifying errors and biases, tracking improvements in performance through model development, and discriminating among models. Adoption of such a system would do much to improve confidence in terrestrial model predictions of climate change impacts and feedbacks.

The energy and water balance of a Eucalyptus plantation in southeast Brazil

The eddy covariance method was used to measure energy and water balance of a plantation of Eucalyptus (grandis x urophylla) hybrids over a 2 year period. The average daily evaporation rates were 5.4 (+/- 2.0) mm day(-1) in summer, but fell to 1.2 (+/- 0.3) mm day(-1) in winter. In contrast, the sensible heat flux was relatively low in summer but dominated the energy balance in winter. Evaporation accounted for 80% and 26% of the available energy, in summer and winter respectively. The annual evaporation was 82% (1124 mm) and 96% (1235 mm) of the annual rainfall recorded during the first and second year, respectively. Daily average canopy and aerodynamic conductance to water vapour were in the summer 51.9 (+/- 38.4) mm s(-1) 84.1 (+/- 25.6) mm s(-1), respectively; and in the winter 6.0 (+/- 10.5) mm s(-1) and 111.6 (+/- 24.6) mm s(-1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Estimating Annual CO2 Flux for Lutjewad Station Using Three Different Gap-Filling Techniques

Long-term measurements of CO2 flux can be obtained using the eddy covariance technique, but these datasets are affected by gaps which hinder the estimation of robust long-term means and annual ecosystem exchanges. We compare results obtained using three gap-fill techniques: multiple regression (MR), multiple imputation (MI), and artificial neural networks (ANNs), applied to a one-year dataset of hourly CO2 flux measurements collected in Lutjewad, over a flat agriculture area near the Wadden Sea dike in the north of the Netherlands. The dataset was separated in two subsets: a learning and a validation set. The performances of gap-filling techniques were analysed by calculating statistical criteria: coefficient of determination (R2), root mean square error (RMSE), mean absolute error (MAE), maximum absolute error (MaxAE), and mean square bias (MSB). The gap-fill accuracy is seasonally dependent, with better results in cold seasons. The highest accuracy is obtained using ANN technique which is also less sensitive to environmental/seasonal conditions. We argue that filling gaps directly on measured CO2 fluxes is more advantageous than the common method of filling gaps on calculated net ecosystem change, because ANN is an empirical method and smaller scatter is expected when gap filling is applied directly to measurements.

Diurnal and seasonal measurements of seawater chemistry, temperature and PAR in Hurricane Hole, U.S. Virgin Islands - 2010-2012

Risk analyses indicate that more than 90% of the world's reefs will be threatened by climate change and local anthropogenic impacts by the year 2030 under "business-as-usual" climate scenarios. Increasing temperatures and solar radiation cause coral bleaching that has resulted in extensive coral mortality. Increasing carbon dioxide reduces seawater pH, slows coral growth, and may cause loss of reef structure. Management strategies include establishment of marine protected areas with environmental conditions that promote reef resiliency. However, few resilient reefs have been identified, and resiliency factors are poorly defined. Here we characterize the first natural, non-reef coral refuge from thermal stress and ocean acidification and identify resiliency factors for mangrove-coral habitats. We measured diurnal and seasonal variations in temperature, salinity, photosynthetically active radiation (PAR), and seawater chemistry; characterized substrate parameters; and examined water circulation patterns in mangrove communities where scleractinian corals are growing attached to and under mangrove prop roots in Hurricane Hole, St. John, US Virgin Islands. Additionally, we inventoried the coral species and quantified incidences of coral bleaching, mortality, and recovery for two major reef-building corals, Colpophyllia natans and Diploria labyrinthiformis, growing in mangrove-shaded and exposed (unshaded) areas. Over 30 species of scleractinian corals were growing in association with mangroves. Corals were thriving in low-light (more than 70% attenuation of incident PAR) from mangrove shading and at higher temperatures than nearby reef tract corals. A higher percentage of C. natans colonies were living shaded by mangroves, and no shaded colonies were bleached. Fewer D. labyrinthiformis colonies were shaded by mangroves, however more unshaded colonies were bleached. A combination of substrate and habitat heterogeneity, proximity of different habitat types, hydrographic conditions, and biological influences on seawater chemistry generate chemical conditions that buffer against ocean acidification. This previously undocumented refuge for corals provides evidence for adaptation of coastal organisms and ecosystem transition due to recent climate change. Identifying and protecting other natural, non-reef coral refuges is critical for sustaining corals and other reef species into the future.

Sea surface fugacity of carbon dioxide measurements in the Indian and Southern Oceans obtained during MINERVE-29/ANTARES-II cruise

The sub-Antarctic zone (SAZ) lies between the subtropical convergence (STC) and the sub-Antarctic front (SAF), and is considered one of the strongest oceanic sinks of atmospheric CO2. The strong sink results from high winds and seasonally low sea surface fugacities of CO2 (fCO2), relative to atmospheric fCO2. The region of the SAZ, and immediately south, is also subject to mode and intermediate water formation, yielding a penetration of anthropogenic CO2 below the mixed layer. A detailed analysis of continuous measurements made during the same season and year, February - March 1993, shows a coherent pattern of fCO2 distributions at the eastern (WOCE/SR3 at about 145°E) and western edges (WOCE/I6 at 30°E) of the Indian sector of the Southern Ocean. A strong CO2 sink develops in the Austral summer (delta fCO2 < - 50 µatm) in both the eastern (110°-150°E) and western regions (20°-90°E). The strong CO2 sink in summer is due to the formation of a shallow seasonal mixed-layer (about 100 m). The CO2 drawdown in the surface water is consistent with biologically mediated drawdown of carbon over summer. In austral winter, surface fCO2 is close to equilibrium with the atmosphere (delta fCO2 ± 5 µatm), and the net CO2 exchange is small compared to summer. The near-equilibrium values in winter are associated with the formation of deep winter mixed-layers (up to 700 m). For years 1992-95, the annual CO2 uptake for the Indian Ocean sector of the sub Antarctic Zone (40°-50°S, 20°-150°E) is estimated to be about 0.4 GtC/yr. Extrapolating this estimate to the entire sub-Antarctic zone suggests the uptake in the circumpolar SAZ is approaching 1 GtC/yr.

Seaweed fails to prevent ocean acidification impact on foraminifera along a shallow-water CO2 gradient

Ocean acidification causes biodiversity loss, alters ecosystems, and may impact food security, as shells of small organisms dissolve easily in corrosive waters. There is a suggestion that photosynthetic organisms could mitigate ocean acidification on a local scale, through seagrass protection or seaweed cultivation, as net ecosystem organic production raises the saturation state of calcium carbonate making seawater less corrosive. Here, we used a natural gradient in calcium carbonate saturation, caused by shallow-water CO2 seeps in the Mediterranean Sea, to assess whether seaweed that is resistant to acidification (Padina pavonica) could prevent adverse effects of acidification on epiphytic foraminifera. We found a reduction in the number of species of foraminifera as calcium carbonate saturation state fell and that the assemblage shifted from one dominated by calcareous species at reference sites (pH 8.19) to one dominated by agglutinated foraminifera at elevated levels of CO2 (pH 7.71). It is expected that ocean acidification will result in changes in foraminiferal assemblage composition and agglutinated forms may become more prevalent. Although Padina did not prevent adverse effects of ocean acidification, high biomass stands of seagrass or seaweed farms might be more successful in protecting epiphytic foraminifera.

Hysteresis between coral reef calcification and the seawater aragonite saturation state

Some predictions of how ocean acidification (OA) will affect coral reefs assume a linear functional relationship between the ambient seawater aragonite saturation state (Omega a) and net ecosystem calcification (NEC). We quantified NEC in a healthy coral reef lagoon in the Great Barrier Reef during different times of the day. Our observations revealed a diel hysteresis pattern in the NEC versus Omega a relationship, with peak NEC rates occurring before the Omega a peak and relatively steady nighttime NEC in spite of variable Omega a. Net ecosystem production had stronger correlations with NEC than light, temperature, nutrients, pH, and Omega a. The observed hysteresis may represent an overlooked challenge for predicting the effects of OA on coral reefs. If widespread, the hysteresis could prevent the use of a linear extrapolation to determine critical Omega a threshold levels required to shift coral reefs from a net calcifying to a net dissolving state.

Freshwater mesocosm experiment 2012 in Vancouver, Canada on effects of increased temperature and food chain length on oxygen fluxes and plankton community structure

Over broad thermal gradients, the effect of temperature on aerobic respiration and photosynthesis rates explains variation in community structure and function. Yet for local communities, temperature dependent trophic interactions may dominate effects of warming. We tested the hypothesis that food chain length modifies the temperature-dependence of ecosystem fluxes and community structure. In a multi-generation aquatic food web experiment, increasing temperature strengthened a trophic cascade, altering the effect of temperature on estimated mass-corrected ecosystem fluxes. Compared to consumer-free and 3-level food chains, grazer-algae (2-level) food chains responded most strongly to the temperature gradient. Temperature altered community structure, shifting species composition and reducing zooplankton density and body size. Still, food chain length did not alter the temperature dependence of net ecosystem fluxes. We conclude that locally, food chain length interacts with temperature to modify community structure, but only temperature, not food chain length influenced net ecosystem fluxes.

Permeable coral reef sediment dissolution driven by elevated pCO2 and pore water advection

Ocean acidification (OA) is expected to drive the transition of coral reef ecosystems from net calcium carbonate (CaCO3) precipitating to net dissolving within the next century. Although permeable sediments represent the largest reservoir of CaCO3 in coral reefs, the dissolution of shallow CaCO3 sands under future pCO2 levels has not been measured under natural conditions. In situ, advective chamber incubations under elevated pCO2 (~800 µatm) shifted the sediments from net precipitating to net dissolving. Pore water advection more than doubled dissolution rates (1.10 g CaCO3/m**2/day) when compared to diffusive conditions (0.42 g CaCO3/m**2 /day). Sediment dissolution could reduce net ecosystem calcification rates of the Heron Island lagoon by 8% within the next century, which is equivalent to a 25% reduction in the global average calcification rate of coral lagoons. The dissolution of CaCO3 sediments needs to be taken into account in order to address how OA will impact the net accretion of coral reefs under future predicted increases in CO2.

A Functional Calvin Cycle Is Not Indispensable for the Light Activation of C4 Phosphoenolpyruvate Carboxylase Kinase and Its Target Enzyme in the Maize Mutant bundle sheath defective2-mutable11

We used a pale-green maize (Zea mays L.) mutant that fails to accumulate ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) to test the working hypothesis that the regulatory phosphorylation of C4 phosphoenolpyruvate carboxylase (PEPC) by its Ca2+-insensitive protein-serine/threonine kinase (PEPC kinase) in the C4 mesophyll cytosol depends on cross-talk with a functional Calvin cycle in the bundle sheath. Wild-type (W22) and bundle sheath defective2-mutable1 (bsd2-m1) seeds were grown in a controlled environment chamber at 100 to 130 μmol m−2 s−1 photosynthetic photon flux density, and leaf tissue was harvested 11 d after sowing, following exposure to various light intensities. Immunoblot analysis showed no major difference in the amount of polypeptide present for several mesophyll- and bundle-sheath-specific photosynthetic enzymes apart from Rubisco, which was either completely absent or very much reduced in the mutant. Similarly, leaf net CO2-exchange analysis and in vitro radiometric Rubisco assays showed that no appreciable carbon fixation was occurring in the mutant. In contrast, the sensitivity of PEPC to malate inhibition in bsd2-m1 leaves decreased significantly with an increase in light intensity, and there was a concomitant increase in PEPC kinase activity, similar to that seen in wild-type leaf tissue. Thus, although bsd2-m1 mutant plants lack an operative Calvin cycle, light activation of PEPC kinase and its target enzyme are not grossly perturbed.

Overexpression of Iron Superoxide Dismutase in Transformed Poplar Modifies the Regulation of Photosynthesis at Low CO2 Partial Pressures or Following Exposure to the Prooxidant Herbicide Methyl Viologen1

Chloroplast-targeted overexpression of an Fe superoxide dismutase (SOD) from Arabidopsis thaliana resulted in substantially increased foliar SOD activities. Ascorbate peroxidase, glutathione reductase, and monodehydroascorbate reductase activities were similar in the leaves from all of the lines, but dehydroascorbate reductase activity was increased in the leaves of the FeSOD transformants relative to untransformed controls. Foliar H2O2, ascorbate, and glutathione contents were comparable in all lines of plants. Irradiance-dependent changes in net CO2 assimilation and chlorophyll a fluorescence quenching parameters were similar in all lines both in air (21% O2) and at low (1%) O2. CO2-response curves for photosynthesis showed similar net CO2-exchange characteristics in all lines. In contrast, values of photochemical quenching declined in leaves from untransformed controls at intercellular CO2 (Ci) values below 200 μL L−1 but remained constant with decreasing Ci in leaves of FeSOD transformants. When the O2 concentration was decreased from 21 to 1%, the effect of FeSOD overexpression on photochemical quenching at limiting Ci was abolished. At high light (1000 μmol m−2 s−1) a progressive decrease in the ratio of variable (Fv) to maximal (Fm) fluorescence was observed with decreasing temperature. At 6oC the high-light-induced decrease in the Fv/Fm ratio was partially prevented by low O2 but values were comparable in all lines. Methyl viologen caused decreased Fv/Fm ratios, but this was less marked in the FeSOD transformants than in the untransformed controls. These observations suggest that the rate of superoxide dismutation limits flux through the Mehler-peroxidase cycle in certain conditions.

The oxygen and carbon dioxide compensation points of C3 plants: possible role in regulating atmospheric oxygen.

The O2 and CO2 compensation points (O2 and CO2) of plants in a closed system depend on the ratio of CO2 and O2 concentrations in air and in the chloroplast and the specificities of ribulose bisphosphate carboxylase/oxygenase (Rubisco). The photosynthetic O2 is defined as the atmospheric O2 level, with a given CO2 level and temperature, at which net O2 exchange is zero. In experiments with C3 plants, the O2 with 220 ppm CO2 is 23% O2; O2 increases to 27% with 350 ppm CO2 and to 35% O2 with 700 ppm CO2. At O2 levels below the O2, CO2 uptake and reduction are accompanied by net O2 evolution. At O2 levels above the O2, net O2 uptake occurs with a reduced rate of CO2 fixation, more carbohydrates are oxidized by photorespiration to products of the C2 oxidative photosynthetic carbon cycle, and plants senesce prematurely. The CO2 increases from 50 ppm CO2 with 21% O2 to 220 ppm with 100% O2. At a low CO2/high O2 ratio that inhibits the carboxylase activity of Rubisco, much malate accumulates, which suggests that the oxygen-insensitive phosphoenolpyruvate carboxylase becomes a significant component of the lower CO2 fixation rate. Because of low global levels of CO2 and a Rubisco specificity that favors the carboxylase activity, relatively rapid changes in the atmospheric CO2 level should control the permissive O2 that could lead to slow changes in the immense O2 pool.

Lake metabolism and the diel oxygen technique: State of the science

Significant improvements have been made in estimating gross primary production (GPP), ecosystem respiration (R), and net ecosystem production (NEP) from diel, “free-water” changes in dissolved oxygen (DO). Here we evaluate some of the assumptions and uncertainties that are still embedded in the technique and provide guidelines on how to estimate reliable metabolic rates from high-frequency sonde data. True whole-system estimates are often not obtained because measurements reflect an unknown zone of influence which varies over space and time. A minimum logging frequency of 30 min was sufficient to capture metabolism at the daily time scale. Higher sampling frequencies capture additional pattern in the DO data, primarily related to physical mixing. Causes behind the often large daily variability are discussed and evaluated for an oligotrophic and a eutrophic lake. Despite a 3-fold higher day-to-day variability in absolute GPP rates in the eutrophic lake, both lakes required at least 3 sonde days per week for GPP estimates to be within 20% of the weekly average. A sensitivity analysis evaluated uncertainties associated with DO measurements, piston velocity (k), and the assumption that daytime R equals nighttime R. In low productivity lakes, uncertainty in DO measurements and piston velocity strongly impacts R but has no effect on GPP or NEP. Lack of accounting for higher R during the day underestimates R and GPP but has no effect on NEP. We finally provide suggestions for future research to improve the technique.