Recurrent Red-tides in the Southampton Water Estuary Caused by the Phototrophic CiliateMesodinium rubrum

In the Southampton Water estuary (southern England, U.K.), red-tides caused by the planktonic, phototrophic ciliateMesodinium rubrum(=Myrionecta rubra) occur during most summers and sometimes in autumn. These events were investigated in detail between 1985 and 1987 and were characterized by levels of chlorophylla(chl a) of over 100 μg l⁻¹, cell numbers ofM. rubrumof over 1×10³ ml⁻¹, oxygen saturations of around 150%, and depleted numbers of macrozooplankton. Initiation of red-water did not appear to be triggered by irradiance or nutrients, but coincided with an increase in temperature and water column stability. This enhanced stability was promoted by increased surface to bottom gradients of both temperature and salinity, and by reduced mixing during neap tides. Development of red-water was accompanied by removal of most of the dissolved NH⁺₄from the water column, whereas some NO⁻₃persisted, presumably maintained by freshwater input. NO⁻₃and NH⁺₄gradually returned to pre-bloom concentrations as the red-water declined in late summer. Maximal biomass ofM. rubrumappeared to be limited by irradiance, and self-shading probably imposed an upper limit of around 300 mg chl a m⁻²within the water column. At the observed levels of chl a, irradiance values within the population maximum between 1 and 3 m depth were only just of the order (≈15 μmol photons m⁻² s⁻¹) required to balance estimated respiratory demands. Oxygen concentration became undersaturated during the late bloom phase, with minimal values of 20–30% saturation recorded in deeper waters; however, despite this and reduced numbers of macrozooplankton, direct deleterious effects on other organisms were not observed.     

Simulating methane emission from a Chinese rice field as influenced by fertilizer and water level

Natural and agricultural wetlands are considered to be the major sources of global atmospheric methane (CH₄). A one‐dimensional model was developed to simulate methane emission and used to examine the influence of various physical processes on the rate of methane emission. Three processes involved in the methane emission are implemented in the model: production, reoxidation and transport. Three transport pathways were considered: diffusion across water–air or soil–air interfaces, ebullition and diffusion through plants. These pathways are influenced by soil properties, plant growth, water‐table conditions, temperature and external inputs (e.g. fertilizer). The model was used to examine the seasonal variation of the methane emission at a rice field in Hunan, China, which was observed during a field experiment for consecutive (early and late) rice seasons in 1992. The observed seasonal variations of methane emission, and role of plants in transporting methane to the atmosphere, are captured by the model simulation. Further model applications were conducted to simulate effects of fertilizer and water‐level condition on the methane emission. The results indicate that unfermented organic fertilizer produces a higher methane emission rate than mineral fertilizer. The simulations with treatments of a deep‐water covering and constant moisture reduced the methane emission. The rice field study provides a framework for further development of the model towards simulations based on spatially distributed variables (e.g. water table, soil temperature and vegetation) at a regional scale. Copyright © 2003 John Wiley & Sons, Ltd.     

Upscaling of elastic properties for large scale geomechanical simulations

Large scale geomechanical simulations are being increasingly used to model the compaction of stress dependent reservoirs, predict the long term integrity of under‐ground radioactive waste disposals, and analyse the viability of hot‐dry rock geothermal sites. These large scale simulations require the definition of homogenous mechanical properties for each geomechanical cell whereas the rock properties are expected to vary at a smaller scale. Therefore, this paper proposes a new methodology that makes possible to define the equivalent mechanical properties of the geomechanical cells using the fine scale information given in the geological model. This methodology is implemented on a synthetic reservoir case and two upscaling procedures providing the effective elastic properties of the Hooke's law are tested. The first upscaling procedure is an analytical method for perfectly stratified rock mass, whereas the second procedure computes lower and upper bounds of the equivalent properties with no assumption on the small scale heterogeneity distribution. Both procedures are applied to one geomechanical cell extracted from the reservoir structure. The results show that the analytical and numerical upscaling procedures provide accurate estimations of the effective parameters. Furthermore, a large scale simulation using the homogenized properties of each geomechanical cell calculated with the analytical method demonstrates that the overall behaviour of the reservoir structure is well reproduced for two different loading cases. Copyright © 2004 John Wiley & Sons, Ltd.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Contribution of storms to groundwater recharge in the semi‐arid region of Karnataka,India

This paper describes the application of environmental isotopes and injected tracer techniques in estimating the contribution of storms as well as annual precipitation to groundwater recharge and its circulation, in the semi‐arid region of Bagepalli, Kolar district, Karnataka. Environmental isotopes ²H, ¹⁸O and ³H were used to study the effect of storms on the hydrological system, and an isotope balance was used to compute the contribution of a storm component to the groundwater. Some of the groundwater samples collected during the post‐storm periods were highly depleted in stable isotope content with higher deuterium excess relative to groundwater from the pre‐storm periods. Significant variation in deuterium excess in groundwater from the same area, collected in two different periods, indicates the different origin of air masses. The estimated recharge component of a storm event of 600 mm to the groundwater was found to be in the range of 117–165 mm. There was no significant variation in environmental tritium content of post‐storm and pre‐storm groundwater, indicating the fast circulation of groundwater in the system. After completion of the environmental isotope work, an injected radiotracer ³H technique was applied to estimate the direct recharge of total precipitation to the groundwater. The estimated recharge to the groundwater is 33 mm of the 550 mm annual precipitation during 1992. Copyright © 2003 John Wiley & Sons, Ltd.     

Estimation of snow water equivalent using microwave radiometry over Arctic first‐year sea ice

The magnitude and spatial distribution of snow on sea ice are both integral components of the ocean–sea‐ice–atmosphere system. Although there exists a number of algorithms to estimate the snow water equivalent (SWE) on terrestrial surfaces, to date there is no precise method to estimate SWE on sea ice. Physical snow properties and in situ microwave radiometry at 19, 37 and 85 GHz, V and H polarization were collected for a 10‐day period over 20 first‐year sea ice sites. We present and compare the in situ physical, electrical and microwave emission properties of snow over smooth Arctic first‐year sea ice for 19 of the 20 sites sampled. Physical processes creating the observed vertical patterns in the physical and electrical properties are discussed. An algorithm is then developed from the relationship between the SWE and the brightness temperature measured at 37 GHz (55°) H polarization and the air temperature. The multiple regression between these variables is able to account for over 90% of the variability in the measured SWE. This algorithm is validated with a small in situ data set collected during the 1999 field experiment. We then compare our data against the NASA snow thickness algorithm, designed as part of the NASA Earth Enterprise Program. The results indicated a lack of agreement between the NASA algorithm and the algorithm developed here. This lack of agreement is attributed to differences in scale between the Special Sensor Microwave/Imager and surface radiometers and to differences in the Antarctic versus Arctic snow physical and electrical properties. Copyright © 2003 John Wiley & Sons, Ltd.     

One hundred years of mountain hydrology in Switzerland by the WSL

In 1903 the Swiss Federal Research Institute WSL started its first forest hydrology measurements with the aim to deliver a sound scientific basis for the implementation of new forest legislation introduced in Switzerland in 1876. This legislation was triggered by several large floods that occurred in Switzerland, for which a major cause was widely seen as the poor condition of forests at that time. Consequently, hydrologic research at WSL first focused on the influence of forests on floods. In the second half of the 20th century, other hydrological issues such as water quality, snow hydrology and sediment transport complemented the hydrologic research at WSL. Some recent results of this work are presented in three papers joining this introductory paper to mark the 100th anniversary of hydrologic research at WSL. Copyright © 2006 John Wiley & Sons, Ltd.     

A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Hysteresis‐based analysis of overland metal transport

Introducing a concept of equivalent mass depth of flow, this study describes the phenomenon of non‐point source pollutant (metal) transport for pavement (or overland) flow in analogy with wave propagation in wide open channels. Hysteretic and normal mass rating curves are developed for runoff rate and mass of 12 dissolved and particulate‐bound metal elements (pollutants) using the rainfall‐runoff and water quality data of the 15 × 20 m² instrumented pavement in Cincinnati, USA. Normal mass rating curves developed for easy computation of pollutant load are found to be of a form similar to Manning's, which is valid for open channel flows. Based on the hysteresis analysis, wave types for dissolution and mixing of particulate‐bound metals are identified. The analysis finds that the second‐order partial‐differential equation normally used for metal transport does not have the efficacy to describe fully the strong non‐linear phenomena such as is described for various metal elements by dynamic waves. In addition, the proportionality concept of the popular SCS‐CN concept is extended for determining the potential maximum metal mass M_p of all the 12 elements transported by a rain storm and related to the antecedent dry period (ADP). For the primary metal zinc element, M_p is found to increase with the ADP and vice versa. Copyright © 2003 John Wiley & Sons, Ltd.