首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
In this study we image crustal structure beneath a magmatic continental rift to understand the interplay between crustal stretching and magmatism during the late stages of continental rifting: the Main Ethiopian Rift (MER). The northern sector of this region marks the transition from continental rifting in the East African Rift to incipient seafloor spreading in the southern Red Sea and western Gulf of Aden. Our local tomographic inversion exploits 172 broad-band instruments covering an area of 250 × 350 km of the rift and adjacent plateaux. The instruments recorded a total of 2139 local earthquakes over a 16-month period. Several synthetic tests show that resolution is good between 12 and 25 km depth (below sea level), but some horizontal velocity smearing is evident along the axis of the Main Ethiopian Rift below 16 km. We present a 3-D P -wave velocity model of the mid-crust and present the first 3-D Vp / Vs model of the region. Our models show high P -wave velocities (6.5 km s−1) beneath the axis of the rift at a depth of 12–25 km. The presence of high Vp / Vs ratios (1.81–1.84) at the same depth range suggest that they are cooled mafic intrusions. The high Vp / Vs values, along with other geophysical evidence, suggest that dyking is pervasive beneath the axis of the rift from the mid-crustal depths to the surface and that some portion of partial melt may exist at lower crustal depths. Although the crustal stretching factor across the Main Ethiopian Rift is ∼1.7, our results indicate that magma intrusion in narrow zones accommodates a large proportion of extensional strain, with similarities to slow-spreading mid-ocean ridge processes.  相似文献   

2.
Rayleigh wave phase velocity maps in southern Africa are obtained at periods from 6 to 40 s using seismic ambient noise tomography applied to data from the Southern Africa Seismic Experiment (SASE) deployed between 1997 and 1999. These phase velocity maps are combined with those from 45 to 143 s period which were determined previously using a two-plane-wave method by Li & Burke. In the period range of overlap (25–40 s), the ambient noise and two-plane-wave methods yield similar phase velocity maps. Dispersion curves from 6 to 143 s period were used to estimate the 3-D shear wave structure of the crust and uppermost mantle on an 1°× 1° grid beneath southern Africa to a depth of about 100 km. Average shear wave velocity in the crust is found to vary from 3.6 km s–1 at 0–10 km depths to 3.86 km s–1 from 20 to 40 km, and velocity anomalies in these layers correlate with known tectonic features. Shear wave velocity in the lower crust is on average low in the Kaapvaal and Zimbabwe cratons and higher in the surrounding Proterozoic terranes, such as the Limpopo and the Namaqua-Natal belts, which suggests that the lower crust underlying the Archean cratons is probably less mafic than beneath the Proterozoic terranes. Crustal thickness estimates agree well with a previous receiver function study of Nair et al. . Archean crust is relatively thin and light and underlain by a fast uppermost mantle, whereas the Proterozoic crust is thick and dense with a slower underlying mantle. These observations are consistent with the southern African Archean cratons having been formed by the accretion of island arcs with the convective removal of the dense lower crust, if the foundering process became less vigorous in arc environments during the Proterozoic.  相似文献   

3.
4.
Summary. Attention has recently been focused on the structure and composition of the lower crust in continental areas. It is generally believed that, except in special circumstances, ductile behaviour below mid-crustal depths precludes the brittle processes that cause earthquakes. The 1984 July 19 earthquake in North Wales occurred at the unexpected depth of 23 km. We report here the location of the larger aftershocks and the relocation of the main shock with respect to one of them. The lower crustal depths of the events are confirmed by tests with a wide range of models. The occurrence of earthquakes at these depths may be related to low heat flow in the region.  相似文献   

5.
A seismic-array study of the continental crust and upper mantle in the Ivrea-Yerbano and Strona-Ceneri zones (northwestern Italy) is presented. A short-period network is used to define crustal P - and S -wave velocity models from earthquakes. The analysis of the seismic-refraction profile LOND of the CROP-ECORS project provided independent information and control on the array-data interpretation.
Apparent-velocity measurements from both local and regional earthquakes, and time-term analysis are used to estimate the velocity in the lower crust and in the upper mantle. The geometry of the upper-lower crust and Moho boundaries is determined from the station delay times.
We have obtained a three-layer crustal seismic model. The P -wave velocity in the upper crust, lower crust and upper mantle is 6.1±0.2 km s−1, 6.5±0.3 km s−1 and 7.8±0.3 km s−1 respectively. Pronounced low-velocity zones in the upper and lower crust are not observed. A clear change in the velocity structure between the upper and lower crust is documented, constraining the petrological interpretation of the Ivrea-type reflective lower continental crust derived from small-scale petrophysical data. Moreover, we found a V P/ V S ratio of 1.69±0.04 for the upper crust and 1.82±0.08 for the lower crust and upper mantle. This is consistent with the structural and petrophysical differences between a compositionally uniform and seismically transparent upper crust and a layered and reflective lower crust. The thickness of the lower crust ranges from about 8 km in front of the Ivrea body (ARVO, Arvonio station) in the northern part of the array to a maximum of about 15 km in the southern part of the array. The lower crust reaches a minimum depth of 5 km below the PROV (Provola) station.  相似文献   

6.
An Mw 5.9 earthquake occurred in the Lake Rukwa rift, Tanzania, on 1994 August 18, and was well recorded by 20 broad-band seismic stations at distances of 160 to 800 km and 21 broad-band stations at teleseismic distances. The regional and teleseismic waveforms have been used to investigate the source characteristics of the main shock, and also to locate aftershocks that occurred within three weeks of the main shock. Teleseismic body-wave modelling yields the following source parameters for the main shock: source depth of 25 ± 2 km, a normal fault orientation, with a horizontal tension axis striking NE-SW and an almost vertical pressure axis (Nodal Plane I: strike 126°–142°, dip 63°–66°, and rake 280°–290°; Nodal Plane II: strike 273°–289°, dip 28°–31°, and rake 235°–245°), a scalar moment of 4.1 times 1017 N m, and a 2 s impulsive source time function. Four of the largest aftershocks also nucleated at depths of 25 km, as deduced from regional sPmp–Pmp times. The nodal planes are broadly consistent with the orientation of both the Lupa and Ufipa faults, which bound the Rukwa rift to the northeast and southwest, respectively. The rupture radius of the main shock, assuming a circular fault, is estimated to be 4 km with a corresponding stress drop of 6.5 MPa. Published estimates of crustal thickness beneath the Rukwa rift indicate that the foci of the main shock and aftershocks lie at least 10 km above the Moho. The presence of lower-crustal seismicity beneath the Rukwa rift suggests that the pre-rift thermal structure of the rifted crust has not been strongly modified by the rifting, at least to depths of 25 km.  相似文献   

7.
Crustal structure and local seismicity in western Anatolia   总被引:4,自引:0,他引:4  
Western Anatolia is one of the most seismically active continental regions in the world and much of it has been undergoing NS-directed extensional deformation since Early Miocene time. In a cooperative study, seismologists from Saint Louis University, USA and Dokuz Eylül University, Turkey, deployed five broad-band and 45 short-period seismic stations in western Anatolia between 2002 November and 2003 October. The present paper uses data collected by this network and the data from five permanent stations operated by the Kandilli Observatory and Earthquake Research Institute to map the hypocentral distribution of local earthquakes and to determine crustal structure of western Anatolia. We obtained a 1-D P -wave crustal velocity model using a generalized scheme for simultaneously obtaining earthquake locations and a crustal velocity model. Our velocity model is characterized by crustal velocities that are significantly lower than average continental values. The low velocities may be associated with high crustal temperatures, a high degree of fracture, or the presence of fluids at high pore pressure in the crust. We located 725 local earthquakes and classified them in three categories. We found that the level of seismic activity in western Anatolia is higher than previously reported. Station delays resulting from the inversion process correlate with near-surface geology and the thickness of sediments throughout the region. The hypocentral distribution of the events indicates that peak seismicity for the region occurs at depths of about 10 km.  相似文献   

8.
Focal mechanisms determined from moment tensor inversion and first motion polarities of the Himalayan Nepal Tibet Seismic Experiment (HIMNT) coupled with previously published solutions show the Himalayan continental collision zone near eastern Nepal is deforming by a variety of styles of deformation. These styles include strike-slip, thrust and normal faulting in the upper and lower crust, but mostly strike-slip faulting near or below the crust–mantle boundary (Moho). One normal faulting earthquake from this experiment accommodates east–west extension beneath the Main Himalayan Thrust of the Lesser Himalaya while three upper crustal normal events on the southern Tibetan Plateau are consistent with east–west extension of the Tibetan crust. Strike-slip earthquakes near the Himalayan Moho at depths >60 km also absorb this continental collision. Shallow plunging P -axes and shallow plunging EW trending T -axes, proxies for the predominant strain orientations, show active shearing at focal depths ∼60–90 km beneath the High Himalaya and southern Tibetan Plateau. Beneath the southern Tibetan Plateau the plunge of the P -axes shift from vertical in the upper crust to mostly horizontal near the crust–mantle boundary, indicating that body forces may play larger role at shallower depths than at deeper depths where plate boundary forces may dominate.  相似文献   

9.
Deep seismic reflection studies in Israel - an update   总被引:1,自引:0,他引:1  
Summary. The results of three deep crustal reflection lines are presently available from Israel. A 90 km line from near the Dead Sea rift to the Mediterranean coast was carried out for deep study. Two other lines in the Mediterranean coastal area were derived by recorrelation of oil exploration lines. The data shows a division between continental inner Israel and the coastal plain. In the first area a reflective lower crust is apparent with transparent upper crust and almost transparent upper mantle. Near the coast, in an area which was previously suggested as underlain by an ancient fossil oceanic crust, strong reflections characterize the uppermost mantle. Comparison between the reflection pattern and previous deep refraction and MT data indicates some agreement away from the coast and lack of correlation in the area of possible fossil oceanic crust near the coast.  相似文献   

10.
The results of deep reflection profiling studies carried out across the palaeo-meso-Proterozoic Delhi Fold Belt (DFB) and the Archaean Bhilwara Gneissic Complex (BGC) in the northwest Indian platform are discussed in this paper. This region is a zone of Proterozoic collision. The collision appears to be responsible for listric faults in the upper crust, which represent the boundaries of the Delhi exposures. In these blocks the lower crust appears to lie NW of the respective surface exposures and the reflectivity pattern does not correspond to the exposed blocks. A fairly reflective lower crust northwest of the DFB exposures appears to be the downward continuation of the DFB upper crust. The poorly reflective lower crust under the exposed DFB may be the westward extension of the BGC upper crust at depth. Thus, the lower crust in this region can be divided into the fairly reflective Marwar Basin (MB)-DFB crust and a poorly reflective BGC crust. Vertically oriented igneous intrusions may have disturbed the lamellar lower-crustal structure of the BGC, resulting in a dome-shaped poorly reflective lower crust whose base, not traceable in the reflection data, may have a maximum depth of about 50 km, as indicated by the gravity modelling.
The DFB appears to be a zone of thick (45-50 km) crust where the lower crust has doubled in width. This has resulted in three Moho reflection bands, two of which are dipping SE from 12.5 to 15.0 s two-way time (TWT) and from 14.5 to 16.0 s TWT. Another band of subhorizontal Moho reflections, at ≈ 12.5 s TWT, may have developed during the crustal perturbations related to a post-Delhi tectonic orogeny. The signatures of the Proterozoic collision, in the form of strong SE-dipping reflections in the lower crust and Moho, have been preserved in the DFB, indicating that the crust here has not undergone any significant ductile deformation since at least after the Delhi rifting event.  相似文献   

11.
Through the investigation of crustal heterogeneities, sedimentary basin architecture and seismic stratigraphy, we demonstrate how a crust‐scale anisotropy controls the initiation of rifting and the subsequent structural and sedimentological evolution of the Mesozoic Gamtoos Basin, southern South Africa. The results demonstrate that the >90‐km‐long Gamtoos Fault established its length very early in its syn‐rift phase (within ~5 Ma of rift initiation) before accruing over 6 s (two‐way‐travel time (TWT)), or >12 km, of displacement without any significant subsequent increase in length. In addition, there is no evidence at the resolution of the data of fault segmentation, isolated depocentres nor of intra‐basin faults progressively coalescing during the syn‐rift interval. The early establishment of length resulted in a rapid transition from a terrestrial depositional environment to anoxic, deep marine conditions. The Gamtoos Fault has a 90° bend in the fault trace that we propose is inherited from the underlying structure. Immediately adjacent to the bend the basin‐fill is significantly deformed and a high‐amplitude (>1.7s TWT) monoclinal fold is observed. Previous workers proposed that the fold was a consequence of a complex interplay between compression and extension. Through a restoration of the basin‐fill deformation we produce a model that suggests that the fold is a consequence of the accommodation of extension by the unusual plan‐view trace of the fault. The evolution of the basin does not conform to current fault growth models and it is proposed that its unusual and complex development can be attributed to the underlying crustal‐scale anisotropy, a fact that is likely to be important in other areas in which crustal stretching is superimposed on heterogeneous continental crust.  相似文献   

12.
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.  相似文献   

13.
We have analysed the fundamental mode of Love and Rayleigh waves generated by 12 earthquakes located in the mid-Atlantic ridge and Jan Mayen fracture zone. Using the multiple filter analysis technique, we isolated the Rayleigh and Love wave group velocities for periods between 10 and 50  s. The surface wave propagation paths were divided into five groups, and average group velocities calculated for each group. The average group velocities were inverted and produced shear wave velocity models that correspond to a quasi-continental oceanic structure in the Greenland–Norwegian Sea region. Although resolution is poor at shallow depth, we obtained crustal thickness values of about 18  km in the Norwegian Sea area and 9  km in the region between Svalbard and Iceland. The abnormally thick crust in the Norwegian Sea area is ascribed to magmatic underplating and the thermal blanketing effect of sedimentary layers. Maximum crustal shear velocities vary between 3.5 and 3.9  km  s−1 for most paths. An average lithospheric thickness of 60  km was observed, which is lower than expected for oceanic-type structure of similar age. We also observed low shear wave velocities in the lower crust and upper mantle. We suggest that high heat flow extending to depths of about 30  km beneath the surface can account for the thin lithosphere and observed low velocities. Anisotropy coefficients of 1–5 per cent in the shallow layers and >7 per cent in the upper mantle point to the existence of polarization anisotropy in the region.  相似文献   

14.
Curie-temperature depth estimation using a self-similar magnetization model   总被引:4,自引:0,他引:4  
The Earth's crust is magnetized down to the Curie-temperature depth at about 10 to 50 km. This limited depth extent of the crustal magnetization is discernible in the power spectra of magnetic maps of South Africa and Central Asia. At short wavelengths, the power increases as rapidly towards longer wavelengths as expected for a self-similar magnetized crust with unlimited depth extent. Above wavelengths of about 100 km the power starts increasing less rapidly, indicating the absence of deep-seated sources. To quantify this effect we derive the theoretical power spectrum due to a slab carved out of a self-similar magnetization distribution. This model power spectrum matches the power spectra of South Africa and Central Asia for a self-similarity parameter of β = 4 and Curie temperature depths of 15 to 20 km.  相似文献   

15.
Interpretation of long‐offset 2D depth‐imaged seismic data suggests that outer continental margins collapse and tilt basinward rapidly as rifting yields to seafloor spreading and thermal subsidence of the margin. This collapse post‐dates rifting and stretching of the crust, but occurs roughly ten times faster than thermal subsidence of young oceanic crust, and thus is tectonic and pre‐dates the ‘drift stage’. We term this middle stage of margin development ‘outer margin collapse’, and it accords with the exhumation stage of other authors. Outer continental margins, already thinned by rifting processes, become hanging walls of crustal‐scale half grabens associated with landward‐dipping shear zones and zones of low‐shear strength magma at the base of the thinned crust. The footwalls of the shear zones comprise serpentinized sub‐continental mantle that commonly becomes exhumed from beneath the embrittled continental margin. At magma‐poor margins, outer continental margins collapse and tilt basinward to depths of about 3 km subsea at the continent–ocean transition, often deeper than the adjacent oceanic crust (accreted later between 2 and 3 km). We use the term ‘collapse’ because of the apparent rapidity of deepening (<3 Myr). Rapid salt deposition, clastic sedimentation (deltaic), or magmatism (magmatic margins) may accompany collapse, with salt thicknesses reaching 5 km and volcanic piles 1525 km. This mechanism of rapid salt deposition allows mega‐salt basins to be deposited on end‐rift unconformities at global sea level, as opposed to deep, air‐filled sub‐sea depressions. Outer marginal collapse is ‘post‐rift’ from the perspective of faulting in the continental crust, but of tectonic, not of thermal, origin. Although this appears to be a global process, the Gulf of Mexico is an excellent example because regional stratigraphic and structural relations indicate that the pre‐salt rift basin was filled to sea level by syn‐rift strata, which helps to calibrate the rate and magnitude of collapse. We examine the role of outer marginal detachments in the formation of East India, southern Brazil and the Gulf of Mexico, and how outer marginal collapse can migrate diachronously along strike, much like the onset of seafloor spreading. We suggest that backstripping estimates of lithospheric thinning (beta factor) at outer continental margins may be excessive because they probably attribute marginal collapse to thermal subsidence.  相似文献   

16.
Summary. The unified seismic exploration program, consisting of 345 km of deep reflection profiling, a 200 km refraction profile, an expanding spread profile and near-surface high resolution reflection meaasurements, revealed a strongly differentiated crust beneath the Black Forest. The highly reflective lower crust contains numerous horizontal and dipping reflectors at depths of 13-14 km down to the crust-mantle boundary (Moho). The Moho appears as a flat horizontal first order discontinuity at a relatively shallow level of 25–27 km above a transparent upper mantle. From modelling of synthetic near-vertical and wide-angle seismograms using the reflectivity method the lower crust is supposed to be composed of laminae with an average thickness of about 100 m and velocity differences of greater than 10% increasing from top to bottom. The upper crust is characterised by mostly dipping reflectors, associated with bivergent underthrusting and accretion tectonics of Variscan age and with extensional faults of Mesozoic age. A bright spot at 9.5 km depth is characterised by low velocity material suggesting a fluid trap. It appears on all of the three profiles in the centre of the intersection region. The upper crust seems to be decoupled from the lowest crust by a relatively transparent zone which is' also identified as a low-velocity zone. This low velocity channel is situated directly above the laminated lower crust. The laminae in the Rhinegraben area are displaced vertically to greater depths indicating an origin before Tertiary rift formation and a subsidence of the whole graben wedge.  相似文献   

17.
The deep seismicity of the Tyrrhenian Sea   总被引:4,自引:0,他引:4  
The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.  相似文献   

18.
Summary. An ocean bottom seismograph survey of the junction of the East Pacific Rise and the Wilkes fracture zone detected only three microearthquakes beneath the rise crest during seven days of recording. In contrast, during the same period 41 events were detected on the fracture zone, all at distances greater than 10 km from the junction. These results suggest that near the rise crest the thin crust can support sufficient stress only to generate infrequent small earthquakes and that most faulting may take place by aseismic slip. At 10 km from the rise axis part of the crest has become competent enough to support stress, resulting in earthquakes probably at depths of up to 5 km below the sea-bed. Gear 5-waves on the seismometer records indicate that a magma chamber, if it exists near the junction, is less than 10 km across.  相似文献   

19.
The Narmada zone in central India is a zone of weakness that separates the region of Vindhyan (Meso-Neoproterozoic) deposition to the north from Gondwana (Permo-Carboniferous–lower Cretaceous) deposits to the south. The reinterpretation of analogue seismic refraction data, acquired during the early 1980s, using 2-D ray-tracing techniques reveals a basement (velocity 5.8–6.0 km s−1 ) topography suggesting that the Narmada zone, bounded by the Narmada North and Narmada South faults is a region of basement uplift. A layer of anomalously high velocity (6.5–6.7 km s−1 ) at depths between 1.5 and 9.0 km appears to be present in the entire region. Within the Narmada zone this layer occurs at shallower depths than outside the Narmada zone. At two places within the Narmada zone this layer is at a depth of about 1.5 km. This layer cannot be considered as the top of the lower crust because in this case it should have produced large positive gravity anomalies at the shallowest parts. Instead, these parts correspond to Bouguer gravity lows. Furthermore, lower crust at such shallow depths has not been reported from any other part of the Indian shield. Therefore, this layer is likely to represent the top of a high-velocity mafic body that has different thicknesses in different places.  相似文献   

20.
The stratigraphic, subsidence and structural history of Orphan Basin, offshore the island of Newfoundland, Canada, is described from well data and tied to a regional seismic grid. This large (400 by 400 km) rifted basin is part of the non‐volcanic rifted margin in the northwest Atlantic Ocean, which had a long and complex rift history spanning Middle Jurassic to Aptian time. The basin is underlain by variably thinned continental crust, locally <10‐km thick. Our work highlights the complex structure, with major upper crustal faults terminating in the mid‐crust, while lower crustal reflectivity suggests ductile flow, perhaps accommodating depth‐dependent extension. We describe three major stratigraphic horizons connected to breakup and the early post‐rift. An Aptian–Albian unconformity appears to mark the end of crustal rifting in the basin, and a second, more subdued Santonian unconformity was also noted atop basement highs and along the proximal margins of the basin. Only minor thermal subsidence occurred between development of these two horizons. The main phase of post‐rift subsidence was delayed until post‐Santonian time, with rapid subsidence culminating in the development of a major flooding surface in base Tertiary time. Conventional models of rifting events predict significant basin thermal subsidence immediately following continental lithospheric breakup. In the Orphan Basin, however, this subsidence was delayed for about 25–30 Myr and requires more thinning of the mantle lithosphere than the crust. Models of the subsidence history suggest that extreme thinning of the lithospheric mantle continued well into the post‐rift period. This is consistent with edge‐driven, small‐scale convective flow in the mantle, which may thin the lithosphere from below. A hot spot may also have been present below the region in Aptian–Albian time.  相似文献   

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

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