Subduction-zone structure and magmatic processes beneath Costa Rica constrained by local earthquake tomography and petrological modelling

A high-quality data set of 3790 earthquakes were simultaneously inverted for hypocentre locations and 3-D P-wave velocities in Costa Rica. Tests with synthetic data and resolution estimates derived from the resolution matrix indicate that the velocity model is well constrained in central Costa Rica to a depth of 70 km; northwestern and southeastern Costa Rica are less well resolved owing to a lack of seismic stations and seismicity. Maximum H2O content and seismic wave speeds of mid-ocean ridge basalt and harzburgite were calculated for metamorphic phase transformations relevant to subduction. Both the 3-D P-wave velocity structure and petrological modelling indicate the existence of low-velocity hydrous oceanic crust in the subducting Cocos Plate beneath central Costa Rica. Intermediate-depth seismicity correlates well with the predicted locations of hydrous metamorphic rocks, suggesting that dehydration plays a key role in generating intermediate-depth earthquakes beneath Costa Rica. Wadati-Benioff zone seismicity beneath central Costa Rica shows a remarkable decrease in maximum depth toward southeastern Costa Rica. The presence of asthenosphere beneath southeastern Costa Rica, which entered through a proposed slab window, may explain the shallowing of seismicity due to increased temperatures and associated shallowing of dehydration of the slab. Tomographic images further constrain the existence of deeply subducted seamounts beneath central Costa Rica. Large, low P-wave velocity areas within the lower crust are imaged beneath the southeasternmost volcanoes in central Costa Rica. These low velocities may represent anomalously hot material or even melt associated with active volcanism in central Costa Rica. Tomographic images and petrological modelling indicate the existence of a shallow, possibly hydrated mantle wedge beneath central Costa Ric

A lithospheric cross-section through the Swiss Alps—I. Thermokinematic modelling of the Neoalpine orogeny

In this paper we develop a forward 2-D thermokinematic model to investigate the Neoalpine 35-0 Ma phase of orogeny along the European Geotraverse (EGT) through the Swiss Alps on a crustal and lithospheric scale. Using a divergence-free kinematic model (div v=0), we define mass displacements, which subsequently serve as input to a transient thermal model. the thermal model uses critically assessed material prorameters and accounts for the depth dependence of the thermal properties in processes such as crustal thickening and mantle-lithospheric subduction. Based on the presentday density pattern of the deep seismic image and estimated exhumation and shortening rates, we derive, in a first modelling step, a mass-displacement field describing the Neoalpine orogeny as a uniform process in time. In a second—thermal—modelling step, this kinematic scenario is further refined by modelling the non-uniform cooling histories of the southern Lepontine in the Penninic domain. For that purpose we adopt lithospheric shortening rates—and consequently exhumation rates—to agree with total Neoalpine shortening, while keeping the geometry of the kinematic model fixed. the resultant thermokinematic model reflects the main characteristics of Neoalpine tectonics, and shows a good overall agreement with combined geological and geophysical data. the asymmetric feature of the present-day tectonic structure along the profile is strongly reflected in the thermal structure of the lithosphere. This demonstrates the need for a kinematic model to investigate the deep-temperature field in active tectonic provinces. For further refinement of the model, the amounts of shortening have to be more precisely estimated, and a higher spatial density in geochronological and metamorphic data is required. Furthermore, surface heat-flow values are, up to now, too uncertain to constrain the predicted surface heat flow. In summary, our results show that we need, in particular, data constraining the horizontal component of the tectonic and thermal evolution. the results of the Neoalpine orogeny modelling demonstrate that the presented thermokinematic procedure yields a good first-order approximation to investigate crustal-scale and lithospheric processes. We conclude. therefore, that the approach presented provides the potential for application not only to continent-continent collision zones, but also to any active tectonic provinc

Three-dimensional P-wave velocity structure on the shallow part of the Central Costa Rican Pacific margin from local earthquake tomography using off- and onshore networks

The Central Costa Rican Pacific margin is characterized by a high-seismicity rate, coincident with the subduction of rough-relief ocean floor and has generated earthquakes with magnitude up to seven in the past. We inverted selected P-wave traveltimes from earthquakes recorded by a combined on- and offshore seismological array deployed during 6 months in the area, simultaneously determining hypocentres and the 3-D tomographic velocity structure on the shallow part of the subduction zone (<70 km). The results reflect the complexity associated to subduction of ocean-floor morphology and the transition from normal to thickened subducting oceanic crust. The subducting slab is imaged as a high-velocity perturbation with a band of low velocities (LVB) on top encompassing the intraslab seismicity deeper than ∼30 km. The LVB is locally thickened by the presence of at least two subducted seamounts beneath the margin wedge. There is a general eastward widening of the LVB over a relatively short distance, closely coinciding with the onset of an inverted forearc basin onshore and the appearance of an aseismic low-velocity anomaly beneath the inner forearc. The latter coincides spatially with an area of the subaerial forearc where differential uplift of blocks has been described, suggesting tectonic underplating of eroded material against the base of the upper plate crust. Alternatively, the low velocities could be induced by an accumulation of upward migrating fluids. Other observed velocity perturbations are attributed to several processes taking place at different depths, such as slab hydration through outer rise faulting, tectonic erosion and slab dehydratio

Strong paleoclimatic legacies in current plant functional diversity patterns across Europe

Numerous studies indicate that environmental changes during the late Quaternary have elicited long‐term disequilibria between species diversity and environment. Despite its importance for ecosystem functioning, the importance of historical environmental conditions as determinants of FD (functional diversity) remains largely unstudied. We quantified the geographic distributions of plant FD (richness and dispersion) across Europe using distribution and functional trait information for 2702 plant species. We then compared the importance of historical and contemporary factors to determine the relevance of past conditions as predictors of current plant FD in Europe. For this, we compared the strength of the relationships between FD with temperature and precipitation stability since the LGM (Last Glacial Maximum), accessibility to LGM refugia, and contemporary environmental conditions (climate, productivity, soil, topography, and land use). Functional richness and dispersion exhibited geographic patterns with strong associations to the environmental history of the region. The effect size of accessibility to LGM refugia and climate stability since the LGM was comparable to that of the contemporary predictors. Both functional richness and dispersion increased with temperature stability since the LGM and accessibility to LGM refugia. Functional richness' geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days, land use heterogeneity, diversity of soil types, and absolute minimum winter temperature. Functional dispersion's geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days and absolute minimum winter temperature. The high explained variance and model support of historical predictors are consistent with the idea that long‐term variability in environmental conditions supplements contemporary factors in shaping FD patterns at continental scales. Given the importance of FD for ecosystem functioning, future climate change may elicit not just short‐term shifts in ecosystem functioning, but also long‐term functional disequilibria

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

Shallow subduction beneath Italy

This paper presents a velocity model of the Italian (central Mediterranean) lithosphere in unprecedented detail. The model is derived by inverting a set of 166,000 Pg and Pn seismic wave arrival times, restricted to the highest-quality data available. The tomographic images reveal the geometry of the subduction-collision system between the European, Adriatic, and Tyrrhenian plates, over a larger volume and with finer resolution than previous studies. We find two arcs of low-Vp anomalies running along the Alps and the Apennines, describing the collision zones of underthrusting continental lithospheres. Our results suggest that in the Apennines, a significant portion of the crust has been subducted below the mountain belt. From the velocity model we can also infer thermal softening of the crustal wedge above the subducting Adriatic plate. In the Tyrrhenian back-arc region, strong and extensive low-Vp anomalies depict upwelling asthenospheric material. The tomographic images also allow us to trace the boundary between the Adriatic and the Tyrrhenian plates at Moho depth, revealing some tears in the Adriatic-Ionian subducting lithosphere. The complex lithospheric structure described by this study is the result of a long evolution; the heterogeneities of continental margins, lithospheric underthrusting, and plate indentation have led to subduction variations, slab tears, and asthenospheric upwelling at the present day. The high-resolution model provided here greatly improves our understanding of the central Mediterranean’s structural puzzle. The results of this study can also shed light on the evolution of other regions experiencing both oceanic and continental subduction

A PROSPECTIVE CONTROLLED STUDY OF LOW BACK SCHOOL IN THE GENERAL POPULATION

There are no data on the efficacy of a back school in primary prevention of back pain in the general population or on the characteristics of the population who volunteers. After announcement in the local press, 494 healthy adults volunteered and paid for a back school course in Switzerland. A total of 371 controls were matched for sex, age, profession, nationality and back pain. A statistically significant decrease in numbers of doctor's visits was found by the participants during the following 6 months compared with the controls. However, there were no significant between-group differences in the four remaining parameters (presence and intensity of back pain, drug intake and sick leave). Three-quarters of participants changed their attitudes after the back school. Volunteering for a back pain prevention programme was associated with the presence of back pain problems. Reasons for volunteering are further discussed. Overall, the results of this study showed that a back school for the general population may not solve the problem of low back pain, but improves self-help in a subgroup of the populatio

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

Reference seismic velocity Earth model for Italy from local source tomography and 30 years of controlled source seismology data

We present here a new high resolution regional P-wave velocity model for the lithosphere beneath the Italian region obtained by including information on the Moho topography, and integrating results from local earthquake tomography with 30 years of CSS data, applying the method of Waldhauser (1996). For the 3D moho map, we extended the crustal model, already available for the Alps by Lippitsch et al., 2003, to the Italian peninsula, Corsica, Sardinia, and Sicily. The tomographic model is obtained by inverting 166,000 Pg and Pn arrival times large part of which have been automatically picked and consistently weighted with an advanced automatic picking system (Aldersons, 2004). The resolution of the obtained velocity model is consistently higher and the grid spacing consistently smaller than in previous tomographic works targeting the same region. We are able to image the complex geometry of this part of the subduction-collision system between the Eurasian and African plates adding important details to the overview derived by the teleseismic tomography. Our results clearly show the plate boundary at Moho level from the Alps to the Southern Apennines and the Calabrian Arc in a volume unresolved in previous studies. The use of global 1D velocity models based on the flat Earth assumption is a pre-requisite to refine and interpret images and seismic responses of the earth obtained with geophysical studies (P and S tomography, surface wave tomography etc). Our model is suitable as a good starting point for a 3D velocity reference model of the crust and upper mantle beneath the Mediterranean area to be extended to the Adriatic Sea and to the Ionian Sea, with benefit for earthquakes location,teleseismic tomography, focal mechanisms and CM