Adaptively parametrized surface wave tomography: methodology and a new model of the European upper mantle

In this study, we aim to close the gap between regional and global traveltime tomography in the context of surface wave tomography of the upper mantle implementing the principle of adaptive parametrization. Observations of seismic surface waves are a very powerful tool to constrain the 3-D structure of the Earth's upper mantle, including its anisotropy, because they sample this volume efficiently due to their sensitivity over a wide depth range along the ray path. On a global scale, surface wave tomography models are often parametrized uniformly, without accounting for inhomogeneities in data coverage and, as a result, in resolution, that are caused by effective under- or overparametrization in many areas. If the local resolving power of seismic data is not taken into account when parametrizing the model, features will be smeared and distorted in tomographic maps, with subsequent misinterpretation. Parametrization density has to change locally, for models to be robustly constrained without losing any accurate information available in the best sampled regions. We have implemented a new algorithm for upper mantle surface wave tomography, based on adaptive-voxel parametrization, with voxel size defined by both the ‘hit count' (number of observations sampling the voxel) and ‘azimuthal coverage' (how well different azimuths with respect to the voxel are covered by the source-station distribution). High image resolution is achieved in regions with dense data coverage, while lower image resolution is kept in regions where data coverage is poorer. This way, parametrization is everywhere tuned to optimal resolution, minimizing both the computational costs, and the non-uniqueness of the solution. The spacing of our global grid is locally as small as ∼50 km. We apply our method to identify a new global model of vertically and horizontally polarized shear velocity, with resolution particularly enhanced in the European lithosphere and upper mantle. We find our new model to resolve lithospheric thickness and radial anisotropy better than earlier results based on the same data. Robust features of our model include, for example, the Trans-European Suture Zone, the Panonnian Basin, thinned lithosphere in the Aegean and Western Mediterranean, possible small-scale mantle upwellings under Iberia and Massif Central, subduction under the Aegean arc and a very deep cratonic root underneath southern Finlan

Three dimensional interface modelling with two-dimensional seismic data: the Alpine crust-mantle boundary

We present a new approach to determine the 3-D topography and lateral continuity of seismic interfaces using 2-D-derived controlled-source seismic reflector data. The aim of the approach is to give the simplest possible structure consistent with all reflector data and error estimates. We define simplicity of seismic interfaces by the degree of interface continuity (i.e. shortest length of offsets) and by the degree of interface roughness (least surface roughness). The method is applied to structural information of the crust—mantle boundary (Moho) obtained from over 250 controlled-source seismic reflection and refraction profiles in the greater Alpine region. The reflected and refracted phases from the Moho interface and their interpretation regarding crustal thickness are reviewed and their reliability weighted. Weights assigned to each reflector element are transformed to depth errors considering Fresnel volumes. The 2-D-derived reflector elements are relocated in space (3-D migration) and interpolation is performed between the observed reflector elements to obtain continuity of model parameters. Interface offsets are introduced only where required according to the principle of simplicity. The resulting 3-D model of the Alpine crust—mantle boundary shows two offsets that divide the interface into a European, an Adriatic and a Ligurian Moho, with the European Moho subducting below the Adriatic Moho, and with the Adriatic Moho underthrusting the Ligurian Moho. Each sub-interface depicts the smoothest possible (i.e. simplest) surface, fitting the reflector data within their assigned errors. The results are consistent with previous studies for those regions with dense and reliable controlled-source seismic data. The newly derived Alpine Moho interface, however, surpasses earlier studies by its lateral extent over an area of about 600 km by 600 km, by quantifying reliability estimates along the interface, and by obeying the principle of being consistently as simple as possibl

Consistent phase picking for regional tomography models: application to the greater Alpine region

The resolution and reliability of tomographic velocity models strongly depends on quality and consistency of available traveltime data. Arrival times routinely picked by network analysts on a day-to-day basis often yield a high level of noise due to mispicks and other inconsistencies, particularly in error assessment. Furthermore, tomographic studies at regional scales require merging of phase picks from several networks. Since a common quality assessment is not usually available for phase data provided by different networks, additional inconsistencies are introduced by the merging process. Considerable improvement in the quality of phase data can only be achieved through complete repicking of seismograms. Considering the amount of data necessary for regional high-resolution tomography, algorithms combining accurate picking with an automated error assessment represent the best tool to derive large suitable data sets. In this work, we present procedures for consistent automated and routine picking of P-wave arrival times at local to regional scales including consistent picking error assessment. Quality-attributed automatic picks are derived from the MPX picking system. The application to earthquakes in the greater Alpine region demonstrates the potential of such a repicking approach. The final data set consists of more than 13 000 high-quality first-arrivals and it is used to derive regional 1-D and preliminary 3-D P-wave models of the greater Alpine region. The comparison with a tomographic model based on routine phase data extracted from the ISC Bulletin illustrates effects on tomographic results due to consistency and reliability of our high-quality data se

Historical Reconstruction of Population Density of the Echinoid Diadema antillarum on Florida Keys Shallow Bank-Barrier Reefs

The 1983–1984 caribbean-wide mass mortality of the once ubiquitous long-spined sea urchin Diadema antillarum Philippi, 1845, is one of several factors considered responsible for coral reef change throughout the region. Unfortunately, there is a paucity of pre-mortality event density data for D. antillarum in the Florida Keys, making it difficult to determine pre-1983 population density levels. results from surveys conducted during 1970–1973 in the lower Florida Keys, in shallow (\u3c12 m) fore reef habitats, yielded relatively abundant and widespread D. antillarum densities in qualitative transects at five reefs prior to the 1983–1984 die-off. In quantitative surveys at one reef, Middle sambo reef in 1972, up to 7.9 individuals m–2 were recorded using quadrats in high-relief spur and groove habitat. A second mortality event in the Florida Keys, beginning in April 1991, again depressed urchin densities that had begun to recover from the 1983–1984 mass mortality. By 1992, D. antillarum densities (\u3c0.01 m–2) were two orders of magnitude lower than pre-die-off estimates (range of 0.07–0.57 m–2 from several spur and groove reefs in the lower Florida Keys) and remained so through 2009. The pre-mortality echinoid density estimates detailed in the Florida Keys provide a baseline to compare with their current population status and should help inform managers about realistic recovery or restoration targets for D. antillarum

Automatic seismic phase picking and consistent observation error assessment: application to the Italian seismicity

Accuracy of seismic phase observation and consistency of timing error assessment define the quality of seismic waves arrival times. High-quality and large data sets are prerequisites for seismic tomography to enhance the resolution of crustal and upper mantle structures. In this paper we present the application of an automated picking system to some 600 000 seismograms of local earthquakes routinely recorded and archived by the Italian national seismic network. The system defines an observation weighting scheme calibrated with a hand-picked data subset and mimics the picking by an expert seismologist. The strength of this automatic picking is that once it is tuned for observation quality assessment, consistency of arrival times is strongly improved and errors are independent of the amount of data to be picked. The application to the Italian local seismicity documents that it is possible to automatically compile a precise, homogeneous and large data set of local earthquake Pg and Pn arrivals with related polarities. We demonstrate that such a data set is suitable for high-precision earthquake location, focal mechanism determination and high-resolution seismic tomograph

Automatic seismic phase picking and consistent observation error assessment: application to Italian seismicity

Accuracy of seismic phase observation and consistency of timing error assessment define the quality of seismic waves arrival times. High-quality and large data sets are prerequisites for seismic tomography to enhance the resolution of crustal and upper mantle structures. In this paperwe present the application of an automated picking system to some 600000 seismograms of local earthquakes routinely recorded and archived by the Italian national seismic network. The system defines an observation weighting scheme calibrated with a hand-picked data subset and mimics the picking by an expert seismologist. The strength of this automatic picking is that once it is tuned for observation quality assessment, consistency of arrival times is strongly improved and errors are independent of the amount of data to be picked. The application to the Italian local seismicity documents that it is possible to automatically compile a precise, homogeneous and large data set of local earthquake Pg and Pn arrivals with related polarities. We demonstrate that such a data set is suitable for high-precision earthquake location, focal mechanism determination and high-resolution seismic tomography