Probabilistic domain decomposition for the solution of the two-dimensional magnetotelluric problem

Probabilistic domain decomposition is proposed as a novel method for solving the two-dimensional Maxwell’s equations as used in the magnetotelluric method. The domain is split into non-overlapping sub-domains and the solution on the sub-domain boundaries is obtained by evaluating the stochastic form of the exact solution of Maxwell’s equations by a Monte-Carlo approach. These sub-domains can be naturally chosen by splitting the sub-surface domain into regions of constant (or at least continuous) conductivity. The solution over each sub-domain is obtained by solving Maxwell’s equations in the strong form. The sub-domain solver used for this purpose is a meshless method resting on radial basis function-based finite differences. The method is demonstrated by solving a number of classical magnetotelluric problems, including the quarter-space problem, the block-in-half-space problem and the triangle-in-half-space problem.     

Long-Term RST Analysis of Anomalous TIR Sequences in Relation with Earthquakes Occurred in Greece in the Period 2004–2013

Real-time integration of multi-parametric observations is expected to accelerate the process toward improved, and operationally more effective, systems for time-Dependent Assessment of Seismic Hazard (t-DASH) and earthquake short-term (from days to weeks) forecast. However, a very preliminary step in this direction is the identification of those parameters (chemical, physical, biological, etc.) whose anomalous variations can be, to some extent, associated with the complex process of preparation for major earthquakes. In this paper one of these parameters (the Earth’s emitted radiation in the Thermal InfraRed spectral region) is considered for its possible correlation with M ≥ 4 earthquakes occurred in Greece in between 2004 and 2013. The Robust Satellite Technique (RST) data analysis approach and Robust Estimator of TIR Anomalies (RETIRA) index were used to preliminarily define, and then to identify, significant sequences of TIR anomalies (SSTAs) in 10 years (2004–2013) of daily TIR images acquired by the Spinning Enhanced Visible and Infrared Imager on board the Meteosat Second Generation satellite. Taking into account the physical models proposed for justifying the existence of a correlation among TIR anomalies and earthquake occurrences, specific validation rules (in line with the ones used by the Collaboratory for the Study of Earthquake Predictability—CSEP—Project) have been defined to drive a retrospective correlation analysis process. The analysis shows that more than 93 % of all identified SSTAs occur in the prefixed space–time window around (M ≥ 4) earthquake's time and location of occurrence with a false positive rate smaller than 7 %. Molchan error diagram analysis shows that such a correlation is far to be achievable by chance notwithstanding the huge amount of missed events due to frequent space/time data gaps produced by the presence of clouds over the scene. Achieved results, and particularly the very low rate of false positives registered on a so long testing period, seems already sufficient (at least) to qualify TIR anomalies (identified by RST approach and RETIRA index) among the parameters to be considered in the framework of a multi-parametric approach to t-DASH.     

Enhancement of prestack diffraction data and attributes using a traveltime decomposition approach

Diffractions not only carry important information about small-scale subsurface structures, they also possess unique properties, which make them a powerful tool for seismic processing and imaging. Since a point diffractor scatters an incoming wave to all directions, a diffraction event implies better illumination than a reflection, because the rays travel through larger parts of the subsurface. Furthermore, unlike the reflection case, in which the emergence location of the reflected wave depends on the source position, in the case of non-Snell scattering, up-going and down-going raypaths are decoupled. Based on this decoupling, we introduce a diffraction traveltime decomposition principle, which establishes a direct connection between zero-offset and finite-offset diffraction wavefield attributes. By making use of this approach, we are able to enhance diffractions and obtain high-quality diffraction wavefield attributes at arbitrary offsets in the prestack domain solely based on zero-offset processing without any further optimization of attributes. We show the accuracy of the method by fitting diffraction traveltimes, and on simple waveform data. Application to complex synthetic data shows the ability of the proposed approach to enhance diffractions and provide high-quality wavefield attributes even in sparsely illuminated regions such as subsalt areas. The promising results reveal a high potential for improved prestack data enhancement and further applications such as efficient diffraction-based finite-offset tomography.     

SEAWAT‐Based Simulation of Axisymmetric Heat Transport

Simulation of heat transport has its applications in geothermal exploitation of aquifers and the analysis of temperature dependent chemical reactions. Under homogeneous conditions and in the absence of a regional hydraulic gradient, groundwater flow and heat transport from or to a well exhibit radial symmetry, and governing equations are reduced by one dimension (1D) which increases computational efficiency importantly. Solute transport codes can simulate heat transport and input parameters may be modified such that the Cartesian geometry can handle radial flow. In this article, SEAWAT is evaluated as simulator for heat transport under radial flow conditions. The 1971, 1D analytical solution of Gelhar and Collins is used to compare axisymmetric transport with retardation (i.e., as a result of thermal equilibrium between fluid and solid) and a large diffusion (conduction). It is shown that an axisymmetric simulation compares well with a fully three dimensional (3D) simulation of an aquifer thermal energy storage systems. The influence of grid discretization, solver parameters, and advection solution is illustrated. Because of the high diffusion to simulate conduction, convergence criterion for heat transport must be set much smaller (10^?10) than for solute transport (10^?6). Grid discretization should be considered carefully, in particular the subdivision of the screen interval. On the other hand, different methods to calculate the pumping or injection rate distribution over different nodes of a multilayer well lead to small differences only.     

Design,simulation, and large‐scale testing of an innovative vibration mitigation device employing essentially nonlinear elastomeric springs

This study proposes an innovative passive vibration mitigation device employing essentially nonlinear elastomeric springs as its most critical component. Essential nonlinearity denotes the absence (or near absence) of a linear component in the stiffness characteristics of these elastomeric springs. These devices were implemented and tested on a large‐scale nine‐story model building structure. The main focus of these devices is to mitigate structural response under impulse‐like and seismic loading when the structure remains elastic. During the design process of the device, numerical simulations, optimizations, and parametric studies of the structure‐device system were performed to obtain stiffness parameters for the devices so that they can maximize the apparent damping of the fundamental mode of the structure. Pyramidal elastomeric springs were employed to physically realize the optimized essentially nonlinear spring components. Component‐level finite element analyses and experiments were conducted to design the nonlinear springs. Finally, shake table tests using impulse‐like and seismic excitation with different loading levels were performed to experimentally evaluate the performance of the device. Experimental results demonstrate that the properly designed devices can mitigate structural vibration responses, including floor acceleration, displacement, and column strain in an effective, rapid, and robust fashion. Comparison between numerical and experimental results verified the computational model of the nonlinear system and provided a comprehensive verification for the proposed device. Copyright © 2014 John Wiley & Sons, Ltd.     

The use of mycorrhiza for eco‐engineering measures in steep alpine environments: effects on soil aggregate formation and fine‐root development

Steep erosion‐prone and vegetation‐free slopes are widespread in alpine areas and are often discussed since they have a high socio‐economic damage potential. We present an eco‐engineering approach to test whether a mycorrhizal inoculum improves the establishment of hedge brush layers and in turn soil structural stability on a steep, coarse‐grained vegetation‐free slope in the eastern Swiss Alps. We established (i) mycorrhizal and (ii) non‐mycorrhizal treated eco‐engineered research plots on a field experimental scale, covering a total area of approximately 1000 m² on an east‐northeast (ENE) exposed slope, where many environmental parameters can be regarded as homogeneous. After a full vegetation period, we quantified soil aggregate stability, the formation of water stable aggregates and the fine‐root development. Our results illustrate that the establishment of brush layers without mycorrhizal inoculum increased aggregate stability significantly. Against our expectation and glasshouse experiments, the addition of mycorrhizal inoculum did not have a statistically significant effect after one vegetation period although it tended to increase aggregate stability. Analogously, root length density (RLD) tended to be higher at the non‐mycorrhizal treated site. Aggregate stability was significantly correlated with RLD. Studies on a bigger field experimental scale are inevitable, complement glasshouse studies and lead to a better understanding for a successful application of sustainable eco‐engineering measures in alpine environments. Based on our results and considering the fact that the response time in natural ecosystems may be slower than in laboratory approaches, we conclude that long‐term field studies are necessary to validate results gained through laboratory experiments. Copyright © 2014 John Wiley & Sons, Ltd.     

The value of multiple earthquake missions: the EEFIT L’Aquila Earthquake experience

In November 2012 EEFIT launched its first ever return mission to an earthquake affected site. The L’Aquila Earthquake site was chosen as this is a recent European event of interest to the UK and European earthquake engineering community. The main aims of this return mission were to document the earthquake recovery process and this paper presents an overview of the post-disaster emergency phase and transition to reconstruction in the Aquila area after the earthquake. It takes an earthquake engineering perspective, highlighting areas mainly of interest to the fields of structural/seismic engineering and reconstruction management. Within the paper, reference is made to published literature, but also to data collected in the field during the return mission that would not otherwise have been available. The paper presents some specific observations and lessons learned from the L’Aquila return mission. However, in light of current international efforts in conducting return missions, the paper ends with some reflections on the value that return missions can provide to the field of earthquake engineering in general, based on the EEFIT L’Aquila experience.     

Estimates of gravity wave drag on Mars: Indication of a possible lower thermospheric wind reversal

Spectral gravity wave parameterization suitable for planetary thermospheres applied to wind and temperature from Mars Climate Database reveals enormously strong drag incompatible with the wind distribution. It points out to a possible wind reversal in the 110-140 km layer similar to the one in the Earth’s lower thermosphere.     

The 1.10- and 1.18-μm nightside windows of Venus observed by SPICAV-IR aboard Venus Express

Observations of the 1.10- and 1.18-μm nightside windows by the SPICAV-IR instrument aboard Venus Express were analyzed to characterize the various sources of gaseous opacity and determine the H₂O mole fraction in the lower atmosphere of Venus. We showed that the line profile model of Afanasenko and Rodin (Afanasenko, T.S., Rodin, A.V. [2007]. Astron. Lett. 33, 203–210) underestimates the CO₂ absorption in the high-wavelength wing of the 1.18-μm window and we derived an empirical lineshape that matches this wing well. An additional continuum opacity is required to reproduce the variation of the 1.10- and 1.18-μm radiances with surface elevation as observed by the VIRTIS-M instrument aboard Venus Express. A constant absorption coefficient of 0.7 ± 0.2 × 10⁻⁹ cm⁻¹ am⁻² best reproduces the observed variation. We compared spectra calculated with different CO₂ and H₂O line lists. We found that the CDSD line list lacks the 5ν₁ + ν₃ series of CO₂ bands, which provide significant opacity in Venus’ deep atmosphere, and we have constructed a composite line list that best reproduces the observations. We also showed for the first time that HDO brings significant absorption at 1140–1190 nm. Using the best representation of the atmospheric opacity we could reach, we retrieved a water vapor mole fraction of ppmv, pertaining to the altitude range 5–25 km. Combined with previous measurements in the 1.74- and 2.3-μm windows, this result provides strong evidence for a uniform H₂O profile below 40 km, in agreement with chemical models.     

Episodic orogeny deduced from coeval sedimentary sequences in the foreland basin and its implication for uplift process of Longmen Mountain, China

Longmen Mountain located at the boundary between the Sichuan Basin and Tibetan Plateau,representing the steepest gradient of any edges of the plateau.Three endmember models of uplift process and mechanism have been proposed,including crustal thickening,crustal flow,and crustal isostatic rebound.Here we use coeval sedimentary sequences in the foreland basin to restraint uplift process and mechanism in the Longmen Mountain.The more than 10,000 m thick Late TriassicQuaternary strata filled in this foreland basin and can be divided into six megasequences that are distinguished as two distinct types.The first type is the wedge-shaped megasequences which are sedimentary response of strong active thrust loading events,characterized by a high rate of subsidence and sediment accumulation,coarsening-upward succession and a dual-sourced sediment supply.This type includes Late Triassic,Late Jurassic to Early Cretaceous and Late Cretaceous to Paleogene megasequences.The second type is the tabular megasequences,characterized by the low rate of subsidence and sediment accumulation,finingupward succession,and a single-sourced sediment supply,which is sedimentary response of isostatic rebound and erosion unloading.This type includes the Early to Middle Jurassic,Middle Cretaceous and Neogene to Quaternary megasequences.Basing on sedimentary,active tectonic,geomorphic evidence,we infer that the direction has been reversed from SSWdirected sinistral strike-slip to NNE-directed dextral strike-slip during 40-3.6 Ma,and since 3.6 Ma,the Longmen Mountain thrust belt belong to times of isostatic rebound and erosional unloading with NNEdirected dextral strike-slip.This suggests that crustal isostatic rebound is a primary driver for uplift and topography of the present Longmen Mountain.The Wenchuan(Ms8.0) earthquake,which ruptured a large thrust fault with NNE-directed dextral strikeslip along the range front,is an active manifestation of this crustal isostatic rebound process with dextral strike-slipping and shortening.This process may be the cause for the Wenchuan Earthquake and the apparent paradox of high relief,little shortening,the relative dearth of historical seismicity in the region.