首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V P values measured for the soils. The downtown part of Baku city is characterized by low VS30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.  相似文献   

2.
Through combining the soil respiration with the main environmental factors under the planting shelterbelt (Populus woodland) and the natural desert vegetation (Tamarix ramosissima+Phragmites communis community and Haloxylon ammodendron community) in the western Junngar Basin, the difference in soil respiration under different land use/land cover types and the responses of soil respiration to temperature and soil moisture were analyzed. Results showed that the rate of soil respiration increased with temperature. During the daytime, the maximum soil respiration rate occurred at 18:00 for the Populus woodland, 12:00 for T. ramosissima+Ph. communis community, and 14:00 for H. ammodendron community, while the minimum rate all occurred at 8:00. The soil respiration, with the maximum rate in June and July and then declining from August, exhibited a similar trend to the near-surface temperature from May to October. During the growing season, the mean soil respiration rates and seasonal variation differed among the land use/land cover types, and followed the order of Populus woodland > T. ramosissima+Ph. communis community > H. ammodendron community. The difference in the soil respiration rate among different land use/land cover types was significant. The soil respiration of Pouplus woodland was significantly correlated with the near-surface temperature and soil temperature at 10 cm depth (P<0.01) in an exponential manner. The soil respiration of T. ramosissima+Ph. communis and H. ammodendron communities were all linearly correlated with the near-surface temperature and soil surface temperature (P<0.01). Based on the near-surface temperature, the calculated Q10 of Populus woodland, T. ramosissima+Ph. communis community and H. ammodendron community were 1.48, 1.59 and 1.63, respectively. The integrated soil respiration of the three land use/land cover types showed a significant correlation with the soil moisture at 0–5 cm, 5–15 cm and 0–15 cm depths (P<0.01). The quadratic model could best describe the relationship between soil respiration and soil moisture at 0–5 cm depth (P<0.01).  相似文献   

3.
Many studies on global climate have forecast major changes in the amounts and spatial patterns of precipitation that may significantly affect temperate grasslands in arid and semi-arid regions. As a part of ChinaFLUX, eddy covariance flux measurements were made at a semi-arid Leymus chinensis steppe in Inner Mongolia, China during 2003–2004 to quantify the response of carbon exchange to environmental changes. Results showed that gross ecosystem production (F GEP) and ecosystem respiration (R eco) of the steppe were significantly depressed by water stress due to lack of precipitation during the growing season. Temperature was the dominant factor affecting F GEP and R eco in 2003, whereas soil moisture imposed a significant influence on both R eco and F GEP in 2004. Under wet conditions, R eco showed an exponentially increasing trend with temperature (Q 10 = 2.0), but an apparent reduction in the value of R eco and its temperature sensitivity were observed during the periods of water stress (Q 10=1.6). Both heat and water stress can cause decrease in F GEP. The seasonality of ecosystem carbon exchange was strongly correlated with the variation of precipitation. With less precipitation in 2003, the steppe sequestrated carbon in June and July, and went into a senescence in early August due to water stress. As compared to 2003, the severe drought during the spring of 2004 delayed the growth of the steppe until late June, and the steppe became a CO2 sink from early July until mid-September, with ample precipitation in August. The semi-arid steppe released a total of 9.7 g C·m?2 from May 16 to the end of September 2003, whereas the net carbon budget during the same period in 2004 was close to zero. Long-term measurements over various grasslands are needed to quantify carbon balance in temperate grasslands.  相似文献   

4.
The study of the Gutenberg-Richter (GR) parameters a and b has been very important to describe and characterize the seismicity over the different seismic provinces around the world. As far as we know, the possible correlation between the GR parameters a and b has not received enough attention. Bayrak et al. reported the a and b values for 27 active seismic regions around the boundaries of the main tectonic plates of the world. From these data, we found that there exists a positive correlation between the a and b parameters (R =?0.85, R2 =?0.72). On the other hand, we made around 150 computer runs of a spring-block model proposed by Olami et al. (Phys Rev Lett 68(8):1244–1247, 1992). This model roughly emulates the interaction between two fault planes and it reaches a self-organized critical state. With these simulations, we also found that the a and b parameters are positively correlated. Motivated by these results, we propose an analytical demonstration that indeed a and b are positively correlated. In addition, we discuss on other possible applications of the spring-block model to actual seismicity and to frictional experiments made with sandpapers.  相似文献   

5.
Line-of-sight magnetograms acquired by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO) and by the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO) for 14 emerging ARs were used to study the derivative of the total unsigned flux–the flux emergence rate, R(t). We found that the emergence regime is not universal: each AR displays a unique emergence process. Nevertheless, two types of the emergence process can be identified. First type is a “regular” emergence with quasi-constant behavior of R(t) during a 1–3 day emergence interval with a rather low magnitude of the flux derivative, Rmax = (0.57 ± 0.22) × 1022 Mx day–1. The second type can be described as “accelerated” emergence with a long interval (>1 day) of the rapidly increasing flux derivative R(t) that result in a rather high magnitude of Rmax= (0.92 ± 0.29) × 1022 Mx day–1, which later changes to a very short (about a one third of day) interval of R(t) = const followed by a monotonous decrease of R(t). The first type events might be associated with emergence of a flux tube with a constant amount of flux that rises through the photosphere with a quasi-constant speed. Such events can be explained by the traditional largescale solar dynamo generating the toroidal flux deep in the convective zone. The second-type events can be interpreted as a signature of sub-surface turbulent dynamo action that generates additional magnetic flux (via turbulent motions) as the magnetic structure makes its way up to the solar surface.  相似文献   

6.
This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.  相似文献   

7.
The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.  相似文献   

8.
Seismic intensity measures (IMs) perform a pivotal role in probabilistic seismic demand modeling. Many studies investigated appropriate IMs for structures without considering soil liquefaction potential. In particular, optimal IMs for probabilistic seismic demand modeling of bridges in liquefied and laterally spreading ground are not comprehensively studied. In this paper, a coupled-bridge-soil-foundation model is adopted to perform an in-depth investigation of optimal IMs among 26 IMs found in the literature. Uncertainties in structural and geotechnical material properties and geometric parameters of bridges are considered in the model to produce comprehensive scenarios. Metrics such as efficiency, practicality, proficiency, sufficiency and hazard computability are assessed for different demand parameters. Moreover, an information theory based approach is adopted to evaluate the relative sufficiency among the studied IMs. Results indicate the superiority of velocity-related IMs compared to acceleration, displacement and time-related ones. In particular, Housner spectrum intensity (HI), spectral acceleration at 2.0 s (S a-20), peak ground velocity (PGV), cumulative absolute velocity (CAV) and its modified version (CAV 5) are the optimal IMs. Conversely, Arias intensity (I a ) and shaking intensity rate (SIR) which are measures often used in liquefaction evaluation or related structural demand assessment demonstrate very low correlations with the demand parameters. Besides, the geometric parameters do not evidently affect the choice of optimal IMs. In addition, the information theory based sufficiency ranking of IMs shows an identical result to that with the correlation measure based on coefficient of determination (R 2). This means that R 2 can be used to preliminarily assess the relative sufficiency of IMs.  相似文献   

9.
In the Tarim Basin, black shale series at the bottom of Cambrian is one of the important marine facies hydrocarbon source rocks. This research focuses on the analysis of the isotope of noble gas of 11 cherts. The R/R a ratio of chert in the Keping area is 0.032–0.319, and 40Ar/36Ar is 338–430. In Quruqtagh the R/R a ratio is 0.44–10.21, and 40Ar/36Ar is 360–765. The R/R a ratio of chert increases with 40Ar/36Ar from the west to the east accordingly. They have evolved from the crust source area to the mantle source area in a direct proportion. Surplus argon 40ArE in chert is in direct proportion to the R/R a ratio, indicating that it has the same origin of excess argon as in fluid and mantle source helium. Comparison of the R/R a ratios between the west and the east shows that the chert in the eastern part formed from the activity system of the bottom hydrothermal venting driven by the mantle source, where the material and energy of crust and mantle had a strong interaction in exchange; whereas in the western part, chert deposited from the floating of hydrothermal plume undersea bottom, which is far away from the centre of activities of the hydrothermal fluid of ocean bottom. In addition, from noble gas isotope composition of chert, it is suggested that the ocean anoxia incident happened at the black shale of the Cambrian bottom probably because of the large-scaled ocean volcanoes and the following hydrothermal activities.  相似文献   

10.
The position of the auroral oval poleward and equatorward boundary projections on the equatorial plane in the nightside MLT sector during magnetically quiet periods (|AL| < 200 nT, |Dst| < 10 nT) has been determined. The oval boundary positions were determined according to the precipitation model developed at Polar Geophysical Institute (http://apm.pgia.ru/). The isotropy of the averaged plasma pressure and the experimentally confirmed balance of pressures during the nighttime have been taken into account. The morphological mapping method has been used to map the oval poleward and equatorward edges without the use of any magnetic field model on the assumption that the condition of magnetostatic equilibrium is valid. Ion pressures at ionospheric altitudes and in the equatorial plane have been compared. It has been shown that the auroral oval equatorward boundary in the midnight sector is localized at geocentric distances of ~7 RE, which is in good agreement with the position of the energetic particle injection boundary in the equatorial plane. The oval poleward edge is localized at the ~10 RE geocentric distance, which is in good agreement with the position of the equatorward boundary of the region with a high turbulence level in the Earth’s magnetosphere plasma sheet.  相似文献   

11.
Motivated by the scientific desire to align observations with quantities of physical interest, we survey how scalar importance functions depend on vertically resolved water vapor. Definitions of importance begin from familiar examples of water mass I m and TOA clear-sky outgoing longwave flux I OLR, in order to establish notation and illustrate graphically how the sensitivity profile or “kernel” depends on whether specific humidity S, relative humidity R, or ln(R) are used as measures of vapor. Then, new results on the sensitivity of convective activity I con to vapor (with implied knock-on effects such as weather prediction skill) are presented. In radiative-convective equilibrium, organized (line-like) convection is much more sensitive to moisture than scattered isotropic convection, but it exists in a drier mean state. The lesson for natural convection may be that organized convection is less susceptible to dryness and can survive and propagate into regions unfavorable for disorganized convection. This counterintuitive interpretive conclusion, with respect to the narrow numerical result behind it, highlights the importance of clarity about what is held constant at what values in sensitivity or susceptibility kernels. Finally, the sensitivities of observable radiance signals I sig for passive remote sensing are considered. While the accuracy of R in the lower free troposphere is crucial for the physical importance scalars, this layer is unfortunately the most difficult to isolate with passive remote sensing: In high emissivity channels, water vapor signals come from too high in the atmosphere (for satellites) or too low (for surface radiometers), while low emissivity channels have poor altitude discrimination and (in the case of satellites) are contaminated by surface emissions. For these reasons, active ranging (LiDAR) is the preferred observing strategy.  相似文献   

12.
The dynamics of the system of field-aligned currents (FACs) and closing ionospheric Pedersen currents is considered with the use of original processing methods and the data from four substorms of different types. The total current system comprises of two parts. One is the well-known substorm current wedge (SCW) system, in which the zonal (westward ) current closes FACs in the R1 zone (region). The component 2 consists of two pairs of meridional currents flowing equatorward and poleward in the R1 region and creating regions R0 and R2 (according to the classification of Iijima and Potemra). It is shown that the total FAC of the disturbed magnetosphere–ionosphere system is dominated by the contribution of component 2, which contradicts the original version of the SCW model but is consistent with new data. The quantitative characteristics of the dawn–dusk asymmetry are determined for the FAC distribution in the ionosphere for each substorm. It is shown that the ratio of the average intensities of FACs in the regions R0 and R2 was IR0/IR2 ≥ 0.4, which contradicts the popular opinion that there are no FACs in the polar cap. In addition, a relatively rare event of the simultaneous start of the substorm explosive phase and the SSC caused by the dynamic impact of the solar wind when the polar cap expands (rather than compresses as usual) is considered.  相似文献   

13.
A technique to estimate mass erosion rate of surface soil during landing of the Apollo Lunar Module (LM) and total mass ejected due to the rocket plume interaction is proposed and tested. The erosion rate is proportional to the product of the second moment of the lofted particle size distribution N(D), and third moment of the normalized soil size distribution S(D), divided by the integral of S(D)?D2/v(D), where D is particle diameter and v(D) is the vertical component of particle velocity. The second moment of N(D) is estimated by optical extinction analysis of the Apollo cockpit video. Because of the similarity between mass erosion rate of soil as measured by optical extinction and rainfall rate as measured by radar reflectivity, traditional NWS radar/rainfall correlation methodology can be applied to the lunar soil case where various S(D) models are assumed corresponding to specific lunar sites.  相似文献   

14.
The consideration of the relation between the daytime and nighttime values of the critical frequency F2, foF2 of the ionospheric F2 layer, started in the previous publication of the authors, is continued. The main regularities in variations in the correlation coefficient R(foF2) characterizing this relation are confirmed using larger statistical material (more ionospheric stations and longer observational series). Long-term trends in the R(foF2) value are found: at all stations the negative value of R(foF2) increases with time after 1980.  相似文献   

15.
Attenuation of P and S waves has been investigated in Alborz and north central part of Iran using the data recorded by two permanent and one temporary networks during October 20, 2009, to December 22, 2010. The dataset consists of 14,000 waveforms from 380 local earthquakes (2 < M L < 5.6). The extended coda normalization method (CNM) was used to estimate quality factor of P (Q P) and S waves (Q S) at seven frequency bands (0.375, 0.75, 1.5, 3, 6, 12, 24 Hz). The Q P and Q S values have been estimated at lapse times from 40 to 100 s. It has been observed that the estimated values of Q P and Q S are time independent; therefore, the mean values of Q P and Q S at different lapse times have been considered. The frequency dependence of quality factor was determined by using a power-law relationship. The frequency-dependent relationship for Q P was estimated in the form of (62 ± 7)f (1.03 ± 0.07) and (48 ± 5)f (0.95 ± 0.07) in Alborz region and North Central Iran, respectively. These relations for Q S for Alborz region and North Central Iran have estimated as (83 ± 8)f (0.99 ± 0.07) and (68 ± 5)f (0.96 ± 0.05), respectively. The observed low Q values could be the results of thermoelastic effects and/or existing fracture. The estimated frequency-dependent relationships are comparable with tectonically active regions.  相似文献   

16.
Stable isotope paleoaltimetry has provided unprecedented insights into the topographic histories of many of the world’s highest mountain ranges. However, on the Tibetan Plateau (TP), stable isotopes from paleosols generally yield much higher paleoaltitudes than those based on fossils. It is therefore essential when attempting to interpret accurately this region’s paleoaltitudes that the empirical calibrations of local stable isotopes and the relations between them are established. Additionally, it is vital that careful estimations be made when estimate how different isotopes sourced from different areas may have been influenced by different controls. We present here 29 hydrogen isotopic values for leaf wax-derived n-alkanes (i.e., δDwax values, and abundance-weighted average δD values of C29 and C31) in surface soils, as well as the δD values of soil water (δDsw) samples (totaling 22) from Mount Longmen (LM), on the eastern TP (altitude ~0.8–4.0 km above sea level (asl), a region climatically affected by the East Asian Monsoon (EAM). We compared our results with published data from Mount Gongga (GG). In addition, 47 river water samples, 55 spring water samples, and the daily and monthly summer precipitation records (from May to October, 2015) from two precipitation observation stations were collected along the GG transect for δD analysis. LM soil δDwax values showed regional differences and responded strongly to altitude, varying from?160‰ to?219‰, with an altitudinal lapse rate (ALR) of?18‰ km?1 (R 2=0.83; p<0.0001; n=29). These δDwax values appeared more enriched than those from the GG transect by ~40‰. We found that both the climate and moisture sources led to the differences observed in soil δDwax values between the LM and GG transects. We found that, as a general rule, ε wax/rw, ε wax/p and ε wax/sw values (i.e., the isotopic fractionation of δDwax corresponding to δDrw, δDp and δDsw) increased with increasing altitude along both the LM and GG transects (up to 34‰and 50‰, respectively). Basing its research on a comparative study of δDwax, δDp, δDrw(δDspringw) and δDsw, this paper discusses the effects of moisture recycling, glacier-fed meltwater, relative humidity (RH), evapotranspiration (ET), vegetation cover, latitude, topography and/or other factors on ε wax/p values. Clearly, if ε wax-p values at higher altitudes are calculated using smaller ε wax-p values from lower altitudes, the calculated paleowaterδDp values are going to be more depleted than the actual δD values, and any paleoaltitude would therefore be overestimated.  相似文献   

17.
Long-term measurement of carbon metabolism of old-growth forests is critical to predict their behaviors and to reduce the uncertainties of carbon accounting under changing climate. Eddy covariance technology was applied to investigate the long-term carbon exchange over a 200 year-old Chinese broad-leaved Korean pine mixed forest in the Changbai Mountains (128°28′E and 42°24′N, Jilin Province, P. R. China) since August 2002. On the data obtained with open-path eddy covariance system and CO2 profile measurement system from Jan. 2003 to Dec. 2004, this paper reports (i) annual and seasonal variation of F NEE, F GPP and R E; (ii) regulation of environmental factors on phase and amplitude of ecosystem CO2 uptake and release Corrections due to storage and friction velocity were applied to the eddy carbon flux.LAI and soil temperature determined the seasonal and annual dynamics of FGPP and RE separately. VPD and air temperature regulated ecosystem photosynthesis at finer scales in growing seasons. Water condition at the root zone exerted a significant influence on ecosystem maintenance carbon metabolism of this forest in winter.The forest was a net sink of atmospheric CO2 and sequestered ?449 g C·m?2 during the study period; ?278 and ?171 gC·m?2 for 2003 and 2004 respectively. F GPP and F RE over 2003 and 2004 were ?1332, ?1294 g C·m?2. and 1054, 1124 g C·m?2 respectively. This study shows that old-growth forest can be a strong net carbon sink of atmospheric CO2.There was significant seasonal and annual variation in carbon metabolism. In winter, there was weak photosynthesis while the ecosystem emitted CO2. Carbon exchanges were active in spring and fall but contributed little to carbon sequestration on an annual scale. The summer is the most significant season as far as ecosystem carbon balance is concerned. The 90 days of summer contributed 66.9, 68.9% of F GPP, and 60.4, 62.1% of R E of the entire year.  相似文献   

18.
Self-centering buckling-restrained braces (SCBRBs) were proposed recently to minimize residual deformation of the braces induced by yielding or buckling. Although earthquake resilience of structures equipped with the SCBRBs can be well achieved using displacement based designs (DBDs), previously proposed DBD procedures generally involve iterations. In this study, a novel direct displacement-based design method with a non-iterative procedure, named RCR DDBD, is proposed and applied to design of steel braced frame structures with SCBRBs. Unlike previously adopted DBD, the yield displacement does not need to be assumed initially in the proposed procedure. Instead, the yield strength and yield displacement are determined directly by the predetermined objective drift (ratio), using the relation of the strength reduction factor (R) and constant-strength inelastic displacement ratio spectra (CR spectra), i.e. the RCR relation. Since the derived RCR relation is independent with the peak ground acceleration of the earthquake records when stiffness and strength degradation are not considered, the proposed procedure can be accurate for any seismic level. The RCR DDBD is supposed to begin with the knowledge of the seismic excitation level (according to the structure category, site classification and owner’s requirements) and the corresponding target drift; the end of the design is to obtain the cross sections of main frame members and all the bracing parameters. The result of two 7-story buildings designed according to the RCR DDBD procedure demonstrates that this procedure can be effective and fairly simple for practical seismic design.  相似文献   

19.
The Xushen gas field, located in the north of Songliao Basin, is a potential giant gas area for China in the future. Its proved reserves have exceeded 1000×108 m3 by the end of 2005. But, the origin of natural gases from the deep strata is still in debating. Epimetamorphic rocks as a potential gas source are widely spreading in the northern basement of Songliao Basin. According to pyrolysis experiments for these rocks in the semi-confined system, gas production and geochemistry of alkane gases are discussed in this paper. The Carboniferous-Permian epimetamorphic rocks were heated from 300°C to 550°C, with temperature interval of 50°C. The gas production was quantified and measured for chemical and carbon isotopic compositions. Results show that δ 13C1 is less than ?20‰, carbon isotope trend of alkane gas is δ 13C1<δ 13C2<δ 13C3 or δ 13C1<δ 13C2>δ 13C3, these features suggest that the gas would be coal-type gas at high-over maturity, not be inorganic gas with reversal trend of gaseous alkanes (δ 13C1>δ 13C2>δ 13C3). These characteristics of carbon isotopes are similar with the natural gas from the basin basement, but disagree with gas from the Xingcheng reservoir. Thus, the mixing gases from the pyrolysis gas with coal-typed gases at high-over maturity or oil-typed gases do not cause the reversal trend of carbon isotopes. The gas generation intensity for epimetamorphic rocks is 3.0×108–23.8×108 m3/km2, corresponding to R o from 2.0% to 3.5% for organic matter.  相似文献   

20.
Microstructure measurements were performed along two sections through the Halmahera Sea and the Ombai Strait and at a station in the deep Banda Sea. Contrasting dissipation rates (??) and vertical eddy diffusivities (K z ) were obtained with depth-averaged ranges of \(\sim [9 \times 10^{-10}-10^{-5}]\) W kg??1 and of \(\sim [1 \times 10^{-5}-2 \times 10^{-3}]\) m2 s??1, respectively. Similarly, turbulence intensity, \(I={\epsilon }/(\nu N^{2})\) with ν the kinematic viscosity and N the buoyancy frequency, was found to vary seven orders of magnitude with values up to \(10^{7}\). These large ranges of variations were correlated with the internal tide energy level, which highlights the contrast between regions close and far from internal tide generations. Finescale parameterizations of ?? induced by the breaking of weakly nonlinear internal waves were only relevant in regions located far from any generation area (“far field”), at the deep Banda Sea station. Closer to generation areas, at the “intermediate field” station of the Halmahera Sea, a modified formulation of MacKinnon and Gregg (2005) was validated for moderately turbulent regimes with 100 < I < 1000. Near generation areas marked by strong turbulent regimes such as “near field” stations within strait and passages, ?? is most adequately inferred from horizontal velocities provided that part of the inertial subrange is resolved, according to Kolmogorov scaling.  相似文献   

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

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