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

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

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

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

5.
A key issue in the design of pile-supported structures on sloping ground is soil–pile interaction, which becomes more complicated in case of dynamic loading. This study aimed to evaluate the effect of slope on the dynamic behavior of pile-supported structures by performing a series of centrifuge tests. Three models were prepared by varying the slope and soil density of dry sand grounds. The mass supported on 3 by 3 group piles was shaken applying sinusoidal wave with various amplitudes. Test results showed that the location of maximum values and distribution shape of the bending moment below the ground surface varied noticeably with the pile position in the slope case. The relationship between the soil resistance and pile deflection (p–yp loops) was carefully evaluated by applying the piecewise cubic spline method to fit the measured bending moment curves along piles. It was found that the shape of the p–yp loops was irregular due to the effect of slope, and immensely influenced by the movement of the unstable zone. In addition, the effect of the pile group in the horizontal case was evaluated by comparing with the previously suggested curves that represent the relationship between the soil resistance and pile–soil relative displacement (p–y curves) to propose the multiplier coefficients.  相似文献   

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

7.
Over the past decade, the H-κ stacking technique of Zhu and Kanamori (J Geophys Res 105:2969–2980, 2000) has become a standard tool to determine the crustal thickness H and the bulk crustal vP/vS ratio κ from teleseismic receiver functions. It is obvious that unfavorable noise conditions as well as a complex 3D velocity structure can severely hamper the interpretation of receiver-function data. However, we observe that ambiguities can even arise from a simple 1D layered velocity structure which raises a high potential for misinterpretations. To analyze the feasibility and basic limitations of the H-κ stacking method, we conduct a series of tests based on synthetic data. The impact of different given elementary parameters, related either to the velocity structure or to the data processing, is evaluated in a series of eight individual tests. We deliberately exclude complications such as 3D structural variations and/or noise to show that even a simple 1D velocity structure, involving, e.g., an additional inter-crustal discontinuity, can have significant consequences for the interpretation of the results. However, our modeling suggests that more complex crustal structures may lead to even less reliable results. Additionally, our tests illustrate that time shifts of the maxima in the H-κ domain due to the superposition and merging of individual phases can lead to significantly overestimated vP/vS ratios. In general, the depth to the Moho (or other discontinuities of interest) is less significantly affected. Our tests indicate the necessity to accurately check delay times derived from the maxima of the H-κ stacks against corresponding phases in the receiver functions. Repeating the stacking with varied weighting factors and filter ranges can help to reduce the ambiguities and to avoid possible misinterpretation.  相似文献   

8.
A 3D, two-time-level, σS-z-σB hybrid-coordinate Marine Science and Numerical Modeling numerical ocean circulation model (HyMOM) is developed in this paper. In HyMOM, the σ coordinate is employed in the surface and bottom regions, and the z coordinate is used in the intermediate layers. This method can overcome problems with vanishing surface cells and minimize the unwanted deviation in representing bottom topography. The connection between the σ and z layers vertically includes an expanded “ghost” method and the linear interpolation. The governing equations in the σS-z-σB hybrid coordinate based on the complete Reynolds-averaged Navier-Stokes equations are derived in detail. The two-level time staggered and Eulerian forward and backward schemes, which are of second-order of accuracy, are adopted for the temporal difference in internal and external mode, respectively. The computation of the baroclinic gradient force is tested in an analytic test problem; the errors for two methods in HyMOM, which are relatively large only in the bottom layers, are obviously smaller than those in the pure σ and z models in almost all of the vertical layers. A quasi-global climatologic numerical experiment is constructed to test the simulation performance of HyMOM. With the monthly mean Levitus climatology data as reference, the HyMOM can improve the simulating accuracy compared with its pure z or σ coordinate implementation.  相似文献   

9.
This study analyzes and compares the P- and S-wave displacement spectra from local earthquakes and explosions of similar magnitudes. We propose a new approach to discrimination between low-magnitude shallow earthquakes and explosions by using ratios of P- to S-wave corner frequencies as a criterion. We have explored 2430 digital records of the Israeli Seismic Network (ISN) from 456 local events (226 earthquakes, 230 quarry blasts, and a few underwater explosions) of magnitudes Md?=?1.4–3.4, which occurred at distances up to 250 km during 2001–2013 years. P-wave and S-wave displacement spectra were computed for all events following Brune’s source model of earthquakes (1970, 1971) and applying the distance correction coefficients (Shapira and Hofstetter, Teconophysics 217:217–226, 1993; Ataeva G, Shapira A, Hofstetter A, J Seismol 19:389-401, 2015), The corner frequencies and moment magnitudes were determined using multiple stations for each event, and then the comparative analysis was performed.The analysis showed that both P-wave and especially S-wave displacement spectra of quarry blasts demonstrate the corner frequencies lower than those obtained from earthquakes of similar magnitudes. A clear separation between earthquake and explosion populations was obtained for ratios of P- to S-wave corner frequency f 0(P)/f 0(S). The ratios were computed for each event with corner frequencies f 0 of P- and S-wave, which were obtained from the measured f 0 I at individual stations, then corrected for distance and finally averaged. We obtained empirically the average estimation of f 0(P)/f 0(S)?=?1.23 for all used earthquakes, and 1.86 for all explosions. We found that the difference in the ratios can be an effective discrimination parameter which does not depend on estimated moment magnitude M w .The new multi-station Corner Frequency Discriminant (CFD) for earthquakes and explosions in Israel was developed based on ratios P- to S-wave corner frequencies f 0(P)/f 0(S), with the empirical threshold value of the ratio for Israel as 1.48.  相似文献   

10.
This study uses macroseismic data and wave equations to solve the problem of ultra long propagation of felt ground motion (over 9000 km from the epicenter) due to the Sea-of-Okhotsk earthquake. We show that the principal mechanism of this phenomenon could be excitation of a previously unknown standing radial wave as a mode of the Earth’s free oscillations, 0S0, due to the superposition of an incident and a reflected spherical P wave in the epicentral area of the Sea-of-Okhotsk earthquake. The standing wave generates slowly attenuating P waves that travel over the earth’s surface that act as carrying waves; when superposed on these, direct body waves acquire the ability to travel over great distances. We show previously unknown parameters of the radial mode 0S0 for the initial phase of earth deformation due to the large deep-focus earthquake. We used data on the Sea-of-Okhotsk and Bolivian earthquakes to show that large deep-focus earthquakes can excite free oscillations of the Earth that are not only recorded by instrumental means, but are also felt by people, with the amplification of the macroseismic effect being directly related to the phenomenon of resonance for multistory buildings.  相似文献   

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

12.
This paper reports a study of the Tolud earthquake sequence; the sequence was a burst of shallow seismicity between November 28 and December 7, 2012; it accompanied the initial phase in the Tolbachik Fissure Eruption of 2012?2013. The largest earthquake (the Tolud earthquake of November 30, 2012, to be referred to as the Tolud Earthquake in what follows, with KS = 11.3, ML = 4.9, MC = 5.4, and MW = 4.8) is one of the five larger seismic events that have been recorded at depths shallower than 10 km beneath the entire Klyuchevskoi Volcanic Cluster in 1961?2015. It was found that the Tolud earthquake sequence was the foreshock–aftershock process of the Tolud Earthquake. This is one of the larger seismicity episodes ever to have occurred in the volcanic areas of Kamchatka. Data of the Kamchatka seismic stations were used to compute some parameters for the Tolud Earthquake and its largest (ML = 4.3) aftershock; the parameters include the source parameters and mechanisms, and the moment magnitudes, since no information on these is available at the world seismological data centers. The focal mechanisms for the Tolud Earthquake and for its aftershock are consistent with seismic ruptures at a tension fault in the rift zone. Instrumental data were used to estimate the intensity of shaking due to the Tolud Earthquake. We discuss the sequence of events that was a signature of the time-dependent seismic and volcanic activity that took place in the Tolbachik zone in late November 2012 and terminated in the Tolud burst of seismicity. Based on the current ideas of the tectonics and magma sources for the Tolbachik volcanic zone, we discuss possible causes of these earthquakes.  相似文献   

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

14.
The characteristics and interplanetary excitation conditions of isolated bursts of Pi2 geomagnetic pulsations observed during the development of magnetospheric substorms (substorm Pi2) and in its absence (nonsubstorm Pi2) on the night side of the Earth are comparatively analyzed. It is shown that, regardless of the local time and season, the amplitude of isolated Pi2 substorm bursts is always higher than that of the nonsubstorm ones, and the periods and duration of the wave packets of substorm Pi2 bursts are less than those of nonsubstorms. Diurnal and seasonal variations in the characteristics of the two groups of Pi2 bursts differ in the form and position of maxima and minima. It is found that the start of excitation of isolated Pi2 bursts, during substorms and in its absence, is controlled by the preferred direction of the interplanetary magnetic field (IMF) vector perpendicular to the Sun–Earth line (angle θxB = arccos(Bx/B) → 90°). It is assumed that isolated Pi2 bursts of both groups are triggered by reorientation of the IMF vector in the ecliptic plane and the plane perpendicular to it ~15 min before their onset. The most likely source of midlatitude isolated Pi2 bursts during substorm development and in its absence are bursty bulk flows (BBFs) in the plasma sheet of the magnetospheric tail, the regularities of which coincide in many respects with the observed features of Pi2 bursts.  相似文献   

15.
The paper addresses estimation of the Hurst exponent for time series of the hourly values of the Dst index for the period from 1957 to 2011. It is found that the Hurst exponent is 0.79–0.94 for yearly intervals and 0.8–1.0 for monthly intervals. Based on R/S graphs, the Dst cycles are identified; they range from 3–4 months to 2.2 years and from 8.5 to 22 years in length. It is shown that a Dst time series can be quite satisfactorily described by an α-stable Levy process.  相似文献   

16.
The results of studying spread F obtained from the DPS-4 ionosonde data at the observatory of the Pushkov Institute of Terrestrial Magnetism, Ionosphere, and Radio Wave Propagation (Moscow) are presented. The methodical questions that arise during the study of a spread F phenomenon in the ionosphere are considered; the current results of terrestrial observations are compared with previously published data and the results of sounding onboard an Earth-satellite vehicle. The automated algorithm for estimation of the intensity of frequency spread F, which was developed by the authors and was successfully verified via comparison of the data of the digisonde DPS-4 and the results of manual processing, is described. The algorithm makes it possible to quantify the intensity of spread F in megahertz (the dFs parameter) and in the number of points (0, 1, 2, 3). The strongest spread (3 points) is shown to be most likely around midnight, while the weakest spread (0 points) is highly likely to occur during the daytime. The diurnal distribution of a 1–2 point spread F in the winter indicates the presence of additional maxima at 0300–0600 UT and 1400–1700 UT, which may appear due to the terminator. Despite the large volume of processed data, we can not definitively state that the appearance of spread F depends on the magnetic activity indices Kp, Dst, and AL, although the values of the dFs frequency spread interval strongly increased both at day and night during the magnetic storm of March 17–22, 2015, especially in the phase of storm recovery on March 20–22.  相似文献   

17.
A new 3D velocity model of the crust and upper mantle in the southeastern (SE) margin of the Tibetan plateau was obtained by joint inversion of body- and surface-wave data. For the body-wave data, we used 7190 events recorded by 102 stations in the SE margin of the Tibetan plateau. The surface-wave data consist of Rayleigh wave phase velocity dispersion curves obtained from ambient noise cross-correlation analysis recorded by a dense array in the SE margin of the Tibetan plateau. The joint inversion clearly improves the v S model because it is constrained by both data types. The results show that at around 10 km depth there are two low-velocity anomalies embedded within three high-velocity bodies along the Longmenshan fault system. These high-velocity bodies correspond well with the Precambrian massifs, and the two located to the northeast of 2013 M S 7.0 Lushan earthquake are associated with high fault slip areas during the 2008 Wenchuan earthquake. The aftershock gap between 2013 Lushan earthquake and 2008 Wenchuan earthquake is associated with low-velocity anomalies, which also acts as a barrier zone for ruptures of two earthquakes. Generally large earthquakes (M ≥ 5) in the region occurring from 2008 to 2015 are located around the high-velocity zones, indicating that they may act as asperities for these large earthquakes. Joint inversion results also clearly show that there exist low-velocity or weak zones in the mid-lower crust, which are not evenly distributed beneath the SE margin of Tibetan plateau.  相似文献   

18.
The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log10E2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log10E2/3, and log10N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log10E2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log10Ni and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log10Ni time series in the different magnitude ranges (M ≥ Mci) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔMi) in the range of low magnitudes are known. It is concluded that using the parameter log10N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.  相似文献   

19.
Analysis of the annual variation of the E-layer critical frequency median foE in the nighttime (22?02 LT) auroral zone by the data of several stations of the Northern Hemisphere has shown the median maximum in winter and minimum in summer, even though the summer contribution of solar radiation to foE is greater. Thus, a new phenomenon was discovered—an foE median winter anomaly in the nighttime auroral zone. Its amplitude (ratio of winter to summer foE figures) can reach 10–15%; however, this anomaly was weakly expressed and statistically insignificant at particular stations located in the auroral zone. The winter anomaly is more distinct for foE avr, the median of the E-layer critical frequency foE caused by the auroral source of atmospheric ionization, i.e., excluding the solar radiation contribution to foE. For foE avr, the amplitude of the winter anomaly can reach 15–20%. Based on the qualitative analysis, it has been found that foE winter anomaly is stipulated by the winter/summer asymmetry of energy flow of accelerated electrons, which induce discrete aurorae in the nighttime auroral zone.  相似文献   

20.
An algorithm is developed for automated detection of the short-period Pc1 geomagnetic pulsations (frequency band f = 0.2–3 Hz) from the continuous time series of digital recording during 1998–2014 at the midlatitude Borok station. A digital catalog with the indication of time intervals of the presence and main morphological characteristics of Pc1 pulsations is created. Based on this catalog, the annual, seasonal, and diurnal dynamics of the midlatitude Pc1 pulsation activity is studied for 1998–2014. It is shown that the annual variation of the Pc1 occurrence has a maximum in 2005, i.e., at the end of the solar cycle decay phase, just as in the previous cycles. It is found that the minimum of the cases of Pc1 occurrence is observed in 2009, i.e., not at the maximum, just was the case in the previous cycles, but during the deep minimum of solar activity, which testifies to the untypical conditions in the magnetosphere during the unusually long minimum of the 23rd cycle. The seasonal variation of the Pc1 occurrence has a summer minimum when the series of Pc1 pulsations occur almost thrice as rarely as in winter. Besides, there are relatively small maxima at equinox. The diurnal behavior of Pc1 pulsations has the maxima in the morning and midnight sectors of the magnetosphere. By the superposed epoch analysis technique it is established that the maximal number of the cases of occurrence of Pc1 pulsations at the Borok observatory is observed on the fourth day after the global geomagnetic disturbances. The statistical distributions of pulsations amplitude and duration are obtained.  相似文献   

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

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