Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion

Nonlocal moment equations allow one to render deterministically optimum predictions of flow in randomly heterogeneous media and to assess predictive uncertainty conditional on measured values of medium properties. We present a geostatistical inverse algorithm for steady-state flow that makes it possible to further condition such predictions and assessments on measured values of hydraulic head (and/or flux). Our algorithm is based on recursive finite-element approximations of exact first and second conditional moment equations. Hydraulic conductivity is parameterized via universal kriging based on unknown values at pilot points and (optionally) measured values at other discrete locations. Optimum unbiased inverse estimates of natural log hydraulic conductivity, head and flux are obtained by minimizing a residual criterion using the Levenberg-Marquardt algorithm. We illustrate the method for superimposed mean uniform and convergent flows in a bounded two-dimensional domain. Our examples illustrate how conductivity and head data act separately or jointly to reduce parameter estimation errors and model predictive uncertainty.This work is supported in part by NSF/ITR Grant EAR-0110289. The first author was additionally supported by scholarships from CONACYT and Instituto de Investigaciones Electricas of Mexico. Additional support was provided by the European Commission under Contract EVK1-CT-1999-00041 (W-SAHaRA-Stochastic Analysis of Well Head Protection and Risk Assessment).     

Diurnal variation of the dayside, ionospheric, mid-latitude trough in the southern hemisphere at 800 km: Model and measurement comparison

Our high latitude ionospheric model predicts the existence of a pronounced “dayside” trough in plasma concentration equatorward of the auroral oval in both the Northern and Southern Hemispheres for solar maximum, winter, and low geomagnetic activity conditions. The trough in the Southern Hemisphere is much deeper than that in the Northern Hemisphere, with the minimum trough density at 800 km being 2 × 10³ cm⁻³ in the Southern Hemisphere and 10⁴ cm⁻³ in the Northern Hemisphere. The dayside trough has a strong longitudinal (diurnal) dependence and appears between 11:00 and 19:00 U.T. in the Southern Hemisphere and between 02:00 and 08:00 U.T. in the Northern Hemisphere. This dayside trough is a result of the auroral oval moving to larger solar zenith angles at those universal times when the magnetic pole is on the antisunward side of the geographic pole. As the auroral ionization source moves to higher geographic latitudes, it leaves a region of declining photoionization on the dayside. For low convection speeds, the ionosphere decays and a dayside trough forms. The trough is deeper in the Southern Hemisphere than in the Northern Hemisphere because of the greater offset between the geomagnetic and geographic poles. Satellite data taken in both the Northern and Southern Hemispheres confirm the gross features of the dayside trough, including its strong longitudinal dependence, its depth, and the asymmetry between the Northern and Southern Hemisphere troughs.     

A smoothed particle hydrodynamics modelling of soil–water mixing and resulting changes in average strength

Soil–water interaction is a pivotal process in many underwater geohazards such as underwater landslides where soil sediments gradually evolve into turbidity currents after interactions with ambient water. Due to the large deformations, multiphase interactions and phase changes this involves, investigations from numerical modelling of the transition process have been limited so far. This study explores a simple numerical replication of such soil–water mixing with respect to changes in average strength using smoothed particle hydrodynamics (SPH). A uniform viscoplastic model is used for both the solid-like and fluid-like SPH particles. The proposed numerical solution scheme is verified by single-phase dam break tests and multiphase simple shear tests. SPH combinations of solid-like and fluid-like particles can replicate the clay–water mixture as long as the liquidity index of the solid-like particles is larger than unity. The proposed numerical scheme is shown to capture key features of an underwater landslide such as hydroplaning, water entrainment and wave generation and thus shows promise as a tool to simulate the whole process of subaquatic geohazards involving solid–fluid transition during mass transport.     

A quantitative analysis on the Egyptian Mediterranean waters

ＡｑｕａｎｔｉｔａｔｉｖｅａｎａｌｙｓｉｓｏｎｔｈｅＥｇｙｐｔｉａｎＭｅｄｉｔｅｒｒａｎｅａｎｗａｔｅｒｓ￥Ｍ．Ａ．ＳａｉｄａｎｄＦ．Ｍ．Ｅｉｄ（ＲｅｃｅｉｖｅｄＳｅｐｔｅｍｂｅｒ２１，１９９３１ａｃｃｆｅｐｔｅｄＪａｎｕａｒｙ１５，１９９４）Ａｂｓ...     

Radiation Products in Processed Ices Relevant to Edgeworth-Kuiper-Belt Objects

Near the inner edge of the Edgeworth-Kuiper Belt (EKB) are Pluto and Charon, which are known to have N₂- and H₂O-dominated surface ices, respectively. Such non-polar and polar ices, and perhaps mixtures of them, also may be present on other trans-Neptunian objects. Pluto, Charon, and all EKB objects reside in a weak, but constant UV-photon and energetic ion radiation environment that drives chemical reactions in their surface ices. Effects of photon and ion processing include changes in ice composition, volatility, spectra, and albedo, and these have been studied in a number of laboratories. This paper focuses on ice processing by ion irradiation and is aimed at understanding the volatiles, ions, and residues that may exist on outer solar system objects. We summarize radiation chemical products of N₂-rich and H₂O-rich ices containing CO or CH₄, including possible volatiles such as alcohols, acids, and bases. Less-volatile products that could accumulate on EKB objects are observed to form in the laboratory from acid-base reactions, reactions promoted by warming, or reactions due to radiation processing of a relatively pure ice (e.g., CO → C₃O₂). New IR spectra are reported for the 1–5 mu;m region, along with band strengths for the stronger features of carbon suboxide, carbonic acid, the ammonium and cyanate ions, polyoxymethylene, and ethylene glycol. These six materials are possible contributors to EKB surfaces, and will be of interest to observers and future missions.     

Chemical and mineralogical characterization of the Mineo (Sicily,Italy) pallasite: A unique sample

The Mineo pallasite is characterized here for the first time. The only 42 g still available worldwide is part of the collection of the Department of Physics and Geology, University of Perugia. A multianalytical approach was used, joining field-emission scanning electron microscopy, Raman analysis, X-ray powder diffraction, electron-probe microanalysis, and laser ablation inductively coupled plasma mass spectrometry. Results highlighted that (1) the Mineo pallasite belongs to the Main Group pallasites; (2) the silicate component is essentially olivine, with no pyroxene component; (3) the olivine chemical composition varies in terms of both iron and trace elements; (4) the metal phase is essentially kamacite with the taenite mainly found in the plessite structure; (5) phosphide phases are present as schreibersite and barringerite. The observed compositional variability in olivines as well as their occurrence as both angular and rounded crystals suggest that the Mineo pallasite could have been derived from a large impact of a differentiated parent body with a larger solid body. The resulting pallasite conglomerate consists of the compositionally different olivines, likely coming from different areas of the same differentiated parent body, and the residual molten Fe-Ni.     

Seasonal temperatures for the past ∼400 years reconstructed from diatom and chironomid assemblages in a high-altitude lake (Lej da la Tscheppa,Switzerland)

We analysed a 42 cm long sediment record from Lej da la Tscheppa, a high-altitude lake (2,616 m a.s.l.) in the Upper Engadine valley (Switzerland) for subfossil diatoms, chironomids and pollen. The chronology of the top 21 cm of the record was established using ²¹⁰Pb analysis using a constant-rate-of-supply model, and validated with ¹³⁷Cs measurements and the content of spheroidal carbonaceous particles. A tentative chronology for the lower part of the core was obtained through extrapolation of the sedimentation rates in the uppermost part of the record. Pollen assemblages in the record reflect regional changes in forestation and land-use patterns in the Upper Engadine valley and show no evidence of significant local human activity in the lake’s catchment. Diatom assemblages record a distinct increase in planktonic taxa since the early 19th century, suggesting a decrease in the duration of ice-cover. In contrast, chironomid assemblages remained stable during a large part of the record. We applied an established chironomid-based July air temperature transfer function and a newly developed diatom-based spring air temperature transfer function to reconstruct past seasonal air temperature changes at Lej da la Tscheppa. The reconstructions indicate a diatom-inferred warming trend in spring temperatures during the past ca. 400 years, whereas chironomid-inferred summer temperatures suggest a slight cooling trend. These biota-based reconstructions are in good agreement with the centennial-scale temperature trend in an independent reconstruction of regional temperatures in the Upper Engadine region based on instrumental records and documentary proxy evidence from the Alps. Our results suggest that, in high-altitude lakes, independent chironomid- and diatom-based seasonal temperature reconstruction is possible and can be successfully used to track seasonal temperature trends.     

Earthen barriers to control lava flows in the 2001 eruption of Mt. Etna

Preceded by four days of intense seismicity and marked ground deformation, a new eruption of Mt. Etna started on 17 July and lasted until 9 August 2001. It produced lava emission and strombolian and phreatomagmatic activity from four different main vents located on a complex fracture system extending from the southeast summit cone for about 4.5 km southwards, from 3000 to 2100 m elevation (a.s.l.). The lava emitted from the lowest vent cut up an important road on the volcano and destroyed other rural roads and a few isolated country houses. Its front descended southwards to about 4 km distance from the villages of Nicolosi and Belpasso. A plan of intervention, including diversion and retaining barriers and possibly lava flow interruption, was prepared but not activated because the flow front stopped as a consequence of a decrease in the effusion rate. Extensive interventions were carried out in order to protect some important tourist facilities of the Sapienza and Mts. Silvestri zones (1900 m elevation) from being destroyed by the lava emitted from vents located at 2700 m and 2550 m elevation. Thirteen earthen barriers (with a maximum length of 370 m, height of 10–12 m, base width of 15 m and volume of 25 000 m³) were built to divert the lava flow away from the facilities towards a path implying considerably less damage. Most of the barriers were oriented diagonally (110–135°) to the direction of the flow. They were made of loose material excavated nearby and worked very nicely, resisting the thrust of the lava without any difficulty. After the interventions carried out on Mt. Etna in 1983 and in 1991–1992, those of 2001 confirm that earthen barriers can be very effective in controlling lava flows.