首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.

Six species of more than 20-year-old desert woody plants in the oasis-desert ecotone were selected for study. The results showed that: (1) in different growing seasons δ13C values of assimilating organ varied between -14%. and -16%. for Haloxylon ammodendron (HA),-14%. - -15%. for Calligonum mongolicum (CM) and -25%. - -28%. for Caragana korshinskii (CK), Nitraria sphaerocarpa (NS) and Hedysarum scoparium (HS). (2) The net photosynthetic rate (Pn) of HA and CM was significantly higher than those of the other species. With the decrease in Pn for the six species, their intercellular CO2 concentration increased, but stomatal limitation value decreased under the intensive light. At the same time, the photochemical efficiency of PS II dropped to different degrees. (3) The CO2 enrichment experiment demonstrated that, Pn of HA and CM increased to different extent under 450 μmol/mol, but their Pn reduced or approximated to the current condition under 650 μmol/mol. Under 450 μmol/mol the efficiency for solar energy utilization of CK and HS significantly reduced and under 650 μmol/mol their respiration rate exceeded photosynthesis rate. It can be concluded that HA and CM have some function of pathway for C4 but the other three species have the function for C3. The decline in their Pn is mainly caused by non-stomatal factors. HA, CM, CK and HS exhibited photoinhibition, which disappeared in a short time. This is a kind of positive readjustment to adapting to the desert environment. HA and CM can adapt to the high CO2 environment, but CK and HS cannot. With the rise in atmospheric CO2 concentration and climate warming, the latter two species in the oasis-desert ecotone may be gradually degraded or even disappear.

  相似文献   

2.
We studied the effects of expected end-of-the-century pCO2 (1000 ppm) on the photosynthetic performance of a coastal marine cyanobacterium Synechococcus sp. PCC7002 during the lag, exponential, and stationary growth phases. Elevated pCO2 significantly stimulated growth, and enhanced the maximum cell density during the stationary phase. Under ambient pCO2 conditions, the lag phase lasted for 6 days, while elevated pCO2 shortened the lag phase to two days and extended the exponential phase by four days. The elevated pCO2 increased photosynthesis levels during the lag and exponential phases, but reduced them during the stationary phase. Moreover, the elevated pCO2 reduced the saturated growth light (Ik) and increased the light utilization efficiency (α) during the exponential and stationary phases, and elevated the phycobilisome:chlorophyll a (Chl a) ratio. Furthermore, the elevated pCO2 reduced the particulate organic carbon (POC):Chl a and particulate organic nitrogen (PON):Chl a ratios during the lag and stationary phases, but enhanced them during the exponential phase. Overall, Synechococcus showed differential physiological responses to elevated pCO2 during different growth phases, thus providing insight into previous studies that focused on only the exponential phase, which may have biased the results relative to the effects of elevated pCO2 in ecology or aquaculture.  相似文献   

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

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

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

6.
The recent seismicity catalogue of metropolitan France Sismicité Instrumentale de l’Hexagone (SI-Hex) covers the period 1962–2009. It is the outcome of a multipartner project conducted between 2010 and 2013. In this catalogue, moment magnitudes (M w) are mainly determined from short-period velocimetric records, the same records as those used by the Laboratoire de Détection Géophysique (LDG) for issuing local magnitudes (M L) since 1962. Two distinct procedures are used, whether M L-LDG is larger or smaller than 4. For M L-LDG >4, M w is computed by fitting the coda-wave amplitude on the raw records. Station corrections and regional properties of coda-wave attenuation are taken into account in the computations. For M L-LDG ≤4, M w is converted from M L-LDG through linear regression rules. In the smallest magnitude range M L-LDG <3.1, special attention is paid to the non-unity slope of the relation between the local magnitudes and M w. All M w determined during the SI-Hex project is calibrated according to reference M w of recent events. As for some small events, no M L-LDG has been determined; local magnitudes issued by other French networks or LDG duration magnitude (M D) are first converted into M L-LDG before applying the conversion rules. This paper shows how the different sources of information and the different magnitude ranges are combined in order to determine an unbiased set of M w for the whole 38,027 events of the catalogue.  相似文献   

7.
An open-path eddy covariance system was set up in Damxung rangeland station to measure the carbon flux from July to October, 2003. The canopy quantum yield (α) of alpine meadow was calculated by the linear function between the net ecosystem carbon dioxide exchange (NEE) and the photosynthetic active radiation (PAR) under low light, and how it was influenced by the temperature was also discussed. Results showed that the canopy α decreased almost linearly with temperature, with the decrease in every 1 °C increase of temperature similar to those measured on leaf level of C3 plant. At the beginning, the decrease of canopy α with temperature was 0.0005 μmol CO2·μmol?1 PAR; while it increased to 0.0008 μmol CO2·μmol?1 PAR in September, showing a rising trend with plant growing stages. Compared with the canopy α calculated with rectangular hyperbola function, the value in the paper was lower. However, the method advanced here has the advantages in examining the relationship between α and the key environmental factors, such as temperature.  相似文献   

8.
The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E y component, dynamic pressure P d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E y value at E y < 12 mV m?1 and are inversely proportional to this value at E y > 12 mV m?1 for the first-type events. On the contrary, the dependence of Dst on E y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.  相似文献   

9.
We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ c and P d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.  相似文献   

10.
Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F 2 region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F 2 region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.  相似文献   

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

12.
Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C v value and to explore its physical meaning experimentally. The result shows that the change of C v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C v value increases firstly and then decreases before the instability, and the instability occurs when C v value decreases to the level before increasing. This indicates that C v value may be a useful parameter for earthquake prediction.  相似文献   

13.
The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N2 and O2 molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O+ on N e ions under winter and summer conditions forms not more than 11% of the N e winter anomaly event in the F 2 layer and topside ionosphere. The model without electronically excited O+ ions reduces the duration of the N e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.  相似文献   

14.
Seasonal metrics and environmental responses to forestry soil surface CO2 emission effluxes among three types of lower subtropical forests were consistently monitored over two years with static chamber-gas chromatograph techniques among three types of lower subtropical forests. Results showed that annual CO2 effluxes (S+L) reached 3942.20, 3422.36 and 2163.02 CO2 g·m?2·a?1, respectively in the monsoon evergreen broadleaf forest, mixed broadleaf-coniferous forest and coniferous forest. All the three types of forests revealed the same characteristics of seasonal changes with the CO2 effluxes peaking throughout June to August. During this peaking period, the effluxes were 35.9%, 38.1% and 40.2% of the total annual effluxes, respectively. The CO2 emission process responding to the environmental factors displayed significantly different patterns in forestry soils of the three types of forests. The coniferous forest (CF) was more sensitive to temperature than the other two types. The Q 10 values were higher, along with greater seasonal variations of the CO2 efflux, indicating that the structurally unique forestry ecosystem has disadvantage against interferences. All the three types of forestry CO2 effluxes showed significant correlation with the soil temperature (T s), soil water content (M s) and air pressure (P a). However, stepwise regression analysis indicated no significant correlation between air pressure and the soil CO2 efflux. With an empirical model to measure soil temperature and water content in 5 cm beneath the soil surface, the CO2 effluxes accounting for 75.7%, 77.8% and 86.5% of the efflux variability respectively in soils of BF, MF and PF were calculated. This model can be better used to evaluate the CO2 emission of soils under water stress and arid or semi-arid conditions.  相似文献   

15.
We propose an algorithm that combines a pre-processing step applied to the a priori state vector prior to retrievals, with the modified damped Newton method (MDNM), to improve convergence. The initial constraint vector pre-processing step updates the initial state vector prior to the retrievals if the algorithm detects that the initial state vector is far from the true state vector in extreme cases where there are CO2 emissions. The MDNM uses the Levenberg-Marquardt parameter γ, which ensures a positive Hessian matrix, and a scale factor α, which adjusts the step size to optimize the stability of the convergence. While the algorithm iteratively searches for an optimized solution using observed spectral radiances, MDNM adjusts parameters γ and α to achieve stable convergence. We present simulated retrieval samples to evaluate the performance of our algorithm and comparing it to existing methods. The standard deviation of our retrievals adding random noise was less than 3.8 ppmv. After pre-processing the initial estimate when it was far from the true value, the CO2 retrieval errors in the boundary layers were within 1.2 ppmv. We tested the MDNM algorithm’s performance using GOSAT L1b data with cloud screening. Our preliminary validations comparing the results to TCCON FTS measurements showed that the average bias was less than 1.8 ppm and the correlation coefficient was approximately 0.88, which was larger than for the GOSAT L2 product.  相似文献   

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

17.
The structure and dynamics of the ionosphere and plasmasphere at high solar activity under quiet geomagnetic conditions of June 2–3, 1979, and January 5–6, 1980, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The plasma drift velocity, determined by comparing the calculated and measured heights of the F 2 layer maximum (hmF2), and the correction of [N2] and [O2], found in the NRLMSISE-00 model, make it possible to coordinate the electron densities (NmF2) calculated at the hmF2 height and the measured anomalous variations in NmF2 over the Argentine Islands ionosonde as well as the calculated and measured NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that, if the interference of the diffusion velocities of O+(4S) and H+ ions is taken into account, the additional heating of plasmaspheric electrons leads to an increase in the flux of O+(4S) ions from the topside ionosphere to lower F 2 layer altitudes, as a result of which an anomalous nighttime increase in NmF2 6, observed on January 6, 1980, over Millstone Hill station, is mainly produced. The second component of the formation of anomalous night-time NmF2 is the plasma drift along the magnetic field caused by the neutral wind, which shifts O+(4S) ions to higher altitudes where the recombination rate of O+(4S) with N2 and O2 is lower and slows down a decrease in NmF2 in the course of time. It has been shown that the influence of electronically excited O+ ions and vibrationally excited N2 and O2 molecules on electron density (N e ) considerably differs under winter and summer conditions. This difference forms significant part of the winter anomaly in N e at heights of the F 2 region and topside ionosphere over Millstone Hill station.  相似文献   

18.
The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b m , N, M 1, M 2) that an observed b m estimate, based on a sample of N magnitudes within an [M 1????≤?ΔM/2,?M 2?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.  相似文献   

19.
We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h Am , which corresponds to the maximum increase, and layer peak height h m F, while h Am is always less than h m F. The difference between h Am and h m F linearly increases with increasing h m F. Whereas the difference is ~38 km for h m F = 280 km, it is ~54 km for h m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.  相似文献   

20.
The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (Ms = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (Ms = 5.6), Andrei-Tas (Ms = 6.1), and Ilin-Tas (Ms = 6.9) earthquakes.  相似文献   

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

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