A data assimilation method for using low-resolution Earth observation data in heterogeneous ecosystems

We present an approach for dealing with coarse-resolution Earth observations (EO) in terrestrial ecosystem data assimilation schemes. The use of coarse-scale observations in ecological data assimilation schemes is complicated by spatial heterogeneity and nonlinear processes in natural ecosystems. If these complications are not appropriately dealt with, then the data assimilation will produce biased results. The “disaggregation” approach that we describe in this paper combines frequent coarse-resolution observations with temporally sparse fine-resolution measurements. We demonstrate the approach using a demonstration data set based on measurements of an Arctic ecosystem. In this example, normalized difference vegetation index observations are assimilated into a “zero-order” model of leaf area index and carbon uptake. The disaggregation approach conserves key ecosystem characteristics regardless of the observation resolution and estimates the carbon uptake to within 1% of the demonstration data set “truth.” Assimilating the same data in the normal manner, but without the disaggregation approach, results in carbon uptake being underestimated by 58% at an observation resolution of 250 m. The disaggregation method allows the combination of multiresolution EO and improves in spatial resolution if observations are located on a grid that shifts from one observation time to the next. Additionally, the approach is not tied to a particular data assimilation scheme, model, or EO product and can cope with complex observation distributions, as it makes no implicit assumptions of normality.

Climate change and coupling of macronutrient cycles along the atmospheric, terrestrial, freshwater and estuarine continuum

This paper provides an introduction to the Special Issue on “Climate Change and Coupling of Macronutrient Cycles along the Atmospheric, Terrestrial, Freshwater and Estuarine Continuum”, dedicated to Colin Neal on his retirement. It is not intended to be a review of this vast subject, but an attempt to synthesize some of the major findings from the 22 contributions to the Special Issue in the context of what is already known. The major research challenges involved in understanding coupled macronutrient cycles in these environmental media are highlighted, and the difficulties of making credible predictions of the effects of climate change are discussed. Of particular concern is the possibility of interactions which will enhance greenhouse gas concentrations and provide positive feedback to global warming.

A depth-integrated 2D coastal and estuarine model with conformal boundary-fitted mesh generation

Details are given of the development and application of a 2D depth-integrated, conformal boundary-fitted, curvilinear model for predicting the depth-mean velocity field and the spatial concentration distribution in estuarine and coastal waters. A numerical method for conformal mesh generation, based on a boundary integral equation formulation, has been developed. By this method a general polygonal region with curved edges can be mapped onto a regular polygonal region with the same number of horizontal and vertical straight edges and a multiply connected region can be mapped onto a regular region with the same connectivity. A stretching transformation on the conformally generated mesh has also been used to provide greater detail where it is needed close to the coast, with larger mesh sizes further offshore, thereby minimizing the computing effort whilst maximizing accuracy. The curvilinear hydrodynamic and solute model has been developed based on a robust rectilinear model. The hydrodynamic equations are approximated using the ADI finite difference scheme with a staggered grid and the solute transport equation is approximated using a modified QUICK scheme. Three numerical examples have been chosen to test the curvilinear model, with an emphasis placed on complex practical applications

Differentiation of neotropical ecosystems by modern soil phytolith assemblages and its implications for palaeoenvironmental and archaeological reconstructions

The interpretation of Neotropical fossil phytolith assemblages for palaeoenvironmental and archaeological reconstructions relies on the development of appropriate modern analogues. We analyzed modern phytolith assemblages from the soils of ten distinctive tropical vegetation communities in eastern lowland Bolivia, ranging from terra firme humid evergreen forest to seasonally-inundated savannah. Results show that broad ecosystems – evergreen tropical forest, semi-deciduous dry tropical forest, and savannah – can be clearly differentiated by examination of their phytolith spectra and the application of Principal Component Analysis (PCA). Differences in phytolith assemblages between particular vegetation communities within each of these ecosystems are more subtle, but can still be identified. Comparison of phytolith assemblages with pollen rain data and stable carbon isotope analyses from the same vegetation plots show that these proxies are not only complementary, but significantly improve taxonomic and ecosystem resolution, and therefore our ability to interpret palaeoenvironmental and archaeological records. Our data underline the utility of phytolith analyses for reconstructing Amazon Holocene vegetation histories and pre-Columbian land use, particularly the high spatial resolution possible with terrestrial soil-based phytolith studies.

Characterisation of Bolivian savanna ecosystems by their modern pollen rain and implications for fossil pollen records

The majority of vegetation reconstructions from the Neotropics are derived from fossil pollen records extracted from lake sediments. However, the interpretation of these records is restricted by limited knowledge of the contemporary relationships between the vegetation and pollen rain of Neotropical ecosystems, especially for more open vegetation such as savannas. This research aims to improve the interpretation of these records by investigating the vegetation and modern pollen rain of different savanna ecosystems in Bolivia using vegetation inventories, artificial pollen traps and surface lake sediments. Two types of savanna were studied, upland savannas (cerrado), occurring on well drained soils, and seasonally-inundated savannas occurring on seasonally water-logged soils. Quantitative vegetation data are used to identify taxa that are floristically important in the different savanna types and to allow modern pollen/vegetation ratios to be calculated. Artificial pollen traps from the upland savanna site are dominated by Moraceae (35%), Poaceae (30%), Alchornea (6%) and Cecropia (4%). The two seasonally-inundated savanna sites are dominated by Moraceae (37%), Poaceae (20%), Alchornea (8%) and Cecropia (7%), and Moraceae (25%), Cyperaceae (22%), Poaceae (19%) and Cecropia (9%), respectively. The modern pollen rain of seasonally-inundated savannas from surface lake sediments is dominated by Cyperaceae (35%), Poaceae (33%), Moraceae (9%) and Asteraceae (5%). Upland and seasonally-flooded savannas were found to be only subtly distinct from each other palynologically. All sites have a high proportion of Moraceae pollen due to effective wind dispersal of this pollen type from areas of evergreen forest close to the study sites. Modern pollen/vegetation ratios show that many key woody plant taxa are absent/under-represented in the modern pollen rain (e.g., Caryocar and Tabebuia). The lower-than-expected percentages of Poaceae pollen, and the scarcity of savanna indicators, in the modern pollen rain of these ecosystems mean that savannas could potentially be overlooked in fossil pollen records without consideration of the full pollen spectrum available.

Differentiation between Neotropical rainforest, dry forest and savannah ecosystems by their modern pollen spectra and implications for the fossil pollen record

Accurate differentiation between tropical forest and savannah ecosystems in the fossil pollen record is hampered by the combination of: i) poor taxonomic resolution in pollen identification, and ii) the high species diversity of many lowland tropical families, i.e. with many different growth forms living in numerous environmental settings. These barriers to interpreting the fossil record hinder our understanding of the past distributions of different Neotropical ecosystems and consequently cloud our knowledge of past climatic, biodiversity and carbon storage patterns. Modern pollen studies facilitate an improved understanding of how ecosystems are represented by the pollen their plants produce and therefore aid interpretation of fossil pollen records. To understand how to differentiate ecosystems palynologically, it is essential that a consistent sampling method is used across ecosystems. However, to date, modern pollen studies from tropical South America have employed a variety of methodologies (e.g. pollen traps, moss polsters, soil samples). In this paper, we present the first modern pollen study from the Neotropics to examine the modern pollen rain from moist evergreen tropical forest (METF), semi-deciduous dry tropical forest (SDTF) and wooded savannah (cerradão) using a consistent sampling methodology (pollen traps). Pollen rain was sampled annually in September for the years 1999–2001 from within permanent vegetation study plots in, or near, the Noel Kempff Mercado National Park (NKMNP), Bolivia. Comparison of the modern pollen rain within these plots with detailed floristic inventories allowed estimates of the relative pollen productivity and dispersal for individual taxa to be made (% pollen/% vegetation or ‘p/v’). The applicability of these data to interpreting fossil records from lake sediments was then explored by comparison with pollen assemblages obtained from five lake surface samples.

A digital ecosystems model of assessment feedback on student learning

The term ecosystem has been used to describe complex interactions between living organisms and the physical world. The principles underlying ecosystems can also be applied to complex human interactions in the digital world. As internet technologies make an increasing contribution to teaching and learning practice in higher education, the principles of digital ecosystems may help us understand how to maximise technology to benefit active, self-regulated learning especially among groups of learners. Here, feedback on student learning is presented within a conceptual digital ecosystems model of learning. Additionally, we have developed a Web 2.0-based system, called ASSET, which incorporates multimedia and social networking features to deliver assessment feedback within the functionality of the digital ecosystems model. Both the digital ecosystems model and the ASSET system are described and their implications for enhancing feedback on student learning are discussed.

The phenology of phytoplankton blooms: ecosystems indicators from remote sensing

Remote sensing offers many advantages in the development of ecosystem indicators for the pelagic zone of the ocean. Particularly suitable in this context are the indicators arising from time series that can be constructed from remotely sensed data. For example, using ocean-colour radiometry, the phenology of phytoplankton blooms can be assessed. Metrics defined in this way show promise as informative indicators for the entire pelagic ecosystem. A simple phytoplankton–substrate model, with forcing dependent on latitude and day number is used to explore the qualitative features of bloom phenology for comparison with the results observed in a suite of 10-year time series of chlorophyll concentration, as assessed by remote sensing, from the Northwest Atlantic Ocean. The model reveals features of the dynamics that might otherwise have been overlooked in evaluation of the observational data.

Terrestrial ecosystems from space: a review of earth observation products for macroecology applications

Aim Earth observation (EO) products are a valuable alternative to spectral vegetation indices. We discuss the availability of EO products for analysing patterns in macroecology, particularly related to vegetation, on a range of spatial and temporal scales. Location Global. Methods We discuss four groups of EO products: land cover/cover change, vegetation structure and ecosystem productivity, fire detection, and digital elevation models. We address important practical issues arising from their use, such as assumptions underlying product generation, product accuracy and product transferability between spatial scales. We investigate the potential of EO products for analysing terrestrial ecosystems. Results Land cover, productivity and fire products are generated from long-term data using standardized algorithms to improve reliability in detecting change of land surfaces. Their global coverage renders them useful for macroecology. Their spatial resolution (e.g. GLOBCOVER vegetation, 300 m; MODIS vegetation and fire, ≥ 500 m; ASTER digital elevation, 30 m) can be a limiting factor. Canopy structure and productivity products are based on physical approaches and thus are independent of biome-specific calibrations. Active fire locations are provided in near-real time, while burnt area products show actual area burnt by fire. EO products can be assimilated into ecosystem models, and their validation information can be employed to calculate uncertainties during subsequent modelling. Main conclusions Owing to their global coverage and long-term continuity, EO end products can significantly advance the field of macroecology. EO products allow analyses of spatial biodiversity, seasonal dynamics of biomass and productivity, and consequences of disturbances on regional to global scales. Remaining drawbacks include inter-operability between products from different sensors and accuracy issues due to differences between assumptions and models underlying the generation of different EO products. Our review explains the nature of EO products and how they relate to particular ecological variables across scales to encourage their wider use in ecological applications.

Millennial climatic fluctuations are key to the structure of last glacial ecosystems

Whereas fossil evidence indicates extensive treeless vegetation and diverse grazing megafauna in Europe and northern Asia during the last glacial, experiments combining vegetation models and climate models have to-date simulated widespread persistence of trees. Resolving this conflict is key to understanding both last glacial ecosystems and extinction of most of the mega-herbivores. Using a dynamic vegetation model (DVM) we explored the implications of the differing climatic conditions generated by a general circulation model (GCM) in “normal” and “hosing” experiments. Whilst the former approximate interstadial conditions, the latter, designed to mimic Heinrich Events, approximate stadial conditions. The “hosing” experiments gave simulated European vegetation much closer in composition to that inferred from fossil evidence than did the “normal” experiments. Given the short duration of interstadials, and the rate at which forest cover expanded during the late-glacial and early Holocene, our results demonstrate the importance of millennial variability in determining the character of last glacial ecosystems.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.