A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.     

3d transition metal L-edge X-ray absorption studies of the dichalcogenides of Fe, Co and Ni

The X-ray absorption spectra of the transition metal L_3,2-edges of the dichalcogenides FeS₂ (pyrite), FeSe₂, FeTe₂, CoS₂, CoSe₂, CoTe₂, NiS₂, NiSe₂ and NiTe₂ have been studied. These spectra are compared with theoretical simulations make using the atomic multiplet approach assuming on octahedral site for the metal atoms. The Ni spectra can be closely simulated using this approach but the Fe and Co spectra show less good fits to the experimental data. The reasons for this and its implications for the use of L-edge spectra in the study of such minerals are discussed in terms of the covalency and the ground states of the compounds; these results indicate the mixing of the ground states increases in the order Ni<Co<Fe.     

Thermomechanical properties of deep argillaceous formations

A difficulty in the interpretation of mechanical and thermo-mechanical tests on specimens drawn from large argillaceous formations is the strong inhomogeneity of void ratio, clay minerals and carbonates content.

In this paper a relationship is developed to link strength and the maximum preconsolidation stress to the initial void ratio and carbonate content. Compressibility is also correlated to carbonates.

Thermal strains in drained and undrained conditions for a Spanish, a belgian and an Italian natural clay are compared. In the elastic state strains are comparable, while in the plastic range thermal strains are highest for the Belgian clay, lower in the Spanish cemented clay and lowest in the Italian clay, very stiff and cemented.     

Solar wind induction in Mercury: Constraints on the formation of a magnetosphere

The discovery of Mercury's magnetosphere by Mariner 10 was surprising since the conventional view of regenerative planetary dynamos had been that the spin requirement would likely have been in excess of the observed spin rate of Mercury. Also internal fluid motions were not expected to be sufficiently large. This paper explores the alternative model of the formation of Mercury's magnetosphere via electromagnetic induction forced by the solar wind. It is shown, however, that the constraints are so severe as to limit severely the applicability of such a model. Although induction is easily observed on the Moon, the modification of the magnetic boundary condition associated with a plasma magnetosphere on Mercury rules out its formation via induction except for interplanetary driving fields which are decreasing in amplitude. That model is explored but retains the difficulty that induced magnetospheres tend to be of small radial and temporal extent compared to that inferred by Ness et al. for Mercury.     

Evidence for downflow following a coronal transient?

On 11 September 1973 a peculiar prominence was observed. The prominence displayed strong ( 50km s^–1) systematic motions toward and away from the observer. The unusual spectrographic appearance of the prominence might have been due to downflowing material lifted into the corona during an earlier coronal transient.The National Center for Atmospheric Research is sponsored by the National Science Foundation.Operated by the Association of Universities for Research in Astronomy, Inc., under contract with the National Science Foundation.     

Stochastic evaluation of simple pairing approaches to reconstruct incomplete rainfall time series

Two-station pairing approaches are routinely used to infill missing information in incomplete rainfall databases. We evaluated the performance of three simple methodologies to reconstruct incomplete time series in presence of variable nonlinear correlation between data pairs. Nonlinearity stems from the statistics describing the marginal peak-over-threshold (POT) values of rainfall events. A Monte Carlo analysis was developed to quantitatively assess expected errors from the use of chronological pairing (CP) with linear and nonlinear regression and frequency pairing (FP). CP is based on a priori selection of regression functions, while FP is based on matching the probability of non-exceedance of an event from one time series with the probability of non-exceedance of a similar event from another time series. We adopted a generalized Pareto (GP) model to describe POT events, and a t-copula algorithm to generate reference nonlinearly correlated pairs of random temporal distributions distributed according with the GP model. The results suggest that the optimal methodology strongly depends on GP statistics. In general, CP seems to provide the lowest errors when GP statistics were similar and correlation became linear; we found that a power-2 function performs well for the selected statistics when the number of missing points is limited. FP outperforms the other methods when POT statistics are different and variables are markedly nonlinearly correlated. Ensemble-based results seem to be supported by the analysis of observed precipitation at two real-world gauge stations.     

Observations of vertical ground accelerations exceeding gravity during the 1997 Umbria-Marche (central Italy) earthquakes

We found extensive evidence that the vertical ground accelerations produced during the largest shock (M = 6.0) of the 1997 Umbria-Marche earthquake sequence exceeded 1g in two areas close to the heavily-damaged villages of Annifo and Colle Croce. This evidence comes from the striking observation of thousands of freshly fractured and broken rocks and stones in these areas. Some of the broken stones lie isolated on soft detritic soil while others had been previously piled up, probably a long time agoto clear the fields for farming. The freshness of the cuts and fractures and the consistency of the observations for thousands of rocks and stones in these areas indicate that these rocks were thrown upwards during the earthquake, with breakage occurring at the time of impact. Ground motion calculations consistent with the static deformation inferred from GPS and interferometry data, show that the broken stones and rocks are found in the zone where the strongest shaking took place during the earthquake and that most of the shaking there was vertical.     

Ambient Changes in Tracer Concentrations from a Multilevel Monitoring System in Basalt

Starting in 2008, a 4‐year tracer study was conducted to evaluate ambient changes in groundwater concentrations of a 1,3,6‐naphthalene trisulfonate tracer that was added to drill water. Samples were collected under open borehole conditions and after installing a multilevel groundwater monitoring system completed with 11 discrete monitoring zones within dense and fractured basalt and sediment layers in the eastern Snake River aquifer. The study was done in cooperation with the U.S. Department of Energy to test whether ambient fracture flow conditions were sufficient to remove the effects of injected drill water prior to sample collection. Results from thief samples indicated that the tracer was present in minor concentrations 28 days after coring, but was not present 6 months after coring or 7 days after reaming the borehole. Results from sampling the multilevel monitoring system indicated that small concentrations of the tracer remained in 5 of 10 zones during some period after installation. All concentrations were several orders of magnitude lower than the initial concentrations in the drill water. The ports that had remnant concentrations of the tracer were either located near sediment layers or were located in dense basalt, which suggests limited groundwater flow near these ports. The ports completed in well‐fractured and vesicular basalt had no detectable concentrations.     

Petrogenesis of voluminous mid-Tertiary ignimbrites of the Sierra Madre Occidental,Chihuahua, Mexico

The mid-Tertiary ignimbrites of the Sierra Madre Occidental of western Mexico constitute the largest continuous rhyolitic province in the world. The rhyolites appear to represent part of a continental magmatic arc that was emplaced when an eastward-dipping subduction zone was located beneath western Mexico.In the Batopilas region of the northern Sierra Madre Occidental the mid-Tertiary Upper Volcanic sequence is composed predominantly of rhyolitic ignimbrites, but volumetrically minor lava flows as mafic as basaltic andesite are also present. The basaltic andesite to rhyolite series is calc-alkalic and contains 1% K₂O at 60% SiO₂. Trace element abundances of a typical ignimbrite with 73% SiO₂ are Sr 225 ppm, Rb 130 ppm, Y 32 ppm, Th 12 ppm, Zr 200 ppm, and Nb 15 ppm. The entire series plots as coherent and continuous trends on variation diagrams involving major and trace elements, and the trends are distinct from those of geographicallyassociated rocks of other suites. We interpret these and other geochemical variations to indicate that the rocks are comagmatic. Mineral chemistry, Sr isotopic data, and REE modelling support this interpretation.Least squares calculations show that the major element variations are consistent with formation of the basaltic andesite to rhyolite series by crystal fractionation of observed phenocryst phases in approximate modal proportions. In addition, calculations modelling the behavior of Sr with the incompatible trace element Th favor a fractional crystallization origin over a crustal anatexis origin for the rock series. The fractionating minerals included plagioclase (> 50%), and lesser amounts of Fe-Ti oxides, pyroxenes, and/or hornblende. The voluminous ignimbrites represent no more than 20% of the original mass of a mantle-derived mafic parental magma.