Relationship between the daytime critical frequencies of the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer and thermospheric characteristics

The relationship between the critical frequency of the F ₂ layer and the atmospheric characteristics has been obtained in a general form. It has been shown that this relation makes it possible to sufficiently accurately describe the daytime values of foF2 while comparing them with the observed monthly median values. Such comparisons were performed, first, for the data of measurements in Irkutsk using the DPS-4 digital ionosonde in 2003–2006 and, second, based on the annual variations in the noon foF2 values at 24 stations of the Northern Hemisphere in 1984. The calculations were performed using the MSIS-86 thermospheric model [Hedin, 1987].     

Equatorial night-time F-region zonal electric fields

Night-time F-region vertical electrodynamic drifts at the magnetic equatorial station, Trivandrum are obtained for a period of 2 years, 1989 and 1990 (corresponding to solar cycle maximum epoch), using ionosonde h^‘F data. The seasonal variation of the vertical drift is found to be associated with the longitudinal gradients of the thermospheric zonal wind. Further, the seasonal variation of the prereversal enhancement of the vertical drift is associated with the time difference between the sunset times of the conjugate E-regions (magnetic field line linked to F-region) which is indicative of the longitudinal gradients of the conductivity (of the E-region). The vertical drifts and the causative zonal electric fields at Trivandrum are compared with those at Jicamarca and F-region zonal electric field models. It is seen that the night-time downward drift (as also the causative westward electric field) at Jicamarca is greater than that at Trivandrum. The prereversal enhancement of the drift is greater at Jicamarca than at Trivandrum during the summer and the equinoxes, whereas during the winter the opposite is the case.     

Day-to-day thermosphere parameter variation as deduced from Millstone Hill incoherent scatter radar observations during March 16-22, 1990 magnetic storm period

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 N_e(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 T_ex 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 [N₂] 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/ N₂ ratio decrease relative to the MSIS-83 model prediction is required to describe the observed N_mF2 decrease on the most disturbed day of March 22, but virtually no change of this ratio is needed for March 21.     

A role of the altitude gradient of the thermospheric wind velocity in ionospheric <Emphasis Type="Italic">F</Emphasis>-region dynamics

By processing the data of vertical ionospheric sounding in Almaty for 2000–2009, we obtained the distributions of the heights of the maximum (h _m F) and bottom (h _bot F) of the F2-layer, incremental changes in its semi-thickness (δh), the characteristic time of losses (τ), and the vertical displacement velocity of the node of the thermospheric wind (V) during the transitional time of the day during nighttime increases in the electron concentration at the layer maximum. The comparison of the measured V and modeled V _m velocities showed a certain discrepancy. The influence of the altitude gradient of the meridional thermospheric wind velocity on the behaviors of h _m F, h _bot F, δh, and τ during nighttime increases in the electron concentration is studied.     

The height of the maximum ionospheric electron density over the MU radar

Ionospheric F₂-layer peak height h_mF₂ variations, as measured over 1986–1995 by the MU radar (34.85°N, 136.1°E) and as calculated with a theoretical model, are discussed. The diurnal variations of the measured peak height for different seasons and levels of solar activity are compared with those estimated from ionosonde M3000F₂ and IRI predictions. Also given are the measured ion drift velocities and meridional neutral winds needed to understand the dynamic behavior of the F₂-layer. It is found that: (1) h_mF₂ is generally higher during periods of the solar maximum than during periods of the solar minimum, and higher in summer than in winter; (2) for the solar maximum, h_mF₂ drops markedly in the morning and in the afternoon, while, for the solar minimum, the h_mF₂ minimum occurs in the morning during summer and usually in the afternoon during winter. In general, the measured h_mF₂ is well reproduced by our model when we use the observed drift velocities and plasma temperatures as inputs. Our modeling study shows that the neutral wind contributes strongly to the diurnal variation of h_mF₂ in winter by lowering the ionization layer by day, particularly for the solar maximum; it also helps to enlarge the day–night difference of h_mF₂ in summer. The northward electromagnetic drifts that usually cancel the neutral wind effect have only a minor effect for the location of the MU radar. Other features of the observed h_mF₂ variations, e.g., the solar maximum–minimum difference, the summer–winter difference, and the morning and afternoon drops, are explained by the basic processes of O⁺ production, loss and diffusion, as influenced by the atomic oxygen concentration and neutral and plasma temperatures.     

Gradient drift instability as a source of spread <Emphasis Type="Italic">F</Emphasis> in the region of large-scale irregularities in the low-latitude topside ionosphere

The expression for the increment of instability and decrement of diffusion damping of gradient drift waves for ionospheric altitudes above the F ₂ layer maximum is obtained. The gradient drift instability is used to interpret the observations of spread F in the region of large-scale horizontal irregularities of the electron density. Two types of such irregularities observed on board the Intercosmos-19 (IC-19) satellite in the region of low latitudes (a peak of the density in the dusk ionosphere and a trough of the density in the dawn ionosphere) are considered. It is shown that the observed gradients of the density and electric field values in the dawn and dusk ionospheric sectors are quite sufficient for the instability development criterion to be satisfied in both considered cases.     

Equatorial ionization anomaly and thermospheric meridional winds during two major storms over Brazilian low latitudes

The equatorial ionosphere responses over Brazil to two intense magnetic storms that occurred during 2001 are investigated. The equatorial ionization anomaly (EIA) and variations in the zonal electric field and meridional winds at different storms phases are studied using data collected by digisondes and GPS receivers. The difference between the F layer peak density (foF2) at an equatorial and a low latitude sites was used to quantify the EIA; while the difference between the true heights (h_F) at the equatorial and an off-equatorial site was used to calculate the magnetic meridional winds. The vertical drift was calculated as dh_F/dt. The results show prompt penetration electric fields causing unusual early morning development of the EIA, and disturbed dynamo electric field producing significant modification in the F region parameters. Variations to different degrees in the vertical drift, the thermospheric meridional winds and the EIA developments were observed depending on the storm phases.     

Daytime F2-layer positive storm effect at middle and lower latitudes

Daytime F2-layer positive storm effects at middle and lower latitudes in the winter thermosphere are analyzed using AE-C, ESRO-4 neutral gas composition data, ground-based ionosonde observations and model calculations. Different longitudinal sectors marked by the storm onset as ‘night-time’ and ‘daytime’ demonstrate different F2-layer positive storm mechanisms. Neutral composition changes in the ‘night-time’ sector with increased [O] and [N₂] absolute concentrations, while (N₂/O)_storm/(N₂/O)_quiet\approx1 at F2-layer heights, are shown to contribute largely to the background N_mF2 increase at lower latitudes lasting during daytime hours. Storm-induced surges of the equatorward wind give rise to an additional N_mF2 increase above this background level. The mid-latitude F2-layer positive storm effect in the ‘daytime’ sector is due to the vertical plasma drift increase, resulting from the interaction of background (poleward) and storm-induced (equatorward) thermospheric winds, but not to changes of [O] and [N₂] concentrations.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ 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 ₂ 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.     

Correlation of vertical electrical sounding and electrical borehole log data for groundwater exploration

We have correlated the longitudinal unit conductance CL obtained from interpreted vertical electrical sounding data with the formation resistivity R_t and the formation resistivity factor F, obtained by carrying out electrical borehole logging. Interpreted geophysical data of eleven soundings and two electrical borehole log records are used for the analysis. The geophysical data used were acquired in a sedimentary basin. The study area is called Lower Maner Basin located in the province of Andhra Pradesh, India. Vertical electrical soundings were carried out using a Schlumberger configuration with half current electrode separation varying from 600–1000 m. For logging the two boreholes, a Widco logger‐model 3200 PLS was used. True formation resistivity R_t was calculated from a resistivity log. Formation resistivity factor F was also calculated at various depths using R_t values. An appreciable inverse relation exists between the correlated parameters. The borehole resistivity R_t and the formation resistivity factor F decrease with the increase in the longitudinal unit conductance CL. We have shown the use of such a relation in computing borehole resistivity R_t and formation resistivity factor F at sites that posses only vertical electrical sounding data, with a fair degree of accuracy. Validation of the correlation is satisfactory. Scope for updating the correlation is discussed. Significance and applications of the relation for exploration of groundwater, namely to update the vertical electrical sounding data interpretation by translating the vertical electrical sounding data into electrical borehole log parameters, to facilitate correlations studies and to estimate the porosity (φ), permeability (K) and water saturation S_w of water bearing zones are discussed.     

Using observations from the GPS and TOPEX satellites to investigate night-time TEC enhancements at mid-latitudes in the southern hemisphere during a low sunspot number period

The state of the ionization of the upper atmosphere at low and mid latitudes in the Australian region has been studied by investigating the total electron content (TEC) obtained by a dual-frequency group path and phase path GPS technique. For the low sunspot number time period of March 1995–February 1996, one week of data centred on the Priority Regular World Day for each month have been used to investigate night-time mid-latitude peaks occurring around midnight in the Australian region. TEC from TOPEX provided additional information related to the formation of the night-time peaks. Although night-time TEC enhancements have been observed previously, there is no general agreement on their origin. From the results of the present study, the development of midnight TEC enhancements coincided with the low latitude processes occurring at around the time of vertical E×B drift velocity reversal. The TOPEX results confirmed that the upward E×B drift velocity reversal and the downward plasma flow from greater heights producing the night-time peaks at mid latitudes are triggered from a common source: the westward electric field.     

Quantitative models of the fallout and dispersal of tephra from volcanic eruption columns

A theoretical model of clast fallout from convective eruption columns has been developed which quantifies how the maximum clast size dispersal is determined by column height and wind strength. An eruption column consists of a buoyant convecting region which rises to a heightH _B where the column density equals that of the atmosphere. AboveH _B the column rises further to a heightH _T due to excess momentum. BetweenH _T andH _B the column is forced laterally into the atmosphere to form an upper umbrella region. Within the eruption column, the vertical and horizontal velocity fields can be calculated from exprimental and theoretical studies and consideration of mass continuity. The centreline vertical velocity falls as a nearly linear function over most of the column's height and the velocity decreases as a gaussian function radially away from the centreline. Both column height and vertical velocity are strong functions of magma discharge rate. From calculations of the velocity field and the terminal fall velocity of clasts, a series of particle support envelopes has been constructed which represents positions where the column vertical velocity and terminal velocity are equal for a clast of specific size and density. The maximum range of a clast is determined in the absence of wind by the maximum width of the clast support envelope.The trajectories of clasts leaving their relevant support envelope at its maximum width have been modelled in columns from 6 to 43 km high with no wind and in a wind field. From these calculations the shapes and areas of maximum grain size contours of the air-fall deposit have been predicted. For the no wind case the theoretical isopleths show good agreement with the Fogo A plinian deposit in the Azores. A diagram has been constructed which plots, for a particular clast size, the maximum range normal to the dispersal axis against the downward range. From the diagram the column height (and hence magma discharge rate) and wind velocity can be determined. Historic plinian eruptions of Santa Maria (1902) and Mount St. Helens (1980) give maximum heights of 34 and 19 km respectively and maximum wind speeds at the tropopause of m/s and 30 m/s respectively. Both estimates are in good agreement with observations. The model has been applied to a number of other plinian deposits, including the ultraplinian phase of theA.D. 180 Taupo eruption in New Zealand which had an estimated column height of 51 km and wind velocity of 27 m/s.     

The ionospheric F2-region at low geomagnetic latitudes during the geomagnetic storms of 22–26 April 1990: Comparison of observed and modeled response

A comparison between the modeled NmF2 and hmF2 and NmF2 and hmF2, which were observed by the Kokubunji, Okinawa, Manila, Vanimo, and Darwin ionospheric sounders and by the middle and upper (MU) atmosphere radar, have been used to study the time-dependent response of the low-latitude ionosphere to geomagnetic forcing during a time series of geomagnetic storms from 22 to 26 April 1990. The reasonable agreement between the model results and data requires the modified equatorial meridional E×B plasma drift, the modified HWM90 wind, and the modified NRLMSISE-00 neutral densities. We found that changes in a flux of plasma into the nighttime equatorial F2-region from higher L-shells to lower L-shells caused by the meridional component of the E×B plasma drift lead to enhancements in NmF2 close to the geomagnetic equator. The equatorward wind-induced plasma drift along magnetic field lines, which cross the Earth equatorward of about 20° geomagnetic latitude in the northern hemisphere and about −19° geomagnetic latitude in the southern hemisphere, contributes to the maintenance of the F2-layer close to the geomagnetic equator. The nighttime weakening of the equatorial zonal electric field (in comparison with that produced by the empirical model of Fejer and Scherliess [Fejer, B.G., Scherliess, L., 1997. Empirical models of storm time equatorial zonal electric fields. J. Geophys. Res. 102, 24047–24056] or Scherliess and Fejer [Scherliess, L., Fejer, B.G., 1999. Radar and satellite global equatorial F region vertical drift model. J. Geophys. Res. 104, 6829–6842) in combination with corrected equatorward nighttime wind-induced plasma drift along magnetic field lines in the both geomagnetic hemispheres are found to be the physical mechanism of the nighttime NmF2 enhancement formation close to the geomagnetic equator over Manila during 22–26 April 1990. The model crest-to-trough ratios of the equatorial anomaly are used to study the relative role of the main mechanisms of the equatorial anomaly suppression for the 22–26 April 1990 geomagnetic storms. During the most part of the studied time period, a total contribution from geomagnetic storm disturbances in the neutral temperature and densities to the equatorial anomaly changes is less than that from meridional neutral winds and variations in the E×B plasma drift. It is shown that the latitudinal positions of the crests are determined by the E×B drift velocity and the neutral wind velocity.     

Converted-wave imaging in anisotropic media: theory and case studies

Common‐conversion‐point binning associated with converted‐wave (C‐wave) processing complicates the task of parameter estimation, especially in anisotropic media. To overcome this problem, we derive new expressions for converted‐wave prestack time migration (PSTM) in anisotropic media and illustrate their applications using both 2D and 3D data examples. The converted‐wave kinematic response in inhomogeneous media with vertical transverse isotropy is separated into two parts: the response in horizontally layered vertical transverse isotrophy media and the response from a point‐scatterer. The former controls the stacking process and the latter controls the process of PSTM. The C‐wave traveltime in horizontally layered vertical transverse isotrophy media is determined by four parameters: the C‐wave stacking velocity V_C2, the vertical and effective velocity ratios γ₀ and γ_eff, and the C‐wave anisotropic parameter χ_eff. These four parameters are referred to as the C‐wave stacking velocity model. In contrast, the C‐wave diffraction time from a point‐scatterer is determined by five parameters: γ₀, V_P2, V_S2, η_eff and ζ_eff, where η_eff and ζ_eff are, respectively, the P‐ and S‐wave anisotropic parameters, and V_P2 and V_S2 are the corresponding stacking velocities. V_P2, V_S2, η_eff and ζ_eff are referred to as the C‐wave PSTM velocity model. There is a one‐to‐one analytical link between the stacking velocity model and the PSTM velocity model. There is also a simple analytical link between the C‐wave stacking velocities V_C2 and the migration velocity V_Cmig, which is in turn linked to V_P2 and V_S2. Based on the above, we have developed an interactive processing scheme to build the stacking and PSTM velocity models and to perform 2D and 3D C‐wave anisotropic PSTM. Real data applications show that the PSTM scheme substantially improves the quality of C‐wave imaging compared with the dip‐moveout scheme, and these improvements have been confirmed by drilling.     

Numerical simulation of the <Emphasis Type="Italic">F</Emphasis>2-layer stratification and appearance of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere: The morphology of the phenomena

This work is devoted to a numerical simulation of the equatorial ionosphere, performed using the GSM TIP model completed with a new block for calculating the electric field. It has been indicated that the usage of the wind system calculated according to the MSIS-90 model makes it possible to reproduce the electromagnetic drift velocities at the equator, the effect of the F2-layer stratification, and the appearance of the F3 layer in the equatorial ionosphere. The calculations performed using the modified GSM TIP model made it possible to detect a maximum in the electron density vertical profile at an altitude of ∼1000 km, formed by H⁺ ions, which we called the G layer. If this layer actually exists, it can be observed during sounding the low-latitude ionosphere from satellites during dark time of day.     

On the propagation of an in-plane shear fault with super-S-wave velocity

The speedv, especially the problem whether super S-wave velocity in the classical model (linear elasticity fracture mechanics) exists, of spontaneous propagation of a shear fault is investigated theoretically. An in-plane shear crack propagating in the crack plane is taken as the model of the shear fault. The results obtained firstly by Kostrov (1975) is extended from sub-Rayleigh wave velocity to super S-wave velocity, and the analytical expression for the stress intensity factorK ₂ in the case ofα>v>β is derived. It is proved that for Poisson mediumK ₂ is positive and real in the velocity range (β, 1.70β). This demonstrates that (β, 1.70β) is the velocity range which fulfils the conditions for spontaneous crack propagation. The existence, convergence and positiveness or negativeness ofK ₂ forv in individual sections are examined, and it is found that for an in-plane shear crack: 1. There are three sections forv, i.e., [0.v _R], (β, 1.70β), andα, respectively, and 2. There are two physically reasonable sections forv, the first is [v _R, β], and the second is [1.70β, α]. These two forbidden sections behave as barriers to fault propagation. The analytical expressions derived in this paper are not only suitable to classical model, but also to the other derivative models (e. g., the slip-weakening model and the renomalization model etc.). The model considered in this paper is more realistic than the static model employed by previous authors. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica, 15, 9–14, 1993.     

Model simulation of the equatorial electrojet in the Peruvian and Philippine sectors

Between 100 and 120 km height at the Earth's magnetic equator, the equatorial electrojet (EEJ) flows as an enhanced eastward current in the daytime E region ionosphere, which can induce a magnetic perturbation on the ground. Calculating the difference between the horizontal components of magnetic perturbation (H) at magnetometers near the equator and about 6–9° away from the equator, ΔH, provides us with information about the strength of the EEJ. The NCAR Thermosphere–Ionosphere–Electrodynamics General Circulation Model (TIE-GCM) is capable of simulating the EEJ current and its magnetic perturbation on the ground. The simulated diurnal, seasonal (March equinox, June solstice, December solstice), and solar activity (F_10.7=80, 140 and 200 units) variations of ΔH in the Peruvian (76°W) and Philippine (121°E) sectors, and the relation of ΔH to the ionospheric vertical drift velocity, are presented in this paper. Results show the diurnal, seasonal and solar activity variations are captured well by the model. Agreements between simulated and observed magnitudes of ΔH and its linear relationship to vertical drift are improved by modifying the standard daytime E region photoionization in the TIE-GCM in order to better simulate observed E region electron densities.     

Modeling of the effect of the“apparent <Emphasis Type="Italic">F</Emphasis>1 layer” observed in the Radar-Progress experiment with the help of the software for automated processing of vertical sounding ionograms

The results of the modeling the anomalous behavior of the F ₁ layer characteristics are presented. The anomaly lies in the jump-like variations in the critical frequency of the layer f ₀ F1 and an abnormally large f ₀ F1 value for the considered season and time. The observed effect is simulated based on the propagation of a large-scale travelling ionospheric disturbance. The quantitative values of the disturbance parameter are obtained and some peculiarities of its propagation are analyzed. The investigation was carried out with the help of the software for automated processing of the vertical sounding ionograms. The experimental data were obtained during the Radar-Progress project.     

Basement structure and yalu river fault belt in Dandong region

Some geophysical surveying works in the northeast part of Dandong, such as shallow shock refracted wave, electrical prospecting, electrical sounding and wave velocity measuring, are introduced in this paper, and the dynamic parameters are calculated. The results show that the basement structure in surveying region is very complex, the overburden thickness of the quaternary period, velocity distribution and dynamic parameters are of regional characteristics. The depth of basement is deep in the north and shallow in the west, the difference between north and west region is about 5–10 m. The south part of Yalu river fault belt is composed of F₁, F₂, F₃, F₄, F₅ fault, their strike direction is NE, we can determine that the F₂ fault is the main one in Yalu river fault belt, and the south part of Yalu river fault belt has no activity since Holocene Epoch. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 282–288, 1993.     

Influence of the plasma zonal E × B drift on the electron concentration in the low-latitude ionospheric F region at the minimum of solar activity near the spring equinox

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.