S-wave splitting intensity analysis and inversion

Analysing S-wave splitting has become a routine step in processing multicomponent data. Typically, this analysis leads to determining the principal directions of a transversely isotropic medium with a horizontal symmetry axis, which is assumed to be responsible for azimuthal anisotropy, and to the time delays between the fast and slow S-waves. These parameters are commonly estimated layer-by-layer from the top. Errors in layer stripping occurring in shallow layers might propagate to deeper layers. We propose a method for S-wave splitting analysis and compensation that consists of inverting interval values of splitting intensity to obtain a model of anisotropic parameters that vary with time and/or depth. Splitting intensity is a robust attribute with respect to structural variations and is commutative, which means that it can be summed along a ray (or throughout a sensitivity kernel volume) and can be linearly related to anisotropic perturbations at depth. Therefore, it is possible to estimate anisotropic properties within a geological formation (e.g. the reservoir) by analysing the differences of splitting intensity measured at the top and at the bottom of the layer. This allows us to avoid layer stripping, in particular, for shallow layers where anisotropic parameters are difficult to estimate due to poor coverage, and it makes S-wave splitting analysis simpler to apply. We demonstrate this method on synthetic and real data. Because the splitting intensity attribute shows usefulness in S-wave splitting analysis in transversely isotropic media, we extend the splitting intensity theory to lower symmetry classes. It enables the characterization of tilted transversely isotropic and tilted orthorhombic media, opening new opportunities for anisotropic model building.     

ROckfall risk MAnagement assessment: the RO.MA. approach

The analysis of risk for vehicles and drivers as a result of rockfall on a road is relevant to design management in geotechnical engineering. This process is very complex due to the large number of parameters involved. In this paper, we discuss risk analysis and management procedures for roads subject to rockfall phenomena. To this aims, we are proposing a quantitative method (the RO.MA. approach). We developed an abacus to define the threshold values of acceptable rockfall risk for a given road. Rockfall risk is calculated using an Event Tree approach and compared with the abacus thresholds to evaluate road safety and the need for additional protective measures to reduce the risk to an acceptable level. The approach was successfully applied at a test site located in Bard, Aosta Valley, north-western Italy.     

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.     

Water table variations in Atacama Desert alluvial fans: discussion of “Evidence of short-term groundwater recharge signal propagation from the Andes to the central Atacama Desert: a singular spectrum analysis approach”*

ABSTRACT

Using a statistical approach, Scheihing attempts to demonstrate the direct influence of recharge events in the Precordillera and the Andean Piedmont on water table variations in downstream alluvial fans. The author “unquestionably” discards the existence of ephemeral recharge processes in alluvial fans although they are frequently impacted by major floods. However, the author does not consider previous publications about recharge processes in arid alluvial fans, the hydraulic setting of the Pampa del Tamarugal Aquifer (PTA) and anthropogenic pressure changes in the PTA. Because the sustainable management of groundwater resources in drylands depends on a good understanding of hydrogeological processes, we propose a thorough reinterpretation of the short- and long-term PTA water table variations addressed by the author. In this comment, we illustrate the limits of the sole use of a statistical approach for characterizing both recharge processes and factors controlling the water table variations in the Atacama Desert.     

Flexural uplift of a lithospheric slab near the Vema transform (Central Atlantic): Timing and mechanisms

The Vema Transverse Ridge (VTR) is a prominent, long and narrow topographic anomaly that runs for over 300 km along a sea floor spreading flow line south of the Vema transform at 11° N in the Atlantic. It rises abruptly about 140 km from the axis of the Mid-Atlantic Ridge (MAR) in 10 Myr old crust and runs continuously up to 25 Myr old crust. It reaches over 3 km above the predicted lithospheric thermal contraction level. It is absent in crust younger than 10 Myr; thus, the uplift of the VTR must have ended roughly 10 Ma. The VTR is interpreted as the exposed edge of a flexured and uplifted slab of oceanic lithosphere that was generated at an 80 km long MAR segment. Based on satellite gravimetry imagery this MAR segment was born roughly 50 Ma and increased its length at an average rate of 1.6 mm/yr. Multibeam data show that the MAR-parallel sea floor fabric south of the VTR shifts its orientation by 5° to 10° clockwise in 11–12 Myr old crust, indicating a change at that time of the orientation of the MAR axis and of the position of the Euler rotation pole. This change caused extension normal to the transform, followed between 12 and 10 Ma by flexure of the edge of the lithospheric slab, uplift of the VTR at a rate of 2 to 4 mm/yr, and exposure of a lithospheric section (Vema Lithospheric Section or VLS) at the northern edge of the slab, parallel to the Vema transform. Ages of pelagic carbonates encrusting ultramafic rocks sampled at the base of the VLS at different distances from the MAR axis suggest that the entire VTR rose vertically as a single block within the active transform offset. A 50 km long portion of the crest of the VTR rose above sea level, subsided, was truncated at sea level and covered by a carbonate platform. Subaerial and submarine erosion has gradually removed material from the top of the VTR and has modified its slopes. Spreading half rate of the crust south of the transform decreased from 17.2 mm/yr between 26 and 19 Ma to 16.9 mm/yr between 19 and 10 Ma, to 13.6 mm/yr from 10 Ma to present. The slowing down of spreading occurred close in time to the change in ridge/transform geometry, suggesting that the two events are related. A numerical model relates lithospheric flexure to extension normal to the transform, suggesting that the extent of the uplift depends on the thickness of the brittle layer, consistent with the observed greater uplift of the older lithosphere along the VTR.     

Effects of pore-scale dispersion,degree of heterogeneity,sampling size,and source volume on the concentration moments of conservative solutes in heterogeneous formations

Pore-scale dispersion (PSD), aquifer heterogeneity, sampling volume, and source size influence solute concentrations of conservative tracers transported in heterogeneous porous formations. In this work, we developed a new set of analytical solutions for the concentration ensemble mean, variance, and coefficient of variation (CV), which consider the effects of all these factors. We developed these models as generalizations of the first-order solutions in the log-conductivity variance of point concentration proposed by [Fiori A, Dagan G. Concentration fluctuations in aquifer transport: a rigorous first-order solution and applications. J Contam Hydrol 2000;45(1–2):139–163]. Our first-order solutions compare well with numerical simulations for small and moderate formation heterogeneity and from small to large sampling and source volumes. However, their performance deteriorates for highly heterogeneous formations. Successively, we used our models to study the interplay among sampler size, source volume, and PSD. Our analysis shows a complex and important interaction among these factors. Additionally, we show that the relative importance of these factors is also a function of plume age, of aquifer heterogeneity, and of the measurement location with respect to the mean plume center of gravity. We found that the concentration moments are chiefly controlled by the sampling volume with pore-scale dispersion playing a minor role at short times and for small source volumes. However, the effect of the source volume cannot be neglected when it is larger than the sampling volume. A different behavior occurs for long periods, which may be relevant for old contaminations, or for small injection volumes. In these cases, PSD causes a significant dilution, which is reflected in the concentration statistics. Additionally, at the center of the mean plume, where high concentrations are most likely to occur, we found that sampling volume and PSD are attenuating mechanisms for both concentration ensemble mean and coefficient of variation, except at very large source and sampler sizes, where the coefficient of variation increases with sampler size and PSD. Formation heterogeneity causes a faster reduction of the ensemble mean concentrations and a larger uncertainty at the center of the mean plume. Therefore, our results highlight the importance of considering the combined effect of formation heterogeneity, exposure volume, PSD, source size, and measurement location in performing risk assessment.     

Anomalously deep earthquakes in northwestern Italy

It is usually assumed that earthquakes in intraplate regions occur in the upper crust, and northwestern Italy is generally assigned to this kind of normal seismicity. In this work, the depth distribution of the events localized in this area by the Istituto Geofisico Geodetico (IGG) seismic network in the period 1991–1997 is analyzed in detail. In particular, the location capability of the network is discussed, adopting as reference quarry blasts (for the epicentral position) and the locations obtained from a dense temporary network (for the depth estimate). Within the so-obtained error limits, the depth distribution of events show a characteristic pattern: while for most of the area covered by the network the well-located seismicity lies within the first 20 km of depth, in a band following the inner arc of the Western Alps, numerous events have anomalously large focal depths, reaching a maximum of 114 km. These depth determinations cannot be attributed to instabilities of the location procedure: different choices of the propagation models used for the hypocentral determination led to very similar depth values, always significantly larger than the standard values for the surrounding areas. A strong correlation has been found between the 3-dimensional distribution of these foci and the P-wave propagation anomalies obtained from tomographic studies, suggesting a direct link between elastic and rheological properties of lower crust and upper mantle in this area.     

Sand transport measurements in Chioggia inlet,Venice lagoon: Theory versus observations

This paper presents results of recent measurements of sand transport made in Chioggia inlet as part of an extensive monitoring programme in the Venetian inlets. Measurements were made in order: (1) to define a relationship between sand transport magnitude and tidal flow; (2) to derive the thresholds for sand transport; (3) to identify the dominant modes of transport; (4) to evaluate the concentration profiles of sand within the benthic boundary layer; (5) to compare bedload transport observations with model predictions using existent bedload formulae; and (6) to produce yearly estimates of bedload transport across the inlet. The vertical distribution of sand in the water column was sampled using modified Helley–Smith bedload samplers at three sites. Transport was found to vary according to the flow and bed grain size, with considerable temporal and spatial variability. A difference of up to three orders of magnitude in transport was observed through the inlet, with higher transport rates measured on the seaward part. The dominant mode of transport in the central inlet was suspension, while bedload was dominant in the mouths. The measured profiles of sand concentration varied with the tidal stage and seabed grain size according to the Rouse parameter (R). R was high at the inlet mouths (1<R<2), indicative of a well-developed bedload layer. The inverse movability number (W_s/U_*) was also higher at these sites and appeared to be grain size dependant. Formulae for bedload transport were tested against field data; stochastic methods such as Einstein–Brown, Engelund–Hansen and Van Rijn produce the best fits. The coupled model SHYFEM-Sedtrans05 appears to simulate well observed transport for most conditions of flow. Long-term bedload predictions indicate a dominant export of sand, with a yearly average of 4500 m³.     

Giant gas bubbles in a rheomorphic vent fill at the Las Cañadas caldera,Tenerife (Canary Islands)

During rheomorphism subsequent to fallout deposition, a portion of the densely welded fallout of the La Grieta Member flowed back into the vent from where it was erupted, while the rest of it flowed down the outer slopes of the Las Cañadas caldera in Tenerife. The welded fallout and conduit-vent structure are physically connected and constitute a rare example of this type of deposits rooted to its feeder conduit and exposed in the caldera wall. The lower part of the vent-filling rheomorphic rocks shows gas bubbles and cavities that increase in size (up to 4 m) down vent. Bubbles are deformed against other bubbles, against the steep vent walls, flattened parallel to the flow foliation planes, and elongated parallel to the flow lineation and flow fold axes. The preservation of such giant bubbles, rather than their formation, seems to be a pretty unique feature of the phonolitic products investigated here and it is likely the result of the combination of factors that acted to preserve, in the surrounding of the glass transition interval, the sealing and the late stage cooling of a pressurized system. In addition, strain drop at the base of the vent-filling rheomorphic flow caused by flow stopping against vertical vent walls may have promoted rapid gas exsolution and the formation of large bubbles.