A control-volume finite-element method for three-dimensional multiphase basin modeling

In this paper we describe a 3D control-volume finite-element method to solve numerically the coupled partial differential equations (PDEs) governing geological processes involved in the evolution of sedimentary basins. These processes include sediment deposition and deformation, hydrocarbon generation, multiphase fluid flow, and heat transfer in deforming porous media.     

The Holocene tephrostratigraphic record of Lake Shkodra (Albania and Montenegro)

Two cores were recovered in the southeastern part of Lake Shkodra (Montenegro and Albania) and sampled for identification of tephra layers. The first core (SK13, 7.8 m long) was recovered from a water depth of 7 m, while the second core (SK19, 5.8 m long) was recovered close to the present‐day shoreline (water depth of 2 m). Magnetic susceptibility investigations show generally low values with some peaks that in some cases are related to tephra layers. Naked‐eye inspection of the cores allowed the identification of four tephra layers in core SK13 and five tephra layers in core SK19. Major element analyses on glass shards and mineral phases allowed correlation of the tephra layers between the two cores, and their attribution to six different Holocene explosive eruptions of southern Italy volcanoes. Two tephra layers have under‐saturated composition of glass shards (foiditic and phonolitic) and were correlated to the AD 472 and the Avellino (ca. 3.9 cal. ka BP) eruptions of Somma‐Vesuvius. One tephra layer has benmoreitic composition and was correlated to the FL eruption of Mount Etna (ca. 3.4 cal. ka BP). The other three tephra layers have trachytic composition and were correlated to Astroni (ca. 4.2 cal. ka BP), Agnano Monte Spina (ca. 4.5 cal. ka BP) and Agnano Pomici Principali (ca. 12.3 cal. ka BP) eruptions of Campi Flegrei. The ages of tephra layers are in broad agreement with eight ¹⁴C accelerator mass spectrometric measurements carried out on plant remains and charcoal from the lake sediments at different depths along the two cores. The recognition of distal tephra layers from Italian volcanoes allowed the physical link of the Holocene archive of Lake Shkodra to other archives located in the central Mediterranean area and the Balkans (i.e. Lake Ohrid). Five of the recognised tephra layers were recognised for the first time in the Balkans area, and this has relevance for volcanic hazard assessment and for ash dispersal forecasting in case of renewed explosive activity from some of the southern Italy volcanoes. Copyright © 2009 John Wiley & Sons, Ltd.     

An integrated procedure to evaluate hydrological models

The growing concern for health‐related problems deriving from pollutants leaching is driving national and international administrations to support the development of tools for evaluating the effects of alternate management scenarios and identifying vulnerable areas. Cropping systems models are powerful tools for evaluating leachates under different environmental, social, and management conditions. As percolating water is the transport vehicle for pollutants transport in soil, a reliable evaluation of water balance models is a fundamental prerequisite for investigating pesticides and nitrate fate. As specific approaches for the evaluation of multi‐layer evolution of state variables are missing, we propose a fuzzy‐based, integrated indicator (I_SWC: 0, best; 1, worst) for a comprehensive evaluation of soil water content (SWC) simulations. We aggregated error metrics with others quantifying the homogeneity of errors across different soil layers, the capability of models to reproduce complex dynamics function of both time and soil depth, and model complexity. We tested I_SWC on a sample dataset where the models CropSyst and CERES‐Wheat were used to simulate SWC for winter wheat systems. I_SWC revealed that, in the explored conditions, the global assessment of the two models' performances allowed identification of CropSyst as the best (average I_SWC = 0·441, with a value of 0·537 obtained by CERES‐Wheat), although each model prevailed for some of the metrics. CropSyst presented the highest accuracy (average agreement module = 0·400), whereas CERES‐Wheat's accuracy was slightly worse, although achieved with a simplified modelling approach (average Akaike Information Criterion = − 230·44), thereby favouring large‐area applicability. The non‐univocal scores achieved by the models for the different metrics support the use of multi‐metric evaluation approaches for quantifying the different aspects of water balance model performances. Copyright © 2010 John Wiley & Sons, Ltd.     

Explanation for the numerous upper floor collapses during the 1985 Mexico city earthquake

Presented herein are the results of a study carried out to investigate the prevalent factors that contributed to the collapse of the upper floors of a large number of buildings during the earthquake of September 19, 1985, in Mexico City. The study involves the review of the available data from such collapses, the identification of common characteristics among the buildings that suffered such type of collapse, the design of a typical 10-storey reinforced concrete frame structure according to the 1976 Mexico City building code, and the non-linear analysis of one of its frames under one of the accelerograms from the earthquake. It is found that in many cases the upper floor collapses might have occurred because the large accelerations and long duration of the earthquake induced in the failed buildings the formation of plastic hinges at their upper columns, and because these plastic hinges in turn induced the lateral instability of the affected storeys and the ones above them. It is also found that in some structures this failure mechanism can occur after plastic hinges are formed in most of their beams, and that it may result from the consequent elongation of their natural periods, their response in a second or third mode and in resonance with the ground motion that excites them, and their vibration in a mode not considered to be the predominant one in their design.     

Reduction seismic response with heavily-damped vibration absorbers

It is shown that two of the damping ratios of certain systems composed of a building and a small attachment in resonance are given by the average of the damping ratios of the two independent components. Based on this fact and the fact that the seismic response of a building can always be reduced by increasing its damping, it is demonstrated that the attachment of a small heavily-damped system in resonance can increase the damping of a building and reduce thus its response to earthquake excitations. Numerical solutions are presented to confirm the demonstration, and recommendations are given to calculate the parameters of such systems.     

Passive seismic control of cable-stayed bridges with damped resonant appendages

A study is conducted to investigate the effectiveness of attaching to cable-stayed bridges resonant appendages with a relatively small mass and a high damping ratio as a means to reduce their response to earthquake excitations. The study is based on a previously developed formulation that shows that the use of these appendages may increase the inherent damping of building structures and, as a result, may reduce their response to seismic disturbances. It includes numerical and experimental tests that are conducted to assess the validity of such a formulation for the case of cable-stayed bridges and the extent to which such appendages can reduce their seismic response. In the numerical study, an actual cable-stayed bridge is modelled with finite elements and analysed with and without the proposed appendages under different earthquake ground motions. Appendages with damping ratios of 10, 15, 20 and 30 per cent and masses that, respectively, represent 0.67, 1.5, 2.7 and 6.0 per cent of the total mass of the bridge are considered. In the experimental test, a 3.7 m long cable-stayed bridge and an appendage consisting of a small mass, a small spring and a small viscous damper are built and the bridge tested, without and with the appendage, on a pair of shaking tables which are set to reproduce ground acceleration records from past earthquakes. The damping ratio of the appendage in this test is 32 per cent and its mass represents about 8 per cent of the total mass of the bridge model. In the numerical test, it is found that the appendages reduce the longitudinal response of the bridge deck up to 88 per cent. Similarly, in the experimental test it is found that the appendage reduces the longitudinal bridge deck response by about 41 per cent. From these studies, it is concluded that the suggested appendages may indeed be effective in reducing the seismic response of cable-stayed bridges as they are for building structures.     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

Ice core evidence for secular variability and 200-year dipolar oscillations in atmospheric circulation over East Antarctica during the Holocene

Two Holocene ice core records from East Antarctica (Vostok and EPICA-Dome C) were analysed for dust concentration and size distribution at a temporal resolution of 1 sample per ~50 years. A series of volcanic markers randomly distributed over the common part of the ice cores (from 9.8 to 3.5 kyear BP) ensures accurate relative dating (±33 years). Dust-size records from the two sites display oscillations structured in cycles with sub-millennial and secular scale frequencies that are apparently asynchronous. The power spectra of the composite sum (Σ) of the two dust-size records display spectral energy mostly for 150- to 500-year periodicities. On the other hand, the 200-year band is common to both records and the 200 year components of the two sites are out-of-phase (100-year lead or lag) over ~5.5 kyear, a phenomenon also reflected by a significant (>99% conf. lev.) band in the power spectra of the composite difference (Δ) of the two size records. During long-range transport, mineral dust originating from the Southern Hemisphere continents is graded to a variable extent depending on the altitude and duration of atmospheric transport. Relatively coarse dust is associated with air mass penetration from the middle–lower troposphere and conversely relatively fine dust with upper troposphere air masses or the influence of subsidence over the Antarctic plateau, a hypothesis already proposed for the changes that occurred during the Last Glacial Maximum to Holocene transition (Delmonte et al. 2004b). Moreover, we assume that the overall fluctuation of air mass advection over Antarctica depends on the meridional pressure gradient with respect to low latitudes, i.e. the Antarctic Oscillation (AAO). We therefore suggest a regional variability in atmospheric circulation over East Antarctica. The 150–500 year power spectrum of the composite (Σ) parameter represents the long term variability of the AAO, imprinted by secular internal oscillations probably related to the southern ocean-climatic system. On the other hand, the Δ dust composite parameter suggests a persistent atmospheric dipole over East Antarctica delivering coarser (finer) dust particles alternatively to Vostok and Dome C regions with a bi-centennial periodicity. Indeed, a seesaw phenomenon in dust size distribution was already observed at three East Antarctic sites during the last deglaciation (Delmonte et al. 2004b) and was interpreted as a progressive reduction of the eccentricity of the polar vortex with respect to the geographic south pole. Interestingly, the Δ parameter shows a pronounced 200-year oscillation mode, throwing new light on the unresolved question of a possible relationship between climate and solar activity.     

Spatial assessment of groundwater use potential for irrigation in Teesta Barrage Project in Bangladesh

In Bangladesh, development of the groundwater resource for irrigation is a vital component of the government’s agricultural strategy to attain food self-sufficiency. Amidst reports of falling groundwater levels in many parts of Bangladesh, the potential of groundwater use has been investigated in the Teesta Barrage Project (TBP) in which large-scale groundwater development for dry-season irrigation has taken place in the recent past. Several techniques and tools have been applied such as the combination of analysis of groundwater hydrographs and mathematical modelling to derive key hydrogeological variables, calculation of net irrigation requirement and the use of geographical information systems. The results show that the economically attractive high-yielding variety (HYV) Boro (dry season) rice cultivation during the groundwater irrigation season may not be sustained in large parts of the project area if the current trends in abstraction are continued. However, due to spatial variation in abstraction, nine thanas (sub-districts)—out of a total of 21 in the project area—may still be able to expand groundwater-irrigated cropland and a groundwater-use potential of 40 mm/year may be created if deep-set shallow tubewells are used by the farmers to abstract groundwater. A structured approach, based on zoning of potential areas, is recommended for groundwater development and use.

---

Resumen En Bangladesh, el desarrollo de las aguas subterráneas para riego es un componente vital del la estrategia agrícola del Gobierno para obtener la auto-suficiencia alimentaria. Un informe intermedio evidencia un descenso de los niveles de aguas subterráneas en muchas partes de Bangladesh, el uso potencial de aguas subterránea ha sido investigado en el Proyecto Tessta Barrage (TBP) en el cual ha tenido lugar el desarrollo de las aguas subterráneas a gran escala para el riego en la estación seca en un pasado reciente. Se ha aplicado algunas técnicas y herramientas, como la combinación del análisis de datos gráficos de aguas subterráneas y modelización matemática para obtener variables hidrogeológicas clave, el cálculo de la necesidad neta de riego y el uso de Sistemas de Información Geográfica. El resultado muestra que el cultivo, económicamente atractivo, de arroz HYV Boro (estación seca) durante la estación de riego con aguas subterráneas puede no ser sostenible en gran parte del área del proyecto si continúan las tendencias actuales en las extracciones. Sin embargo, debido a la variación espacial en la extracción, nueve thanas (subdistritos)—de un total de 21 en el área del Proyecto—pueden todavía sufrir un aumento de la tierra cultivable regada con aguas subterráneas y que se puede asumir un uso potencial de agua subterránea de 40 mm/estación si los granjeros utilizan sondeos someros instalados en fosas para la extracción. Se recomienda una aproximación estructurada para el desarrollo y uso de las aguas subterráneas, basada en la zonificación de áreas potenciales.

---

Résumé Au Bangladesh, l’exploitation des ressources en eau souterraines pour l’irrigation est une composante vitale de la stratégie agricole gouvernementale pour atteindre l’auto-suffisance alimentaire. Parmi les baisses de niveaux piézométriques rapportées sur de nombreuses régions du Bangladesh, le potentiel d’utilisation des eaux souterraines a été étudié sur le Projet de Barrage de Teesta (TBP), où, dans un passé récent, les eaux souterraines ont été massivement exploitées pour l’irrigation en saison sèche. De nombreux outils et techniques ont été utilisés, comme l’analyse combinée des chroniques piézométriques et des modèles mathématiques pour dégager les variables hydrogéologiques clés, le calcul des besoins nets pour l’irrigation, et l’utilisation des systèmes d’information géographique. Les résultats démontrent que la culture en saison d’irrigation de la variété de riz HYV Boro (saison sèche), économiquement attrayante, n’est potentiellement pas viable sur de nombreux secteurs de la zone étudiée si les prélèvements perdurent selon la tendance actuelle. Cependant, du fait de la variation spatiale des prélèvements, neuf thanas (sous-districts)—sur un total de 21 sur le secteur d’étude—pourraient encore augmenter les surfaces irriguées, et une disponibilité en eau de 40 mm/saison pourrait être générée si les agriculteurs utilisaient des puits peu profonds pour exploiter les ressources souterraines. Une approche structurée, basée sur le zonage des secteurs potentiels, est recommandée pour l’exploitation et l’utilisation des eaux souterraines.