ALBERT

Description: The Tonga-Kermadec forearc is deforming in response to on-going subduction of the Pacific Plate beneath the Indo-Australian Plate. Previous research has focussed on the structural development of the forearc where large bathymetric features such as the Hikurangi Plateau and Louisville Ridge seamount chain are being subducted. Consequently, knowledge of the ‘background’ forearc in regions of normal plate convergence is limited. We report on an ~250-km-long multichannel seismic reflection profile that was shot perpendicular to the Tonga-Kermadec trench at ~28°S to determine the lateral and temporal variations in the structure, stratigraphy and deformation of the Kermadec forearc resulting solely from Pacific Plate subduction. Interpretation of the seismic profile, in conjunction with regional swath bathymetry data, shows that the Pacific Plate exhibits horst and graben structures that accommodate bending-induced extensional stresses, generated as the trenchward dip of the crust increases. Trench infill is also much thicker than expected at 1 km which, we propose, results from increased sediment flux into and along the trench. Pervasive normal faulting of the mid-trench slope most likely accommodates the majority of the observed forearc extension in response to basal subduction erosion, and a structural high is located between the mid- and upper-trench slopes. We interpret this high as representing a dense and most likely structurally robust region of crust lying beneath this region. Sediment of the upper-trench slope documents depositional hiatuses and on-going uplift of the arc. Strong along-arc currents appear to erode the Kermadec volcanic arc and distribute this sediment to the surrounding basins, while currents over the forearc redistribute deposits as sediment waves. Minor uplift of the transitional Kermadec forearc, observed just to the north of the profile, appears to relate to an underlying structural trend as well as subduction of the Louisville Ridge seamount chain 250 km to the north. Relative uplift of the Kermadec arc is observed from changes in the tilt of upper-trench slope deposits and extensional faulting of the basement immediately surrounding the Louisville Ridge.

Description: New methods are required to combine the information contained in the passive electrical and seismic signals to detect, localize and monitor hydromechanical disturbances in porous media. We propose a field experiment showing how passive seismic and electrical data can be combined together to detect a preferential flow path associated with internal erosion in a Earth dam. Continuous passive seismic and electrical (self-potential) monitoring data were recorded during a 7-d full-scale levee (earthen embankment) failure test, conducted in Booneschans, Netherlands in 2012. Spatially coherent acoustic emissions events and the development of a self-potential anomaly, associated with induced concentrated seepage and internal erosion phenomena, were identified and imaged near the downstream toe of the embankment, in an area that subsequently developed a series of concentrated water flows and sand boils, and where liquefaction of the embankment toe eventually developed. We present a new 4-D grid-search algorithm for acoustic emissions localization in both time and space, and the application of the localization results to add spatially varying constraints to time-lapse 3-D modelling of self-potential data in the terms of source current localization. Seismic signal localization results are utilized to build a set of time-invariant yet spatially varying model weights used for the inversion of the self-potential data. Results from the combination of these two passive techniques show results that are more consistent in terms of focused ground water flow with respect to visual observation on the embankment. This approach to geophysical monitoring of earthen embankments provides an improved approach for early detection and imaging of the development of embankment defects associated with concentrated seepage and internal erosion phenomena. The same approach can be used to detect various types of hydromechanical disturbances at larger scales.

Description: To estimate the seismic hazard, the geometry (dip, length and orientation) and the dynamics (type of displacements and amplitude) of the faults in the area of interest need to be understood. In this paper, in addition to geomorphologic observations, we present the results of two ground penetrating radar (GPR) campaigns conducted in 2010 and 2011 along the Emeelt fault in the vicinity of Ulaanbaatar, capital of Mongolia, located in an intracontinental region with low deformation rate that induces long recurrence time between large earthquakes. As the geomorphology induced by the fault activity has been highly smoothed by erosion processes since the last event, the fault location and geometry is difficult to determine precisely. However, by using GPR first, a non-destructive and fast investigation, the fault and the sedimentary deposits near the surface can be characterized and the results can be used for the choice of trench location. GPR was performed with a 50 MHz antenna over 2-D lines and with a 500 MHz antenna for pseudo-3-D surveys. The 500 MHz GPR profiles show a good consistency with the trench observations, dug next to the pseudo-3-D surveys. The 3-D 500 MHz GPR imaging of a palaeochannel crossed by the fault allowed us to estimate its lateral displacement to be about 2 m. This is consistent with a right lateral strike-slip displacement induced by an earthquake around magnitude 7 or several around magnitude 6. The 2-D 50 MHz profiles, recorded perpendicular to the fault, show a strong reflection dipping to the NE, which corresponds to the fault plane. Those profiles provided complementary information on the fault such as its location at shallow depth, its dip angle (from 23° to 35°) and define its lateral extension.

Description: The method presented here aims to assess the tsunami threat very rapidly after the occurrence of a large earthquake, using as input the parameters of the seismic source, and an approach based on Green's summation. We show that the main weakness of the approach (the need to consider only linear shallow water propagation) is largely compensated by the advantages in terms of computing performance and independence with respect to pre-computed scenarios. To test the approach and to illustrate its implementation in a real environment, we focus on the Sea of Oman, a tsunamigenic area characterized by Makran subduction zone which detailed structure is partially unknown and where secondary tsunami sources must also be taken into account, both for hazard studies and warning purposes. The potential source area is partitioned into a grid of unity water sources. A shallow water (SW) numerical model is used to pre-compute the corresponding empirical Green's functions on several points of interest located on the coasts of Iran, Pakistan and Oman. The comparison between Green's summation and the direct SW computation using the full resolution of the bathymetric grid shows that the accuracy is good enough for practical applications.

Description: We present a scale- and parameter-adaptive method to pre-condition the gradient of the parameters to be inverted in time-domain 2-D elastic full-waveform inversion (FWI). The proposed technique, which relies on a change of variables of the model parameters, allows to balance the value of the gradient of the Lamé parameters and density throughout the model in each step of the multiscale inversion. The main difference compared to existing gradient pre-conditioners is that the variables are automatically selected based on a least-squares minimization criteria of the gradient weight, which corresponds to the product of the gradient by a power of the parameter to be inverted. Based on numerical tests made with (1) a modified version of the Marmousi-2 model, and (2) a high-velocity and density local anomaly model, we illustrate that the value of the power helps to balance the gradient throughout the model. In addition, we show that a particular value exists for each parameter that optimizes the inversion results in terms of accuracy and efficiency. For the two models, the optimal power is ~2.0–2.5 and ~1.5 for the first and second Lamé parameters, respectively; and between 3 and 6, depending on the inverted frequency, for density. These power values provide the fastest and most accurate inversion results for the three parameters in the framework of multiscale and multishooting FWI using three different optimization schemes.