A geostatistical analysis of soil, vegetation, and image data characterizing land surface variation

The elucidation of spatial variation in the landscape can indicate potential wildlife habitats or breeding sites for vectors, such as ticks or mosquitoes, which cause a range of diseases. Information from remotely sensed data could aid the delineation of vegetation distribution on the ground in areas where local knowledge is limited. The data from digital images are often difficult to interpret because of pixel-to-pixel variation, that is, noise, and complex variation at more than one spatial scale. Landsat Thematic Mapper Plus (ETM+) and Satellite Pour l'Observation de La Terre (SPOT) image data were analyzed for an area close to Douna in Mali, West Africa. The variograms of the normalized difference vegetation index (NDVI) from both types of image data were nested. The parameters of the nested variogram function from the Landsat ETM+ data were used to design the sampling for a ground survey of soil and vegetation data. Variograms of the soil and vegetation data showed that their variation was anisotropic and their scales of variation were similar to those of NDVI from the SPOT data. The short- and long-range components of variation in the SPOT data were filtered out separately by factorial kriging. The map of the short-range component appears to represent the patterns of vegetation and associated shallow slopes and drainage channels of the tiger bush system. The map of the long-range component also appeared to relate to broader patterns in the tiger bush and to gentle undulations in the topography. The results suggest that the types of image data analyzed in this study could be used to identify areas with more moisture in semiarid regions that could support wildlife and also be potential vector breeding sites.

Probabilistic simulations of crop yield over western India using the DEMETER seasonal hindcast ensembles

Process-based integrated modelling of weather and crop yield over large areas is becoming an important research topic. The production of the DEMETER ensemble hindcasts of weather allows this work to be carried out in a probabilistic framework. In this study, ensembles of crop yield (groundnut, Arachis hypogaea L.) were produced for 10 2.5 degrees x 2.5 degrees grid cells in western India using the DEMETER ensembles and the general large-area model (GLAM) for annual crops. Four key issues are addressed by this study. First, crop model calibration methods for use with weather ensemble data are assessed. Calibration using yield ensembles was more successful than calibration using reanalysis data (the European Centre for Medium-Range Weather Forecasts 40-yr reanalysis, ERA40). Secondly, the potential for probabilistic forecasting of crop failure is examined. The hindcasts show skill in the prediction of crop failure, with more severe failures being more predictable. Thirdly, the use of yield ensemble means to predict interannual variability in crop yield is examined and their skill assessed relative to baseline simulations using ERA40. The accuracy of multi-model yield ensemble means is equal to or greater than the accuracy using ERA40. Fourthly, the impact of two key uncertainties, sowing window and spatial scale, is briefly examined. The impact of uncertainty in the sowing window is greater with ERA40 than with the multi-model yield ensemble mean. Subgrid heterogeneity affects model accuracy: where correlations are low on the grid scale, they may be significantly positive on the subgrid scale. The implications of the results of this study for yield forecasting on seasonal time-scales are as follows. (i) There is the potential for probabilistic forecasting of crop failure (defined by a threshold yield value); forecasting of yield terciles shows less potential. (ii) Any improvement in the skill of climate models has the potential to translate into improved deterministic yield prediction. (iii) Whilst model input uncertainties are important, uncertainty in the sowing window may not require specific modelling. The implications of the results of this study for yield forecasting on multidecadal (climate change) time-scales are as follows. (i) The skill in the ensemble mean suggests that the perturbation, within uncertainty bounds, of crop and climate parameters, could potentially average out some of the errors associated with mean yield prediction. (ii) For a given technology trend, decadal fluctuations in the yield-gap parameter used by GLAM may be relatively small, implying some predictability on those time-scales.

Climate change and coastal waterbird populations – past declines and future impacts

Considerable attention has been given to the impact of climate change on avian populations over the last decade. In this paper we examine two issues with respect to coastal bird populations in the UK: (1) is there any evidence that current populations are declining due to climate change, and (2) how might we predict the response of populations in the future? We review the cause of population decline in two species associated with saltmarsh habitats. The abundance of Common Redshank Tringa totanus breeding on saltmarsh declined by about 23% between the mid-1980s and mid-1990s, but the decline appears to have been caused by an increase in grazing pressure. The number of Twite Carduelis flavirostris wintering on the coast of East Anglia has declined dramatically over recent decades; there is evidence linking this decline with habitat loss but a causal role for climate change is unclear. These examples illustrate that climate change could be having population-level impacts now, but also show that it is dangerous to become too narrowly focused on single issues affecting coastal birds. Making predictions about how populations might respond to future climate change depends on an adequate understanding of important ecological processes at an appropriate spatial scale. We illustrate this with recent work conducted on the Icelandic population of Black-tailed Godwits Limosa limosa islandica that shows large-scale regulatory processes. Most predictive models to date have focused on local populations (single estuary or a group of neighbouring estuaries). We discuss the role such models might play in risk assessment, and the need for them to be linked to larger-scale ecological processes. We argue that future work needs to focus on spatial scale issues and on linking physical models of coastal environments with important ecological processes.

Climate change and coastal waterbird populations - past declines and future impacts

Considerable attention has been given to the impact of climate change on avian populations over the last decade. In this paper we examine two issues with respect to coastal bird populations in the UK: (1) is there any evidence that current populations are declining due to climate change, and (2) how might we predict the response of populations in the future? We review the cause of population decline in two species associated with saltmarsh habitats. The abundance of Common Redshank Tringa totanus breeding on saltmarsh declined by about 23% between the mid-1980s and mid-1990s, but the decline appears to have been caused by an increase in grazing pressure. The number of Twite Carduelis flavirostris wintering on the coast of East Anglia has declined dramatically over recent decades; there is evidence linking this decline with habitat loss but a causal role for climate change is unclear. These examples illustrate that climate change could be having population-level impacts now, but also show that it is dangerous to become too narrowly focused on single issues affecting coastal birds. Making predictions about how populations might respond to future climate change depends on an adequate understanding of important ecological processes at an appropriate spatial scale. We illustrate this with recent work conducted on the Icelandic population of Black-tailed Godwits Limosa limosa islandica that shows large-scale regulatory processes. Most predictive models to date have focused on local populations (single estuary or a group of neighbouring estuaries). We discuss the role such models might play in risk assessment, and the need for them to be linked to larger-scale ecological processes. We argue that future work needs to focus on spatial scale issues and on linking physical models of coastal environments with important ecological processes.

The relationship between convection and sea surface temperature on intraseasonal timescales

The relationship between tropical convection, surface fluxes, and sea surface temperature (SST) on intraseasonal timescales has been examined as part of an investigation of the possibility that the intraseasonal oscillation is a coupled atmosphere–ocean phenomenon. The unique feature of this study is that 15 yr of data and the whole region from the Indian Ocean to the Pacific Ocean have been analyzed using lag-correlation analysis and compositing techniques. A coherent relationship between convection, surface fluxes, and SST has been found on intraseasonal timescales in the Indian Ocean, Maritime Continent, and west Pacific regions of the Tropics. Prior to the maximum in convection, there are positive shortwave and latent heat flux anomalies into the surface, followed by warm SST anomalies about 10 days before the convective maximum. Coincident with the convective maximum, there is a minimum in the shortwave flux, followed by a cooling due to increased evaporation associated with enhanced westerly wind stress, leading to negative SST anomalies about 10 days after the convection. The relationships are robust from year to year, including both phases of the El Niño–Southern Oscillation (ENSO) although the eastward extent of the region over which the relationship holds varies with the phase of ENSO, consistent with the variations in the eastward extent of the warm pool and westerly winds. The spatial scale of the anomalies is about 60° longitude, consistent with the scale of the intraseasonal oscillation. The spatial and temporal characteristics of the surface flux and SST perturbations are consistent with the surface flux variations forcing the ocean, and the magnitudes of the anomalies are consistent with mixed-layer depths appropriate to the Indian Ocean and west Pacific

The observed relationships between wheat and climate in China

Recent changes in climate have had a measurable impact on crop yield in China. The objective of this study is to investigate how climate variability affects wheat yield in China at different spatial scales. First the response of wheat yield to the climate at the provincial level from 1978 to 1995 for China was analysed. Wheat yield variability was only correlated with climate variability in some regions of China. At the provincial level, the variability of precipitation had a negative impact on wheat yield in parts of southeast China, but the seasonal mean temperature had a negative impact on wheat yield in only a few provinces, where significant variability in precipitation explained about 23–60% of yield variability, and temperature variability accounted for 37–41% of yield variability from 1978 to 1995. The correlation between wheat yield and climate for the whole of China from 1985 to 2000 was investigated at five spatial scales using climate data. The Climate Research Unit (CRU) and National Centers for Environmental Prediction (NCEP) proportions of the grid cells with a significant yield–precipitation correlation declined progressively from 14.6% at 0.5° to 0% at 5° scale. In contrast, the proportion of grid cells significant for the yield–temperature correlation increased progressively from 1.9% at 0.5° scale to 16% at 5° scale. This indicates that the variability of precipitation has a higher association with wheat yield at small scales (0.5°, 2°/2.5°) than at larger scales (4°/5.0°); but wheat yield has a good association with temperature at all levels of aggregation. The precipitation variable at the smaller scales (0.5°, 2°/2.5°) is a dominant factor in determining inter-annual wheat yield variability more so than at the larger scales (4°/5°). We conclude that in the current climate the relationship between wheat yield and each of precipitation and temperature becomes weaker and stronger, respectively, with an increase in spatial scale.

An ocean modelling and assimilation guide to using GOCE geoid products

We review the procedures and challenges that must be considered when using geoid data derived from the Gravity and steady-state Ocean Circulation Explorer (GOCE) mission in order to constrain the circulation and water mass representation in an ocean 5 general circulation model. It covers the combination of the geoid information with timemean sea level information derived from satellite altimeter data, to construct a mean dynamic topography (MDT), and considers how this complements the time-varying sea level anomaly, also available from the satellite altimeter. We particularly consider the compatibility of these different fields in their spatial scale content, their temporal rep10 resentation, and in their error covariances. These considerations are very important when the resulting data are to be used to estimate ocean circulation and its corresponding errors. We describe the further steps needed for assimilating the resulting dynamic topography information into an ocean circulation model using three different operational fore15 casting and data assimilation systems. We look at methods used for assimilating altimeter anomaly data in the absence of a suitable geoid, and then discuss different approaches which have been tried for assimilating the additional geoid information. We review the problems that have been encountered and the lessons learned in order the help future users. Finally we present some results from the use of GRACE geoid in20 formation in the operational oceanography community and discuss the future potential gains that may be obtained from a new GOCE geoid.

The general observation period 2007 within the priority program on quantitative precipitation forecasting: concept and first results

In the year 2007 a General Observation Period (GOP) has been performed within the German Priority Program on Quantitative Precipitation Forecasting (PQP). By optimizing the use of existing instrumentation a large data set of in-situ and remote sensing instruments with special focus on water cycle variables was gathered over the full year cycle. The area of interest covered central Europe with increasing focus towards the Black Forest where the Convective and Orographically-induced Precipitation Study (COPS) took place from June to August 2007. Thus the GOP includes a variety of precipitation systems in order to relate the COPS results to a larger spatial scale. For a timely use of the data, forecasts of the numerical weather prediction models COSMO-EU and COSMO-DE of the German Meteorological Service were tailored to match the observations and perform model evaluation in a near real-time environment. The ultimate goal is to identify and distinguish between different kinds of model deficits and to improve process understanding.

Using learning networks to understand complex systems: a case study of biological, geophysical and social research in the Amazon

Developing high-quality scientific research will be most effective if research communities with diverse skills and interests are able to share information and knowledge, are aware of the major challenges across disciplines, and can exploit economies of scale to provide robust answers and better inform policy. We evaluate opportunities and challenges facing the development of a more interactive research environment by developing an interdisciplinary synthesis of research on a single geographic region. We focus on the Amazon as it is of enormous regional and global environmental importance and faces a highly uncertain future. To take stock of existing knowledge and provide a framework for analysis we present a set of mini-reviews from fourteen different areas of research, encompassing taxonomy, biodiversity, biogeography, vegetation dynamics, landscape ecology, earth-atmosphere interactions, ecosystem processes, fire, deforestation dynamics, hydrology, hunting, conservation planning, livelihoods, and payments for ecosystem services. Each review highlights the current state of knowledge and identifies research priorities, including major challenges and opportunities. We show that while substantial progress is being made across many areas of scientific research, our understanding of specific issues is often dependent on knowledge from other disciplines. Accelerating the acquisition of reliable and contextualized knowledge about the fate of complex pristine and modified ecosystems is partly dependent on our ability to exploit economies of scale in shared resources and technical expertise, recognise and make explicit interconnections and feedbacks among sub-disciplines, increase the temporal and spatial scale of existing studies, and improve the dissemination of scientific findings to policy makers and society at large. Enhancing interaction among research efforts is vital if we are to make the most of limited funds and overcome the challenges posed by addressing large-scale interdisciplinary questions. Bringing together a diverse scientific community with a single geographic focus can help increase awareness of research questions both within and among disciplines, and reveal the opportunities that may exist for advancing acquisition of reliable knowledge. This approach could be useful for a variety of globally important scientific questions.

Local management and landscape drivers of pollination and biological control services in a Kenyan agro-ecosystem

Arthropods that have a direct impact on crop production (i.e. pests, natural enemies and pollinators) can be influenced by both local farm management and the context within which the fields occur in the wider landscape. However, the contributions and spatial scales at which these drivers operate and interact are not fully understood, particularly in the developing world. The impact of both local management and landscape context on insect pollinators and natural enemy communities and on their capacity to deliver related ecosystem services to an economically important tropical crop, pigeonpea was investigated. The study was conducted in nine paired farms across a gradient of increasing distance to semi-native vegetation in Kibwezi, Kenya. Results show that proximity of fields to semi-native habitats negatively affected pollinator and chewing insect abundance. Within fields, pesticide use was a key negative predictor of pollinator, pest and foliar active predator abundance. On the contrary, fertilizer application significantly enhanced pollinator and both chewing and sucking insect pest abundance. At a 1 km spatial scale of fields, there were significant negative effects of the number of semi-native habitat patches within fields dominated by mass flowering pigeonpea on pollinators abundance. For service provision, a significant decline in fruit set when insects were excluded from flowers was recorded. This study reveals the interconnections of pollinators, predators and pests with pigeonpea crop. For sustainable yields and to conserve high densities of both pollinators and predators of pests within pigeonpea landscapes, it is crucial to target the adoption of less disruptive farm management practices such as reducing pesticide and fertilizer inputs.

A complete atomistic model of molten polyethylene from neutron scattering data: a new methodology for polymer structure

A new approach to the study of the local organization in amorphous polymer materials is presented. The method couples neutron diffraction experiments that explore the structure on the spatial scale 1–20 Å with the reverse Monte Carlo fitting procedure to predict structures that accurately represent the experimental scattering results over the whole momentum transfer range explored. Molecular mechanics and molecular dynamics techniques are also used to produce atomistic models independently from any experimental input, thereby providing a test of the viability of the reverse Monte Carlo method in generating realistic models for amorphous polymeric systems. An analysis of the obtained models in terms of single chain properties and of orientational correlations between chain segments is presented. We show the viability of the method with data from molten polyethylene. The analysis derives a model with average C-C and C-H bond lengths of 1.55 Å and 1.1 Å respectively, average backbone valence angle of 112, a torsional angle distribution characterized by a fraction of trans conformers of 0.67 and, finally, a weak interchain orientational correlation at around 4 Å.

Extracting quantitative structural parameters for disordered polymers from neutron scattering data

The organization of non-crystalline polymeric materials at a local level, namely on a spatial scale between a few and 100 a, is still unclear in many respects. The determination of the local structure in terms of the configuration and conformation of the polymer chain and of the packing characteristics of the chain in the bulk material represents a challenging problem. Data from wide-angle diffraction experiments are very difficult to interpret due to the very large amount of information that they carry, that is the large number of correlations present in the diffraction patterns.We describe new approaches that permit a detailed analysis of the complex neutron diffraction patterns characterizing polymer melts and glasses. The coupling of different computer modelling strategies with neutron scattering data over a wide Q range allows the extraction of detailed quantitative information on the structural arrangements of the materials of interest. Proceeding from modelling routes as diverse as force field calculations, single-chain modelling and reverse Monte Carlo, we show the successes and pitfalls of each approach in describing model systems, which illustrate the need to attack the data analysis problem simultaneously from several fronts.

Sea surface temperature from a geostationary satellite by optimal estimation

Optimal estimation (OE) is applied as a technique for retrieving sea surface temperature (SST) from thermal imagery obtained by the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) on Meteosat 9. OE requires simulation of observations as part of the retrieval process, and this is done here using numerical weather prediction fields and a fast radiative transfer model. Bias correction of the simulated brightness temperatures (BTs) is found to be a necessary step before retrieval, and is achieved by filtered averaging of simulations minus observations over a time period of 20 days and spatial scale of 2.5° in latitude and longitude. Throughout this study, BT observations are clear-sky averages over cells of size 0.5° in latitude and longitude. Results for the OE SST are compared to results using a traditional non-linear retrieval algorithm (“NLSST”), both validated against a set of 30108 night-time matches with drifting buoy observations. For the OE SST the mean difference with respect to drifter SSTs is − 0.01 K and the standard deviation is 0.47 K, compared to − 0.38 K and 0.70 K respectively for the NLSST algorithm. Perhaps more importantly, systematic biases in NLSST with respect to geographical location, atmospheric water vapour and satellite zenith angle are greatly reduced for the OE SST. However, the OE SST is calculated to have a lower sensitivity of retrieved SST to true SST variations than the NLSST. This feature would be a disadvantage for observing SST fronts and diurnal variability, and raises questions as to how best to exploit OE techniques at SEVIRI's full spatial resolution.

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.

The population dynamics of disease on short and long time-scales

Observational evidence is scarce concerning the distribution of plant pathogen population sizes or densities as a function of time-scale or spatial scale. For wild pathosystems we can only get indirect evidence from evolutionary patterns and the consequences of biological invasions.We have little or no evidence bearing on extermination of hosts by pathogens, or successful escape of a host from a pathogen. Evidence over the last couple of centuries from crops suggest that the abundance of particular pathogens in the spectrum affecting a given host can vary hugely on decadal timescales. However, this may be an artefact of domestication and intensive cultivation. Host-pathogen dynamics can be formulated mathematically fairly easily–for example as SIR-type differential equation or difference equation models, and this has been the (successful) focus of recent work in crops. “Long-term” is then discussed in terms of the time taken to relax from a perturbation to the asymptotic state. However, both host and pathogen dynamics are driven by environmental factors as well as their mutual interactions, and both host and pathogen co-evolve, and evolve in response to external factors. We have virtually no information about the importance and natural role of higher trophic levels (hyperpathogens) and competitors, but they could also induce long-scale fluctuations in the abundance of pathogens on particular hosts. In wild pathosystems the host distribution cannot be modelled as either a uniform density or even a uniform distribution of fields (which could then be treated as individuals). Patterns of short term density-dependence and the detail of host distribution are therefore critical to long-term dynamics. Host density distributions are not usually scale-free, but are rarely uniform or clearly structured on a single scale. In a (multiply structured) metapopulation with coevolution and external disturbances it could well be the case that the time required to attain equilibrium (if it exists) based on conditions stable over a specified time-scale is longer than that time-scale. Alternatively, local equilibria may be reached fairly rapidly following perturbations but the meta-population equilibrium be attained very slowly. In either case, meta-stability on various time-scales is a more relevant than equilibrium concepts in explaining observed patterns.