Developing service-oriented applications in a grid environment

The Open-source Project for a Network Data Access Protocol (OPeNDAP) software framework has evolved over the last 10 years to become a robust, high performance, service oriented architecture for the access and transport of scientific data from a broad variety of disciplines, over the Internet. Starting with version 4.0 of the server release, Hyrax, has at its core, the Back-End Server (BES). The BES offers ease of programming, extensibility and reliability allowing added functionality especially related to high-performance data Grids, that goes beyond its original goal of distributed data access. We present the fundamentals of the BES server as a component based architecture, as well as our experiences using the BES to distribute data processing in Grid-oriented intensive parallel computational tasks.     

Effects of CO2 on benthic biota: An in situ benthic chamber experiment in Storfjorden (Norway)

Carbon capture and storage (CCS) methods, either sub-seabed or in ocean depths, introduces risk of CO₂ leakage and subsequent interaction with the ecosystem. It is therefore important to obtain information on possible effects of CO₂. In situ CO₂ exposure experiments were carried out twice for 10 days during 2005 using a Benthic Chamber system at 400 m depth in Storfjorden, Norway. pCO₂ in the water above the sediment in the chambers was controlled at approximately 500, 5000 and 20,000 μatm, respectively. This article describes the experiment and the results from measured the biological responses within the chamber sediments. The results show effects of elevated CO₂ concentrations on biological processes such as increased nanobenthos density. Methane production and sulphate reduction was enhanced in the approximately 5000 μatm chamber.     

The use of reference models from a priori data to guide 2D inversion of electrical resistivity tomography data

Clay-rich till plains cover much of the UK. Such sites are attractive locations for landfills, since clay aquitards lower the risk of landfill leachate entering groundwater. However, such tills often contain sand and gravel bodies that can act as leachate flow routes. Such bodies may not be detected by conventional site investigation techniques such as drilling boreholes and trial pitting. A method of guided inversion, where a priori data are used to construct structural reference models for use in inverting electrical resistivity tomography data, was proposed as a tool to improve the detection of sand and gravel bodies within clay-rich till sequences.
Following a successful 2D guided inversion synthetic modelling study, a field study was undertaken. Wenner 2D electrical resistivity tomography lines, resistivity cone penetrometry bores and electromagnetic induction ground resistivity data were collected over a site on the East Yorkshire coast, England, where sand and gravel lenses were known to exist from cliff exposures. A number of equally valid geoelectrical models were constructed using the electromagnetic and resistivity cone data. These were used as structural reference models in the inversion of the resistivity tomography data. Blind inversion using an homogenous reference model was also carried out for comparison.
It was shown for the first time that the best solution model produced by 2D inversion of one data set with a range of structural reference models could be determined by using the l ₂ model misfit between the solution models and associated reference models (reference misfit) as a proxy for the l ₂ misfit between the solution models and the synthetic model or 'best-guess' geoelectrical model (true misfit). The 2D methodology developed here is applicable in clay-rich till plains containing sand and gravel bodies throughout the UK.     

Estimating overland flow erosion capacity using unit stream power

Soil erosion caused by water flow is a complex problem. Both empirical and physically based approaches were used for the estimation of surface erosion rates. Their applications are mainly limited to experimental areas or laboratory studies. The maximum sediment concentration overland flow can carry is not considered in most of the existing surface erosion models. The lack of erosion capacity limitation may cause over estimations of sediment concentration. A correlation analysis is used in this study to determine significant factors that impact surface erosion capacity. The result shows that the unit stream power is the most dominant factor for overland flow erosion which is consistent with experimental data. A bounded regression formula is used to reflect the limits that sediment concentration cannot be less than zero nor greater than a maximum value. The coefficients used in the model are calibrated using published laboratory data. The computed results agree with laboratory data very well. A one dimensional overland flow diffusive wave model is used in conjunction with the developed soil erosion equation to simulate field experimental results. This study concludes that the non-linear regression method using unit stream power as the dominant factor performs well for estimating overland flow erosion capacity.     

Effects of Faraday rotation observed in filter magnetograph data

In this paper we analyze the effects of Faraday rotation on the azimuth of a transverse magnetic field as determined from the linear polarization in the inverse Zeeman effect. Observations of a simple sunspot were obtained with the Marshall Space Flight Center's vector magnetograph over the wavelength interval of 170 mÅ redward of line center of the Fe i 5250.22 Å spectral line to 170 mÅ to the blue, in steps of 10 mÅ. These data were analyzed to produce the variation of the azimuth as a function of wavelength at each pixel over the field of view of the sunspot. At selected locations in the sunspot, curves of the observed variation of azimuth with wavelength were compared with model calculations for the azimuth at each wavelength as derived from the inverse Zeeman effect modified by Faraday rotation. From these comparisons we derived the maximum amount of rotation as functions of both the magnitude and inclination of the sunspot's field. These results show that Faraday rotation of the azimuth will be a significant problem in observations taken near the center of a spectral line for fields as low as 1200 G and inclinations of the field in the range 20–80 deg. Conversely, they show that measurements taken in the wing of a spectral line are relatively free of the effects of Faraday rotation.     

Copper(II) humate mobility in kaolinite soil

An investigation of the influence of humate on the mobility of copper(II) ions in a kaolinite soil using leaching tests and electrokinetic experiments is reported. The data are interpreted in terms of humate–copper–clay interactions and humate electrical charge. Humate is mostly immobile below pH8 but is more mobile in alkaline conditions (sorption to kaolinite reduces its mobility in neutral conditions). Copper humate complexes are mobile in both acidic and alkaline conditions, but not in neutral conditions where they are sorbed. The dissolved copper humate complexes that form in acidic conditions are positively charged. The net effect of humate is to increase cupric ion mobility in kaolinite soil, especially in alkaline conditions.     

Exploring the sensitivity of coastal inundation modelling to DEM vertical error

As sea level is projected to rise throughout the twenty-first century due to climate change, there is a need to ensure that sea level rise (SLR) models accurately and defensibly represent future flood inundation levels to allow for effective coastal zone management. Digital elevation models (DEMs) are integral to SLR modelling, but are subject to error, including in their vertical resolution. Error in DEMs leads to uncertainty in the output of SLR inundation models, which if not considered, may result in poor coastal management decisions. However, DEM error is not usually described in detail by DEM suppliers; commonly only the RMSE is reported. This research explores the impact of stated vertical error in delineating zones of inundation in two locations along the Devon, United Kingdom, coastline (Exe and Otter Estuaries). We explore the consequences of needing to make assumptions about the distribution of error in the absence of detailed error data using a 1 m, publically available composite DEM with a maximum RMSE of 0.15 m, typical of recent LiDAR-derived DEMs. We compare uncertainty using two methods (i) the NOAA inundation uncertainty mapping method which assumes a normal distribution of error and (ii) a hydrologically correct bathtub method where the DEM is uniformly perturbed between the upper and lower bounds of a 95% linear error in 500 Monte Carlo Simulations (HBM+MCS). The NOAA method produced a broader zone of uncertainty (an increase of 134.9% on the HBM+MCS method), which is particularly evident in the flatter topography of the upper estuaries. The HBM+MCS method generates a narrower band of uncertainty for these flatter areas, but very similar extents where shorelines are steeper. The differences in inundation extents produced by the methods relate to a number of underpinning assumptions, and particularly, how the stated RMSE is interpreted and used to represent error in a practical sense. Unlike the NOAA method, the HBM+MCS model is computationally intensive, depending on the areas under consideration and the number of iterations. We therefore used the HBM+ MCS method to derive a regression relationship between elevation and inundation probability for the Exe Estuary. We then apply this to the adjacent Otter Estuary and show that it can defensibly reproduce zones of inundation uncertainty, avoiding the computationally intensive step of the HBM+MCS. The equation-derived zone of uncertainty was 112.1% larger than the HBM+MCS method, compared to the NOAA method which produced an uncertain area 423.9% larger. Each approach has advantages and disadvantages and requires value judgements to be made. Their use underscores the need for transparency in assumptions and communications of outputs. We urge DEM publishers to move beyond provision of a generalised RMSE and provide more detailed estimates of spatial error and complete metadata, including locations of ground control points and associated land cover.