首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
A mathematical model of the middle and high latitude ionosphere   总被引：5，自引：0，他引：5  
  相似文献   

2.
  相似文献   

3.
The Global Coupled Ionosphere–Thermosphere-Electrodynamics Model developed at Institute of Geology and Geophysics, Chinese Academy of Sciences (GCITEM-IGGCAS), is introduced in this paper. This new model self-consistently calculates the time-dependent three-dimensional (3-D) structures of the main thermospheric and ionospheric parameters in the height range from 90 to 600 km, including neutral number density of major species O2, N2, and O and minor species N(2D), N(4S), NO, He and Ar; ion number densities of O+ ,O2+, N2+, NO+, N+ and electron; neutral, electron and ion temperature; and neutral wind vectors. The mid- and low-latitude electric fields can also be self-consistently calculated. GCITEM-IGGCAS is a full 3-D code with 5° latitude by 7.5° longitude cells in a spherical geographical coordinate system, which bases on an altitude grid. We show two simulations in this paper: a March Equinox one and a June Solstice one, and compare their simulation results to MSIS00 and IRI2000 empirical model. GCITEM-IGGCAS can reproduce the main features of the thermosphere and ionosphere in both cases.  相似文献   

4.
Near-earth plasma parameters were calculated using a global numerical self-consistent and time-dependent model of the thermosphere, ionosphere and protonosphere (GSM TIP). The model results are compared with experimental data of different origin, mainly EISCAT measurements and simultaneous satellite data (Ne and ion composition). Model runs with varying inputs of auroral FAC distributions, temperature of vibrationally excited nitrogen and photoelectron energy escape fluxes are used to make adjustments to the observations. The satellite data are obtained onboard Active and its subsatellite Magion –2 when they passed nearby the EISCAT station around 0325 and 1540 UT on 31 July 1990 at a height of about 2000 and 2200 km, respectively. A strong geomagnetic disturbance was observed two days before the period under study. Numerical calculations were performed with consideration of vibrationally excited nitrogen molecules for high solar-activity conditions. The results show good agreement between the incoherent-scatter radar measurements (Ne, Te, Ti) and model calculations, taking into account the excited molecular nitrogen reaction rates. The comparison of model results of the thermospheric neutral wind shows finally a good agreement with the HWM93 empirical wind model.  相似文献   

5.
The period of interplanetary, geomagnetic and solar disturbances of September 7–15, 2005, is characterized by two sharp increases of solar wind velocity to 1000 km/s and great Dst variation of the geomagnetic field (~140 nT). The time variations of theoretical and experimental geomagnetic thresholds observed during this strong geomagnetic storm, their connection with solar wind parameters and the Dst index, and the features of latitudinal behavior of geomagnetic thresholds at particular times of the storm were studied. The theoretical geomagnetic thresholds were calculated with cosmic ray particle tracing in the magnetic field of the disturbed magnetosphere described by Ts01 model. The experimental geomagnetic thresholds were specified by spectrographic global survey according to the data of cosmic ray registration by the global station network.  相似文献   

6.
The three-dimensional nonstationary theoretical model of the concentrations and temperatures of electrons and ions in the ionospheric F region and plasmasphere at low and middle latitudes is used to study variations in the concentration NmF2 and height hmF2 of the ionospheric F2 layer under the action of the plasma zonal drift in the direction geomagnetic west-geomagnetic east perpendicularly to the electric E and geomagnetic B fields. The calculated and measured values of NmF2 and hmF2 for 16 ionospheric sounding stations during the quiet geomagnetic period on March 28–29, 1964 at low solar activity are compared. This comparison made it possible to correct the input parameters of the model: [O] from the NRLMSISE-00 model and the meridional component of the neutral wind velocity from the HWW90 model. It is shown that the nighttime NmF2 values decrease up to twice at low solar activity in the low-latitude ionosphere, and the hmF2 values change by up to 16 km, if the plasma zonal E×B drift is not taken into account. Under the daytime conditions, the influence of the plasma zonal E×B drift on NmF2 can be neglected.  相似文献   

7.
Variations of the upper boundary of the ionosphere (UBI) are investigated based on three sources of information: (i) ionosonde-derived parameters: critical frequency foF2, propagation factor M3000F2, and sub-peak thickness of the bottomside electron density profile; (ii) total electron content (TEC) observations from signals of the Global Positioning System (GPS) satellites; (iii) model electron densities of the International Reference Ionosphere (IRI*) extended towards the plasmasphere. The ionospheric slab thickness is calculated as ratio of TEC to the F2 layer peak electron density, NmF2, representing a measure of thickness of electron density profile in the bottomside and topside ionosphere eliminating the plasmaspheric slab thickness of GPS-TEC with the IRI* code. The ratio of slab thickness to the real thickness in the topside ionosphere is deduced making use of a similar ratio in the bottomside ionosphere with a weight Rw. Model weight Rw is represented as a superposition of the base-functions of local time, geomagnetic latitude, solar and magnetic activity. The time-space variations of domain of convergence of the ionosphere and plasmasphere differ from an average value of UBI at ∼1000 km over the earth. Analysis for quiet monthly average conditions and during the storms (September 2002, October–November 2003, November 2004) has shown shrinking UBI altitude at daytime to 400 km. The upper ionosphere height is increased by night with an ‘ionospheric tail’ which expands from 1000 km to more than 2000 km over the earth under quiet and disturbed space weather. These effects are interposed on a trend of increasing UBI height with solar activity when both the critical frequency foF2 and the peak height hmF2 are growing during the solar cycle.  相似文献   

8.
Numerical calculations of the thermospheric and ionospheric parameters above EISCAT are presented for quiet geomagnetic conditions in summer. The Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) was used. The numerical results were obtained both with a self-consistent calculation of the electric fields of magnetospheric and dynamo-action origin and with the magnetospheric electric fields only. It was found that the dynamo-electric field has some effect on the ionospheric convection pattern during quiet geomagnetic conditions. It has a marked effect mainly on the zonal neutral wind component above EISCAT (±20m/s at 140 km altitude). We have studied the effects of various field-aligned current (FAC) distributions on thermosphere/ionosphere parameters and we show that a qualitative agreement can be obtained with region-I and -II FAC zones at 75° and 65° geomagnetic latitude, respectively. The maximum FAC intensities have been assumed at 03–21 MLT for both regions with peak values of 2.5 × 10–7 Am–2 (region I) and 1.25 × 10–7 A m–2 (region II). These results are in agreement with statistical potential distribution and FAC models constructed by use of EISCAT data. The lack of decreased electron density in the night-time sector as observed by the EISCAT radar was found to be due to the spatial distribution of ionospheric convection resulting from electric fields of magnetospheric origin.  相似文献   

9.
A self-consistent method for day-time F2-region modelling was applied to the analysis of Millstone Hill incoherent scatter observations during the storm period of March 16-22, 1990. The method allows us to calculate in a self-consistent way neutral composition, temperature and meridional wind as well as the ionized species height distribution. Theoretically calculated Ne(h) profiles fit the observed daytime ones with great accuracy in the whole range of heights above 150 km for both quiet and disturbed days. The overall increase in Tex by 270 K from March 16 to March 22 reflects the increase of solar activity level during the period in question. A 30% decrease in [O] and a twofold increase in [A^] are calculated for the disturbed day of March 22 relative to quiet time prestorm conditions. Only a small reaction to the first geomagnetic disturbance on March 18 and the initial phase of the second storm on March 20 was found in [O] and [N2] variations. The meridional neutral wind inferred from plasma vertical drift clearly demonstrates the dependence on the geomagnetic activity level being more equatorward on disturbed days. Small positive F2-layer storm effects on March 18 and 20 are totally attributed to the decrease in the northward neutral wind but not to changes in neutral composition. A moderate (by a factor of 1.5) O/ N2 ratio decrease relative to the MSIS-83 model prediction is required to describe the observed NmF2 decrease on the most disturbed day of March 22, but virtually no change of this ratio is needed for March 21.  相似文献   

10.
The total electron content (TEC) is a key ionospheric parameter for various space weather applications. Over the last decade an extensive database of TEC measurements has become available from both space- and ground-based observations, and these measurements have established the general morphology of the global TEC distributions. In particular, the TOPEX TEC measurements have shown strong longitudinal variations of TEC in addition to the observed day-to-day variabilities. To better understand the observed TEC variations and to better guide its modeling, we have studied the sensitivity of quiet-time TEC to the following key atmospheric and ionospheric parameters: neutral density, neutral wind, plasma temperatures, plasmaspheric flux, and the O+–O collision frequency. These parameters are often only roughly known and can cause large uncertainties in model results. For this study, we have developed a numerical mid-latitude ionospheric model, which solves the momentum and continuity equations for the O+ density and a simplified set of equations for the H+ density. To obtain TEC, the calculated ion densities have been integrated from the bottom altitude (100 km) to the altitude of the TOPEX satellite (1336 km). Our study shows that during the day the neutral wind and the neutral composition have the most important effect on TEC. In particular, the zonal component of the neutral wind can have a large effect on TEC in the southern hemisphere where the magnetic declination angle is large. During the night, most of the above-mentioned parameters can play a significant role in the TEC morphology, except for the plasma temperature, which has only a small effect on TEC. Finally, the TEC varies roughly linearly with respect to all of the parameters except for the neutral wind.  相似文献   

11.
The F-region peak electron densities NmF2 measured during daytime quiet geomagnetic conditions at low solar activity on January 22, 2008, April 8, 1997, July 12, 1986, and October 26, 1995, are compared. Ionospheric parameters are measured by the ionosonde and incoherent scatter radar at Millstone Hill and calculated with the use of a 1D nonstationary ionosphere–plasmasphere model of number densities and temperatures of electrons and ions at middle geomagnetic latitudes. The formation of the semiannual anomaly of the midlatitudinal NmF2 under daytime quiet geomagnetic conditions at low solar activity is studied. The study shows that the semiannual NmF2 anomaly occurs due to the total impact of three main causes: seasonal variations in the velocity of plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity; seasonal variations in the composition and temperature of the neutral atmosphere; and the dependence of the solar zenith angle on a number of the day in the year at the same solar local time.  相似文献   

12.
The F2-layer peak density, NmF2, and peak altitude, hmF2, which were observed by 12 ionospheric sounders during the 20 September 1964 geomagnetically quiet time period at solar minimum are compared with those calculated by the three-dimensional time-dependent theoretical model of the Earth's low and middle latitude ionosphere and plasmasphere. The modeled NmF2 are also compared with those measured during the geomagnetically quiet time periods of 12–15, 18–21, and 26 September 1964 to take into account observed day-to-day ionospheric variability. Major features of the data are reproduced by the model if the corrected HWM90 neutral wind is used. The changes in NmF2 due to the zonal E×B plasma drift are found to be less than 20% in the daytime low latitude ionosphere. The model, which does not take into account the zonal E×B plasma drift, underestimates night-time NmF2 up to the maximum factor of 2 at low geomagnetic latitudes. The night-time increase of NmF2 caused by the zonal E×B plasma drift is less pronounced at −20° and 20° geomagnetic latitudes in comparison with that between −10° and 10° geomagnetic latitude. The longitude dependence of the calculated night-time low latitude influence of the zonal E×B plasma drift on NmF2 is explained in terms of the longitudinal asymmetry in B (the eccentric magnetic dipole is displaced from the Earth's center and the Earth's eccentric tilted magnetic dipole moment is inclined with respect to the Earth's rotational axis), and the variations of the wind induced plasma drift and the meridional E×B plasma drift in geomagnetic longitude. The difference between the hmF2 values calculated by including the effect of zonal E×B drift and that obtained when it is excluded does not exceed 19 km in the low latitude ionosphere. Over the geomagnetic equator the zonal E×B plasma drift produces the maximum increase in the electron density by a factor of 1.06–1.48 and 1.05–1.30 at 700 and 1000 km altitude, respectively, and this increase is not significant above about 1500 km. Changes in the vertical electron content, VEC, caused by the zonal E×B plasma do not exceed 16% during the day, while the value of the night-time VEC is increased up to a factor of 1.4 due to this drift. The maximum effects of the zonal E×B plasma drift on the night-time electron density derived from the model results corresponding to solar minimum and maximum are quite comparable.  相似文献   

13.
The Atmospheric Ultraviolet Radiance Integrated Code (AURIC) is a software package developed by Computational Physics, Inc. (CPI) under the sponsorship of the Air Force Phillips Laboratory/Geophysics Directorate (PL/GP) (currently the Air Force Research Laboratory) for middle and upper atmospheric radiance modeling from the far ultraviolet to the near infrared. The AURIC has been considered as a general model for the radiative transfer simulation of airglow. Based on the theory of MODerate resolution atmospheric TRANsmission (MODTRAN), the AURIC extends calculation to altitudes above 100 km and the wavelength down to 80 nm. A package of AURIC v1.2 was released in 2002, which can be used for single-point simulation from 1947 to 1999. It means that the model is not suitable for atmospheric simulation of large datasets or for atmospheric parameters retrieval from amount of satellite measurements. In this paper, AURIC v1.2 is upgraded to AURIC-2012 based on MATLAB with improvements for modules of the Geomagnetic Parameter (GEOPARM), Atmosphere neutral composition (ATOMS), and Ionospheric electron density (IONOS). The improved AURIC can be used for global automatic airglow simulation and also for automatic retrieval of atmospheric compositions from satellite global observations, such as O/N2 and electron density etc. Besides, the model supplies possibilities for further improvement of airglow radiative mechanism and for substitution of other modules. Based on the AURIC-2012, Limb Column Emission Intensity (L-CEI) and Volume Emission Rate (VER) are calculated. For validation, the results were compared with measurements of the Global Ultraviolet Imager (GUVI) and TIMED Doppler Interferometer (TIDI), respectively. The averaged relative errors of L-CEI and VER at peak altitude are both within 20%. Finally, L-CEI varying with latitude, altitude, solar activity, and geomagnetic activity is simulated, and the distribution characteristics of the simulation and their influencing factors are analyzed subsequently.  相似文献   

14.
Formation mechanism of the spring–autumn asymmetry of the F2-layer peak electron number density of the midlatitudinal ionosphere, NmF2, under daytime quiet geomagnetic conditions at low solar activity are studied. We used the ionospheric parameters measured by the ionosonde and incoherent scatter radar at Millstone Hill on March 3, 2007, March 29, 2007, September 12, 2007, and September 18, 1984. The altitudinal profiles of the electron density and temperature were calculated for the studied conditions using a one-dimensional, nonstationary, ionosphere–plasmasphere theoretical model for middle geomagnetic latitudes. The study has shown that there are two main factors contributing to the formation of the observed spring–autumn asymmetry of NmF2: first, the spring–autumn variations of the plasma drift along the geomagnetic field due to the corresponding variations in the components of the neutral wind velocity, and, second, the difference between the composition of the neutral atmosphere under the spring and autumn conditions at the same values of the universal time and the ionospheric F2-layer peak altitude. The seasonal variations of the rate of O+(4S) ion production, which are associated with chemical reactions with the participation of the electronically excited ions of atomic oxygen, does not significantly affect the studied NmF2 asymmetry. The difference in the degree of influence of O+(4S) ion reactions with vibrationally excited N2 and O2 on NmF2 under spring and autumn conditions does not significantly change the spring–autumn asymmetry of NmF2.  相似文献   

15.
The influence of atmospheric solar radiation absorption on the photodetachment, dissociative photodetachment, and photodissociation rate coefficients (photodestruction rate coefficients) of O?, Cl?, O2 ?, O3 ?, OH?, NO2 ?, NO3 ?, O4 ?, OH?(H2O), CO3 ?, CO4 ?, ONOO?, HCO3 ?, CO3 ?(H2O), NO3 ?(H2O), O2 +(H2O), O4 +, N4 +, NO+(H2O), NO+(H2O)2, H+(H2O) n for n = 2–4, NO+(N2), and NO+(CO2) at D-region altitudes of the ionosphere is studied. A numerical one-dimensional time-dependent neutral atmospheric composition model has been developed to estimate this influence. The model simulations are carried out for the geomagnetically quiet time period of 15 October 1998 at moderate solar activity over the Boulder ozonesonde. If the solar zenith angle is not more than 90° then the strongest influence of atmospheric solar radiation absorption on photodestruction of ions is found for photodissociation of CO4 ? ions when CO3 ? ions are formed. It follows from the calculations that decreases in the photodestruction rate coefficients of ions under consideration caused by this influence are less than 2 % at 70 km altitude and above this altitude if the solar zenith angle does not exceed 90°.  相似文献   

16.
It has been clearly established that there is a substantial outflow of ionospheric plasma from the Earth's ionosphere in both the northern and southern polar regions. The outflow consists of both light thermal ions (H+ and He+) and an array of energized ions (NO+, O2+, N2+, O+, N+, He+, and H+). If the outflow is driven by thermal pressure gradients in the ionosphere, the outflow is called the “classical” polar wind. On the other hand, if the outflow is driven by energization processes either in the auroral oval or at high altitudes in the polar cap, the outflow is called the “generalized” polar wind. In both cases, the field-aligned outflow occurs in conjunction with magnetospheric convection, which causes the plasma to drift into and out of the sunlit hemisphere, cusp, polar cap, nocturnal auroral oval, and main trough. Because the field-aligned and horizontal motion are both important, three-dimensional (3-D) time-dependent models of the ionosphere–polar wind system are needed to properly describe the flow. Also, as the plasma executes field-aligned and horizontal motion, charge exchange reactions of H+ and O+ with the background neutrals (H and O) act to produce low-energy neutrals that flow in all directions (the neutral polar wind). This review presents recent simulations of the “global” ionosphere–polar wind system, including the classical, generalized, and neutral polar winds. The emphasis is on displaying the 3-D and dynamical character of the polar wind.  相似文献   

17.
Low frequency stochastic variations of the geomagnetic AE-index characterized by 1/fb-like power spectrum (where f is a frequency) are studied. Based on the analysis of experimental data we show that the Bz-component of IMF, velocity of solar wind plasma, and the coupling function of Akasofu are insufficient factors to explain these behaviors of the AE-index together with the 1/fb fluctuations of geomagnetic intensity. The effect of self-organized criticality (SOC) is proposed as an internal mechanism to generate 1/fb fluctuations in the magnetosphere. It is suggested that localized spatially current instabilities, developing in the magnetospheric tail at the initial substorm phase can be considered as SOC avalanches or dynamic clusters, superposition of which leads to the 1/fb fluctuations of macroscopic characteristics in the system. Using the sandpile model of SOC, we undertake numerical modeling of space-localized and global disturbances of magnetospheric current layer. Qualitative conformity between the disturbed dynamics of self-organized critical state of the model and the main phases of real magnetospheric substorm development is demonstrated. It is also shown that power spectrum of sandpile model fluctuations controlled by real solar wind parameters reproduces all distinctive spectral features of the AE fluctuations.  相似文献   

18.
The spatial distributions of electric fields and currents in the Earth’s atmosphere are calculated. Electric potential distributions typical of substorms and quiet geomagnetic conditions are specified in the ionosphere. The Earth is treated as a perfect conductor. The atmosphere is considered as a spherical layer with a given height dependence of electrical conductivity. With the chosen conductivity model and an ionospheric potential of 300 kV with respect to the Earth, the electric field near the ground is vertical and reaches 110 Vm−1. With the 60-kV potential difference in the polar cap of the ionosphere, the electric field disturbances with a vertical component of up to 13 V m−1 can occur in the atmosphere. These disturbances are maximal near the ground. If the horizontal scales of field nonuniformity are over 100 km, the vertical component of the electric field near the ground can be calculated with the one-dimensional model. The field and current distributions in the upper atmosphere can be obtained only from the three-dimensional model. The numerical method for solving electrical conductivity problems makes it possible to take into account conductivity inhomogeneities and the ground relief.  相似文献   

19.
The time variations in the CR geomagnetic cutoff rigidity and their relation to the interplanetary parameters and the Dst index during a strong magnetic storm of November 18–24, 2003, have been analyzed. The Tsyganenko (Ts03) model of a strongly disturbed magnetosphere [Tsyganenko, 2002a, 2002b; Tsyganenko et al., 2003] have been used to calculate effective geomagnetic thresholds with the help of the method for tracing CR particle trajectories in the magnetospheric magnetic field. The geomagnetic thresholds have been calculated using the method of global spectrographic survey (GSS), based on the data from the global network of CR stations, and the results have been compared with the effective geomagnetic cutoff rigidities. The daily anisotropy of effective geomagnetic thresholds during the Dst variation minimum have been estimated. The relation of the theoretical and experimental geomagnetic thresholds, obtained using the GSS method, to the interplanetary parameters and Dst variation is analyzed. The Dst variations, IMF B z , and solar wind density are most clearly defined in the geomagnetic thresholds during this storm. The correlation between B y and experimental geomagnetic thresholds is higher than such a correlation between this parameter and theoretical thresholds by a factor 2–3, which suggests that a real dawn-dusk asymmetry during this storm was stronger than such an asymmetry represented by the Ts03 model.  相似文献   

20.
This study compares the OV1-10 satellite measurements of the integral airglow intensities at 630 nm in the SAR arc regions observed in the northern and southern hemisphere as a conjugate phenomenon, with the model results obtained using the time-dependent one-dimensional mathematical model of the Earth ionosphere and plasmasphere (the IZMIRAN model) during the geomagnetic storm of the period 15–17 February 1967. The major enhancements to the IZMIRAN model developed in this study are the inclusion of He+ ions (three major ions: O+ H+ and He+ and three ion temperatures), the updated photochemistry and energy balance equations for ions and electrons, the diffusion of NO+ and O+2 ions and O(1D) and the revised electron cooling rates arising from their collisions with unexcited N2, O2 molecules and N2 molecules at the first vibrational level. The updated model includes the option to use the models of the Boltzmann or non-Boltzmann distributions of vibrationally excited molecular nitrogen. Deviations from the Boltzmann distribution for the first five vibrational levels of N2 were calculated. The calculated distribution is highly non-Boltzmann at vibrational levels v > 2 and leads to a decrease in the calculated electron density and integral intensity at 630 nm in the northern and southern hemispheres in comparison with the electron density and integral intensity calculated using the Boltzmann vibrational distribution of N2. It is found that the intensity at 630 nm is very sensitive to the oxygen number densities. Good agreement between the modeled and measured intensities is obtained provided that at all altitudes of the southern hemisphere a reduction of about factor 1.35 in MSIS-86 atomic oxygen densities is included in the IZMIRAN model with the non-Boltzm-ann vibrational distribution of N2. The effect of using of the O(1D) diffusion results in the decrease of 4–6% in the calculated integral intensity of the northern hemisphere and 7–13% in the calculated integral intensity of the southern hemisphere. It is found that the modeled intensities of the southern hemisphere are more sensitive to the assumed values of the rate coefficients of O+(4S) ions with vibrationally excited nitrogen molecules and quenching of O+(2D) by atomic oxygen than the modeled intensities of the northern hemisphere.  相似文献   

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

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