The value of subsidence data in ground water model calibration

The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data.     

Assessment of chemical dispersant effectiveness in a wave tank under regular non-breaking and breaking wave conditions

Current chemical dispersant effectiveness tests for product selection are commonly performed with bench-scale testing apparatus. However, for the assessment of oil dispersant effectiveness under real sea state conditions, test protocols are required to have hydrodynamic conditions closer to the natural environment, including transport and dilution effects. To achieve this goal, Fisheries and Oceans Canada and the US Environmental Protection Agency (EPA) designed and constructed a wave tank system to study chemical dispersant effectiveness under controlled mixing energy conditions (regular non-breaking, spilling breaking, and plunging breaking waves). Quantification of oil dispersant effectiveness was based on observed changes in dispersed oil concentrations and oil-droplet size distribution. The study results quantitatively demonstrated that total dispersed oil concentration and breakup kinetics of oil droplets in the water column were strongly dependent on the presence of chemical dispersants and the influence of breaking waves. These data on the effectiveness of dispersants as a function of sea state will have significant implications in the drafting of future operational guidelines for dispersant use at sea.     

Characterization of Damage in Sandstones along the Mojave Section of the San Andreas Fault: Implications for the Shallow Extent of Damage Generation

Following theoretical calculations that suggest shallow generation of rock damage during an earthquake rupture, we measure the degree of fracture damage in young sedimentary rocks from the Juniper Hills Formation (JHF) that were displaced 21 km along the Mojave section of San Andreas Fault (SAF) and were not exhumed significantly during their displacement. In exposures adjacent to the fault, the JHF typically displays original sedimentary fabrics and little evidence of bulk shear strain at the mesoscopic scale. The formation is, however, pervasively fractured at the microscopic scale over a zone that is about a 100 m wide on the southwest side of the SAF near Little Rock. The abundance of open fractures, the poor consolidation, and the shallow inferred burial depth imply that the damage was generated close to the surface of the Earth. The spatial correlation of this damage with a seismically active trace of the SAF suggests that it was generated by SAF slip events that by assumption were of a seismic nature throughout the displacement history of the JHF. Thus the JHF provides a very shallow upper bound for the generation of brittle damage in a seismic fault zone. The fracture fabric is characterized by preferred orientations of fractures that split grains between contact points and is consistent with overall deformation under directed compression. However, the available results cannot be used to distinguish between proposed off-fault damage mechanisms. Fracture orientations are compatible with a maximum compressive stress oriented at a high angle to the fault at about 10 m, and at a lower, more variable angle farther away from the fault. The fracture distribution and fabric are consistent with observations made of the microscale damage characteristics of the Hungry Valley Formation in the northwestern section of the SAF in the Mojave, and with previous observations of exhumed, ancestral strands of the SAF.     

Kinematic transformations for planar multi‐directional pseudodynamic testing

The pseudodynamic (PSD) test method imposes command displacements to a test structure for a given time step. The measured restoring forces and displaced position achieved in the test structure are then used to integrate the equations of motion to determine the command displacements for the next time step. Multi‐directional displacements of the test structure can introduce error in the measured restoring forces and displaced position. The subsequently determined command displacements will not be correct unless the effects of the multi‐directional displacements are considered. This paper presents two approaches for correcting kinematic errors in planar multi‐directional PSD testing, where the test structure is loaded through a rigid loading block. The first approach, referred to as the incremental kinematic transformation method, employs linear displacement transformations within each time step. The second method, referred to as the total kinematic transformation method, is based on accurate nonlinear displacement transformations. Using three displacement sensors and the trigonometric law of cosines, this second method enables the simultaneous nonlinear equations that express the motion of the loading block to be solved without using iteration. The formulation and example applications for each method are given. Results from numerical simulations and laboratory experiments show that the total transformation method maintains accuracy, while the incremental transformation method may accumulate error if the incremental rotation of the loading block is not small over the time step. A procedure for estimating the incremental error in the incremental kinematic transformation method is presented as a means to predict and possibly control the error. Copyright © 2009 John Wiley & Sons, Ltd.     

Time–space continuity of daily maps of fractional snow cover and albedo from MODIS

Using reflectance values from the seven MODIS “land” bands with 250 or 500 m resolution, along with a corresponding cloud product, we estimate the fraction of each 500 m pixel that snow covers, along with the albedo of that snow. The daily products have data gaps and errors because of cloud cover and sensor viewing geometry. Rather than make users interpolate and filter these patchy daily maps without completely understanding the retrieval algorithm and instrument properties, we use the daily time series to improve the estimate of the measured snow properties for a particular day. We use a combination of noise filtering, snow/cloud discrimination, and interpolation and smoothing to produce our best estimate of the daily snow cover and albedo. We consider two modes: one is the “predictive” mode, whereby we estimate the snow-covered area and albedo on that day using only the data up to that day; the other is the “retrospective” mode, whereby we reconstruct the history of the snow properties for a previous period.     

The last 40 ka tephrostratigraphic record of Lake Ohrid,Albania and Macedonia: a very distal archive for ash dispersal from Italian volcanoes

A 1075 cm long core (Lz1120) was recovered in the south-eastern part of the Lake Ohrid (Republics of Macedonia and Albania) and sampled for identification of tephra layers. Magnetic susceptibility investigations show rather high magnetic values throughout the core, with peaks unrelated to the occurrence of tephra layers but instead to the relative abundance of detrital magnetic minerals in the sediment. Naked-eye inspection of the core allowed us to identify of two tephra layers, at 896–897 cm and 1070–1075 cm. Laboratory inspection of the grain-size fraction > 125 μm allowed for the identification of a third cryptotephra at 310–315 cm. Major element analyses on glass shards of the tephra layers at 896–897 cm and 1070–1075 cm show a trachytic composition, and indicate a correlation with the regionally dispersed Y-3 and Y-5 tephra layers, dated at ca 30 and 39 cal ka BP. The cryptotephra at 310–315 cm has a mugearitic–benmoreitic composition, and was correlated with the FL eruption of Mt. Etna, dated at 3370 ± 70 cal yr BP. These ages are in agreement with five ¹⁴C AMS measurements carried out on plant remains and macrofossils from the lake sediments at different depths along the core.     

Fe,Mn, P and S Speciation in Sediments from the Capivara Hydroelectric Dam Lake (Brazil) as an Indicator of Anthropogenic Influences

The conservation of large water resources is essential for the preservation of human life. The quantification and, more importantly, the speciation of chemical substances that indicate the presence of anthropogenic contamination in water resources are of great importance. This paper presents the results of analysis for the determination of organic, inorganic and total phosphorus, pseudo‐sulfur, and iron and manganese, in five fractions, in water sediments collected from the Capivara Hydroelectric Power Plant, Brazil. A study on the seasonal variation of these parameters was conducted, with data having been collected in the winter and in the summer, at two sites along the dam, 5 km apart, close to the city of Primeiro de Maio. Phosphorous was found in sediments and adjacent soil in the organic form (OP), and was used as an indicator of anthropogenic influence on the reservoir banks. Speciation of potentially toxic Mn showed that it is present in the exchangeable fraction of the 0–5 cm depth layer (sediment/water interface), making its transfer to the water column possible. Results from this study showed that domestic and industrial effluent treatment measures are needed for the preservation of the quality of aquatic environments.     

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.     

The November 2002 eruption of Piton de la Fournaise,Réunion: tracking the pre-eruptive thermal evolution of magma using melt inclusions

The November 2002 eruption of Piton de la Fournaise in the Indian Ocean was typical of the activity of the volcano from 1999 to 2006 in terms of duration and volume of magma ejected. The first magma erupted was a basaltic liquid with a small proportion of olivine phenocrysts (Fo₈₁) that contain small numbers of melt inclusions. In subsequent flows, olivine crystals were more abundant and richer in Mg (Fo_83–84). These crystals contain numerous melt and fluid inclusions, healed fractures, and dislocation features such as kink bands. The major element composition of melt inclusions in this later olivine (Fo_83–84) is out of equilibrium with that of its host as a result of extensive post-entrapment crystallization and Fe²⁺ loss by diffusion during cooling. Melt inclusions in Fo₈₁ olivine are also chemically out of equilibrium with their hosts but to a lesser degree. Using olivine–melt geothermometry, we determined that melt inclusions in Fo₈₁ olivine were trapped at lower temperature (1,182 ± 1°C) than inclusions in Fo_83–84 olivine (1,199–1,227°C). This methodology was also used to estimate eruption temperatures. The November 2002 melt inclusion compositions suggest that they were at temperatures between 1,070°C and 1,133°C immediately before eruption and quenching. This relatively wide temperature range may reflect the fact that most of the melt inclusions were from olivine in lava samples and therefore likely underwent minor but variable amounts of post-eruptive crystallization and Fe²⁺ loss by diffusion due to their relatively slow cooling on the surface. In contrast, melt inclusions in tephra samples from past major eruptions yielded a narrower range of higher eruption temperatures (1,163–1,181°C). The melt inclusion data presented here and in earlier publications are consistent with a model of magma recharge from depth during major eruptions, followed by storage, cooling, and crystallization at shallow levels prior to expulsion during events similar in magnitude to the relatively small November 2002 eruption.