共查询到20条相似文献，搜索用时 31 毫秒
1.
Elastic finite element models are applied to investigate the effects of topography and medium heterogeneities on the surface deformation and the gravity field produced by volcanic pressure sources. Changes in the gravity field cannot be interpreted only in terms of gain of mass disregarding the ground deformation of the rocks surrounding the source. Contributions to gravity changes depend also on surface and subsurface mass redistribution driven by dilation of the volcanic source. Both ground deformation and gravity changes were firstly evaluated by solving a coupled axisymmetric problem to estimate the effects of topography and medium heterogeneities. Numerical results show significant discrepancies in the ground deformation and gravity field compared to those predicted by analytical solutions, which disregard topography, elastic heterogeneities and density subsurface structures. With this in mind, we reviewed the expected gravity changes accompanying the 1993–1997 inflation phase on Mt Etna by setting up a fully 3-D finite element model in which we used the real topography, to include the geometry, and seismic tomography, to infer the crustal heterogeneities. The inflation phase was clearly detected by different geodetic techniques (EDM, GPS, SAR and levelling data) that showed a uniform expansion of the overall volcano edifice. When the gravity data are integrated with ground deformation data and a coupled FEM modelling was solved, a mass intrusion could have occurred at depth to justify both ground deformation and gravity observations. 相似文献
2.
Forward modelling receiver functions for crustal structure beneath station TBZ (Trabzon, Turkey) 总被引：4，自引：0，他引：4
Özcan Çak&#;r Murat Erduran Hakan Ç&#;nar & Ahmet Y&#;lmaztürk 《Geophysical Journal International》2000,140(2):341-356
We use teleseismic three-component digital data from the Trabzon, Turkey broadband seismic station TBZ to model the crustal structure by the receiver function method. The station is located at a structural transition from continental northeastern Anatolia to the oceanic Black Sea basin. Rocks in the region are of volcanic origin covered by young sediments. By forward modelling the radial receiver functions, we construct 1-D crustal shear velocity models that include a lower crustal low-velocity zone, indicating a partial melt mechanism which may be the source of surfacing magmatic rocks and regional volcanism. Within the top 5 km, velocities increase sharply from about 1.5 to 3.5 km s−1 . Such near-surface low velocities are caused by sedimentation, extending from the Black Sea basin. Velocities at around 20 km depth have mantle-like values (about 4.25 km s−1 ), which easily correlate to magmatic rocks cropping out on the surface. At 25 km depth there is a thin low-velocity layer of about 4.0 km s−1 . The average Moho velocity is about 4.6 km s−1 , and its depth changes from 32 to 40 km. Arrivals on the tangential components indicate that the Moho discontinuity dips approximately southwards, in agreement with the crustal thickening to the south. We searched for the solution of receiver functions around the regional surface wave group velocity inversion results, which helped alleviate the multiple solution problem frequently encountered in receiver function modelling.
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region. 相似文献
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region. 相似文献
3.
We have developed an inversion formuialion for velocity and attenuation structure using seismic attributes, including envelope amplitude, instantaneous frequency and arrival times of selected seismic phases. We refer to this approach as AFT inversion for amplitude, (instantaneous) frequency and time. Complex trace analysis is used to extract the different seismic attributes. The instantaneous frequency data are converted to t * using a matching procedure that approximately removes the effects of the source spectra. To invert for structure, ray-perturbation methods are used to compute the sensitivity of the seismic attributes to variations in the model. An iterative inversion procedure is then performed from smooth to less smooth models that progressively incorporates the shorter-wavelength components of the model. To illustrate the method, seismic attributes are extracted from seismic-refraction data of the Ouachita PASSCAL experiment and used to invert for shallow crustal velocity and attenuation structure. Although amplitude data are sensitive to model roughness, the inverted velocity and attenuation models were required by the data to maintain a relatively smooth character. The amplitude and t * data were needed, along with the traveltimes, at each step of the inversion in order to fit all the seismic attributes at the final iteration. 相似文献
4.
The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect. 相似文献
5.
Summary. Moment tensor inversion methods can be applied with success in the determination of source properties of simple earthquakes. However, these methods utilize the assumption of a point source, which is inadequate for modelling many complicated, shallow earthquakes. For complex earthquakes, an inversion using finite faulting models is desirable but the number of parameters involved requires that a good starting model be found or that independent constraints be placed on some of the parameters. A method is presented for low-pass filtering both the data and Green's functions, passing only signals with wavelengths greater than the dimension of the entire fault. The filter tends to smooth complications in the waveforms and allows application of the point source moment tensor inversion. This method is applied to body waves from the 1978 Thessaloniki, Greece, earthquake, the 1971 San Fernando earthquake and to a multiple-point source synthetic model of the San Fernando event. For the Thessaloniki event, although a multiple-source mechanism has been suggested, inversion results before and after filtering were essentially identical, indicating that a point source mechanism is sufficient in modelling the long-period, teleseismic body waves. In the case of the San Fernando earthquake, the point source Green's functions were incapable of simultaneously modelling the P - and SH -waves. Inversion of P -waves alone resulted in extreme parameter resolution problems, but allowed constraint in one axis of the moment tensor and suggested an overall source time function. Inversion of a synthetic San Fernando data set yielded similar results, but allowed an investigation of the shortcomings of the method under controlled circumstances. Although the results may require substantial interpretation, the method presented represents a simple first step in the analysis of complex earthquakes. 相似文献
6.
Summary. In order to examine the development of the oceanic crust in the neighbourhood of a slowly spreading ridge, a seismic refraction experiment was carried out at 59° 30'N on the Reykjanes Ridge. Three 120 km long overlapped split profiles were shot parallel to the trend of the ridge, on the eastern flank, and recorded on up to five recording sonobuoys. The profiles were at distances of 0, 30 and 90km from the ridge axis, corresponding to approximate crustal ages of 0, 3 and 9 Myr. Data from the main profiles were supplemented by using a large chamber air gun during recovery of the buoys.
The analysis of the data combined standard travel-time interpretation, the 'tau' method of systematic travel-time inversion and detailed amplitude modelling using the Reflectivity Method to calculate synthetic seismograms. Detailed velocity-depth models were constructed for each of the profiles.
There is no indication of a significant magma chamber at the ridge crest, although a slight velocity inversion in layer 3 suggests a zone of elevated temperature. Away from the crest there was a slight positive velocity gradient in layer 3. Layer 2 was most effectively modelled by a region of varying velocity gradients, which thinned with age and the transition to layer 3 is marked by a sharp change in velocity gradient. The transition to mantle velocities is also best modelled by a high velocity gradient rather than an interface.
Although some lateral variation in properties is apparent along the profiles, the lateral velocity gradients were sufficiently weak to allow an effective analysis in terms of laterally uniform models. 相似文献
The analysis of the data combined standard travel-time interpretation, the 'tau' method of systematic travel-time inversion and detailed amplitude modelling using the Reflectivity Method to calculate synthetic seismograms. Detailed velocity-depth models were constructed for each of the profiles.
There is no indication of a significant magma chamber at the ridge crest, although a slight velocity inversion in layer 3 suggests a zone of elevated temperature. Away from the crest there was a slight positive velocity gradient in layer 3. Layer 2 was most effectively modelled by a region of varying velocity gradients, which thinned with age and the transition to layer 3 is marked by a sharp change in velocity gradient. The transition to mantle velocities is also best modelled by a high velocity gradient rather than an interface.
Although some lateral variation in properties is apparent along the profiles, the lateral velocity gradients were sufficiently weak to allow an effective analysis in terms of laterally uniform models. 相似文献
7.
R. S. Spratt 《Geophysical Journal International》1982,71(1):173-186
Summary. The 1973 Hawaii earthquake occurred north of Hilo, at a depth of 40 to 50km. The location was beneath the east flank of Mauna Kea, a volcano dormant historically, but active within the last 4000 yr. Aftershocks were restricted to a depth of 55–35km. The event and its aftershock sequence are located in an area not normally associated with the seismicity of the Mauna Loa and Kilauea calderas. The earthquake was a double event, the epicentres trending NE-SW. The events were of similar size and faulting mechanism. The fault plane solutions obtained by seismic waveform analysis are a strike-slip fault striking EW and dipping 55° S, the auxiliary plane a NS vertical plane with a faulting plunge of 35°. The axis of maximum compressive stress is aligned with the direction of the gravity gradient associated with the island of Hawaii. The fault plane striking EW parallels a surface feature, the Mauna Kea east rift zone. The earthquakes were clearly not associated with volcanic activity normally associated with Mauna Loa and Kilauea and may indicate a deep seated prelude to a resumption of activity at Mauna Kea. 相似文献
8.
Peter M. Shearer 《Geophysical Journal International》1993,115(3):878-904
Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies. 相似文献
9.
About 50 000 P and S arrival times and 25 000 values of t * recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P , S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp , Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting. 相似文献
10.
We present a neural network approach to invert surface wave data for a global model of crustal thickness with corresponding uncertainties. We model the a posteriori probability distribution of Moho depth as a mixture of Gaussians and let the various parameters of the mixture model be given by the outputs of a conventional neural network. We show how such a network can be trained on a set of random samples to give a continuous approximation to the inverse relation in a compact and computationally efficient form. The trained networks are applied to real data consisting of fundamental mode Love and Rayleigh phase and group velocity maps. For each inversion, performed on a 2°× 2° grid globally, we obtain the a posteriori probability distribution of Moho depth. From this distribution any desired statistic such as mean and variance can be computed. The obtained results are compared with current knowledge of crustal structure. Generally our results are in good agreement with other crustal models. However in certain regions such as central Africa and the backarc of the Rocky Mountains we observe a thinner crust than the other models propose. We also see evidence for thickening of oceanic crust with increasing age. In applications, characterized by repeated inversion of similar data, the neural network approach proves to be very efficient. In particular, the speed of the individual inversions and the possibility of modelling the whole a posteriori probability distribution of the model parameters make neural networks a promising tool in seismic tomography. 相似文献
11.
We compare 3-D upper mantle anisotropic structures beneath the North American continent obtained using standard and improved crustal corrections in the framework of Non-linear Asymptotic Coupling Theory (NACT) applied to long period three component fundamental and higher mode surface waveform data. Our improved approach to correct for crustal structure in high-resolution regional waveform tomographic models goes beyond the linear perturbation approximation, and is therefore more accurate in accounting for large variations in Moho topography within short distances as observed, for instance, at ocean–continent margins. This improved methodology decomposes the shallow-layer correction into a linear and non-linear part and makes use of 1-D sensitivity kernels defined according to local tectonic structure, both for the forward computation and for the computation of sensitivity kernels for inversion. The comparison of the 3-D upper mantle anisotropic structures derived using the standard and improved crustal correction approaches shows that the model norm is not strongly affected. However, significant variations are observed in the retrieved 3-D perturbations. The largest differences in the velocity models are present below 250 km depth and not in the uppermost mantle, as would be expected. We suggest that inaccurate crustal corrections preferentially map into the least constrained part of the model and therefore accurate corrections for shallow-layer structure are essential to improve our knowledge of parts of the upper mantle where our data have the smallest sensitivity. 相似文献
12.
Seismic tomographic imaging of P- and S-waves velocity perturbations in the upper mantle beneath Iran 总被引：2，自引：0，他引：2
The inverse tomography method has been used to study the P - and S -waves velocity structure of the crust and upper mantle underneath Iran. The method, based on the principle of source–receiver reciprocity, allows for tomographic studies of regions with sparse distribution of seismic stations if the region has sufficient seismicity. The arrival times of body waves from earthquakes in the study area as reported in the ISC catalogue (1964–1996) at all available epicentral distances are used for calculation of residual arrival times. Prior to inversion we have relocated hypocentres based on a 1-D spherical earth's model taking into account variable crustal thickness and surface topography. During the inversion seismic sources are further relocated simultaneously with the calculation of velocity perturbations. With a series of synthetic tests we demonstrate the power of the algorithm and the data to reconstruct introduced anomalies using the ray paths of the real data set and taking into account the measurement errors and outliers. The velocity anomalies show that the crust and upper mantle beneath the Iranian Plateau comprises a low velocity domain between the Arabian Plate and the Caspian Block. This is in agreement with global tomographic models, and also tectonic models, in which active Iranian plateau is trapped between the stable Turan plate in the north and the Arabian shield in the south. Our results show clear evidence of the mainly aseismic subduction of the oceanic crust of the Oman Sea underneath the Iranian Plateau. However, along the Zagros suture zone, the subduction pattern is more complex than at Makran where the collision of the two plates is highly seismic. 相似文献
13.
Massif Central (France): new constraints on the geodynamical evolution from teleseismic tomography 总被引：2，自引：0，他引：2
M. Granet G. Stoll J. Dorel U. Achauer G. Poupinet K. Fuchs 《Geophysical Journal International》1995,121(1):33-48
The Massif Central, the most significant geomorphological unit of the Hercynian belt in France, is characterized by graben structures which are part of the European Cenozoic Rift System (ECRIS) and also by distinct volcanic episodes, the most recent dated at 20 Ma to 4000 years BP. In order to study the lithosphere-asthenosphere system beneath this volcanic area, we performed a teleseismic field experiment.
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase. 相似文献
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase. 相似文献
14.
Several years of broad-band teleseismic data from the GRSN stations have been analysed for crustal structure using P -to- S converted waves at the crustal discontinuities. An inversion technique was developed which applies the Thomson-Haskell formalism for plane waves without slowness integration. The main phases observed are Moho conversions, their multiples in the crust, and conversions at the base of the sediments. The crustal thickness derived from these data is in good agreement with results from other studies. For the Gräfenberg stations, we have made a more detailed comparison of our model with a previously published model obtained from refraction seismic experiments. The refraction seismic model contains boundaries with strong velocity contrasts and a significant low-velocity zone, resulting in teleseismic waveforms that are too complicated as compared to the observed simple waveforms. The comparison suggests that a significant low-velocity zone is not required and that internal crustal boundaries are rather smooth. 相似文献
15.
Upper crustal shear velocity models from higher mode Rayleigh wave dispersion in Scotland 总被引：1，自引：0，他引：1
Summary. Group velocities for first and second higher mode Rayleigh waves, in the frequency range 0.8–4.8 Hz, generated from a local earthquake of magnitude 3.7 M L in western Scotland, are measured at stations along the 1974 LISPB line. These provide detailed information about the crustal structure west of the line. The data divide the region into seven apparently homogeneous provinces. Averaged higher mode velocity dispersion curves for each province are analysed simultaneously using a linearized inversion technique, yielding regionalized shear velocity profiles down to a depth of 17 km into the upper crust. Shear wave velocity is between 3.0 and 3.4 km s−1 in the upper 2 km, with a slow increase to around 3.8 km s−1 . P -wave models computed using these results agree with profiles from the LISPB and LUST refraction experiments. 相似文献
16.
Volcano seismology often deals with rather shallow seismic sources and seismic stations deployed in their near field. The complex stratigraphy on volcanoes and near‐field source effects have a strong impact on the seismic wavefield, complicating the interpretation techniques that are usually employed in earthquake seismology. In addition, as most volcanoes have a pronounced topography, the interference of the seismic wavefield with the stress‐free surface results in severe waveform perturbations that affect seismic interpretation methods. In this study we deal predominantly with the surface effects, but take into account the impact of a typical volcano stratigraphy as well as near‐field source effects. We derive a correction term for plane seismic waves and a plane‐free surface such that for smooth topographies the effect of the free surface can be totally removed. Seismo‐volcanic sources radiate energy in a broad frequency range with a correspondingly wide range of different Fresnel zones. A 2‐D boundary element method is employed to study how the size of the Fresnel zone is dependent on source depth, dominant wavelength and topography in order to estimate the limits of the plane wave approximation. This approximation remains valid if the dominant wavelength does not exceed twice the source depth. Further aspects of this study concern particle motion analysis to locate point sources and the influence of the stratigraphy on particle motions. Furthermore, the deployment strategy of seismic instruments on volcanoes, as well as the direct interpretation of the broad‐band waveforms in terms of pressure fluctuations in the volcanic plumbing system, are discussed. 相似文献
17.
Terry C. Wallace 《Geophysical Journal International》1986,85(2):349-363
Summary . Seismograms recorded at regional distances (2°–12°) are quite complicated due to the waveguide nature of the crust. Generalized ray theory can be used to model the body waves in this distance range but a very large number of rays is required. Here I present a series of approximations to streamline generalized ray theory for the waveguide problem. If a layer over a half-space is used for the structure, then the de Hoop contour for a given ray is most strongly dependent on the fastest velocity of any leg of the ray. This results in analytic approximations to locate the contour. Each ray has two body wave arrivals (a headwave and a reflected arrival) so the displacement response of the ray need only be evaluated at a few points in time about the two arrival times and interpolated in between. A change in structure (increasing crustal thickness or Pn velocity) most strongly affects the relative timing of the headwave and the reflected arrival, so it is possible to 'stretch' or 'squeeze' the waveform of a representative model to simulate a whole suite of models.
Also discussed is the applicability of a single layer over a half-space structure for modelling the observed regional distance waveforms for shallow earthquakes. At periods greater than a few seconds crustal layering can be replaced by a single layer having the appropriate average velocities. Lateral variations in crustal thickness with scale lengths of less than about 100 km can also be modelled with a simple horizontal layer of appropriate average thickness. 相似文献
Also discussed is the applicability of a single layer over a half-space structure for modelling the observed regional distance waveforms for shallow earthquakes. At periods greater than a few seconds crustal layering can be replaced by a single layer having the appropriate average velocities. Lateral variations in crustal thickness with scale lengths of less than about 100 km can also be modelled with a simple horizontal layer of appropriate average thickness. 相似文献
18.
Analysis of the crustal velocity structure of the British Isles using teleseismic receiver functions
J. P. Tomlinson P. Denton P. K. H. Maguire D. C. Booth 《Geophysical Journal International》2006,167(1):223-237
The onshore crustal and upper mantle velocity structure of the British Isles has been investigated by teleseismic receiver function analysis. The results of the study augment the dense offshore and sparse onshore models of the velocity structure beneath the area. In total almost 1500 receiver functions have been analysed, which have been calculated using teleseismic data from 34 broadband and short-period, three-component seismic recording instruments. The crustal structure has primarily been investigated using 1-D grid search and forward modelling techniques, returning crustal thicknesses, bulk crustal Vp / Vs ratio and velocity-depth models. H −κ stacking reveals crustal thicknesses between 25 and 36 km and Vp / Vs ratios between 1.6 and 1.9. The crustal thicknesses correlate with the results of previous seismic reflection and refraction profiles to within ±2 km. The significant exceptions are the stations close to the Iapetus Suture where the receiver function crustal thicknesses are up to 5 km less than the seismic refraction Moho. This mismatch could be linked to the presence of underplated magmatic material at the base of the crust. 1-D forward modelling has revealed subcrustal structures in northern Scotland. These correlate with results from other UK receiver function studies, and correspond with the Flannan and W-reflectors. The structures are truncated or pinch out before they reach the Midland Valley of Scotland. The isolated subcrustal structure at station GIM on the Isle of Man may be related to the closure of the Iapetus Ocean. 相似文献
19.
Timing of depth‐dependent lithosphere stretching on the S. Lofoten rifted margin offshore mid‐Norway: pre‐breakup or post‐breakup? 总被引：1，自引：0，他引：1
Depth‐dependent stretching, in which whole‐crustal and whole‐lithosphere extension is significantly greater than upper‐crustal extension, has been observed at both non‐volcanic and volcanic rifted continental margins. A key question is whether depth‐dependent stretching occurs during pre‐breakup rifting or during sea‐floor spreading initiation and early sea‐floor spreading. Analysis of post‐breakup thermal subsidence and upper‐crustal faulting show that depth‐dependent lithosphere stretching occurs on the outer part of the Norwegian volcanic rifted margin. For the southern Lofoten margin, large breakup lithosphere β stretching factors approaching infinity are required within 100 km of the continent–ocean boundary to restore Lower Eocene sediments and flood basalt surfaces (～54 Ma) to interpreted sub‐aerial depositional environments at sea level as indicated by well data. For the same region, the upper crust shows no significant Palaeocene and Late Cretaceous faulting preceding breakup with upper‐crustal β stretching factors <1.05. Further north on the Lofoten margin, reverse modelling of post‐breakup subsidence with a β stretching factor of infinity predicts palaeo‐bathymetries of ～1500 m to the west of the Utrøst Ridge and fails to restore Lower Eocene sediments and flood basalt tops to sea level at ～54 Ma. If these horizons were deposited in a sub‐aerial depositional environment, as indicated by well data to the south, an additional subsidence event younger than 54 Ma is required compatible with lower‐crustal thinning during sea‐floor spreading initiation. For the northern Vøring margin, breakup lithosphere β stretching factors of ～2.5 are required to restore Lower Eocene sediments and basalts to sea level at deposition, while Palaeocene and Late Cretaceous upper‐crustal β stretching factors for the same region are < 1.1. The absence of significant Palaeocene and late Cretaceous extension on the southern Lofoten and northern Vøring margins prior to continental breakup supports the hypothesis that depth‐dependent stretching of rifted margin lithosphere occurs during sea‐floor spreading initiation or early sea‐floor spreading rather than during pre‐breakup rifting. 相似文献
20.
Germn Bayona Mauricio Baquero Catalina Ramírez Manuela Tabares Ana M. Salazar Giovanny Nova Edward Duarte Andrs Pardo Angelo Plata Carlos Jaramillo Guillermo Rodríguez Victor Caballero Agustín Cardona Camilo Montes Sebastin Gmez Marulanda Andrs L. Crdenas‐Rozo 《Basin Research》2021,33(1):809-845
The onset of deformation in the northern Andes is overprinted by subsequent stages of basin deformation, complicating the examination of competing models illustrating potential location of earliest synorogenic basins and uplifts. To establish the width of the earliest northern Andean orogen, we carried out field mapping, palynological dating, sedimentary, stratigraphic and provenance analyses in Campanian to lower Eocene units exposed in the northern Eastern Cordillera of Colombia (Cocuy region) and compare the results with coeval succession in adjacent basins. The onset of deformation is recorded in earliest Maastrichtian time, as terrigenous detritus arrived into the basin marking the end of chemical precipitation and the onset of clastic deposition produced by the uplift of a western source area dominated by shaly Cretaceous rocks. Disconformable contacts within the upper Maastrichtian to middle Palaeocene succession document increasing supply of quartzose sandy detritus from Cretaceous quartzose rocks exposed in eastern source areas. The continued unroofing of both source areas produced a rapid shift in depositional environments from shallow marine in Maastrichtian to fluvial‐lacustrine systems during the Palaeocene‐early Eocene. Supply of immature Jurassic sandstones from nearby western uplifts, together with localized plutonic and volcanic Cretaceous rocks, caused a shift in Palaeocene sandstones composition from quartzarenites to litharenites. Supply of detrital sandy fragments, unstable heavy minerals and Cretaceous to Ordovician detrital zircons, were derived from nearby uplifted blocks and from SW fluvial systems within the synorogenic basin, instead of distal basement rocks. The presence of volcanic rock fragments and 51–59 Ma volcanic zircons constrain magmatism within the basin. The Maastrichtian–Palaeocene sequence studied here documents crustal deformation that correlates with coeval deformation farther south in Ecuador and Peru. Slab flattening of the subducting Caribbean plate produced a wider orogen (>400 km) with a continental magmatic arc and intra‐basinal deformation and magmatism. 相似文献