ALBERT

Description: Remote Sensing Letters, Volume 3, Issue 8, Page 697-706, 20 December 2012.

Description: SUMMARY In the absence of seismological measurements, observations of the topography and gravity fields of solid planets are the primary constraints on their internal structure. To compute the synthetic geoid and topography induced by the dynamics of planetary interiors, we introduce a 3-D numerical tool describing mantle convection beneath an elastic lithosphere. Although the energy conservation is treated in the whole spherical domain, the deformation aspect is solved using a hybrid technique (finite volume method for the viscous flow, spectral method for elastic deformation). The mechanical coupling is achieved via the imposition of the traction at the surface of the viscous flow as a basal boundary condition for the elastic deformation. We present both response functions and full thermal convection cases computed with our new method for planetary bodies of varying dimensions: the filtering effect of the lithosphere on the dynamic topography and geoid is specific for each planetary body, justifying the importance of such a tool. Furthermore, since our approach specifically focuses on the mechanical coupling at the base of the lithosphere, it will permit future, more elaborate, rheological treatments. It also enables to discriminate between the radial and tangential components of the viscous traction. The latter is found to have a significant influence on the elastic deformation. The effect on geoid is prominent. More specifically, while a thin elastic lithosphere is usually considered to play little role on the dynamic topography and geoid of Venus, a ∼35 per cent reduction is obtained for geoid height in the numerical example we propose. On a planet with thicker elastic lithosphere such as Mars, the consequence of this filtering effect is to rule out the possibility of a dynamical support for the Tharsis Rise, even for the lowest admissible values of elastic thickness in this region.

Description: SUMMARY The Central Iran plateau appears aseismic during the last few millenniums based on instrumental and historical seismic records. Nevertheless, it is sliced by several strike-slip faults that are hundreds of kilometres long. These faults display along-strike, horizontal offsets of intermittent gullies that suggest the occurrence of earthquakes in the Holocene. Establishing this is crucial for accurately assessing the regional seismic hazard. The first palaeoseismic study performed on the 200-km long, NS striking Anar fault shows that this right-lateral fault hosted three large ( M w ≈ 7) earthquakes during the Holocene or possibly Uppermost Pleistocene for the older one. These three seismic events are recorded within a sedimentary succession, which is not older than 15 ka, suggesting an average recurrence of at most 5 ka. The six optically stimulated luminescence ages available provide additional constraints and allow estimating that the three earthquakes have occurred within the following time intervals: 4.4 ± 0.8, 6.8 ± 1 and 9.8 ± 2 ka. The preferred age of the more recent event, ranging between 3600 and 5200 yr, suggests that the fault is approaching the end of its seismic cycle and the city of Anar could be under the threat of a destructive earthquake in the near future. In addition, our results confirm a previous minimum slip rate estimate of 0.8 ± 0.1 mm yr −1 for the Anar fault indicating that the westernmost prominent right-lateral faults of the Central Iran plateau are characterized by slip rates close to 1 mm yr −1 . These faults, which have repeatedly produced destructive earthquakes with large magnitudes and long recurrence interval of several thousands of years during the Holocene, show that the Central Iran plateau does not behave totally as a rigid block and that its moderate internal deformation is nonetheless responsible for a significant seismic hazard.

Description: SUMMARY Passive seismic interferometry is a new promising methodology for seismic exploration. Interferometry allows information about the subsurface structure to be extracted from ambient seismic noise. In this study, we apply the cross-correlation technique to approximately 25 hr of recordings of ambient seismic noise at the Ketzin experimental CO 2 storage site, Germany. Common source gathers were generated from the ambient noise for all available receivers along two seismic lines by cross-correlation of noise records. This methodology isolates the interstation Green's functions that can be directly compared to active source gathers. We show that the retrieved response includes surface waves, refracted waves and reflected waves. We use the dispersive behaviour of the retrieved surface waves to infer geological properties in the shallow subsurface and perform passive seismic imaging of the subsurface structure by processing the retrieved reflected waves.

Description: SUMMARY Normal mode observations play an important role in studying broad-scale lateral variations in the Earth. Such studies require the calculation of accurate synthetic spectra in realistic earth models, and this remains a computationally challenging problem. Here, we describe a new implementation of the direct solution method for calculating normal mode spectra in laterally heterogeneous earth models. In this iterative direct solution method , the mode-coupling equations are solved in the frequency-domain using the preconditioned biconjugate gradient algorithm, and the time-domain solution is recovered using a numerical inverse Fourier transform. A number of example calculations are presented to demonstrate the accuracy and efficiency of the method for performing large ‘full coupling’ calculations as compared to methods based on matrix diagonalization and the traditional direct solution method.

Description: SUMMARY We report on a unique set of infrasound observations from a single earthquake, the 2011 January 3 Circleville earthquake ( M w 4.7, depth of 8 km), which was recorded by nine infrasound arrays in Utah. Based on an analysis of the signal arrival times and backazimuths at each array, we find that the infrasound arrivals at six arrays can be associated to the same source and that the source location is consistent with the earthquake epicentre. Results of propagation modelling indicate that the lack of associated arrivals at the remaining three arrays is due to path effects. Based on these findings we form the working hypothesis that the infrasound is generated by body waves causing the epicentral region to pump the atmosphere, akin to a baffled piston. To test this hypothesis, we have developed a numerical seismoacoustic model to simulate the generation of epicentral infrasound from earthquakes. We model the generation of seismic waves using a 3-D finite difference algorithm that accounts for the earthquake moment tensor, source time function, depth and local geology. The resultant acceleration–time histories on a 2-D grid at the surface then provide the initial conditions for modelling the near-field infrasonic pressure wave using the Rayleigh integral. Finally, we propagate the near-field source pressure through the Ground-to-Space atmospheric model using a time-domain Parabolic Equation technique. By comparing the resultant predictions with the six epicentral infrasound observations from the 2011 January 3, Circleville earthquake, we show that the observations agree well with our predictions. The predicted and observed amplitudes are within a factor of 2 (on average, the synthetic amplitudes are a factor of 1.6 larger than the observed amplitudes). In addition, arrivals are predicted at all six arrays where signals are observed, and importantly not predicted at the remaining three arrays. Durations are typically predicted to within a factor of 2, and in some cases much better. These results suggest that measured infrasound from the Circleville earthquake is consistent with the generation of infrasound from body waves in the epicentral region.

Description: SUMMARY We present new dedicated core surface field models spanning the decade from 2000.0 to 2010.0. These models, called gufm-sat , are based on CHAMP, Ørsted and SAC-C satellite observations along with annual differences of processed observatory monthly means. A spatial parametrization of spherical harmonics up to degree and order 24 and a temporal parametrization of sixth-order B-splines with 0.25 yr knot spacing is employed. Models were constructed by minimizing an absolute deviation measure of misfit along with measures of spatial and temporal complexity at the core surface. We investigate traditional quadratic or maximum entropy regularization in space, and second or third time derivative regularization in time. Entropy regularization allows the construction of models with approximately constant spectral slope at the core surface, avoiding both the divergence characteristic of the crustal field and the unrealistic rapid decay typical of quadratic regularization at degrees above 12. We describe in detail aspects of the models that are relevant to core dynamics. Secular variation and secular acceleration are found to be of lower amplitude under the Pacific hemisphere where the core field is weaker. Rapid field evolution is observed under the eastern Indian Ocean associated with the growth and drift of an intense low latitude flux patch. We also find that the present axial dipole decay arises from a combination of subtle changes in the southern hemisphere field morphology.

Description: SUMMARY We present a novel frequency-domain inverse solution to recover the 3-D electrical conductivity distribution in the mantle. The solution is based on analysis of local C -responses. It exploits an iterative gradient-type method—limited-memory quasi-Newton method—for minimizing the penalty function consisting of data misfit and regularization terms. The integral equation code is used as a forward engine to calculate responses and data misfit gradients during inversion. An adjoint approach is implemented to compute misfit gradients efficiently. Further improvements in computational load come from parallelizing the scheme with respect to frequencies, and from setting the most time-consuming part of the forward calculations—calculation of Green’s tensors—apart from the inversion loop. Convergence, performance, and accuracy of our 3-D inverse solution are demonstrated with a synthetic numerical example. A companion paper applies the strategy set forth here to real data.

Description: Remote Sensing Letters, Volume 3, Issue 8, Page 687-696, 20 December 2012.

Description: Marine Geodesy, Volume 35, Issue 1, Page 49-65, January-March 2012.