Recharge estimation for transient ground water modeling

Reliable ground water models require both an accurate physical representation of the system and appropriate boundary conditions. While physical attributes are generally considered static, boundary conditions, such as ground water recharge rates, can be highly variable in both space and time. A practical methodology incorporating the hydrologic model HELP3 in conjunction with a geographic information system was developed to generate a physically based and highly detailed recharge boundary condition for ground water modeling. The approach uses daily precipitation and temperature records in addition to land use/land cover and soils data. The importance of the method in transient ground water modeling is demonstrated by applying it to a MODFLOW modeling study in New Jersey. In addition to improved model calibration, the results from the study clearly indicate the importance of using a physically based and highly detailed recharge boundary condition in ground water quality modeling, where the detailed knowledge of the evolution of the ground water flowpaths is imperative. The simulated water table is within 0.5 m of the observed values using the method, while the water levels can differ by as much as 2 m using uniform recharge conditions. The results also show that the combination of temperature and precipitation plays an important role in the amount and timing of recharge in cooler climates. A sensitivity analysis further reveals that increasing the leaf area index, the evaporative zone depth, or the curve number in the model will result in decreased recharge rates over time, with the curve number having the greatest impact.     

A Lost Realm in the Internal Domains of the Betic-Rif Orogen (Spain and Morocco): Evidence from Conglomerates and Consequences for Alpine Geodynamic Evolution

The Malaguide-Ghomaride Complex is capped by Upper Oligocene-Aquitanian clastic deposits postdating early Alpine orogenesis but predating the main tectonic-metamorphic evolution, end of nappe emplacement, unroofing, and exhumation of the metamorphic units of the Betic-Rif Orogen. Two conglomerate intervals within these deposits are characterized by clasts of sedimentary, epimetamorphic, and mafic volcanic rocks derived from Malaguide-Ghomaride units and by clasts of acidic magmatic and orthogneissic rocks of unknown provenance, here studied. Magmatic rocks originated from late-Variscan two-mica cordierite-bearing granitoids and, subordinately, from aplitic dikes. Orthogneisses derive from similar plutonic rocks but are affected by an Alpine metamorphic overprint evolving from greenschist (T=510&j0;-530 degrees C and P=5-6 kbar) to low-temperature amphibolite facies (T>550&j0;C and P<3 kbar). Such a plutonic rock suite is unknown in any Betic-Rif unit or in the basement of the Alboran Sea, and the metamorphic evolution in the orthogneisses is different from (and older than) that of Alpujarride-Sebtide rocks to which they were formerly ascribed. Magmatic and metamorphic rocks very similar to those studied characterize the basements of some Kabylia and Calabria-Peloritani units. Therefore, the source area is a currently lost continental-crust realm of Calabria-Peloritani-Kabylia type, located to the ESE of the Malaguide-Ghomaride Domain and affected by a pre-latest Oligocene Alpine metamorphism. Increasingly active tectonics transformed this realm into rising areas from which erosion fed small subsiding synorogenic basins formed on the Malaguide-Ghomaride Complex. This provenance analysis demonstrates that all these domains constituted a single continental-crust block until Aquitanian-Burdigalian times, before its dispersal around nascent western Mediterranean basins.     

Development and validation of a metallic haunch seismic retrofit solution for existing under‐designed RC frame buildings

The feasibility and efficiency of a seismic retrofit solution for existing reinforced concrete frame systems, designed before the introduction of modern seismic‐oriented design codes in the mid 1970s, is conceptually presented and experimentally investigated. A diagonal metallic haunch system is introduced at the beam–column connections to protect the joint panel zone from extensive damage and brittle shear mechanisms, while inverting the hierarchy of strength within the beam–column subassemblies and forming a plastic hinge in the beam. A complete step‐by‐step design procedure is suggested for the proposed retrofit strategy to achieve the desired reversal of strength hierarchy. Analytical formulations of the internal force flow at the beam–column‐joint level are derived for the retrofitted joints. The study is particularly focused on exterior beam–column joints, since it is recognized that they are the most vulnerable, due to their lack of a reliable joint shear transfer mechanism. Results from an experimental program carried out to validate the concept and the design procedure are also presented. The program consisted of quasi‐static cyclic tests on four exterior, ? scaled, beam–column joint subassemblies, typical of pre‐1970 construction practice using plain round bars with end‐hooks, with limited joint transverse reinforcement and detailed without capacity design considerations. The first (control specimen) emulated the as‐built connection while the three others incorporated the proposed retrofitted configurations. The experimental results demonstrated the effectiveness of the proposed solution for upgrading non‐seismically designed RC frames and also confirmed the applicability of the proposed design procedure and of the analytical derivations. Copyright © 2006 John Wiley & Sons, Ltd.     

Decoupling of the coupled poroelastic equations for quasistatic flow in deformable porous media containing two immiscible fluids

The study of the poroelastic behavior of sedimentary materials containing two immiscible fluids in response to either applied stress or pore pressure change in a quasistatic limit, i.e., negligible second time-derivatives, is of great importance to many hydrogelogical problems, e.g., land subsidence caused by withdrawal of subsurface fluids. The poroelasticity models developed for analyzing these problems feature partial differential equations that are coupled in the terms describing viscous damping and strain field. To determine closed-form analytical solutions for induced volumetric strain (dilatation) of the solid framework and its interaction with fluid flows, the choice of normal coordinates whose transformation can be performed to decouple these poroelastic equations is highly desirable. In this paper, we show that normal coordinates for decoupling these equations are real-valued and equal to three different linear combinations of the dilatations of the solid and the fluids (or equivalently, three different linear combinations of two individual fluid pressures and solid dilatation). In contrast to fully saturated porous media, it is found that the viscous damping effect must be represented in normal coordinates in the presence of the second fluid. The resulting decoupled equations representing independent motional modes are a Laplace equation and two diffusion equations, which can be solved analytically under a variety of initial and boundary conditions. Thus, after inverse transformation of normal coordinates is performed, the closed-form analytical solutions for induced solid volumetric strain and excess pore fluid pressures can be obtained simultaneously from our decoupled partial differential equations.     

The origin of the fumaroles of La Solfatara (Campi Flegrei, South Italy)

The analysis of gaseous compositions from Solfatara (Campi Flegrei, South Italy) fumaroles since the early 1980s, clearly reveals a double thermobarometric signature. A first signature at temperatures of about 360 °C was inferred by methane-based chemical-isotopic geoindicators and by the H₂/Ar geothermometer. These high temperatures, close to the critical point of water, are representative of a deep zone where magmatic gases flash the hydrothermal liquid, forming a gas plume. A second signature was found to be at around 200-240 °C. At these temperatures, the kinetically fast reactive species (H₂ and CO) re-equilibrate in a pure vapor phase during the rise of the plume. A combination of these observations with an original interpretation of the oxygen isotopic composition of the two dominant species, i.e. H₂O and CO₂, shed light on the origin of fumarolic fluids by showing that effluents are mixture between fluids degassed from a magma body and the vapor generated at about 360 °C by the vaporization of hydrothermal liquids. A typical ‘andesitic’ water type (δD ∼ −20‰, δ¹⁸O ∼10‰) and a CO₂-rich composition (X_CO2∼0.4) has been inferred for the magmatic fluids, while for the hydrothermal component a meteoric origin and a CO₂ fugacity fixed by fluid-rock reaction at high temperatures have been estimated. In the time the fraction of magmatic fluids in the fumaroles increased (up to ∼0.5) at each seismic and ground uplift crisis (bradyseism) which occurred at Campi Flegrei, suggesting that bradyseismic crises are triggered by periodic injections of CO₂-rich magmatic fluids at the bottom of the hydrothermal system.     

Geophysical and geodetic observations in the messina straits

The present note is intended as a contribution to the clarification of the collision process, generally assumed to be active between the African and the Eurasian plates in the area of the Messina Straits.For this purpose detailed analyses are made of the surveys carried out in the Straits of Messina in order to investigate the 1908 earthquake. The magnitude and intensity are established; limits are set for the linear dimension of the fault and of the dislocation that gave rise to the earthquake and an estimate is then given of these parameters. The results are then set out of the investigation into the focal mechanism of the earthquake of 1908. The results are also studied of geodetic surveys carried out at the time for the purpose of measuring the variations in height that occurred during the earthquake.A study is also made of the results of geodetic surveys (some of which are still under way) designed to measure the horizontal movements of Sicily with respect to Calabria in order to study locally the continental collision process in the area of the Messina Straits.It is found that the old and recent geodetic measurements as well as the fault-plane solution of the old and recent earthquakes of that area and the slips associated with those earthquakes indicate a motion towards the north of Sicily with respect to Calabria.     

Compositions of Formosan basalts and aspects of their petrogeneses

Twenty-three basalts from northwestern Formosa were analyzed for Si, Ti, Al, Fe, Mg, Ca, Na, K, Rb, Sr, Hf, Co, Sc, Cr, Th, and 7 REE. Formosan alkalic basalts are generally similar to most alkalic basalts, and Formosan tholeiites are somewhat similar to island arc and continental tholeiites in terms of trends in K/Rb, Rb/Sr, K/Sr, and Ca/Sr ratios. Compared to most submarine tholeiites, Formosan tholeiites are enriched in Rb relative to K and Sr, and enriched in Sr relative to K and Ca. Both types of Formosan basalts show, relative to chondritic proportions, enrichments in light REE which smoothly decrease through the heavy REE. Their REE patterns are similar to those observed in Dogo-Oki alkalic basalts, Hawaiian basalts, and Deccan Plateau basalt, but differ strikingly from those of Japanese tholeiites and oceanic ridge basalts. These relationships and contrasts suggest that Formosan basalts originated in a plagioclase-free portion of the upper mantle and at greater depths that did Japanese tholeiites and oceanic ridge basalts.