首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 717 毫秒
1.
Shear wave splitting analyses have been carried out using teleseismic data from broad-band seismograph stations deployed at temporary and permanent locations in Dronning Maud Land (DML), Antarctica. In most cases, the observed anisotropy can be related to major tectonic events that formed the present-day Antarctic continent. We rule out an anisotropic contribution from recent asthenospheric flow. At the Russian base Novolazarevskaya near the coast in central DML, waveform inversion suggests a two-layer model where the fast direction of the upper layer is oriented parallel to Archean fabrics in the lithosphere, whereas the anisotropy of the lower layer is interpreted to have been created during the Jurassic Gondwana break-up. Recordings at the South African base Sanae IV, however, show enigmatic results. For narrow backazimuthal segments, splitting parameters show strong variations together with a multitude of isotropic measurements, indicative of complex scattering that cannot be explained by simple one- or two-layer anisotropic models. In the interior of the continent, the data are consistent with single-layer anisotropy, but show significant spatial variations in splitting parameters. A set of temporary stations across the Heimefront shear zone in western DML yield splitting directions that we interpret as frozen anisotropy from Proterozoic assembly of the craton. An abrupt change in fast axis direction appears to mark a suture between the Grunehogna craton, a fragment of the Kalahari–Kaapvaal craton in southern Africa and the Mesoproterozoic Maudheim Province.  相似文献   

2.
Summary. Reduced Pn travel times from the Archaean Pilbara Craton of north-west Australia show a strong correlation with azimuth, which could be used as evidence of anisotropy. However, the azimuthal correlation could also be explained by a southerly dip of between 1 and 2° on the crust–mantle boundary, although the models from several reversed seismic profiles across the craton suggest a smaller dip.
A time-term analysis of the Pn date yielded several models. The preferred solution, in which the dip on the crust–mantle boundary is similar to that in the models from the reversed profiles, has approximately 2 per cent anisotropy in the uppermost mantle, with the direction of maximum velocity 30° east of north. One possible cause of the anisotropy is that olivine crystals were aligned by syntectonic recrystallization and/or power law creep in the tensional environment caused at the base of the lithosphere by flexure during loading of the lithosphere by the strata of the Hamersley Basin which overlies the Pilbara Craton.
A seismic discontinuity occurs about 15 km below the crust–mantle boundary under the craton. A qualitative analysis of all available seismic data suggests that the velocity below the boundary is probably also anisotropic, with the direction of maximum velocity between north and 40° west of north. The direction of minimum velocity below the sub-Moho boundary correlates loosely with the direction of basement lineaments in the Proterozoic Capricorn Orogenic Belt to the south of the craton, suggesting that the anisotropy under the boundary may be younger than that immediately under the crust/mantle boundary. This is consistent with the notion that the Archaean lithosphere was thinner than the present lithosphere.  相似文献   

3.
Upper-mantle flow beneath French Polynesia from shear wave splitting   总被引：1，自引：0，他引：1  
Upper-mantle flow beneath the South Pacific is investigated by analysing shear wave splitting parameters at eight permanent long-period and broad-band seismic stations and 10 broad-band stations deployed in French Polynesia from 2001 to 2005 in the framework of the Polynesian Lithosphere and Upper Mantle Experiment (PLUME). Despite the small number of events and the rather poor backazimuthal coverage due to the geographical distribution of the natural seismicity, upper-mantle seismic anisotropy has been detected at all stations except at Tahiti where two permanent stations with 15 yr of data show an apparent isotropy. The median value of fast polarization azimuths (N67.5°W) is parallel to the present Pacific absolute plate motion direction in French Polynesia (APM: N67°W). This suggests that the observed SKS fast polarization directions result mainly from olivine crystal preferred orientations produced by deformation in the sublithospheric mantle due to viscous entrainment by the moving Pacific Plate and preserved in the lithosphere as the plate cools. However, analysis of individual measurements highlights variations of splitting parameters with event backazimuth that imply an actual upper-mantle structure more complex than a single anisotropic layer with horizontal fast axis. A forward approach shows that a two-layer structure of anisotropy beneath French Polynesia better explains the splitting observations than a single anisotropic layer. Second-order variations in the measurements may also indicate the presence of small-scale lateral heterogeneities. The influence of plumes or fracture zones within the studied area does not appear to dominate the large-scale anisotropy pattern but may explain these second-order splitting variations across the network.  相似文献   

4.
High noise levels hamper teleseismic shear wave splitting measurements, which bandpass filtering does not always help. To investigate how robust splitting measurements are to noise, we analysed a set of synthetic records with known splitting parameters and added fixed levels of noise. In the presence of weak anisotropy, single-waveform splitting measurements are unreliable when operating with noisy data sets. A practical rule in terms of S/N ratio and splitting delay time parameters is that splitting is confidently detectable at S/N > 8, regardless of the wave's original polarization orientation. However, for the evidence of weak anisotropy to be detectable and measurable at an S/N value of 4, the backazimuth separation of the phases from the fast polarization direction needs to be higher than 20°. Stacks of individual measurements consistently yield reliable results down to S/N values of 4. Applying stacking to data from DSB (Dublin, Ireland), the fast polarization direction φ and lag time δt are 58° and 0.95  s. This orientation reflects surface trends of deformation in the area, as found elsewhere in the UK. Our result thus reinforces the proposed model that the detected anisotropy in the British Isles originates from lithospheric coherent deformation preserved from the last main tectonic episode.  相似文献   

5.
We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.  相似文献   

6.
Data from 90 permanent broad-band stations spread over central and eastern Europe were analysed using Ps receiver functions to study the crustal and upper-mantle structure down to the mantle transition zone. Receiver functions provide valuable information on structural features, which are important for the resolution of European lithospheric dynamics. Moho depths vary from less than 25 km in extensional areas in central Europe to more than 50 km at stations in eastern Europe (Craton) and beneath the Alpine–Carpathian belt. A very shallow Moho depth can be observed at stations in the Upper Rhine Graben area ( ca. 25 km), whereas, for example, stations in the SW Bohemian Massif show a significantly deeper Moho interface at a depth of 38 km. Vp / Vs ratios vary between 1.60 and 1.96, and show no clear correlation to the major tectonic units, thus probably representing local variations in crustal composition. Delayed arrivals of converted phases from the mantle transition zone are observed at many stations in central Europe, whereas stations in the cratonic area show earlier arrivals compared with those calculated from the IASP91 Earth reference model. Differential delay times between the P410s and P660s phases indicate a thickened mantle transition zone beneath the eastern Alps, the Carpathians and the northern Balkan peninsula, whereas the transition zone thickness in eastern and central Europe agrees with the IASP91 value. The thickening of the mantle transition zone beneath the eastern Alps and the Carpathians could be caused by cold, deeply subducted oceanic slabs.  相似文献   

7.
As a baseline measurement for understanding the Himalayan–Tibetan orogen, a product of continent–continent collision between India and Eurasia, we analyse digital seismic data in order to constrain the seismic anisotropy of the Indian shield. Based on spatially sparse data that are currently available in the public domain, there is little shear-wave birefringence for SKS phases under the Indian shield, even though it is part of a fast-moving plate in the hotspot frame of reference. If most of the northern Indian mantle has little transverse anisotropy, the onset of significant anisotropy under Tibet marks the northern terminus of intact Indian lithosphere that is thrusting under the Himalayan–Tibetan orogen. Beyond this terminus, tectonic fabric such as that associated with the deforming lithospheric mantle of Eurasia must be present in the upper mantle. Along the profile from Yadong to Golmud, the only profile in Tibet where a number of shear-wave birefringence data are available, the amount of birefringence shows two marked increases, near 30° and 33°N, between which a local high in Bouguer gravity anomaly is observed. Such a correlation between patterns of shear-wave birefringence and gravity anomalies is explained by the juxtaposition of Indian lithosphere against the overlying Eurasian lithosphere: while the Eurasian lithospheric mantle appears only to the north of 30°N, the Indian lithospheric mantle extends northwards to near 33°N.  相似文献   

8.
Ultra‐large rift basins, which may represent palaeo‐propagating rift tips ahead of continental rupture, provide an opportunity to study the processes that cause continental lithosphere thinning and rupture at an intermediate stage. One such rift basin is the Faroe‐Shetland Basin (FSB) on the north‐east Atlantic margin. To determine the mode and timing of thinning of the FSB, we have quantified apparent upper crustal β‐factors (stretching factors) from fault heaves and apparent whole‐lithosphere β‐factors by flexural backstripping and decompaction. These observations are compared with models of rift basin formation to determine the mode and timing of thinning of the FSB. We find that the Late Jurassic to Late Palaeocene (pre‐Atlantic) history of the FSB can be explained by a Jurassic to Cretaceous depth‐uniform lithosphere thinning event with a β‐factor of ~1.3 followed by a Late Palaeocene transient regional uplift of 450–550 m. However, post‐Palaeocene subsidence in the FSB of more than 1.9 km indicates that a Palaeocene rift with a β‐factor of more than 1.4 occurred, but there is only minor Palaeocene or post‐Palaeocene faulting (upper crustal β‐factors of less than 1.1). The subsidence is too localized within the FSB to be caused by a regional mantle anomaly. To resolve the β‐factor discrepancy, we propose that the lithospheric mantle and lower crust experienced a greater degree of thinning than the upper crust. Syn‐breakup volcanism within the FSB suggests that depth‐dependent thinning was synchronous with continental breakup at the adjacent Faroes and Møre margins. We suggest that depth‐dependent continental lithospheric thinning can result from small‐scale convection that thins the lithosphere along multiple offset axes prior to continental rupture, leaving a failed breakup basin once seafloor spreading begins. This study provides insight into the structure and formation of a generic global class of ultra‐large rift basins formed by failed continental breakup.  相似文献   

9.
Shear wave splitting measurements from S arrivals of local earthquakes recorded at the Incorporated Research Institutions for Seismology (IRIS) broadband sensor SNZO are used to determine a basic anisotropic structure for the subduction zone in the Wellington region. With the use of high-frequency filters, fast anisotropic polarization ( φ ) and splitting time ( δt ) measurements typical of crustal anisotropy are evident, but the larger splitting expected from the mantle is often not resolved. The small splitting seen agrees well with the results of previous studies concerning shallow crustal anisotropy. With the use of lower-frequency filters, measurements more consistent with mantle anisotropy are made. Anisotropy of 4.4 ± 0.9 per cent with a fast polarization of 29° ± 38° is calculated for the subducting slab, from 20 to 70  km depth. Using this result in addition to the results of previous studies, a model is proposed. The model requires a frequency-dependent anisotropy of less than 1.4 per cent when measured with a period of ~2  s to be present in the sub-slab mantle.
Separate from this population, a band of events in northern Cook Strait with an 86° ± 10° fast polarization is seen. This is at about 40° from the strike of the Hikurangi margin, and suggests a source of shear strain 40° removed from that found in the majority of the region. The cause of this is probably a deformation in the subducting slab in this region, as it moves towards a greater incline to the south.  相似文献   

10.
The diagenetic history of the Ediacaran sedimentary rocks in the East European Craton (EEC) over the area extending from Arkhangelsk (Russia) in the north to Podolia (Ukraine) in the south was revealed by means of the XRD characterization and K–Ar dating of clay fractions, mudstone porosity measurements and organic geochemistry investigations. Mudstone porosity measurements produced direct evidence of shallow maximum burial of the Ediacaran sediments on the craton (Russia, Lithuania, Belarus, Volyn), not exceeding 1.5 km, and much deeper burial at the cratonic margin, in Podolia and Poland. In general, illitization of smectite and biomarker indices indicates more advanced diagenesis at the cratonic margin. K–Ar dating of authigenic illite–smectite and aluminoceladonite revealed the Palaeozoic age of mineral diagenesis (ca. 450–300 Ma) both on the craton and its margin, with older ages generally observed in the north. When the maximum palaeotemperatures were evaluated from illite–smectite and biomarkers, based on the calibrations from the conventional burial diagenetic sections, a major mismatch was detected for the cratonic area: 100°C–130°C from illite––smectite and tens of oC lower from the lipid biomarkers. This diagenetic pattern was interpreted as the result of short‐lasting (in ky scale) pulses of potassium‐bearing hot fluids migrating from the Caledonian and Variscan orogens deep in the craton interior, effectively promoting illitization in porous rocks without altering the organic matter. Analogous short pulses of fluids were responsible for numerous diagenetic phenomena, including Mississippi Valley‐Type ore deposits, in the American Midwest, in front of the Appalachians. K–Ar dating indicates that the entire Proterozoic sedimentary cover of the Great Unconformity on the EEC remained untouched by measureable post‐sedimentary changes until the early Palaeozoic, thus for over 1000 My, which is an unprecedented finding.  相似文献   

11.
We infer the lithospheric structure in eastern Turkey using teleseismic and regional events recorded by 29 broad-band stations from the Eastern Turkey Seismic Experiment (ETSE). We combine the surface wave group velocities (Rayleigh and Love) with telesesimic receiver functions to jointly invert for the S -wave velocity structure, Moho depth and mantle-lid (lithospheric mantle) thickness. We also estimated the transverse anisotropy due to Love and Rayleigh velocity discrepancies. We found anomalously low shear wave velocities underneath the Anatolian Plateau. Average crustal thickness is 36 km in the Arabian Plate, 44 km in Anatolian Block and 48 km in the Anatolian Plateau. We observe very low shear wave velocities at the crustal portion (30–38 km) of the northeastern part of the Anatolian Plateau. The lithospheric mantle thickness is either not thick enough to resolve it or it is completely removed underneath the Anatolian Plateau. The shear velocities and anisotropy down to 100 km depth suggest that the average lithosphere–asthenosphere boundary in the Arabian Plate is about 90 and 70 km in Anatolian block. Adding the surface waves to the receiver functions is necessary to constrain the trade-off between velocity and the thickness. We find slower velocities than with the receiver function data alone. The study reveals three different lithospheric structures in eastern Turkey: the Anatolian plateau (east of Karliova Triple Junction), the Anatolian block and the northernmost portion of the Arabian plate. The boundary of lithospheric structure differences coincides with the major tectonic boundaries.  相似文献   

12.
A spectral method is employed to study the response to surface loads of a Maxwell earth including lateral viscosity variations. In particular, we focus on the effects of lithospheric cratons on the long-wavelength time-dependent displacement field for simple earth models. The viscosity contrast of the craton with respect to the surrounding mantle is kept fixed, whereas its thickness is allowed to vary. We show that the long-wavelength vertical displacement is not greatly affected by the presence of a lithospheric craton, while the tangential displacement is severely modified for the case of a homogeneous mantle. With increasing harmonic degree and thickness of the craton, the load-deformation coefficients deviate from those pertaining to a homogeneous mantle with a viscosity of 1021 Pa s. These deviations are particularly enhanced on timescales larger than a few hundred years. These findings indicate that the interpretation of the viscosity structure of the mantle inferred from postglacial rebound signatures based on radially stratified models is affected by the presence of lateral viscosity variations.  相似文献   

13.
Classical models of lithosphere thinning predict deep synrift basins covered by wider and thinner post‐rift deposits. However, synextensional uplift and/or erosion of the crust are widely documented in nature (e.g. the Base Cretaceous unconformity of the NE Atlantic), and generally the post‐rift deposits dominate basins fills. Accordingly, several basin models focus on this discrepancy between observations and the classical approach. These models either involve differential thinning, where the mantle thins more than the crust thereby increasing average temperature of the lithosphere, or focus on the effect of metamorphic reactions, showing that such reactions decrease the density of lithospheric rocks. Both approaches result in less synrift subsidence and increased post‐rift subsidence. The synextensional uplift in these two approaches happens only for special cases, that is for a case of initially thin crust, specific mineral assemblage of the lithospheric mantle or extensive differential thinning of the lithosphere. Here, we analyse the effects of shear heating and tectonic underpressure on the evolution of sedimentary basins. In simple 1D models, we test the implications of various mechanisms in regard to uplift, subsidence, density variations and thermal history. Our numerical experiments show that tectonic underpressure during lithospheric thinning combined with pressure‐dependent density is a widely applicable mechanism for synextensional uplift. Mineral phase transitions in the subcrustal lithosphere amplify the effect of underpressure and may result in more than 1 km of synextensional erosion. Additional heat from shear heating, especially combined with mineral phase transitions and differential thinning of the lithosphere, greatly decreases the amount of synrift deposits.  相似文献   

14.
Amplitude measurements of the transverse component of SKS waves, the so-called splitting intensity, can be used to formulate a non-linear inverse problem to image the 3-D variations of upper mantle anisotropy. Assuming transverse isotropy (or hexagonal symmetry), one can parametrize anisotropy by two anisotropic parameters and two angles describing the orientation of the symmetry axis. These can also be written as two collinear pseudo-vectors. The tomographic process consists of retrieving the spatial distribution of these pseudo-vectors, and thus resembles surface wave vectorial tomography. Spatial resolution results from the sensitivity of low-frequency SKS waves to seismic anisotropy off the ray path. The expressions for the 3-D sensitivity kernels for splitting intensity are derived, including the near-field contributions, and validated by comparison with a full wave equation solution based upon the finite element method. These sensitivity kernels are valid for any orientation of the symmetry axis, and thus generalize previous results that were only valid for a horizontal symmetry axis. It is shown that both lateral and vertical subwavelength variations of anisotropy can be retrieved with a dense array of broad-band stations, even in the case of vertically propagating SKS waves.  相似文献   

15.
Summary. Magnetotelluric soundings (MT) were conducted along the 14th parallel on the eastern border of the West African craton in the Republic of Niger.
This profile includes eight sites an average of 30 km apart. We determined the regional structure of electric conductivity and its relation to the various geological belts covered. This study took place within a 15–500 s period band. Two-dimensional modelling suggests that, in the sites located on the sedimentary basin and the mobile belt, there exists a conductive layer in the upper mantle at a depth of 80 km. Within the craton we were unable to prove the existence of this conductive layer.
There also exists another conductive layer at the crust—upper mantle boundary at a depth of 30 km, but this seems to disappear in the cratonic belt.
A significant electric discontinuity is present between the mobile belt and the sedimentary basin, due to a variation in resistivity in the substratum and a thickening in the surface cover.  相似文献   

16.
Finite deformation during fluid flow   总被引：6，自引：0，他引：6  
Summary. Typical upper mantle circulations obtained by solving Stokes' equation produce finite deformations which differ in important ways from those produced by pure or simple shear. Finite strain, defined by the ratio of the long to the short axis of the deformation ellipse, in most cases shows a steady increase with superimposed oscillations. Similarity solutions for the flow near plate boundaries demonstrate that the observed seismic anisotropy in the oceanic lithosphere can be produced by the finite deformation beneath the ridge axes. The same mechanism should give rise to strong anisotropy in the mantle above sinking slabs. Such anisotropy has not yet been detected, perhaps because the observed high velocities have been attributed to thermal effects. Convection in the mantle remote from plate boundaries produces complicated deformation which varies rapidly with position and will therefore be difficult to map seismically. The fabrics of nodules in lavas and kimberlites suggest that large strains can occur in the mantle under stresses which are too small to produce dislocation movement. The large and complicated finite deformation produced by the convective circulation in the mantle also affects closed geochemical systems, and leads to thorough mixing of any convecting region.  相似文献   

17.
Summary. It is known that flow in the mantle can produce preferred orientation in olivine crystals with seismic anisotropy as a consequence. Flow in the subcrustal lithosphere is unlikely because of the high viscosity. Lenses of high temperature and low-viscosity ( anomalous mantle ) are located under the crust in many tectonically active regions, and viscous flow can easily arise in such material resulting in seismic anisotropy. After cooling, such anomalous mantle acquires high viscosity and becomes incorporated into the lithospheric layer preserving the anisotropy produced by the flows which existed previously. The interaction of the stresses with cracks in the upper crust can be one of the causes of anisotropy in this layer.  相似文献   

18.
We describe a waveform modelling technique and demonstrate its application to determine the crust- and upper-mantle velocity structure beneath Africa. Our technique uses a parallelized reflectivity method to compute synthetic seismograms and fits the observed waveforms by a global optimization technique based on a Very Fast Simulated Annealing (VFSA). We match the S , Sp, SsPmP and shear-coupled PL phases in seismograms of deep (200–800 km), moderate-to-large magnitude (5.5–7.0) earthquakes recorded teleseismically at permanent broad-band seismic stations in Africa. Using our technique we produce P - and S -wave velocity models of crust and upper mantle beneath Africa. Additionally, our use of the shear-coupled PL phase, wherever observed, improves the constraints for lower crust- and upper-mantle velocity structure beneath the corresponding seismic stations. Our technique retains the advantages of receiver function methods, uses a different part of the seismogram, is sensitive to both P - and S -wave velocities directly, and obtains helpful constraints in model parameters in the vicinity of the Moho. The resulting range of crustal thicknesses beneath Africa (21–46 km) indicates that the crust is thicker in south Africa, thinner in east Africa and intermediate in north and west Africa. Crustal P - (4.7–8 km s−1) and S -wave velocities (2.5–4.7  km s−1) obtained in this study show that in some parts of the models, these are slower in east Africa and faster in north, west and south Africa. Anomalous crustal low-velocity zones are also observed in the models for seismic stations in the cratonic regions of north, west and south Africa. Overall, the results of our study are consistent with earlier models and regional tectonics of Africa.  相似文献   

19.
李栋梁  邵鹏程  王慧 《中国沙漠》2013,33(5):1511-1519
利用1951—2009年中国662个气象台站逐日降水资料和NCEP/NCAR逐日再分析资料,分析了东亚副热带夏季风北边缘地理位置的年际、年代际变化以及夏季风影响北边缘的时间演变特征。结果表明:基于过程透雨量定义的夏季风北边缘能很好地反映边缘带降水特点,东亚副热带夏季风北边缘带大致呈东北-西南走向,在110°E处,北边缘带南界大致位于36°N,北界大致位于41°N。边缘带宽度与其南界纬度具有很好的负相关,且有逐年变宽的趋势。边缘带存在显著的准3年和9年的周期特征。夏季风平均在6月28日开始影响北边缘带,9月21日南撤离开北边缘带。110°E是夏季风边缘带中南风分量最大的经度,与偏南年相比,偏北年夏季风影响边缘带时间开始早,结束迟。  相似文献   

20.
A 3-D P -velocity map of the crust and upper mantle beneath the southeastern part of India has been reconstructed through the inversion of teleseismic traveltimes. Salient geological features in the study region include the Archean Dharwar Craton and Eastern Ghat metamorphic belt (EGMB), and the Proterozoic Cuddapah and Godavari basins. The Krishna–Godavari basin, on the eastern coastal margin, evolved in response to the Indo–Antarctica breakup. A 24-station temporary network provided 1161 traveltimes, which were used to model 3-D P -velocity variation. The velocity model accounts of 80 per cent of the observed data variance. The velocity picture to a depth of 120 km shows two patterns: a high velocity beneath the interior domain (Dharwar craton and Cuddapah basin), and a lower velocity beneath the eastern margin region (EGMB and coastal basin). Across the array velocity variations of 7–10 per cent in the crust (0–40 km) and 3–5 per cent in the uppermost mantle (40–120 km) are observed. At deeper levels (120–210 km) the upper-mantle velocity differences are insignificant among different geological units. The presence of such a low velocity along the eastern margin suggests significantly thin lithosphere (<100 km) beneath it compared to a thick lithosphere (>200 km) beneath the eastern Dharwar craton. Such lithospheric thinning could be a consequence of Indo–Antarctica break-up.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号