A BME solution of the inverse problem for saturated groundwater flow

In most real-world hydrogeologic situations, natural heterogeneity and measurement errors introduce major sources of uncertainty in the solution of the inverse problem. The Bayesian Maximum Entropy (BME) method of modern geostatistics offers an efficient solution to the inverse problem by first assimilating various physical knowledge bases (hydrologic laws, water table elevation data, uncertain hydraulic resistivity measurements, etc.) and then producing robust estimates of the subsurface variables across space. We present specific methods for implementing the BME conceptual framework to solve an inverse problem involving Darcys law for subsurface flow. We illustrate one of these methods in the case of a synthetic one-dimensional case study concerned with the estimation of hydraulic resistivity conditioned on soft data and hydraulic head measurements. The BME framework processes the physical knowledge contained in Darcys law and generates accurate estimates of hydraulic resistivity across space. The optimal distribution of hard and soft data needed to minimize the associated estimation error at a specified sampling cost is determined. This work was supported by grants from the National Institute of Environmental Health Sciences (Grant no. 5 P42 ES05948 and P30ES10126), the National Aeronautics and Space Administration (Grant no. 60-00RFQ041), the Army Research Office (Grant no. DAAG55-98-1-0289), and the National Science Foundation under Agreement No. DMS-0112069.     

Cosmic‐ray exposure history of two Frontier Mountain H‐chondrite showers from spallation and neutron‐capture products

Abstract— We measured the concentrations of ¹⁰Be, ²⁶Al, ³⁶Cl, ⁴¹Ca and ¹⁴C in the metal and/or stone fractions of 27 Antarctic chondrites from Frontier Mountain (FRO), including two large H‐chondrite showers. To estimate the pre‐atmospheric size of the two showers, we determined the contribution of neutron‐capture produced ³⁶Cl (half‐life = 3.01 times 10⁵ years) and ⁴¹Ca (1.04 times 10⁵ years) in the stone fraction. The measured activities of neutron‐capture ³⁶Cl and ⁴¹Ca, as well as spallation produced ¹⁰Be and ²⁶Al, were compared with Monte Carlo‐based model calculations. The largest shower, FRO 90174, includes eight fragments with an average terrestrial age of (100 ± 30) × 10³ years; the neutron‐capture saturation activities extend to 27 dpm/kg stone for ³⁶Cl and 19 dpm/kg stone for ⁴¹Ca. The concentrations of spallation produced ¹⁰Be, ²⁶Al and ³⁶Cl constrain the radius (R) to 80–100 cm, while the neutron‐capture ⁴¹Ca activities indicate that the samples originated from the outer 25 cm. With a pre‐atmospheric radius of 80–100 cm, FRO 90174 is among the largest of the Antarctic stony meteorites. The large pre‐atmospheric size supports our hypothesis that at least 50 of the ～150 classified H5/H6‐chondrites from the Frontier Mountain stranding area belong to this single fall; this hypothesis does not entirely account for the high H/L ratio at Frontier Mountain. The smaller shower, FRO 90001, includes four fragments with an average terrestrial age of (40 ± 10) × 10³ years; they contain small contributions of neutron‐capture ³⁶Cl, but no excess of ⁴¹Ca. FRO 90001 experienced a complex exposure history with high shielding conditions in the first stage (150 < R < 300 cm) and much lower shielding in the second stage (R < 30 cm), the latter starting ～1.0 million years (Ma) ago. Based on the measured ¹⁰Be/²¹Ne and ²⁶Al/²¹Ne ratios, the cosmic‐ray exposure ages of the two showers are 7.2 ± 0.5 Ma for FRO 90174 and 8 ± 1 Ma for FRO 90001. These ages coincide with the well‐established H‐chondrite peak and corroborate the observation that the exposure age distribution of FRO H‐chondrites is similar to that of non‐Antarctic falls. In addition, we found that corrections for neutron‐capture ³⁶Ar (from decay of ³⁶Cl) result in concordant ²¹Ne and ³⁸Ar exposure ages.     

Modeling the effects of riparian planting strategies on stream temperature: Increasing suitable habitat for endangered Formosan Landlocked Salmon in Shei‐Pa National Park,Taiwan

Reducing or stabilizing the stream temperature of ChiChiaWan Creek is a crucial work for Formosan Landlocked Salmon because ChiChiaWan Creek is the only one habitat for this endangered species. Planting trees in the riparian zone would be one of the alternatives. The purpose of this study was to evaluate the effects of several planting strategies on daily maximum stream temperature along the river. The results showed the effective vegetative shading angles should be more than 50° along ChiChiaWan Creek to reduce the direct solar radiation heating effectively. Upstream planting with 70° vegetative shading angle could be the most effective way among all the scenarios. However, this planting strategy could not improve the worst situations in summer because of the large solar elevation angles. The upstream planting in ChiChiaWan Creek was strongly recommended because the canopies could be easier to extend to totally cover the narrow width of river producing the most effective shades. Practicing the upstream planting with 90° vegetative shading angle can increase more than 1 km long suitable habitats for the endangered Salmon in summer. Alternatively, the west‐side planting scenario was the second effective way for temperature reduction. Our result provided a useful suggestion for the authorities in charge of saving the Formosan Landlocked Salmon, particularly under the stress of global warming. Copyright © 2011 John Wiley & Sons, Ltd.     

Site effects at Euroseistest—I. Determination of the valley structure and confrontation of observations with 1D analysis

This paper describes the process of construction of the 2D model of Volvi's geological structure and results of empirical and theoretical approaches to the evaluation of site response at Euroseistest. The construction of the 2D model is based on a re-interpretation of the available geophysical and geotechnical data in an effort to improve the definition of the subsoil structure at Euroseistest in terms of the most important parameters needed to model site response. The results of this re-interpretation are compared with a previous published 2D model of the same alluvial valley. Different analysis of the measurements and different criteria in the synthesis of data have led to a different model, even if both studies had access to the same field measurements. This underscores the fact that a model results of an interpretation and is not uniquely determined by the data, no matter how detailed they are. The well known subsoil structure opened the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from a detailed analysis of two events in frequency and time domains and 1D numerical modeling. The study of site response shows the important effect of the lateral variations on the ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In the case of Volvi's graben, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects in this basin.     

EISCAT/CRRES observations: nightside ionospheric ion outflow and oxygen-rich substorm injections

We present combined observations made near midnight by the EISCAT radar, all-sky cameras and the combined released and radiation efects satellite (CRRES) shortly before and during a substorm. In particular, we study a discrete, equatorward-drifting auroral arc, seen several degrees poleward of the onset region. The arc passes through the field-aligned beam of the EISCAT radar and is seen to be associated with a considerable upflow of ionospheric plasma. During the substorm, the CRRES satellite observed two major injections, 17 min apart, the second of which was dominated by O⁺ ions. We show that the observed are was in a suitable location in both latitude and MLT to have fed O⁺ ions into the second injection and that the upward flux of ions associated with it was sufficient to explain the observed injection. We interpret these data as showing that arcs in the nightside plasma-sheet boundary layer could be the source of O⁺ ions energised by a dipolarisation of the mid- and near-Earth tail, as opposed to ions ejected from the dayside ionosphere in the cleft ion fountain.     

Asymptotic Rayleigh instantaneous unit hydrograph

The instantaneous unit hydrograph for a channel network under general linear routing and conditioned on the network magnitude,N, tends asymptotically, asN grows large, to a Rayleigh probability density function. This behavior is identical to that of the width function of the network, and is proven under the assumption that the network link configuration is topologically random and the link hydraulic and geometric properties are independent and identically distributed random variables. The asymptotic distribution depends only on a scale factor, , where is a mean link wave travel time.     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper