Post‐stack stratigraphic inversion workflow applied to carbon dioxide storage: application to the saline aquifer of Sleipner field‡

In the Norwegian North Sea, the Sleipner field produces gas with a high CO₂ content. For environmental reasons, since 1996, more than 11 Mt of this carbon dioxide (CO₂) have been injected in the Utsira Sand saline aquifer located above the hydrocarbon reservoir. A series of seven 3D seismic surveys were recorded to monitor the CO₂ plume evolution. With this case study, time‐lapse seismics have been shown to be successful in mapping the spread of CO₂ over the past decade and to ensure the integrity of the overburden. Stratigraphic inversion of seismic data is currently used in the petroleum industry for quantitative reservoir characterization and enhanced oil recovery. Now it may also be used to evaluate the expansion of a CO₂ plume in an underground reservoir. The aim of this study is to estimate the P‐wave impedances via a Bayesian model‐based stratigraphic inversion. We have focused our study on the 1994 vintage before CO₂ injection and the 2006 vintage carried out after a CO₂ injection of 8.4 Mt. In spite of some difficulties due to the lack of time‐lapse well log data on the interest area, the full application of our inversion workflow allowed us to obtain, for the first time to our knowledge, 3D impedance cubes including the Utsira Sand. These results can be used to better characterize the spreading of CO₂ in a reservoir. With the post‐stack inversion workflow applied to CO₂ storage, we point out the importance of the a priori model and the issue to obtain coherent results between sequential inversions of different seismic vintages. The stacking velocity workflow that yields the migration model and the a priori model, specific to each vintage, can induce a slight inconsistency in the results.     

Strain-induced magma degassing: insights from simple-shear experiments on bubble bearing melts

Experiments have been performed to determine the effect of deformation on degassing of bubble-bearing melts. Cylindrical specimens of phonolitic composition, initial water content of 1.5 wt.% and 2 vol.% bubbles, have been deformed in simple-shear (torsional configuration) in an internally heated Paterson-type pressure vessel at temperatures of 798–848 K, 100–180 MPa confining pressure and different final strains. Micro-structural analyses of the samples before and after deformation have been performed in two and three dimensions using optical microscopy, a nanotomography machine and synchrotron tomography. The water content of the glasses before and after deformation has been measured using Fourier Transform Infrared Spectroscopy (FTIR). In samples strained up to a total of γ ∼ 2 the bubbles record accurately the total strain, whereas at higher strains (γ ∼ 10) the bubbles become very flattened and elongate in the direction of shear. The residual water content of the glasses remains constant up to a strain of γ ∼ 2 and then decreases to about 0.2 wt.% at γ ∼ 10. Results show that strain enhances bubble coalescence and degassing even at low bubble volume-fractions. Noticeably, deformation produced a strongly water under-saturated melt. This suggests that degassing may occur at great depths in the volcanic conduit and may force the magma to become super-cooled early during ascent to the Earth’s surface potentially contributing to the genesis of obsidian.     

Source location of long period seismicity at Volcàn de Colima,México

This paper presents an analysis of seismicity associated with the volcanic activity of Volcàn de Colima (México) and recorded in the period November 2005–April 2006 during a field survey by the Istituto Nazionale di Geofisica e Vulcanologia (INGV)–Osservatorio Vesuviano, the Observatorio Vulcanologico de Colima of Colima University and the Instituto Andaluz de Geofisica, University of Granada. Three different types of volcanic earthquakes have been identified on the basis of their spectral properties: Type A (0.3–1 Hz), Type B (1–5 Hz) and Type C (3–4 Hz). Results of polarization analysis applied to Type A events show a predominance of radial motion, indicating that the wavefield comprises compressional waves (P) and shear waves polarized in the vertical plane (SV), while the signal always begins with a negative polarity. Type A, B and C earthquakes have been located using both a flat layered model and a 3D model including topography. Hypocentre distributions indicate that the source of Type A signals is very shallow and confined to a small volume lying about 1 km below the crater. In contrast, the source of Type B and C events is significantly deeper, with most hypocentres located in a volume of about 1 km³ centred at 2.5–3 km depth. A cluster analysis based on the cross-correlation among the waveforms of different events recorded at the same station was applied to Type A earthquakes. Only two clusters, which include only a small percentage of events were found, indicating that earthquake families were uncommon during the period of our survey.     

Decadal variability and scales of the sea surface structure in the northern Ionian

The variability and scales of the sea surface structure of the northern Ionian Sea from January 1993 to December 2007 were studied by means of altimeter remotely-sensed weekly Sea Level Anomaly (SLA) objective maps. Variability in the sea surface structure was addressed by means of empirical orthogonal function (EOF) analysis and, assuming an exponential correlation model, scales of the SLA field were quantified as e-folding distances of the SLA autocorrelation function. The variability in the sea surface structure, described by the first three EOFs, which cumulatively explain 60.3% of the data set variance, is characterized by a large-scale structure with variability on a time scale of ∼10-13 years and, on shorter scales, an eddy system with variability on an annual scale. The variability in the large-scale structure describes an overturning of the SLA field, which took place in 1997, and determines a reversal of the geostrophic upper-layer circulation. As the large-scale circulation transition takes place, time-dependent spectral analysis of EOF coefficients shows a redistribution of the spectral energy from inter-annual to semi-annual and monthly components. Spatial scales display variability on an annual and inter-annual time scale. On the annual time scale, variability in spatial scales is characterized by longer values in summer-fall and shorter in winter-spring. Inter-annual variability in spatial scales is demonstrated by a remarkable drop in the values during fall in the period 1998-2000. We propose an explanation of the variability in horizontal scales in terms of the redistribution of water masses and related modifications of the vertical structure of the water column associated with different regimes of the basin-scale circulation.     

Field observations from the Aquila, Italy earthquake of April 6, 2009

On April 6, 2009 an earthquake of magnitude 6.2 (M_w) struck the Abbruzzo region of Italy causing widespread damage to buildings in the city of L??Aquila and surrounding areas. This paper summarizes field observations made by the Earthquake Engineering Field Investigation Team (EEFIT) after the event. The paper presents an overview of seismological and geotechnical aspects of the earthquake as well as a summary of the observed damage to buildings and infrastructure. A brief overview of the earthquake casualties is also reported.     

High-temperature vesuvianite: crystal chemistry and surface considerations

A multi-methodical approach has been applied for characterizing the bulk and surface crystal chemical features of a high-temperature vesuvianite crystal from skarns of Mount Somma-Vesuvius Volcano (Naples, Italy). Vesuvianite belongs to the space group P4/nnc with unit cell parameters a = 15.633(1) Å, c = 11.834(1) Å and chemical formula (Ca_18.858 Na_0.028 Ba_0.004 K_0.006 Sr_0.005 □_0.098)_19.000 (Al_8.813 Ti_0.037 Mg_2.954 Mn_0.008 Fe_0.114²⁺ Fe_1.375³⁺ Cr_0.008 B_0.202)_13.511 Si_18.000(O_0.261 F_0.940 OH_7.799)_9.000. Structure refinement, which converges at R = 0.0328, demonstrates a strong positional disorder down the fourfold axes, indicating that the Y1 site is split into two positions (Y1A and Y1B) alternatively occupied. However, because of X4 proximity to Y1B and Y1A, X4 cannot be occupied if Y1B or Y1A are. Overall Y1 occupancy (Y1A + Y1B) reaches approximately 0.5, as common in vesuvianite and occupancy of Y1B site is extremely limited. Moreover, T1 position, limitedly occupied, accommodates the excess of cations generally related to Y position. A small quantity (0.202 apfu) of boron is sited at the T2 site that, like T1, is poorly occupied. The determination of the amount of each element on the (100) vesuvianite surface, obtained through X-ray photoelectron spectroscopy high-resolution spectra in the region of the Si_2p, Al_2p, Mg_1s, and Ca_2p core levels, evidences that a greater amount of aluminum and a smaller amount of calcium characterize the surface with respect to the bulk. Although both of these features require further investigation, we may consider the Al increase can be related to preferential orientation of Al-rich sites on the (100) plane. Furthermore, the surface structure of vesuvianite suggests that Al, Ca, and Mg cations maintain coordination features at the surface similar to the bulk. Silica, however, while presenting fourfold coordination, shows also a [1]-fold small coordinated component at binding energy 99.85 eV, due to broken Si–O bonds at the surface. The presence of eight- and nine-coordinated Ca cations is suggested by a large XPS feature resulting from the merging of Ca_2p3/2 and Ca_2p1/2 peaks at 348.45 and 352.05 eV, respectively.     

Identifying sources of chlorinated aliphatic hydrocarbons in a residential area in Italy using the integral pumping test method

The results of integral pumping tests (IPTs) performed in the city of Fabriano, Italy, are presented. The IPT methodology was developed by the European Union project INCORE, as a tool for groundwater investigation and source localization in contaminated areas. This methodology consists of a multiple-well pumping test in which the wells are positioned along a control plane downstream of suspected contaminant source zones and perpendicular to the mean groundwater flow direction. During the pumping, concentration time series of target contaminants are measured. In Fabriano, two control planes were realized to identify a chlorinated aliphatic hydrocarbon plume, to estimate the mass fluxes and draw up a ranked list of the main contamination sources. A numerical flow model was implemented to support the IPT design and to interpret the results. This study revealed low-level trichloroethylene contamination (concentration below 8 μg/l), tetrachloroethylene contamination (mean concentration up to 500 μg/l) and a mass flow rate of about 300 g/day. Through the application of the IPT method, the mean contaminant concentrations, the spatial distribution of concentration values along the control planes, and the total contaminant mass flow rates were evaluated, and the investigation area was reduced for further and deeper investigation activities.     

Effect of tidal fluctuations on transient dispersion of simulated contaminant concentrations in coastal aquifers

Variable-density groundwater models require extensive computational resources, particularly for simulations representing short-term hydrologic variability such as tidal fluctuations. Saltwater-intrusion models usually neglect tidal fluctuations and this may introduce errors in simulated concentrations. The effects of tides on simulated concentrations in a coastal aquifer were assessed. Three analyses are reported: in the first, simulations with and without tides were compared for three different dispersivity values. Tides do not significantly affect the transfer of a hypothetical contaminant into the ocean; however, the concentration difference between tidal and non-tidal simulations could be as much as 15%. In the second analysis, the dispersivity value for the model without tides was increased in a zone near the ocean boundary. By slightly increasing dispersivity in this zone, the maximum concentration difference between the simulations with and without tides was reduced to as low as 7%. In the last analysis, an apparent dispersivity value was calculated for each model cell using the simulated velocity variations from the model with tides. Use of apparent dispersivity values in models with a constant ocean boundary seems to provide a reasonable approach for approximating tidal effects in simulations where explicit representation of tidal fluctuations is not feasible.     

Multi-scale European Soil Information System (MEUSIS): a multi-scale method to derive soil indicators

The Multi-scale Soil Information System (MEUSIS) can be a suitable framework for building a nested system of soil data that could facilitate interoperability through a common coordinate reference system, a unique grid coding database, a set of detailed and standardized metadata, and an open exchangeable format. In the context of INSPIRE Directive, MEUSIS may be implemented as a system facilitating the update of existing soil information and accelerating the harmonization of various soil information systems. In environmental data like the soil one, it is common to generalize accurate data obtained at the field to coarser scales using either the pedotransfer rules or knowledge of experts or even some statistical solutions which combine single values of spatially distributed data. The most common statistical process for generalization is averaging the values within the study area. In this paper, we do not present a simple averaging of numerical values without any further processed information. The upscaling process is accompanied with significant statistical analysis in order to demonstrate the method suitability. The coarser resolution nested grids cells (10??? 10?km) represent broad regions where the calculated soil property (e.g., organic carbon) can be accurately upscaled. Multi-scaled approaches are urgently required to integrate different disciplines (such as Statistics) and provide a meta-model platform to improve current mechanistic modeling frameworks, request new collected data, and identify critical research questions. Past papers have described in detail the upscaling methodology while our present approach is to demonstrate an important application of this methodology accompanied with statistical evidence.     

Realised added value in dynamical downscaling of Australian climate change

Climate Dynamics - Coarse resolution global climate models (GCMs) cannot resolve fine-scale drivers of regional climate, which is the scale where climate adaptation decisions are made. Regional...