首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M w and energy magnitude M e is outlined and critically discussed. The formulas for M w and M e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M 0 and radiated seismic energy E S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M w and M e are linked via the parameter Θ = log(E S/M 0), and the formula for M e can be written as M e = M w + (Θ + 4.7)/1.5. This relationship directly links M e with M w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M w and M e can significantly differ. In fact, Θ is assumed to be constant when calculating M w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V R and seismic wave radiation efficiency η R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M e may sometimes differ by more than one magnitude unit from M w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M e is more suitable than M w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E S and M e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E S in characterizing the seismic source, the use of M e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E S estimations in order to allow M e to be extensively used as an important complement to M w in common seismological practice and its applications.  相似文献   

2.
The static closed chamber technique is used in the study on the CH4 and N2O fluxes from the soils of primeval Abies fabri forest, the succession Abies fabri forest and the clear-cut areas of mid-aged Abies fabri forest in the Gongga Mountain from May 1998 to September 1999. The results indicate the following: (i) The forest soil serves as the source of atmospheric N2O at the three measurement sites, while the fluxes of CH4 are all negative, and soil is the sink of atmospheric CH4. The comparative relations of N2O emissions between the three sites are expressed as primeval Abies fabri forest < clear-cut areas < succession Abies fabri forest, and those of CH4 consumption fluxes are primeval Abies fabri forest < succession Abies fabri forest < clear-cut areas, (ii) Significant seasonal variations of N2O emission at various sites were observed, and two emission peaks of N2O occurr during summer (July—August) and spring (February—March), whereas N2O emission is relatively low in winter and spring (mid March—April). Seasonal variations of CH4 consumption at each measurement site fluctuate drastically with unclear regularities. Generally, CH4 consumption fluxes of succession Abies fabri forest and clear-cut areas are higher from mid May to late July but lower in the rest of sampling time, while the CH4 flux keeps a relatively high value even up to September in primeval Abies fabri forest. In contrast to primeval Abies fabri forest, the CH4 absorbabilities of succession Abies fabri forest and clear-cut areas of mid-aged Abies fabri forest are weaker. Particularly, the absorbability of the clear-cut areas is even weaker as compared with the other two sites, for the deforestation reduces the soil absorbability of atmospheric CH4. (iii) Evident diurnal variation regularity exists in the N2O emissions of primeval Abies fabri forest, and there is a statistic positive correlation between the fluxes of N2O and air temperature (R=0.95, n=11, α· 0.01), and also the soil temperature of 5-cm layer (R=0.81, n=11, α> 0.01), whereas the CH4 diurnal variation regularities are unclear and have no significant correlation with the soil temperature of 5-cm layer and air temperature.  相似文献   

3.
(Fe4Cr4Ni)9C4 is a metal carbide mineral formed by combination of Fe, Cr and Ni with C. It occurs in a chromite deposit in the Luobusha ophiolite, Tibet. Based on the determination of its crystal structure, the empirical formula is (Fe4.12Cr3.84Ni0.96)8.92C3.70 and the simplified formula is (Fe4Cr4Ni)4C9. The mineral is hexagonal with a = 1.38392(2) nm, c = 0.44690(9) nm, pace group P63 m c, Z=6 and the calculated specific gravity Dx = 7.089 g/cm3. Fe, Cr and Ni occupy different crystallographic sites and their coordination numbers are approximately 12, forming an alternate stacking sequence of flat and puckered layers along the c axis. Some metallic atoms have a defect structure. The interatomic distances of Fe, Cr and Ni are 0.2525-0.2666 nm, and the distances between Fe, Cr, Ni and C are 0.1893-0.2169 nm. The coordination number of carbon is 6. It occurs in interstices of the metallic atoms Fe, Cr and Ni to form trigonalprismatically coordinated polyhedra. These coordination polyhedra are linked with each other via shared corners or shared edges into a new type of metal carbide structure.  相似文献   

4.
On the basis of Parry’s method (1986), an improved method was established to determine the molar volume (Vm) and compositions (X) of the NaCl-H2O-CO2 (NHC) system inclusion. To use this method, the determination of Vm-X only requires three microthermometric data of a NHC inclusion: partial homog-enization temperature (Th ,CO2), salinity (S) and total homogenization temperature (Th). Theoretically, four associated equations are needed containing four unknown parameters: X CO2, XNaCl, Vm and F (volume fraction of CO2 phase in total inclusion when occurring partial homogenization). When they are released, the Vm-X are determined. The former three equations, only correlated with Th ,CO2, S and F, have simplified expressions:XCO2=f1(Th,CO2,S,F),XNaCl=f2(Th,CO2,S,F),Vm=f3(Th,CO2,S,F). The last one is the thermodynamic relationship of X CO2, XNaCl, Vm and Th:f4(XCO2,XNaCl,Vm,Th)=0.Since the above four associated equations are complicated, it is necessary to adopt iterative technique to release them. The technique can be described by:(i) Freely input a F value (0≤F≤1),with Th ,CO2 and S, into the former three equations. As a result,X CO 2,XNaCl and the molar volume value recorded as Vm1 are derived. (ii) Input the X CO2 and XNaCl gotten in the step above into the last equation, and another molar volume value recorded as Vm2 is determined. (iii) If Vm1 is unequal to Vm2, the calculation will be restarted from “(i)”. The iteration is completed until Vm1 is equal to Vm2, which means that the four associated equations are released. Compared to Parry’s (1986) solution method, the improved method is more convenient to use, as well as more accurate to determine X CO 2. It is available for a NHC inlusion whose partial homogenization temperature is higher than clatherate melting temperature and there are no solid salt crystals in the inclusion at parital homogenization.  相似文献   

5.
Source rock extracts and crude oils from the Songliao Basin were analyzed by high-temperature gas chromatography (HTGC), gas chromatography-mass spectrometry (HTGC-MS) and gas chromatography-isotope ratio-mass spectrometry (GC-IRMS), for high molecular-weight alkanes. The distributions of n-alkanes in the Nenjiang Formation extracts are in the C14―C63 range; a bimodal distribution occurs in the C-21 and C21―40 regions. The C30―C37 n-alkanes are accompanied by C29―C35 hopanes, whereas the high molecular-weight C45―C47 n-alkanes co-occur with abundant isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The high δ 13C values of the n-alkanes and the microscopic maceral compositions indicate a highly diversified organic source input for the Nenjiang Formation source rocks, ranging from aquatic plants, blue alge-bacteria, to land plant material. In contrast, n-alkanes in the rock extracts of the Qingshankou Formation are characterized by a single modal distribution, with relatively low abundances of C29―C35 hopanes, but high molecular-weight isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The relatively low δ 13C values of C22―C44 n-alkanes and organic material compositions indicate that the source rocks in the Qingshankou Formation contain dominantly type I algal organic matter. The relative abundance of C 40 compounds in source rocks changes little at low maturity stage, but decreases drastically at higher maturity levels, with a concurrent reduction in the odd/even carbon predominance. In crude oils, in contrast, the relative abundance of C 40 compounds appears to relate closely with the oil source and oil viscosity.  相似文献   

6.
This paper provides a generic equation for the evaluation of the maximum earthquake magnitude mmax for a given seismogenic zone or entire region. The equation is capable of generating solutions in different forms, depending on the assumptions of the statistical distribution model and/or the available information regarding past seismicity. It includes the cases (i) when earthquake magnitudes are distributed according to the doubly-truncated Gutenberg-Richter relation, (ii) when the empirical magnitude distribution deviates moderately from the Gutenberg-Richter relation, and (iii) when no specific type of magnitude distribution is assumed. Both synthetic, Monte-Carlo simulated seismic event catalogues, and actual data from Southern California, are used to demonstrate the procedures given for the evaluation of mmax.The three estimates of mmax for Southern California, obtained by the three procedures mentioned above, are respectively: 8.32 ± 0.43, 8.31 ± 0.42 and 8.34 ± 0.45. All three estimates are nearly identical, although higher than the value 7.99 obtained by Field et al. (1999). In general, since the third procedure is non-parametric and does not require specification of the functional form of the magnitude distribution, its estimate of the maximum earthquake magnitude mmax is considered more reliable than the other two which are based on the Gutenberg-Richter relation.  相似文献   

7.
We use 576 earthquakes of magnitude, M w, 3.3 to 6.8 that occurred within the region 33° N–42.5° N, 19° E–30° E in the time period 1969 to 2007 to investigate the stability of the relation between moment magnitude, M w, and local magnitude, M L, for earthquakes in Greece and the surrounding regions. We compare M w to M L as reported in the monthly bulletins of the National Observatory of Athens (NOA) and to M L as reported in the bulletins of the Seismological Station of the Aristotle University of Thessaloniki. All earthquakes have been analyzed through regional or teleseismic waveform inversion, to obtain M w, and have measured maximum trace amplitudes on the Wood–Anderson seismograph in Athens, which has been in operation since 1964. We show that the Athens Wood–Anderson seismograph performance has changed through time, affecting the computed by NOA M L by at least 0.1 magnitude units. Specifically, since the beginning of 1996, its east–west component has been recording systematically much larger amplitudes compared to the north–south component. From the comparison between M w and M L reported by Thessaloniki, we also show that the performance of the sensors has changed several times through time, affecting the calculated M L’s. We propose scaling relations to convert the M L values reported from the two centers to M w. The procedures followed here can be applied to other regions as well to examine the stability of magnitude calculations through time.  相似文献   

8.
This paper presents a robust H∞ output feedback control approach for structural systems with uncertainties in model parameters by using available acceleration measurements and proposes conditions for the existence of such a robust output feedback controller. The uncertainties of structural stiffness, damping and mass parameters are assumed to be norm-bounded. The proposed control approach is formulated within the framework of linear matrix inequalities, for which existing convex optimization techniques, such as the LMI toolbox in MATLAB, can be used effectively and conveniently. To illustrate the effectiveness of the proposed robust H∞ strategy, a six-story building was subjected both to the 1940 El Centro earthquake record and to a suddenly applied Kanai-Tajimi filtered white noise random excitation. The results show that the proposed robust H∞ controller provides satisfactory results with or without variation of the structural stiffness, damping and mass parameters.  相似文献   

9.
We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q p ), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q p with an interface at 660 km depth and take Q p constant with period. Our data set is too small to solve for a possible frequency dependence of Q p . The upper mantle Q 1 is 476 [299–1176] and the lower mantle Q 2 is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q p –1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135.  相似文献   

10.
The work presents statistical methods for estimating the distribution parameters of rare, strong earthquakes. Using the two main theorems of extreme value theory (EVT), the distribution of T-maximum (the maximum magnitude over the time period T). Two methods for estimating the parameters of this distribution are proposed using the Generalized Pareto Distribution (GPD) and the General Extreme Value Distribution (GEV). In addition, the that allow the determination of the distribution of the T-maximum for an arbitrary value of T are proposed. The approach being used clarifies the nature of the instability of the widely accepted M max parameter. In the work, instead of unstable values of the M max parameter, the robust parameter Q T (q), the q level quantile for the distribution of the T-maximum, is proposed to be used. The described method has been applied to the Harvard Catalogue of Seismic Moments of 1977–2006 and to the Magnitude Catalogue for Fennoscandia in 1900–2005. Moreover, the estimates of parameters of the corresponding GPD and GEV distributions, in particular, the most interesting shape parameter and the values of the M max and Q T (q) parameters are given.  相似文献   

11.
We present an uncertainty analysis of ecological process parameters and CO2 flux components (R eco, NEE and gross ecosystem exchange (GEE)) derived from 3 years’ continuous eddy covariance measurements of CO2 fluxes at subtropical evergreen coniferous plantation, Qianyanzhou of ChinaFlux. Daily-differencing approach was used to analyze the random error of CO2 fluxes measurements and bootstrapping method was used to quantify the uncertainties of three CO2 flux components. In addition, we evaluated different models and optimization methods in influencing estimation of key parameters and CO2 flux components. The results show that: (1) Random flux error more closely follows a double-exponential (Laplace), rather than a normal (Gaussian) distribution. (2) Different optimization methods result in different estimates of model parameters. Uncertainties of parameters estimated by the maximum likelihood estimation (MLE) are lower than those derived from ordinary least square method (OLS). (3) The differences between simulated Reco, NEE and GEE derived from MLE and those derived from OLS are 12.18% (176 g C·m−2·a−1), 34.33% (79 g C·m−2·a−1) and 5.4% (92 g C·m−2·a−1). However, for a given parameter optimization method, a temperature-dependent model (T_model) and the models derived from a temperature and water-dependent model (TW_model) are 1.31% (17.8 g C·m−2·a−1), 2.1% (5.7 g C·m−2·a−1), and 0.26% (4.3 g C·m−2·a−1), respectively, which suggested that the optimization methods are more important than the ecological models in influencing uncertainty in estimated carbon fluxes. (4) The relative uncertainty of CO2 flux derived from OLS is higher than that from MLE, and the uncertainty is related to timescale, that is, the larger the timescale, the smaller the uncertainty. The relative uncertainties of Reco, NEE and GEE are 4%−8%, 7%−22% and 2%−4% respectively at annual timescale. Supported by the National Natural Science Foundation of China (Grant No. 30570347), Innovative Research International Partnership Project of the Chinese Academy of Sciences (Grant No. CXTD-Z2005-1) and National Basic Research Program of China (Grant No. 2002CB412502)  相似文献   

12.
Soil H2 and CO2 surveys were carried out along seven active faults and around the aftershock region of the 2000 Tottori-ken Seibu earthquake in Japan. Diffuse CO2 effluxes were also measured along one fault and around the 2000 aftershock region. The results show highly variable H2 concentration in space and time and it seems that the maximum H2 concentration at each active fault correlates with fault activity as exemplified by the time of the latest big earthquakes. Even though observed H2 concentrations in four faults were markedly lower than those collected previously in the latter half of the 1970s, it is evident that the higher H2 concentrations in this study are due to the addition of the fault gases. Comparing the chemical composition of trapped gases (H2: 5–20% and CO2/H2: 0.5–12) in fractured rocks of drill cores bored at the Nojima fault, a soil gas sample with the highest H2 concentration showed large amounts of the trapped fault gas, diluted with atmospheric component. The profile experiment across a fracture zone at the Yamasaki fault showed higher H2 concentrations and lower CO2/H2 ratios as was observed in soil gas from the fracture zone. A few days after the 2000 Tottori-kei Seibu earthquake, no CO2 effluxes related to the occurrence of earthquakes were observed at the aftershock region. However, only above the epicenter zone, relatively high H2 concentrations in soil gases were observed.  相似文献   

13.
Variations in carbon isotopic ratios (δ 13C) of C3 plants and distribution of C4 plants were investigated along an altitudinal transect on the eastern slope of Mount Gongga, and the environmental effects on them were discussed. It is shown that plants with C4 photosynthetic pathway mainly occur at altitudes below 2100 m a.s.l., suggesting that the low summer temperature is responsible for the distributional pattern. In addition, δ 13C of C3 plants increases with elevation at the region above 2000 m a.s.l. with the characteristics of humid climate, and the increase rate in δ 13C for C3 plants is about 1.3% per kilometer. Temperature determines the altitudinal trend of δ 13C.  相似文献   

14.
The solubility measurements of WO3 in NaCl-H2O system with 4.0% by weight salt were carried out to study the supercritical phenomena of solubility. These experiments were carried out using rapid-quench pressure vessels, with quantitative Na2WO4-2H2O, HCl, NaCl and H2O in Pt capsule to determine the solubility of WO3 from supersaturated solution. The pressure was kept at 34 MPa, which was near the critical pressure 31.4 MPa and the temperatures varied from 250°C to 550°C. The experimental results indicate that the solubility of WO3 has the features of supercritical phenomenon in critical region and is sensitive to the change of system temperature and solvent density. The temperature and pressure in this experiment are close to those of ore-forming fluids in many tungsten deposits, so the features of supercritical geofluids can be helpful to investigating the ore-forming mechanism of tungsten deposit.  相似文献   

15.
Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy   总被引:1,自引:0,他引:1  
At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO2–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO2 discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d–1. The steam associated with the CO2 output, computed assuming that the original H2O/CO2 ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×1012 J d–1 (17 MW). The location of the CO2 and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO2 degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO2 flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara  相似文献   

16.
The present-day state of studies of vibrationally excited N2 and O2 in the upper atmosphere and the mechanisms by which these molecules affect the electron density and temperature in the ionosphere and neutral densities in the upper atmosphere are considered.  相似文献   

17.
A new modified magnitude scale M S (20R) is elaborated. It permits us to extend the teleseismic magnitude scale M S (20) to the regional epicenter distances. The data set used in this study contains digital records at 12 seismic stations of 392 earthquakes that occured in the northwest Pacific Ocean in the period of 1993–2008. The new scale is based on amplitudes of surface waves of a narrow range of the periods (16–25 s) close to the period of 20 s, for distances of 80–3000 km. The digital Butterworth filter is used for processing. On the basis of the found regional features concerning distance dependence for seismic wave attenuation, all the stations of the region have been subdivided into two groups, namely, “continental” and “island-arc.” For each group of stations, its own calibration function is proposed. Individual station corrections are used to compensate for the local features.  相似文献   

18.
Earthquake surface rupture is the result of transformation from crustal elastic strain accumulation to permanent tectonic deformation. The surface rupture zone produced by the 2001 Kunlunshan earthquake (M w 7.8) on the Kusaihu segment of the Kunlun fault extends over 426 km. It consists of three relatively independent surface rupture sections: the western strike-slip section, the middle transtensional section and the eastern strike-slip section. Hence this implies that the Kunlunshan earthquake is composed of three earthquake rupturing events, i.e. the M w =6.8, M w =6.2 and M w ⩽=7.8 events, respectively. The M w =7.8 earthquake, along the eastern section, is the main shock of the Kunlunshan earthquake, further decomposed into four rupturing subevents. Field measurements indicate that the width of a single surface break on different sections ranges from several meters to 15 m, with a maximum value of less than 30 m. The width of the surface rupture zone that consists of en echelon breaks depends on its geometric structures, especially the stepover width of the secondary surface rupture zones in en echelon, displaying a basic feature of deformation localization. Consistency between the Quaternary geologic slip rate, the GPS-monitored strain rate and the localization of the surface ruptures of the 2001 Kunlunshan earthquake may indicate that the tectonic deformation between the Bayan Har block and Qilian-Qaidam block in the northern Tibetan Plateau is characterized by strike-slip faulting along the limited width of the Kunlun fault, while the blocks themselves on both sides of the Kunlun fault are characterized by block motion. The localization of earthquake surface rupture zone is of great significance to determine the width of the fault-surface-rupture hazard zone, along which direct destruction will be caused by co-seismic surface rupturing along a strike-slip fault, that should be considered before the major engineering project, residental buildings and life line construction. Supported by the National Natural Science Foundation of China (Grant No. 40474037) and the National Basic Research Program of China (Grant No. 2004CB418401)  相似文献   

19.
The characteristics of dayside auroras during the large (16–24 nT) positive values of the IMF B z component, observed on January 14, 1988, during the interaction between the Earth’s magnetosphere and the body of the interplanetary magnetic cloud, have been studied based on the optical observations on Heiss Island. A wide band of diffuse red luminosity with an intensity of 1–2 kilorayleigh (kR) was observed during 6 h in the interval 1030–1630 MLT at latitudes higher than 75° CGL. Rayed auroral arcs, the brightness of which in the 557.7 nm emission sharply increased to 3–7 kR in the postnoon sector immediately after the polarity reversal of the IMF B y component from positive to negative, were continuously registered within the band. Bright auroral arcs were observed at the equatorward edge of red luminosity. It has been found out that the red auroral intensity increases and the band equatorward boundary shifts to lower latitudes with increasing solar wind dynamic pressure. However, a direct proportional dependence of the variations in the auroral features on the dynamic pressure variations has not been found. It has been concluded that the source of bright discrete auroras is located in the region of the low-latitude boundary layer (LLBL) on closed geomagnetic field lines. The estimated LLBL thickness is ∼3 R e . It has been concluded that the intensity of the dayside red band depends on the solar wind plasma density, whereas the position of the position equatorward boundary depends on the dynamic pressure value and its variations.  相似文献   

20.
The problem of determining trends in thermospheric dynamics parameters (horizontal winds) based on analysis of trends in various combinations of ionospheric F 2-layer parameters is formulated. The previous attempts of the authors in this direction are briefly described. It is shown that all studied parameters lead to the same result: after the “boundary date” (approximately 1980) a systematic change in these parameters (a long-term trend) is observed, this fact manifesting changes in the dynamical regime of the thermosphere because of cooling and contraction of the entire middle and upper atmosphere. The results of a search for trends in the hmF2 height for the moment (T(ss) + 2 h) are described. These trends are found higher than the hmF2 trends obtained earlier by various authors analyzing the hmF2 behavior at fixed moments of local time.  相似文献   

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

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