首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
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.  相似文献   

2.
We use data from the Chile Argentina Geophysical Experiment (CHARGE) broad-band seismic deployment to refine past observations of the geometry and deformation within the subducting slab in the South American subduction zone between 30°S and 36°S. This region contains a zone of flat slab subduction where the subducting Nazca Plate flattens at a depth of ∼100 km and extends ∼300 km eastward before continuing its descent into the mantle. We use a grid-search multiple-event earthquake relocation technique to relocate 1098 events within the subducting slab and generate contours of the Wadati-Benioff zone. These contours reflect slab geometries from previous studies of intermediate-depth seismicity in this region with some small but important deviations. Our hypocentres indicate that the shallowest portion of the flat slab is associated with the inferred location of the subducting Juan Fernández Ridge at 31°S and that the slab deepens both to the south and the north of this region. We have also determined first motion focal mechanisms for ∼180 of the slab earthquakes. The subhorizontal T -axis solutions for these events are almost entirely consistent with a slab pull interpretation, especially when compared to our newly inferred slab geometry. Deviations of T -axes from the direction of slab dip may be explained with a gap within the subducting slab below 150 km in the vicinity of the transition from flat to normal subducting geometry around 33°S.  相似文献   

3.
Summary. We present the results of a systematic study of events with M s > 6 in northern Chile (20–33°S), for the period between 1963 and 1971. Medium to large earthquakes near the coast of this region are of three types: (1) Interplate events at the interface between the downgoing slab and the overriding South American plate. These events can be very large reaching magnitudes greater than 8. (2) Intra-plate earthquakes 20–30 km inside the downgoing slab. They have fault mechanisms indicating extension along the dip of the slab and may have magnitudes up to 7.5. (3) Less frequent, M s∼ 6 events that occur near the top of the downgoing slab and have thrust mechanisms with an almost horizontal E-W compressional axis. This type of mechanism is very different from that of the events of type 1 which are due to shallow dipping reverse faulting. There is a rotation of about 30° of the compressional axis in the vertical plane between events of types (1) and (3). Three groups of events near 32.5°, 25.5° and 21°s were studied in detail. Depth and mechanisms were redetermined by P -wave modelling and relative locations were obtained by a master event technique. Near 32.5°S, only events of types 1 and 2 were found in the time period of this study. At the two other sites, the three types of events were identified. This shows clearly that there are compressive stresses at the top of the slab and extension at the centre, a situation which is usually found in the areas where a double Benioff-zone has been identified in the seismicity.  相似文献   

4.
Summary. Over 80 earthquakes, exclusively from the Hindukush focal region, which were recorded at the Gauribidanur seismic array (GBA) have been used in this study. These events have similar epicentral distances and a narrow azimuthal range from GBA but varying focal depths from 10 to 240 km. A fault plane dipping steeply (75°) in the north-west direction and striking N 66° E has been investigated on the basis of the spatial distribution of earthquakes in two vertical planes through 68° E and 32° N. Short period P -wave recordings up to 30 s were processed using the adaptive cross-correlation filtering technique. Slowness and azimuthal anomalies were obtained for first arrivals. These anomalies show positive as well as negative bias and are attributed to a steep velocity gradient in the upper mantle between the 400–700 km depth range where the seismic rays have their maximum penetration. Relative time residuals between the stations of GBA owe their origin very near to the surface beneath the array. A search of the signals across the array revealed that most of the events occurring at shallower depths had complex signatures as compared to the deeper events. The structure near the source region, complicated source functions and the scattering confined to the crust—upper mantle near source are mainly responsible for the complexity of the Hindukush earthquakes as the transmission zone of the ray tubes from turning point to the recording station is practically the same.  相似文献   

5.
Summary. Bulletins of the International Seismological Centre (ISC) show very large residuals, up to 15 s early, for arrivals from events in the Tonga–Kermadec subduction zone to the New Zealand network of seismometers. The very early arrivals are confined to events south of about 22°S, and shallower than about 350 km. The waveforms show two distinct phases: an early, emergent, first phase with energy in the high-frequency band 2–10 Hz, and a distinct second phase, containing lower frequency energy, arriving at about the time predicted by JB tables.
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects.  相似文献   

6.
The North Canterbury region marks the transition from Pacific plate subduction to continental collision in the South Island of New Zealand. Details of the seismicity, structure and tectonics of this region have been revealed by an 11-week microearthquake survey using 24 portable digital seismographs. Arrival time data from a well-recorded subset of microearthquakes have been combined with those from three explosions at the corners of the microearthquake network in a simultaneous inversion for both hypocentres and velocity structure. The velocity structure is consistent with the crust in North Canterbury being an extension of the converging Chatham Rise. The crust is about 27 km thick, and consists of an 11 km thick seismic upper crust and 7 km thick seismic lower crust, with the middle part of the crust being relatively aseismic. Seismic velocities are consistent with the upper and middle crust being composed of greywacke and schist respectively, while several lines of evidence suggest that the lower crust is the lower part of the old oceanic crust on which the overlying rocks were originally deposited.
The distribution of relocated earthquakes deeper than 15 km indicates that the seismic lower crust changes dip markedly near 43S. To the south-west it is subhorizontal, while to the north-east it dips north-west at about 10. Fault-plane solutions for these earthquakes also change near 43S. For events to the south, P -axes trend approximately normal to the plate boundary (reflecting continental collision), while for events to the north, T -axes are aligned down the dip of the subducted plate (reflecting slab pull). While lithospheric subduction is continuous across the transition, it is not clear whether the lower crust near 43S is flexed or torn.  相似文献   

7.
An analysis of the Zihuatanejo, Mexico, earthquake of 1994 December 10 ( M = 6.6), based on teleseismic and near-source data, shows that it was a normal-faulting, intermediate-depth ( H = 50 ± 5 km) event. It was located about 30 km inland, within the subducted Cocos plate. The preferred fault plane has an azimuth of 130°, a dip of 79° and a rake of −86°. The rupture consisted of two subevents which were separated in time by about 2 s, with the second subevent occurring downdip of the first. The measured stress drop was relatively high, requiring a Δσ of about a kilobar to explain the high-frequency level of the near-source spectra. A rough estimate of the thickness of the seismogenic part of the oceanic lithosphere below Zihuatanejo, based on the depth and the rupture extent of this event, is 40 km.
This event and the Oaxaca earthquake of 1931 January 15 ( M = 7.8) are the two significant normal-faulting, intermediate-depth shocks whose epicentres are closest to the coast. Both of these earthquakes were preceded by several large to great shallow, low-angle thrust earthquakes, occurring updip. The observations in other subduction zones show just the opposite: normal-faulting events precede, not succeed, updip, thrust shocks. Indeed, the thrust events, soon after their occurrence, are expected to cause compression in the slab, thus inhibiting the occurrence of normal-faulting events. To explain the occurrence of the Zihuatanejo earthquake, we note that the Cocos plate, after an initial shallow-angle subduction, unbends and becomes subhorizontal. In the region of the unbending, the bottom of the slab is in horizontal extension. We speculate that the large updip seismic slip during shallow, low-angle thrust events increases the buckling of the slab, resulting in an incremental tensional stress at the bottom of the slab and causing normal-faulting earthquakes. This explanation may also hold for the 1931 Oaxaca event.  相似文献   

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

9.
Summary. Multichannel seismic reflection sections recorded across Vancouver Island have revealed two extensive zones of deep seismic reflections that dip gently to the northeast, and a number of moderate northeasterly dipping reflections that can be traced to the surface where major faults are exposed. Based on an integrated interpretation of these data with information from gravity, heat flow, seismicity, seismic refraction, magnetotelluric and geological studies it is concluded that the lower zone of gently dipping reflections is due to underplated oceanic sediments and igneous rocks associated with the current subduction of the Juan de Fuca plate, and that the upper zone represents a similar sequence of accreted rocks associated with an earlier episode of subduction. The high density/high velocity material between the two reflection zones is either an underplated slab of oceanic lithosphere or an imbricated package of mafic rocks. Reprocessing of data from two of the seismic lines has produced a remarkable image of the terrane bounding Leech River fault, with its dip undulating from >60° near the surface to 20° at 3 km depth and ∼38° at 6 km depth.  相似文献   

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

11.
Summary. The tectonics of the region at the convergence of the Tyrrhenian and Ionian basins is very complicated. The seismicity is particularly interesting in the area 38°–39° 30'N, 14° 30'–16° 30'E, where both deep events occur with compressive mechanism and crustal events with predominant normal faulting in excess of magnitude 7.0. The tectonic origin of this apparently contradictory behaviour can be traced to a remanent downgoing slab in the southern Tyrrhenian basin which is now inactive but is still undergoing considerable horizontal bending. Sedimentological studies suggest the bending to be accompanied by a marked uplift of Calabira reaching its highest rates in the area of the Strait of Messina. We performed a spirit levelling campaign in the Strait and found that a differential uplift of the inner lands with respect to coastlines took place in the period 1970–82. he question regarding the possibilities of either an absolute uplift of inner lands or a subsidence of coastlines can only be resolved after detailed tidal studies; the available tidal data seem to support the uplift hypothesis but they are substantially incomplete. The apparent pattern of vertical deformations seems to be a marked steady uplift with episodic subsidence periods in connection with large normal faulting earthquakes. Furthermore, our measurements and the results of a recently completed trilateration campaign tend to confirm the hypothesis of the existence of an angular graben in the Strait and a south–north right lateral slip of Sicily with respect to Calabria.  相似文献   

12.
Summary. Travel times of seismic waves from teleseismically recorded events that occurred within the Adak Island region of the Central Aleutian Islands are used to constrain both the velocity structure of the region and the locations of events. It is found that the P-wave velocity within the descending lithosphere is approximately 7–8 per cent faster than that in the surrounding mantle for a slab length of 300–400 km.
The velocity structure so determined is then used as a model for the relocation of events recorded only locally. It is found that features of the Benioff zone seen in locations made with a spherically symmetric model can be seriously in error. In particular, an apparent increase in dip of the zone with depth disappears when the locations are made with the inhomogeneous model. The new locally determined hypocentres are now also more consistent with the locations based on teleseismic data.  相似文献   

13.
The migration of teleseismic receiver functions yields high-resolution images of the crustal structure of western Crete. Data were collected during two field campaigns in 1996 and 1997 by networks of six and 47 short-period three-component seismic stations, respectively. A total of 1288 seismograms from 97 teleseismic events were restituted to true ground displacement within a period range from 0.5 to 7 s. The application of a noise-adaptive deconvolution filter and a new polarization analysis technique helped to overcome problems with local coda and noise conditions. The computation and migration of receiver functions results in images of local crustal structures with unprecedented spatial resolution for this region. The crust under Crete consists of a continental top layer of 15–20 km thickness above a 20–30 km thick subducted fossil accretionary wedge with a characteristic en echelon fault sequence. The downgoing oceanic Moho lies at a depth of 40–60 km and shows a topography or undulation with an amplitude of several kilometres. As a consequence of slab depth and distribution of local seismicity, the Mediterranean Ridge is interpreted as the recent accretionary wedge.  相似文献   

14.
Data recorded by the Italian Telemetered Seismic Network (ITSN) of the Istituto Nazionale di Geofisica (ING) have been widely used in recent years to image slab structures and to find evidence for active processes along the Italian Peninsula. However, the use of seismic data for geostructural purposes may be affected by the well-known trade-off between earthquake location and seismic-velocity parameters. Furthermore, the confidence ellipse predicted by standard procedures may be inadequate for the representation of the probable error of a computed localization. This paper evaluates the probable errors on the hypocentre determinations of the seismic events recorded by the ITSN, using a Monte Carlo method.
We compute synthetic arrival times using a 1-D velocity model appropriate as an average for the Italian area. The hypocentres used are all those recorded by the ITSN during the period January 1992 to March 1994 (1972 events). Station locations are those of the current ITSN configuration. The synthetic arrival times are perturbed with a Gaussian distribution of errors and input to ING's standard hypocentral location procedure, but using crustal velocities differing by 10 per cent from those used to generate them. Each simulation is repeated at least 30 times. Average absolute shifts of hypocentres are assessed in grid cells of linear dimension 33 km covering the whole Italian region.
For regions within the ITSN, shifts are typically 5–10 km in location and up to 20 km in depth. However, for offshore and coastal regions, they are much greater: 50 km or more in both location and depth (far exceeding the equivalent uncertainties quoted by ING bulletins). Possible consequences of this are highlighted by producing a cross-section of subcrustal hypocentres from the Adriatic to the Tyrrhenian Sea, where the large uncertainty in depth precludes any confident interpretation of dipping tectonic features.  相似文献   

15.
Signature of remnant slabs in the North Pacific from P-wave tomography   总被引：1，自引：0，他引：1  
A 3-D ray-tracing technique was used in a global tomographic inversion in order to obtain tomographic images of the North Pacific. The data reported by the Geophysical Survey of Russia (1955–1997) were used together with the catalogues of the International Seismological Center (1964–1991) and the US Geological Survey National Earthquake Information Center (1991–1998), and the recompiled catalogue was reprocessed. The final data set, used for following the inversion, contained 523 430 summary ray paths. The whole of the Earth's mantle was parametrized by cells of 2° × 2° and 19 layers. The large and sparse system of observation equations was solved using an iterative LSQR algorithm.
A subhorizontal high-velocity anomaly is revealed just above the 660 km discontinuity beneath the Aleutian subduction zone. This high-velocity feature is observed at latitudes of up to ~70°N and is interpreted as a remnant of the subducted Kula plate, which disappeared through ridge subduction at about 48 Ma. A further positive velocity perturbation feature can be identified beneath the Chukotka peninsula and Okhotsk Sea, extending from ~300 to ~660 km depth and then either extending further down to ~800 km (Chukotka) or deflecting along the 660 km discontinuity (Okhotsk Sea). This high-velocity anomaly is interpreted as a remnant slab of the Okhotsk plate accreted to Siberia at ~55 Ma.  相似文献   

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

17.
The highest intermediate depth moment release rates in Indonesia occur in the slab beneath the largely submerged segment of the Banda arc in the Banda Sea to the east of Roma, termed the Damar Zone. The most active, western-part of this zone is characterized by downdip extension, with moment release rates (∼1018 Nm yr–1 per 50 km strike length) implying the slab is stretching at ∼10−14 s−1 consistent with near complete slab decoupling across the 100–200 km depth range. Differential vertical stretching along the length of the Damar Zone is consistent with a slab rupture front at ∼100–200 km depth beneath Roma propagating eastwards at ∼100 km Myr–1. Complexities in the slab deformation field are revealed by a narrow zone of anomalous in-plane P -axis trends beneath Damar, where subhorizontal constriction suggests extreme stress concentrations ∼100 km ahead of the slab rupture front. Such stress concentrations may explain the anomalously deep ocean gateways in this region, in which case ongoing slab rupture may have played a key role in modulating the Indonesian throughflow in the Banda Sea over the last few million years.  相似文献   

18.
Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps.  相似文献   

19.
Data from apatite fission track analysis are presented for 20 outcrop samples collected in the southern Adelaide Fold Belt, South Australia. Interpretation of these data, with the aid of numerical models which allow inference of multiphase cooling histories, indicate three discrete cooling events that are likely to correlate with sedimentation events in surrounding depositional settings. An event beginning some time after 85 Ma (Late Cretaceous) was characterized by cooling throughout the study area from temperatures of roughly 50 to 70 °C. An event beginning at 300–270 Ma (Late Palaeozoic) was characterized by cooling from temperatures >120 °C in all areas except for the Mount Lofty Ranges and Murray Bridge region, where peak temperatures were only 95–115 °C prior to Palaeozoic cooling. Some samples from these subregions of relatively cool Late Palaeozoic temperatures also retain evidence for even earlier cooling from temperatures >120 °C, beginning prior to 350 Ma. We interpret the post 85-Ma event as the consequence of regional exhumation from a depth of 1.0–1.6 km. The Late Palaeozoic event (300–270 Ma) is interpreted as cooling associated with the termination of the Alice Springs Orogeny, while cooling prior to 350 Ma probably represents the final stages of Early Middle Palaeozoic unroofing of the southern Adelaide Fold Belt.
The results highlight the importance of regional, episodic postorogenic exhumation of Palaeozoic fold belts, where – in some cases – conventional methods have erroneously suggested relatively long-term stability.  相似文献   

20.
Seismic anisotropy within the uppermost mantle of southern Germany   总被引：1，自引：0，他引：1  
This paper presents an updated interpretation of seismic anisotropy within the uppermost mantle of southern Germany. The dense network of reversed and crossing refraction profiles in this area made it possible to observe almost 900 traveltimes of the Pn phase that could be effectively used in a time-term analysis to determine horizontal velocity distribution immediately below the Moho. For 12 crossing profiles, amplitude ratios of the Pn phase compared to the dominant crustal phase were utilized to resolve azimuthally dependent velocity gradients with depth. A P -wave anisotropy of 3–4 per cent in a horizontal plane immediately below the Moho at a depth of 30 km, increasing to 11 per cent at a depth of 40 km, was determined. For the axis of the highest velocity of about 8.03 km s−1 at a depth of 30 km a direction of N31°F was obtained. The azimuthal dependence of the observed Pn amplitude is explained by an azimuth-dependent sub-Moho velocity gradient decreasing from 0.06 s−1 in the fast direction to 0 s−1 in the slow direction of horizontal P -wave velocity. From the seismic results in this study a petrological model suggesting a change of modal composition and percentage of oriented olivine with depth was derived.  相似文献   

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

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