iSERVO: Implementing the International Solid Earth Research Virtual Observatory by Integrating Computational Grid and Geographical Information Web Services

We describe the goals and initial implementation of the International Solid Earth Virtual Observatory (iSERVO). This system is built using a Web Services approach to Grid computing infrastructure and is accessed via a component-based Web portal user interface. We describe our implementations of services used by this system, including Geographical Information System (GIS)-based data grid services for accessing remote data repositories and job management services for controlling multiple execution steps. iSERVO is an example of a larger trend to build globally scalable scientific computing infrastructures using the Service Oriented Architecture approach. Adoption of this approach raises a number of research challenges in millisecond-latency message systems suitable for internet-enabled scientific applications. We review our research in these areas.     

Evaluation of an upscaled model for DNAPL dissolution kinetics in heterogeneous aquifers

Estimates of contaminant fluxes from DNAPL sources as a function of time and DNAPL mass reduction are important to assess the long-term sustainability and costs of monitored natural attenuation and to determine the benefits of partial source removal. We investigate the accuracy of the upscaled mass transfer function (MTF) proposed by Parker and Park [Parker JC, Park E. Modeling field-scale dense nonaqueous phase liquid dissolution kinetics in heterogeneous aquifers. WRR 2004;40:W05109] to describe field-scale dissolved phase fluxes from DNAPL sources for a range of scenarios generated using high-resolution 3-D numerical simulations of DNAPL infiltration and long-term dissolved phase transport. The results indicate the upscaled MTF is capable of accurately describing field-scale DNAPL dissolution rates as a function of time. For finger-dominated source regions, an empirical mass depletion exponent in the MTF takes on values greater than one which results in predicted mass flux rates that decrease continuously with diminishing DNAPL mass over time. Lens-dominated regions exhibit depletion exponents less than one, which results in more step-function like mass flux versus time behavior. Mass fluxes from DNAPL sources exhibiting both lens- and finger-dominated subregions were less accurately described by the simple MTF, but were well described by a dual-continuum model of the same form for each subregion. The practicality of calibrating a dual-continuum model will likely depend on the feasibility of obtaining spatially resolved field measurements of contaminant fluxes or concentrations associated with the subregions using multilevel sampling or some other means.     

Water quality in Lemoa, Guatemala

A variety of water supply sources are used in the rural Guatemalan Highlands. Formal municipal systems are rare, but villagers frequently form local water committees that construct spring-based supply systems. Recent work has indicated that groundwater may be a viable water supply source. The water quality from these two sources was characterized as well as water from other common sources including surface water and precipitation collection systems. Typically, all of the water sources contained unacceptable levels of coliform bacteria except for groundwater from a drilled well. Water monitoring indicated that not only did the water contain coliforms at the actual sources, but the infrastructure used to transmit the water to points of use may also be a source of coliforms. A cost comparison indicates that groundwater may be a cost-competitive, higher quality alternative to traditional spring-fed water systems.     

Differential supply of autochthonous organic carbon and nitrogen to the microbial loop in the Delaware Estuary

Using stable isotope tracer techniques in 4-h bottle incubations, the importance of organic matter transfer from phytoplankton to heterotrophic bacteria (bacteria) has been re-evaluated in the Delaware Estuary, considering carbon (C) and nitrogen (N) cycles separately. The hypothesis is that the transfer of C and N from phytoplankton to bacteria varies both temporally and spatially along estuarine gradients in response to variation in factors such as terrestrial organic C supply, inorganic N speciation and concentrations, and extracellular release of dissolved organic matter by phytoplankton. The percentage of autochthonous dissolved organic C being assimilated by bacteria varied between 3% and 10% of primary production and was not related to the rate of primary production. The transfer of N was considerably more variable when compared to C transfer, averaging ca. 20% of phytoplankton N assimilation; individual experiments yielded rates as high as 50%. Unlike C, autochthonous dissolved organic N transfer appears to vary with the magnitude of primary production, and its assimilation by bacteria accounted for 0–56% of the total measured bacterial N uptake. The results highlight the importance of separate consideration of C and N elemental cycles in evaluating sources of organic matter to the estuarine microbial loop.     

Reducing the bias of multitaper spectrum estimates

The power spectral density of geophysical signals provides information about the processes that generated them. We present a new approach to determine power spectra based on Thomson's multitaper analysis method. Our method reduces the bias due to the curvature of the spectrum close to the frequency of interest. Even while maintaining the same resolution bandwidth, bias is reduced in areas where the power spectrum is significantly quadratic. No additional sidelobe leakage is introduced. In addition, our methodology reliably estimates the derivatives (slope and curvature) of the spectrum. The extra information gleaned from the signal is useful for parameter estimation or to compare different signals.     

Enhancing the human benefits of flood warnings

This article evaluates some of the factors which limit the human benefits of hazard warnings, with specific reference to flood warnings, and we conclude by suggesting ways of enhancing these benefits. We focus mainly upon the economic benefits generated by flood damage savings by households that warnings facilitate; health effects of flooding and flood warnings; and the effects of warnings on loss of life and physical injury. Our results, based partly upon surveys of flooded households, reveal that economic benefits are currently more limited than we previously thought, but that for several reasons these benefits are likely to be under-estimated. We argue that the intangible benefits to public health, safety and security must also be taken into account in decisions about investment in flood warnings. In England and Wales, the public’s response to flood warnings is currently low and is a key benefit-limiting factor which could begin to undermine a recent major shift in national flood risk management policy towards a more people-centred, portfolio approach in which changing human behaviour is viewed as important. Using a trans-disciplinary approach, we discuss the evidence and literature surrounding this poor response, and suggest a number of ways in which the issue may be addressed in future.     

Modelling potential impacts of bottom trawl fisheries on soft sediment biogeochemistry in the North Sea

Bottom trawling causes physical disturbance to sediments particularly in shelf areas. The disturbance due to trawling is most significant in deeper areas with softer sediments where levels of natural disturbance due to wave and tidal action are low. In heavily fished areas, trawls may impact the same area of seabed more than four times per year. A single pass of a beam trawl, the heaviest gear routinely used in shelf sea fisheries, can kill 5–65% of the resident fauna and mix the top few cm of sediment. We expect that sediment community function, carbon mineralisation and biogeochemical fluxes will be strongly affected by trawling activity because the physical effects of trawling are equivalent to those of an extreme bioturbator, and yet, unlike bioturbating macrofauna, trawling does not directly contribute to community metabolism. We used an existing box-model of a generalised soft sediment system to examine the effects of trawling disturbance on carbon mineralisation and chemical concentrations. We contrasted the effects of a natural scenario, where bioturbation is a function of macrobenthos biomass, with an anthropogenic impact scenario where physical disturbance results from trawling rather than the action of bioturbating macrofauna. Simulation results suggest that the effects of low levels of trawling disturbance will be similar to those of natural bioturbators but that high levels of trawling disturbance prevent the modelled system from reaching equilibrium due to large carbon fluxes between oxic and anoxic carbon compartments. The presence of macrobenthos in the natural disturbance scenario allowed sediment chemical storage and fluxes to reach equilibrium. This is because the macrobenthos are important carbon consumers in the system whose presence reduces the magnitude of available carbon fluxes. In soft sediment systems, where the level physical disturbance due to waves and tides is low, model results suggest that intensive trawling disturbance could cause large fluctuations in benthic chemical fluxes and storage.     

Radial variation of differential rotation in the solar electron corona

Long-lived brightness structures in the solar electron corona persist over many solar rotation periods and permit an observational determination of coronal magnetic tracer rotation as a function of latitude and height in the solar atmosphere. For observations over 1964–1976 spanning solar cycle 20, we compare the latitude dependence of rotation at two heights in the corona. Comparison of rotation rates from East and West limbs and from independent computational procedures is used to estimate uncertainty. Time-averaged rotation rates based on three methods of analysis demonstrate that, on average, coronal differential rotation decreases with height from 1.125 to 1.5 R _S. The observed radial variation of differential rotation implies a scale height of approximately 0.7 R _S for coronal differential rotation.Model calculations for a simple MHD loop show that magnetic connections between high and low latitudes may produce the observed radial variations of magnetic tracer rotation. If the observed tracer rotation represents the rotation of open magnetic field lines as well as that of closed loops, the small scale height for differential rotation suggests that the rotation of solar magnetic fields at the base of the solar wind may be only weakly latitude dependent. If, instead, closed loops account completely for the radial gradients of rotation, outward extrapolation of electron coronal rotation may not describe magnetic field rotation at the solar wind source. Inward extrapolations of observed rotation rates suggest that magnetic field and plasma are coupled a few hundredths of a solar radius beneath the photosphere.     

Basic properties of magnetic flux tubes and restrictions on theories of solar activity

It is shown that the mean longitudinal field in a magnetic flux tube is reduced, rather than enhanced, by twisting the tube to form a rope. It is shown that there is no magnetohydrostatic equilibrium when one twisted rope is wound around another. Instead there is rapid line cutting (neutral point annihilation). It is shown that the twisting increases, and the field strength decreases, along a flux tube extending upward through a stratified atmosphere.These facts are at variance with Piddington's recent suggestion that solar activity is to be understood as the result of flux tubes which are enormously concentrated by twisting, which consist of several twisted ropes wound around each other, and which came untwisted where they emerge through the photosphere.This work was supported in part by the National Aeronautics and Space Administration under Grant NGL 14-001-001.     

Upward and downward continuation as inverse problems

Summary . The formalism of Backus & Gilbert is applied to the problems of upward and downward continuation of harmonic functions. We first treat downward continuation of a two-dimensional field to a level surface everywhere below the observation locations; the calculation of resolving widths and solution estimates is a straightforward application of Backus—Gilbert theory. The extension to the downward continuation of a three-dimensional field uses a delta criterion giving resolving areas rather than widths. A feature not encountered in conventional Backus—Gilbert problems is the requirement of an additional constraint to guarantee the existence of the resolution integrals. Finally, we consider upward continuation of a two-dimensional field to a level above all observations. We find that solution estimates must be weighted averages of the field not only on this level, but also on a line passing between the observations and sources. Weighting on the lower line may be traded off against resolution on the upper level.