Uncertainty in the estimation of potential evapotranspiration under climate change.

21st century climate change is projected to result in an intensification of the global hydrological cycle, but there is substantial uncertainty in how this will impact freshwater availability. A relatively overlooked aspect of this uncertainty pertains to how different methods of estimating potential evapotranspiration (PET) respond to changing climate. Here we investigate the global response of six different PET methods to a 2 °C rise in global mean temperature. All methods suggest an increase in PET associated with a warming climate. However, differences in PET climate change signal of over 100% are found between methods. Analysis of a precipitation/PET aridity index and regional water surplus indicates that for certain regions and GCMs, choice of PET method can actually determine the direction of projections of future water resources. As such, method dependence of the PET climate change signal is an important source of uncertainty in projections of future freshwater availability.

Microbial and biochemical soil quality indicators and their potential for differentiating areas under contrasting agricultural management regimes

The aim of this study was to examine interrelationships between functional biochemical and microbial indicators of soil quality, and their suitability to differentiate areas under contrasting agricultural management regimes. The study included five 0.8 ha areas on a sandy-loam soil which had received contrasting fertility and cropping regimes over a 5 year period. These were organically managed vegetable, vegetable -cereal and arable rotations, an organically managed grass clover ley, and a conventional cereal rotation. The organic areas had been converted from conventional cereal production 5 years prior to the start of the study. All of the biochemical analyses, including light fraction organic matter (LFOM) C and N, labile organic N (LON), dissolved organic N and water-soluble carbohydrates showed significant differences between the areas, although the nature of the relationships between the areas varied between the different parameters, and were not related to differences in total soil organic matter content. The clearest differences were seen in LFOM C and N and LON, which were higher in the organic arable area relative to the other areas. In the case of the biological parameters, there were differences between the areas for biomass-N, ATP, chitin content, and the ratios of ATP: biomass and basal respiration: biomass. For these parameters, the precise relationships between the areas varied. However, relative to the conventionally managed area, areas under organic management generally had lower biomass-N and higher ATP contents. Arbuscular mycorrhizal fungus colonization potential was extremely low in the conventional area relative to the organic areas. Further, metabolic diversity and microbial community level physiological profiles, determined by analysis of microbial community metabolism using Biolog GN plates and the activities of eight key nutrient cycling enzymes, grouped the organic areas together, but separated them from the conventional area. We conclude that microbial parameters are more effective and consistent indicators of management induced changes to soil quality than biochemical parameters, and that a variety of biochemical and microbial analyses should be used when considering the impact of management on soil quality. (C) 2004 Elsevier Ltd. All rights reserved.

A comparison of diatom phosphorus transfer functions and export coefficient models as tools for reconstructing lake nutrient histories

1. We compared the baseline phosphorus (P) concentrations inferred by diatom-P transfer functions and export coefficient models at 62 lakes in Great Britain to assess whether the techniques produce similar estimates of historical nutrient status. 2. There was a strong linear relationship between the two sets of values over the whole total P (TP) gradient (2-200 mu g TP L-1). However, a systematic bias was observed with the diatom model producing the higher values in 46 lakes (of which values differed by more than 10 mu g TP L-1 in 21). The export coefficient model gave the higher values in 10 lakes (of which the values differed by more than 10 mu g TP L-1 in only 4). 3. The difference between baseline and present-day TP concentrations was calculated to compare the extent of eutrophication inferred by the two sets of model output. There was generally poor agreement between the amounts of change estimated by the two approaches. The discrepancy in both the baseline values and the degree of change inferred by the models was greatest in the shallow and more productive sites. 4. Both approaches were applied to two lakes in the English Lake District where long-term P data exist, to assess how well the models track measured P concentrations since approximately 1850. There was good agreement between the pre-enrichment TP concentrations generated by the models. The diatom model paralleled the steeper rise in maximum soluble reactive P (SRP) more closely than the gradual increase in annual mean TP in both lakes. The export coefficient model produced a closer fit to observed annual mean TP concentrations for both sites, tracking the changes in total external nutrient loading. 5. A combined approach is recommended, with the diatom model employed to reflect the nature and timing of the in-lake response to changes in nutrient loading, and the export coefficient model used to establish the origins and extent of changes in the external load and to assess potential reduction in loading under different management scenarios. 6. However, caution must be exercised when applying these models to shallow lakes where the export coefficient model TP estimate will not include internal P loading from lake sediments and where the diatom TP inferences may over-estimate TP concentrations because of the high abundance of benthic taxa, many of which are poor indicators of trophic state.

Phosphorus dynamics and export in streams draining micro-catchments: Development of empirical models

Annual total phosphorus (TP) export data from 108 European micro-catchments were analyzed against descriptive catchment data on climate (runoff), soil types, catchment size, and land use. The best possible empirical model developed included runoff, proportion of agricultural land and catchment size as explanatory variables but with a low explanation of the variance in the dataset (R-2 = 0.37). Improved country specific empirical models could be developed in some cases. The best example was from Norway where an analysis of TP-export data from 12 predominantly agricultural micro-catchments revealed a relationship explaining 96% of the variance in TP-export. The explanatory variables were in this case soil-P status (P-AL), proportion of organic soil, and the export of suspended sediment. Another example is from Denmark where an empirical model was established for the basic annual average TP-export from 24 catchments with percentage sandy soils, percentage organic soils, runoff, and application of phosphorus in fertilizer and animal manure as explanatory variables (R-2 = 0.97).

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.

Utilizing chromophoric dissolved organic matter measurements to derive export and reactivity of dissolved organic carbon exported to the Arctic Ocean: A case study of the Yukon River, Alaska

The quality and quantity of dissolved organic matter (DOM) exported by Arctic rivers is known to vary with hydrology and this exported material plays a fundamental role in the biogeochemical cycling of carbon at high latitudes. We highlight the potential of optical measurements to examine DOM quality across the hydrograph in Arctic rivers. Furthermore, we establish chromophoric DOM (CDOM) relationships to dissolved organic carbon (DOC) and lignin phenols in the Yukon River and model DOC and lignin loads from CDOM measurements, the former in excellent agreement with long-term DOC monitoring data. Intensive sampling across the historically under-sampled spring flush period highlights the importance of this time for total export of DOC and particularly lignin. Calculated riverine DOC loads to the Arctic Ocean show an increase from previous estimates, especially when new higher discharge data are incorporated. Increased DOC loads indicate decreased residence times for terrigenous DOM in the Arctic Ocean with important implications for the reactivity and export of this material to the Atlantic Ocean.

A sensitivity analysis of soil moisture retrieval from the t-w microwave emission model

The potential of the τ-ω model for retrieving the volumetric moisture content of bare and vegetated soil from dual polarisation passive microwave data acquired at single and multiple angles is tested. Measurement error and several additional sources of uncertainty will affect the theoretical retrieval accuracy. These include uncertainty in the soil temperature, the vegetation structure and consequently its microwave singlescattering albedo, and uncertainty in soil microwave emissivity based on its roughness. To test the effects of these uncertainties for simple homogeneous scenes, we attempt to retrieve soil moisture from a number of simulated microwave brightness temperature datasets generated using the τ-ω model. The uncertainties for each influence are estimated and applied to curves generated for typical scenarios, and an inverse model used to retrieve the soil moisture content, vegetation optical depth and soil temperature. The effect of each influence on the theoretical soil moisture retrieval limit is explored, the likelihood of each sensor configuration meeting user requirements is assessed, and the most effective means of improving moisture retrieval indicated.

The omega equation and potential vorticity

Two fundamental perspectives on the dynamics of midlatitude weather systems are provided by potential vorticity (PV) and the omega equation. The aim of this paper is to investigate the link between the two perspectives, which has so far received very little attention in the meteorological literature. It also aims to give a quantitative basis for discussion of quasi-geostrophic vertical motion in terms of components associated with system movement, maintaining a constant thermal structure, and with the development of that structure. The former links with the isentropic relative-flow analysis technique. Viewed in a moving frame of reference, the measured development of a system depends on the velocity of that frame of reference. The requirement that the development should be a minimum provides a quantitative method for determining the optimum system velocity. The component of vertical velocity associated with development is shown to satisfy an omega equation with forcing determined from the relative advection of interior PV and boundary temperature. The analysis carries through in the presence of diabatic heating provided the omega equation forcing is based on the interior PV and boundary thermal tendencies, including the heating effect. The analysis is shown to be possible also at the level of the semi-geostrophic approximation. The analysis technique is applied to a number of idealized problems that can be considered to be building blocks for midlatitude synoptic-scale dynamics. They focus on the influences of interior PV, boundary temperature, an interior boundary, baroclinic instability associated with two boundaries, and also diabatic heating. In each case, insights yielded by the new perspective are sought into the dynamical behaviour, especially that related to vertical motion. Copyright © 2003 Royal Meteorological Society

The potential to narrow uncertainty in projections of regional precipitation change

We separate and quantify the sources of uncertainty in projections of regional (*2,500 km) precipitation changes for the twenty-first century using the CMIP3 multi-model ensemble, allowing a direct comparison with a similar analysis for regional temperature changes. For decadal means of seasonal mean precipitation, internal variability is the dominant uncertainty for predictions of the first decade everywhere, and for many regions until the third decade ahead. Model uncertainty is generally the dominant source of uncertainty for longer lead times. Scenario uncertainty is found to be small or negligible for all regions and lead times, apart from close to the poles at the end of the century. For the global mean, model uncertainty dominates at all lead times. The signal-to-noise ratio (S/N) of the precipitation projections is highest at the poles but less than 1 almost everywhere else, and is far lower than for temperature projections. In particular, the tropics have the highest S/N for temperature, but the lowest for precipitation. We also estimate a ‘potential S/N’ by assuming that model uncertainty could be reduced to zero, and show that, for regional precipitation, the gains in S/N are fairly modest, especially for predictions of the next few decades. This finding suggests that adaptation decisions will need to be made in the context of high uncertainty concerning regional changes in precipitation. The potential to narrow uncertainty in regional temperature projections is far greater. These conclusions on S/N are for the current generation of models; the real signal may be larger or smaller than the CMIP3 multi-model mean. Also note that the S/N for extreme precipitation, which is more relevant for many climate impacts, may be larger than for the seasonal mean precipitation considered here.

Excess open solar magnetic flux from satellite data: 1. Analysis of the third perihelion Ulysses pass

We use the third perihelion pass by the Ulysses spacecraft to illustrate and investigate the “flux excess” effect, whereby open solar flux estimates from spacecraft increase with increasing heliocentric distance. We analyze the potential effects of small-scale structure in the heliospheric field (giving fluctuations in the radial component on timescales smaller than 1 h) and kinematic time-of-flight effects of longitudinal structure in the solar wind flow. We show that the flux excess is explained by neither very small-scale structure (timescales < 1 h) nor by the kinematic “bunching effect” on spacecraft sampling. The observed flux excesses is, however, well explained by the kinematic effect of larger-scale (>1 day) solar wind speed variations on the frozen-in heliospheric field. We show that averaging over an interval T (that is long enough to eliminate structure originating in the heliosphere yet small enough to avoid cancelling opposite polarity radial field that originates from genuine sector structure in the coronal source field) is only an approximately valid way of allowing for these effects and does not adequately explain or account for differences between the streamer belt and the polar coronal holes.

A baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change

This report forms part of a larger research programme on 'Reinterpreting the Urban-Rural Continuum', which conceptualises and investigates current knowledge and research gaps concerning 'the role that ecosystems services play in the livelihoods of the poor in regions undergoing rapid change'. The report aims to conduct a baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. The appraisal is conducted at three spatial scales: global, regional (four consortia areas), and meso scale (case studies within the four regions). At all three scales of analysis water resources form the interweaving theme because water provides a vital provisioning service for people, supports all other ecosystem processes and because water resources are forecast to be severely affected under climate change scenarios. This report, combined with an Endnote library of over 1100 scientific papers, provides an annotated bibliography of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. After an introductory, section, Section 2 of the report defines water-related ecosystem services and how these are affected by human activities. Current knowledge and research gaps are then explored in relation to global scale climate and related hydrological changes (e.g. floods, droughts, flow regimes) (section 3). The report then discusses the impacts of climate changes on the ESPA regions, emphasising potential responses of biomes to the combined effects of climate change and human activities (particularly land use and management), and how these effects coupled with water store and flow regime manipulation by humans may affect the functioning of catchments and their ecosystem services (section 4). Finally, at the meso-scale, case studies are presented from within the ESPA regions to illustrate the close coupling of human activities and catchment performance in the context of environmental change (section 5). At the end of each section, research needs are identified and justified. These research needs are then amalgamated in section 6.

A lightweight technique for assessing risks in requirements analysis

A simple and practical technique for assessing the risks, that is, the potential for error, and consequent loss, in software system development, acquired during a requirements engineering phase is described. The technique uses a goal-based requirements analysis as a framework to identify and rate a set of key issues in order to arrive at estimates of the feasibility and adequacy of the requirements. The technique is illustrated and how it has been applied to a real systems development project is shown. How problems in this project could have been identified earlier is shown, thereby avoiding costly additional work and unhappy users.

Time series analysis of lava flux

We have applied time series analytical techniques to the flux of lava from an extrusive eruption. Tilt data acting as a proxy for flux are used in a case study of the May–August 1997 period of the eruption at Soufrière Hills Volcano, Montserrat. We justify the use of such a proxy by simple calibratory arguments. Three techniques of time series analysis are employed: spectral, spectrogram and wavelet methods. In addition to the well-known ~9-hour periodicity shown by these data, a previously unknown periodic flux variability is revealed by the wavelet analysis as a 3-day cycle of frequency modulation during June–July 1997, though the physical mechanism responsible is not clear. Such time series analysis has potential for other lava flux proxies at other types of volcanoes.

Horizontal potential vorticity dipoles on the convective storm scale

The structure and dynamics of potential vorticity (PV) anomalies generated by convective storms is investigated both theoretically and in a numerical model case study. Linear theory suggests that if the storm-induced heating is on a sufficiently small scale (relative to the Rossby radius of deformation), and the environment contains moderate vertical wind shear (of order 1 m s(-1) km(-1)), then the dominant mode of a diabatically generated PV anomaly is a horizontally oriented dipole. The horizontal dipoles are typically of O(10 PVU), compared with the O(1 PVU) vertical dipoles that have been studied extensively throughout the literature. Furthermore, the horizontal PV dipoles are realized almost entirely as relative vorticity anomalies (on a time-scale of the order of tens of minutes after the heating has been turned on). The analysis of horizontal PV dipoles offers a new perspective on the vorticity dynamics of individual convective cells, implying that moist processes play a role in the maintenance of vertical vorticity in the convective storm environment.

Analysis of the eyes formed in simulated tropical cyclones and polar lows

The structure and size of the eyes generated in numerically simulated tropical cyclones and polar lows have been studied. A primitive-equation numerical model simulated systems in which the structures of the eyes formed were consistent with available observations. Whilst the tropical cyclone eyes generated were usually rapidly rotating, it appeared impossible for an eye formed in a system with a polar environment to develop this type of structure. The polar low eyes were found to be unable to warm through the subsidence of air with high values of potential temperature, as the environment was approximately statically neutral. Factors affecting the size of the eye were investigated through a series of controlled experiments. In mature tropical cyclone systems the size of the eye was insensitive to small changes in initial conditions, surface friction and latent and sensible heating from the ocean. In contrast, the eye size was strongly dependent on these parameters in the mature polar lows. Consistent with the findings, a mechanism is proposed in which the size of the eye in simulated polar lows is controlled by the strength of subsidence within it.