首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 390 毫秒
1.
The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO2 flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO2 flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO2 leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol?1·m?2·s?1) flux data during windy conditions (u* > 0.2 m·s?1), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C02 exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO2·μmol?1 photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO2·m?2·s?1. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m?2mon?1, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as ?43.2±29.6 gC·m?2·mon?1. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as ?563.0 and ?441.2 gC·m?2·a?1 respectively, accounting for about 32% of GPP.  相似文献   

2.
Particle hygroscopicity plays a key role in understanding the mechanisms of haze formation and particle optical properties. The present study developed a method for predicting the effective hygroscopic parameter k and the water content of PM_(2.5) on the basis of the k-K?hler theory and bulk chemical components of PM_(2.5). Our study demonstrated that the effective hygroscopic parameter can be estimated using the PM_(2.5) mass concentration, water-soluble ions, and total water-soluble carbon. By combining the estimated k and ambient relative humidity, the water content of PM_(2.5) can be further estimated. As an example, the k and water content of PM_(2.5) in Beijing were estimated utilizing the method proposed in this study. The annual average value of k of PM_(2.5) in Beijing was 0.25±0.09, the maximum k value 0.26±0.08 appeared in summer, and the seasonal variation is insignificant. The PM_(2.5) water content was determined by both the PM_(2.5) hygroscopicity and the ambient relative humidity(RH). The annual average mass ratio of water content and PM_(2.5) was 0.18±0.20, and the maximum value 0.31±0.25 appeared in summer. Based on the estimated water content of PM_(2.5) in Beijing, the relationship between the PM_(2.5) water content and RH was parameterized as: m(%)=0.03+(5.73×10~(-8)) ×RH~(3.72).This parametric formula helps to characterize the relationship between the PM_(2.5) mass concentration and atmospheric visibility.  相似文献   

3.
The structure and theory of a tiltmeter with a liquid used as an inertial body and damper and a float with a string converter fitted with reliable remote system of arrestment used as acceleration indicator are described. The device’s calculated error is about ±4 × 10−5 (±2 × 10−10) g in the range of ±4 (±2 × 10−5) g. In the arrestment state, a second sensitive system of a capacitive type is turned on and the device starts to register seismic accelerations.  相似文献   

4.
The structure and theory of a tiltmeter, in which liquid is used as an inertial body and a damper, and a float with a string converter equipped with a reliable remote system of arrestment is employed as an acceleration indicator are described. The calculated error of the device is ±4 × 10−5 (±2 × 10−10)g in the range of ±4(±2 × 10−5)g. A second sensitive system of a capacitive type is turned on in the arrestment state, and the device registers seismic accelerations.  相似文献   

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.
Seismic strain and b value are used to quantify seismic potential in the Zagros region (Iran). Small b values (0.69 and 0.69) are related to large seismic moment rates (9.96×1017 and 4.12×1017) in southern zones of the Zagros, indicating more frequent large earthquakes. Medium to large b values (0.72 and 0.92) are related to small seismic moment rates (2.94×1016 and 6.80×1016) in middle zones of the Zagros, indicating less frequent large earthquakes. Small b value (0.64) is related to medium seismic moment rate (1.38×1017) in middle to northern zone of the Zagros, indicating frequent large earthquakes. Large b value (0.87) is related to large seismic moment rate (2.29×1017) in northwestern zone, indicating more frequent large earthquakes. Recurrence intervals of large earthquakes (M > 6) are short in southern (10 and 14 years) and northwestern (13 years) zones, while the recurrence intervals are long in the middle (46 and 114 years) and middle to northern (25 years) zones.  相似文献   

7.
Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R 1.0 as the geometrical spreading function in the NMSZ, the Q P and Q S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q P?=?(115.80?±?1.36) f (0.495?±?0.129) and Q S?=?(161.34?±?1.73) f (0.613?±?0.067), respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q S/Q P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.  相似文献   

8.
In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q S (f)?=?(147?±?16)f 0.71?±?0.02 and (126?±?12)f 0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).  相似文献   

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

10.
To describe the cumulative distribution of the number of lithospheric plates over areas S, the dependence N(≥S) ～ S ?0.33 was proposed by Bird [2003]. Based on dimension considerations, the dependence N(≥S) ～ (?/S)1/3, where ? is the generation rate of kinetic energy of convection in the mantle estimated at 10?11 m2/s3, is proposed. The analogy of plate formation with developed hydrodynamic turbulence and other processes involving an energy input into the system and its dissipation is considered. Simple experiments on random partitioning of a surface into polygons gave their cumulative distributions over areas resembling those observed for lithospheric plates. This has led to the conclusion that the plate distribution pattern is characteristic of the random partitioning of surfaces.  相似文献   

11.
Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M L ?≤?4.4. Up to five seismograph stations were utilised to determine Q c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q 0 and α were determined for the attenuation relation Q c(f)?=?Q 0 f α . The result was Q 0?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t s???t 0) (time from origin time t 0 to the start of coda analysis window is 1.9 times the S-travel time, t s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q c was not calculated. The derived Q c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q S was Q 0?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q c result.  相似文献   

12.
Most-probable-number (MPN) dilution series were used to enumerate and isolate bacteria from bulk water, suspended aggregates, the oxic layer, and the oxic–anoxic transition zone of the sediment of a tidal flat ecosystem in the southern North Sea. The heterotrophic aerobic bacteria were able to grow on agar-agar, alginate, cellulose, chitin, dried and ground Fucus vesiculosus, Marine Broth 2216, palmitate, and starch. MPN counts of bulk water and aggregate samples ranged between 0.18?×?101 and 1.1?×?106 cells per milliliter and those of the sediment surface and the transition zone between 0.8?×?101 and 5.1?×?107 cells per gram dry weight. Marine Broth and F. vesiculosus yielded the highest values of all substrates tested and corresponded to 2.3–32% of 4,6-diamidinophenyl indole cell counts. Strains of seven phylogenetic classes were obtained: Actinobacteria, Bacilli, α- and γ-Proteobacteria, Sphingobacteria, Flavobacteria, and Planctomycetacia. Only with agar-agar as substrate could organisms of all seven classes be isolated.  相似文献   

13.
Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s???1 (for Rayleigh waves) and 8 ×10???3 s???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.  相似文献   

14.
Antarctic Intermediate Water (AAIW) formation constitutes an important mechanism for the export of macronutrients out of the Southern Ocean that fuels primary production in low latitudes. We used quality-controlled gridded data from five hydrographic cruises between 1990 and 2014 to examine decadal variability in nutrients and dissolved inorganic carbon (DIC) in the AAIW (neutral density range 27 < γ n <? 27.4) along the Prime Meridian. Significant positive trends were found in DIC (0.70 ± 0.4 μmol kg??1 year??1) and nitrate (0.08 ± 0.06 μ mol kg??1 year??1) along with decreasing trends in temperature (??0.015 ± 0.01°C year??1) and salinity (??0.003 ± 0.002 year??1) in the AAIW. Accompanying this is an increase in apparent oxygen utilization (AOU, 0.16 ± 0.07 μ mol kg??1 year??1). We estimated that 75% of the DIC change has an anthropogenic origin. The remainder of the trends support a scenario of a strengthening of the upper-ocean overturning circulation in the Atlantic sector of the Southern Ocean in response to the positive trend in the Southern Annular Mode. A decrease in net primary productivity (more nutrients unutilized) in the source waters of the AAIW could have contributed as well but cannot fully explain all observed changes.  相似文献   

15.
Based on the light-photosynthesis response measurement at leaf level, combined with over-and under-canopy eddy covariance measurements, research on photosynthetic characteristics of single trees and forest canopy was conducted. The relationship between light intensity and photosynthetic rates for leaves and canopy can be well fitted by a non-rectangular hyperbola model. Mongolian oak presented a high light compensation point, L cp (28 μmol·m?2·s?1), a light saturation point L sp (>1800 μmol·m?2·s?1), and a maximal net photosynthetic rate P max (9.96 μmol·m?2·s?1), which suggest that it is a typical heliophilous plant. Mono maple presented the highest apparent quantum efficiency α (0.066) but the lowest, L cp (16 μmol·m?2·s?1), L sp (≈800 μmol·m?2·s?1), and P max (4.51 μmol·m?2·s?1), which suggest that it is heliophilous plant. Korean pine showed the lowest α value but a higher P max, which suggest that it is a semi-heliophilous plant. At the canopy level, the values of both α and P max approached the upper limit of reported values in temperate forests, while L cp was within the lower limit. Canopy photosynthetic characteristics were well consistent with those of leaves. Both showed a high ability to photosynthesize. However, environmental stresses, especially high vapor pressure deficits, could significantly reduce the photosynthetic ability of leaves and canopy.  相似文献   

16.
The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M0?=?1.12×1017 Nm, M w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω?2 source model) by attenuation parameters Q(f)?=?500f0.9, κ?=?0.04 s, and fmax?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.  相似文献   

17.
Forty-six mining-induced seismic events with moment magnitude between ?1.2 and 2.1 that possibly caused damage were studied. The events occurred between 2008 and 2013 at mining level 850–1350 m in the Kiirunavaara Mine (Sweden). Hypocenter locations were refined using from 6 to 130 sensors at distances of up to 1400 m. The source parameters of the events were re-estimated using spectral analysis with a standard Brune model (slope ?2). The radiated energy for the studied events varied from 4.7 × 10?1 to 3.8 × 107 J, the source radii from 4 to 110 m, the apparent stress from 6.2 × 102 to 1.1 × 106 Pa, energy ratio (E s/E p) from 1.2 to 126, and apparent volume from 1.8 × 103 to 1.1 × 107 m3. 90% of the events were located in the footwall, close to the ore contact. The events were classified as shear/fault slip (FS) or non-shear (NS) based on the E s/E p ratio (>10 or <10). Out of 46 events 15 events were classified as NS located almost in the whole range between 840 and 1360 m, including many events below the production. The rest 31 FS events were concentrated mostly around the production levels and slightly below them. The relationships between some source parameters and seismic moment/moment magnitude showed dependence on the type of the source mechanism. The energy and the apparent stress were found to be three times larger for FS events than for NS events.  相似文献   

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

19.
The concentration of rock-forming elements, the static magnetic susceptibility κ, spectra of electron paramagnetic resonance, and their relative intensities I are studied in samples from a borehole drilled in Cenozoic sedimentary deposits of southern Western Siberia. All measured values experience appreciable irregular variations with depth. A linear dependence exists between κ and I within the range of their medium and large values; κ and I have maximum values in the same sample, and κmax = 1920 × 10?6SI, κmin = 210 × 10?6 SI, and κav = 630 × 10?6 SI. The magnetic properties of the samples are controlled by Fe2+ ions present in clastic material and by microphases (clusters) with Fe3+ ions of the goethite and lepidocrocite type present in the cement. The theoretically possible magnetic susceptibility of the Fe2+ ion system (provided that all iron exists in this form) is quite comparable with κmin but, even with very high concentrations of Fe2+, does not reach half of κav: (154 < κ(Fe2+) < 254) × 10?6 SI. Anomalously high values of κ are due to a large number of clusters with Fe3+ ions if structural units FeOOH do not dissociate and the interaction of the clusters with hydroxides of aluminum and precipitation medium impedes the process of their coagulation. Otherwise, the cluster sizes gradually increase, an antiferromagnetic structure develops in clusters, and the magnetic susceptibility decreases.  相似文献   

20.
High-biomass red tides occur frequently in some semi-enclosed bays of Hong Kong where ambient nutrients are not high enough to support such a high phytoplankton biomass. These high-biomass red tides release massive inorganic nutrients into local waters during their collapse. We hypothesized that the massive inorganic nutrients released from the collapse of red tides would fuel growth of other phytoplankton species. This could influence phytoplankton species composition. We tested the hypothesis using a red tide event caused by Mesodinium rubrum (M. rubrum) in a semi-enclosed bay, Port Shelter. The red tide patch had a cell density as high as 5.0×105 cells L?1, and high chlorophyll a (63.71 μg L?1). Ambient inorganic nutrients (nitrate: \(\rm{NO}_3^-\), ammonium: \(\rm{NH}_4^+\), phosphate: \(\rm{PO}_4^{3-}\), silicate: \(\rm{SiO}_4^{3-}\)) were low both in the red tide patch and the non-red-tide patch (clear waters outside the red tide patch). Nutrient addition experiments were conducted by adding all the inorganic nutrients to water samples from the two patches followed by incubation for 9 days. The results showed that the addition of inorganic nutrients did not sustain high M. rubrum cell density, which collapsed after day 1, and did not drive M. rubrum in the non-red-tide patch sample to the same high-cell density in the red tide patch sample. This confirmed that nutrients were not the driving factor for the formation of this red tide event, or for its collapse. The death of M. rubrum after day 1 released high concentrations of \(\rm{NO}_3^-\), \(\rm{PO}_4^{3-}\), \(\rm{SiO}_4^{3-}\), \(\rm{NH}_4^+\), and urea. Bacterial abundance and heterotrophic activity increased, reaching the highest on day 3 or 4, and decreased as cell density of M. rubrum declined. The released nutrients stimulated growth of diatoms, such as Chaetoceros affinis var. circinalis, Thalassiothrix frauenfeldii, and Nitzschia sp., particularly with additions of \(\rm{SiO}_4^{3-}\) treatments, and other species. These results demonstrated that initiation of M. rubrum red tides in the bay was not directly driven by nutrients. However, the massive inorganic nutrients released from the collapse of the red tide could induce a second bloom in low-ambient nutrient water, influencing phytoplankton species composition.  相似文献   

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

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