The use of GPS arrays in detecting shock-acoustic waves generated during rocket launchings

This paper is concerned with the parameters of shock-acoustic waves (SAW) generated during rocket launchings. We have developed the interferometric method for determining SAW parameters (including angular characteristics of the wave vector, and the SAW phase velocity, as well as the direction towards the source) using GPS-arrays. Contrary to the conventional radio-probing techniques, the proposed method provides an estimate of SAW parameters without a priori information about the site and time of a rocket launching. The application of the method is illustrated by a case study of ionospheric effects from launchings of rockets PROTON, SOYUZ and SPACE SHUTTLE from Baikonur and Kennedy Space Center cosmodromes in 1998–2000. In spite of a difference of rocket characteristics, the ionospheric response for all launchings had the character of an N-wave corresponding to the form of a shock wave. The SAW period T is 270–360 s, and the amplitude exceeds the standard deviation of total electron content background fluctuations in this range of periods under quiet and moderate geomagnetic conditions by factors of 2–5 as a minimum. The angle of elevation of the SAW wave vector varies from 30° to 60°, and the SAW phase velocity (900–1200 m/s) approaches the sound velocity at heights of the ionospheric F-region maximum. The position of the SAW source, inferred by neglecting refraction corrections, corresponds to the segment of the rockets path at a distance no less than 200–900 km from the launch pad, and to the rocket flying altitude no less than 100 km. Our data are consistent with the existing view that SAW are generated during a nearly horizontal flight of the rocket with its engine in operation in the acceleration segment of the path at 100–130 km altitudes in the lower atmosphere.     

The Pn wave velocities and the relief of Moho in the Tibetan Plateau

Eleven PASSCAL broadband digital seismic stations were employed in the Tibetan Plateau by the Sino-US team from September, 1991 to June, 1992. Seven of them were distributed along the Qinghai-Tibet highway, others in Maqin and Yushu of Qinghai Province, Linzhi and Xigatze of Tibet. The data included 31 local earthquakes recorded by these stations from July, 1991 to January, 1992. Considering the characters of digital data, we identified the seismic phases carefully in several methods using SAC softwares (Seismic Analysis Code) in SUN workstation. We compared the seismic phases after converting the seismograms of the single stations to the seismic profiles; analyzed the first arrivals of P waves in the incident planes by rotating 3 component seismic records; identified the seismic phases from the particle motion pictures. The Pn apparent velocities were calculated in the distance range of 230–1200 km from Linzhi earthquakes, western Changtang earthquakes, Xitieshan and Gonghe earthquakes and the earthquakes in India. The results show that the Pn velocities change slightly in the Tibetan Plateau (8.0–8.1 km/s). These values near the velocities at the uppermost mantle of the normal continents. The Moho undulation in the Tibetan Plateau was also studied by using Pn data by the time-term method. The primary results indicate that the Moho beneath the Tibetan Plateau is flat, and its depths are 67–70 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 593–600, 1992.     

First in situ measurement of electric field fluctuations during strong spread F in the Indian zone

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip 14°N) along with other experiments, as a part of equatorial spread F (ESF) campaign, to study the nature of irregularities in electric field and electron density. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of vertical and horizontal electric field fluctuations are presented here. Scale sizes of electric field fluctuations were measured in the vertical direction only. Strong ESF irregularities were observed in three regions, viz., 160/190 km, 210/257 km and 290/330 km. Some of the valley region vertical electric field irregularities (at 165 km and 168 km), in the intermediate-scale size range, observed during this flight, show spectral peak at kilometer scales and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of a new type. Scale sizes of vertical electric field fluctuations showed a decrease with increasing altitude. The most prominent scales were of the order of a few kilometers around 170 km and a few hundred meters around 310 km. Spectra of intermediate-scale vertical electric field fluctuations below the base of the F region (210/257 km) showed a tendency to become slightly flatter (spectral index n = –2.1 ± 0.7) as compared to the valley region (n = –3.6 ± 0.8) and the region below the F peak (n = –2.8 ± 0.5). Correlation analysis of the electron density and vertical electric field fluctuations suggests the presence of a sheared flow of current in 160/330 km region.     

Shallow Seismic Velocity Structure of the Betic Cordillera (Southern Spain) from Modelling of Rayleigh Wave Dispersion

A detailed dispersion analysis of Rayleigh waves generated by local earthquakes and occasionally by blasts that occurred in southern Spain, was undertaken to obtain the shear-wave velocity structure of the region at shallow depth. Our database includes seismograms generated by 35 seismic events that were recorded by 15 single-component short-period stations from 1990 to 1995. All these events have focal depths less than 10 km and body-wave magnitudes between 3.0 and 4.0, and they were all recorded at distances between 40 and 300 km from the epicentre. We analysed a total of 90 source-station Rayleigh-wave paths. The collected data were processed by standard digital filtering techniques to obtain Rayleigh-wave group-velocity dispersion measurements. The path-averaged group velocities vary from 1.12 to 2.25 km/s within the 1.0-6.0 s period interval. Then, using a stochastic inversion approach we obtained 1-D shear-wave velocity–depth models across the study area, which were resolved to a depth of circa 5 km. The inverted shear-wave velocities range approximately between 1.0 and 3.8 km/s with a standard deviation range of 0.05–0.16 km/s, and show significant variations from region to region. These results were combined to produce 3-D images via volumetric modelling and data visualization. We present images that show different shear velocity patterns for the Betic Cordillera. Looking at the velocity distribution at various depths and at vertical sections, we discuss of the study area in terms of subsurface structure and S-wave velocity distribution (low velocity channels, basement depth, etc.) at very shallow depths (0–5 km). Our results characterize the region sufficiently and lead to a correlation of shear-wave velocity with the different geological units features.     

The phase velocities of Rayleigh waves and the lateral variation of lithospheric structure in Tibetan Plateau

According to a Sino-U. S. joint project, eleven broadband digital PASSCAL seismometers had been deployed inside the Tibetan Plateau, of which 7 stations were on the profile from Lhasa to Golmud and other 4 stations situated at Maxin, Yushu, Xigatze and Linzhi. Dispersions and phase velocities of the Rayleigh surface waves (10s–120s) were obtained on five paths distributed in the different blocks of Tibetan Plateau. Inversions of the S-wave velocity structures in Songpan-Ganzi block, Qiang-Tang block, Lhasa block and the faulted rift zone were obtained from the dispersion data. The results show that significant lateral variation of the S-wave velocity structures among the different blocks exists. The path from Wenquan to Xigatze (abbreviated as Wndo-Xiga) passes through the rift-zone of Yadong-Anduo. The phase velocities of Rayleigh waves from 10s to 100s on this path are significantly higher than that on other paths. The calculated mean crustal velocity on this path is 3.8 km/s, much greater than that on other paths, where mean crustal velocities of 3.4–3.5 km/s are usually observed. Low velocity zones with different thicknesses and velocities are observed in the middle-lower crust for different paths. Songpan-Ganzi block, located in the northern part of Tibetan Plateau is characterized by a thinner crust of 65 km thick and a prominent low velocity zone in the upper mantle. The low velocity zone with a velocity of 4.2 km/s is located at a depth form 115 km to 175 km. While in other blocks, no low velocity zone in the upper mantle is observed. The value of Sn in Songpan-Ganzi is calculated to be 4.5 km/s, while those in Qiang-Tang and Lhasa blocks are about 4.6 km/s. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 566–573, 1992.     

Geomagnetic field effects of the Chelyabinsk meteoroid

An analysis was conducted of time variations in geomagnetic field components on the day of the Chelyabinsk meteorite event (February 15, 2013) and on control days (February 12 and 16, 2013). The analysis uses the data collected by magnetic observatories in Novosibirsk, Almaty, Kyiv, and Lviv. The distance R from the explosion site to the observatories varies in the range 1200–2700 km. The flyby and explosion of the Chelyabinsk cosmic body is found to have been accompanied by variations mainly in the horizontal component of the geomagnetic field. The variations are quasi-periodic with a period of 30–40 min, an amplitude of 0.5–2 nT for R ≈ 2700?1200 km, respectively, and a duration of 2–3 h. The horizontal velocity of the geomagnetic field disturbances is close to 260–370 m/s. A theoretical model of wave disturbances is proposed. According to the model, wave disturbances in the geomagnetic field are caused (a) by the motion of the gravity wave generated in the atmosphere by the falling space body and (b) by traveling ionospheric disturbances, which modulate the ionospheric current at dynamo altitudes. The calculated amplitudes of the wave disturbances are 0.6–1.8 nT for R ≈ 2700?1200 km, respectively. The estimates are in good agreement with the observational data. Disturbances in the geomagnetic field level (geomagnetic pulsations) in the period range 1–1000 s are negligible (less than 1 nT).     

Reversed polarity patches at the CMB and geomagnetic field reversal

The International Geomagnetic Reference Field models (IGRF) for 1900–2000 are used to calculate the geomagnetic field distribution in the Earth’ interior from the ground surface to the core-mantle boundary (CMB) under the assumption of insulated mantle. Four reversed polarity patches, as one of the most important features of the CMB field, are revealed. Two patches with +Z polarity (downward) at the southern African and the southern American regions stand out against the background of ™Z polarity (upward) in the southern hemisphere, and two patches of ™Z polarity at the North Polar and the northern Pacific regions stand out against the +Z background in the northern hemisphere. During the 1900–2000 period the southern African (SAF) patch has quickly drifted westward at a speed of 0.20–.3° /a; meanwhile its area has expanded 5 times, and the magnetic flux crossing the area has intensified 30 times. On the other hand, other three patches show little if any change during this 100-year period. Extending upward, each of the reversed polarity patches at the CMB forms a chimney-shaped “reversed polarity column” in the mantle with the bottom at the CMB. The height of the SAF column has grown rapidly from 200km in 1900 to 900km in 2000. If the column grows steadily at the same rate in the future, its top will reach to the ground surface in 600–700 years. And then a reversed polarity patch will be observed at the Earth’s surface, which will be an indicator of the beginning of a magnetic field reversal. On the basis of this study, one can describe the process of a geomagnetic polarity reversal, the polarity reversal may be observed firstly in one or several local regions; then the areas of these regions expand, and at the same time, other new reversed polarity regions may appear. Thus several poles may exist during a polarity reversal.     

Analytical comparison of geomagnetic total field between Sino-US stations

Using observational data of geomagnetic total intensity from 13 stations in the Beijing-Tianjin region, 3 stations in the western Yunnan region of China, and 6 stations in California of U. S. A., the daily variations and their spectra of geomagnetic total intensity were analyzed and compared. The results show that the morphology, the range and spectrum of daily variations in geomagnetic total intensity are basically the same within the local extent of 100–200 km and are different in the large extent of 500 km. The latitude factor of the daily variation range of geomagnetic total intensity is about 1–2 nT/degree within the latitude extent of 25°–40°. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, 83–89, 1992. This work is supported by the State Seismological Bureau and the Chinese Joint Seismological Science Foundation, and U.S. Geological Survey.     

On the onset and meridional propagation of the ionospheric F2-region response to geomagnetic storms

The meridional propagation velocities of the ionospheric F2-region response to 268 geomagnetic storms are calculated. Ionospheric vertical sounding data of 1 h time resolution from several stations located in a longitude sector approximately centred along the great circle that contains both the geomagnetic poles and the geographic poles are used.Most meridional propagation velocities from high to low latitudes are less than 600 m/s. The smaller velocities are typical of global neutral meridional wind circulation and the larger are representative of traveling atmospheric disturbances.Simultaneous disturbances at several locations are more frequent during positive phases than during negative phases. Negative phase meridional propagation velocities associated with meridional neutral winds are less frequent in the southern hemisphere when compared with corresponding velocities observed in the northern hemisphere. This may be related to the fact that the distance between the geomagnetic pole and the equator is smaller in the northern hemisphere.Most negative phase onsets are within the 06–10 LT interval. For middle geomagnetic latitudes a “forbidden time interval” between 11 and 14 LT is present. The positive phase onsets show the “dusk effect”.     

Elastic-anisotropic properties of garnet granulites from the lower crust of the belomorian mobile belt: Results of experimental study

The elastic-anisotropic properties of fifteen samples of garnet granulites from xenoliths in the diatreme on Elovyi Island, in Kandalaksha Bay in the White Sea are investigated. The PT-conditions of the formation of these rocks correspond to a depth of 25–40 km, which is the current depth of the lower crust in the region. According to acoustopolariscopy, the elastic anisotropy of these deep rocks is relatively low. The average seismic P and S velocities corresponding to the in-situ PT-conditions are 6.5–6.7 km/s and 3.7–3.8 km/s, respectively. It is shown that the garnet content in the deep-seated rocks can be estimated from the seismic data. According to the obtained results, the formation of the structural and textural pattern of the lower crustal rocks in the region was poorly affected by the paleotectonic factors.     

Parameters of traveling ionospheric disturbances estimated from satellite beacon observations in low-earth orbit

Based on beacon sounding the ionosphere using coherent signals of low-orbiting navigation satellites, the following parameters of medium-scale traveling ionospheric disturbances have been determined: the effective thickness of the atmospheric waveguide and the altitude of its longitudinal axis, the horizontal spatial period, the velocity, the maximum amplitudes of disturbances, and the inclination of the disturbance phase front. It has been found that the horizontal spatial period of studied disturbances increases with increasing distance to an initiating impact location and its delay. In this case the minimum value of the horizontal spatial period always exceeds ∼150 km, and the disturbance amplitude increases with increasing horizontal period and distance from the initiating disturbance location and with decreasing altitude of the main ionization maximum. All disturbance parameters are independent of the initiating impact nature. It has been found that disturbances with relative amplitudes of 0.1–0.7 are often observed. Disturbances that are registered more often travel southward, and their absolute velocities are 7–60 m/s.     

CPO-induced seismic anisotropy in UHP eclogites

Ultrahigh-pressure (UHP) eclogites often show strong plastic deformation and anisotropy of seismic properties. We report in this paper the seismic velocity and anisotropy of eclogite calculated from the crystallographic preferred orientations (CPOs) of constituent minerals (garnet, omphacite, quartz and rutile) and single crystal elastic properties. We also compared the calculated results with the measured results in similar eclogites. Our results suggest that (1) Except that garnet is a seismically quasi-isotropic mineral, omphacite, quartz, coesite and rutile all have strong seismic anisotropies (AVp = 23.0%―40.9%, Max. AVs = 18.5%―47.1%). They are the major sources for anisotropy in eclogite. The average seismic velocities are fast in garnet and rutile, moderate in omphacite and coesite, and slow in quartz. (2) The deformed eclogites have the maximum Vp (8.33―8.75 km/s) approximately parallel to foliation and lineation, the minimum Vp (8.25―8.62 km/s) approximately normal to foliation and lineation and the Vp anisotropies of 1.0―1.7%. Their Vs are 4.93―4.97 km/s. The corresponding maximum anisotropies (0.73%―1.78%) of Vs are at 45° to both foliation and lineation and the minimum anisotropies at positions normal to lineation on the foliation plane. The Vs1 polarization planes are approximately parallel to foliation. The mean Vp and Vs of eclogite under UHP peak metamorphism conditions (P = 3―5 GPa, T = 900―1100℃) are estimated to be 3.4%―7.2% and 6.3%―12.1% higher than those at ambient pressure and temperature conditions, respectively. (3) Omphacite component dominates the anisotropy of eclogite while garnet component reduces the anisotropy and increases the seismic velocities. Quartz component has a small effect on the anisotropy but reduces the seismic velocities of eclogite. The effect of rutile component is negligible on seismic properties of eclogite due to its trivial volume fraction. (4) The increase of volume fraction of omphacite in eclogite will reduce the seismic velocities and increase the anisotropy. Omphacitite has seismic velocities reduced by 6%―8% and anisotropies increased to 3%―4% compared to those of garnetite. Our results suggest that the seismic properties calculated with single crystal elastic properties and CPOs are equivalent to those measured in laboratory. Moreover, it provides insights into the mineral physical interpretations of eclogite seismic properties.     

Inversion of seismic and gravity data for the composition and core sizes of the Moon

We model the internal structure of the Moon, initially homogeneous and later differentiated due to partial melting. The chemical composition and the internal structure of the Moon are retrieved by the Monte-Carlo inversion of the gravity (the mass and the moment of inertia), seismic (compressional and shear velocities), and petrological (balance equations) data. For the computation of phase equilibrium relations and physical properties, we have used a method of minimization of the Gibbs free energy combined with a Mie-Gr@uneisen equation of state within the CaO-FeO-MgO-Al₂O₃-SiO₂ system. The lunar models with a different degree of constraints on the solution are considered. For all models, the geophysically and geochemically permissible ranges of seismic velocities and concentrations in three mantle zones and the sizes of Fe-10%S core are estimated. The lunar mantle is chemically stratified; different mantle zones, where orthopyroxene is the dominant phase, have different concentrations of FeO, Al₂O₃, and CaO. The silicate portion of the Moon (crust + mantle) may contain 3.5–5.5% Al₂O₃ and 10.5–12.5% FeO. The chemical boundary between the middle and the lower mantle lies at a depth of 620–750 km. The lunar models with and without a chemical boundary at a depth of 250–300 km are both possible. The main parameters of the crust, the mantle, and the core of the Moon are estimated. At the depths of the lower mantle, the P and S velocities range from 7.88 to 8.10 km/s and from 4.40 to 4.55 km/s, respectively. The radius of a Fe-10%S core is 340 ± 30 km.     

Crustal Structure of the Southern Korean Peninsula from Seismic Waves Generated by Large Explosions in 2002 and 2004

In order to investigate the velocity structure of the southern part of the Korean peninsula, seismic refraction profiles were obtained along a 294-km WNW-ESE line and a 335-km NNW-SSE line in 2002 and 2004, respectively. Seismic waves were generated by detonating 500–1000 kg explosives in drill holes at depths of 80–150 m. The seismic signals were recorded by portable seismometers at nominal intervals of 1.5–1.7 km. Separate velocity tomograms were derived from first arrival times using a series expansion method of travel-time inversion. The raypaths indicate several mid-crust interfaces including those at approximate depths of 2–3, 15–17, and 22 km. The Moho discontinuity with refraction velocity of 7.8 to 8.4 km/s has a maximum depth of 37–39 km under the southern central portion of the peninsula. The Moho becomes shallower as the Yellow Sea and the East Sea are approached on the west and east coasts of the peninsula, respectively. The depth of the 7.6 km/s velocity contour varies from 29.4 km to 36.5 km. The discrepancy in depth between the seismological Moho and the interpreted critically refracting interface may result from the presence of a gradual transition between the crust and mantle. The velocity tomograms show particular crustal structures including (1) the existence of an over 70-km wide low-velocity zone centered at 6–7 km depth under the Okchon fold belt and Ryeongnam massif, (2) existence of high-velocity materials under the Gyeongsang basin, and (3) the downward extension of the Yeongdong fault to depths greater than 10 km.     

Investigation of the broadband ELF turbulence by observations of the FAST satellite

Scaling properties of variable electric fields in the topside ionosphere have been investigated on scales s from ∼30 m to 2 km by FAST electric field observations with sample rate of 512 s⁻¹, in sixteen events of the broadband ELF turbulence. It is shown that down to scales of a few hundred meters, the power of turbulent electric fluctuations is a power law, ∼s ^α. Scaling index α derived from the slope of logarithmic diagrams (LD) constructed by the discrete wavelet transform of data can be estimated as α = 2.2 ± 0.3, which is close to α estimate earlier reported for scales 1–30 km by electric field observations of the Dynamics Explorer 2 satellite. The behavior of α index is analyzed near the scale of the order of electron inertial length λ_e = c/ω₀ (ω₀ being the electron plasma frequency). At altitudes considered (700–2500 km), λ_e makes 100–900 m. We demonstrate that at scales ≤λ_e, a decrease of LD slope and deviation from the power law are typically observed. As pointed out in the discussion, this feature cannot be identified as a transition to the diffusion range, where dissipation of the turbulence occurs.     

Patterns of tidal flow asymmetry on shoreface-connected ridge topography off Spiekeroog Island,German Bight

Summary Asymmetry in spring-neap, flood-ebb, and acceleration-deceleration near-bottom tidal current parameters is assessed on the basis of fair-weather measurements at nineteen locations off Spiekeroog Island, German Bight. The study is aimed at a better understanding of flow-substrate interactions and, by implication, the development of shoreface-conneeted ridges which occur in 8–18 m water depth off the East Frisian coast. The ridges are spaced 1–1.5 km apart, have a relief of 3–5 m, widths of 1–2 km, and flank steepnesses of < 1‡. Flood-current peak-velocity asymmetry on each ridge was found to increase from the trough across the crest towards the seaward flank at both spring and neap tide. Tidal excursion was flood asymmetric at spring tide, but the asymmetry diminished with offshore distance. Flood and ebb accelerating currents had higher mean velocities than the decelerating counterparts at spring tide. The mean velocity of decelerating currents tended to be flood-asymmetric at neap tide. Asymmetry was also noted between flood and ebb current directions, on the one hand, and between spring tide and neap tide cross-ridge flow angle (0–44‡) of these currents, on the other. The observed asymmetries in tidal flow parameters are considered critical to ridge growth and maintenance, and most probably result from a combination of overtides as well as geometric and inertial effects due to the ridge topography.

---

Tidestromasymmetrien der Zungenrifftopographie vor Spiekeroog

Zusammenfassung Basierend auf Daten von 19 Me\stationen vor Spiekeroog wurden die Asymmetrien bodennaher Tidestromparameter bei Spring-und Nipptide, Ebbe und Flut sowie Tidebeschleunigung/-verz?gerung bei ruhigem Wetter untersucht. Ziel der Studie ist ein besseres Verst?ndnis der Wechselwirkung zwischen Str?mung und Boden und somit der Entwicklung der Zungenriffe vor der ostfriesischen Küste in 8 bis 18 m Tiefe. Die 1–1,5 km auseinanderliegenden Riffe sind 3–5 m hoch, 1–2 km breit und haben eine Neigung von < 1‡. Bei maximalem Flutstrom nahm die Asymmetrie auf den Zungenriffen von der Mulde über den Scheitel zum seew?rtigen Riffhang hin zu, und zwar bei Spring-und bei Nipptide. Die Gezeitenauslenkung war bei Springtide asymmetrisch, jedoch nahm die Asymmetrie mit zunehmender Entfernung von der Küste ab. Die Flut-und Ebbstr?me hatten bei Springtide h?here mittlere Geschwindigkeiten als die Str?me bei Nipptide. Die mittlere Geschwindigkeit der Str?me bei Nipptide war eher flutasymmetrisch. Asymmetrie wurde auch zwischen den Richtungen bei Flut-und Ebbstrom einerseits und dem Str?mungswinkel (0–44‡) über die Riffzunge bei Nipp-und Springtide andererseits beobachtet. Die Asymmetrien der Tidestrompara-meter sind offenbar für Wachstum und Erhalt der Riffzungen von Bedeutung. Sie sind wahrscheinlich auf Obertiden sowie geometrie-und inertialbedingte Effekte der Riffzungentopographie zurückzuführen.

     

Gradients and phase velocities of ULF geomagnetic disturbances used to determine the source of an impending strong earthquake

Results of studying the behavior of the vectors of gradients and phase velocities of ULF geomagnetic disturbances (F < 1 Hz) in the Japan seismic zone are presented. The gradient and phase velocity vectors along the Earth’s surface were determined using data of the group of three high-sensitivity three-component magnetovariation stations (MVC-3DS) located at triangle vertices at a small (～5 km) distance from one another (magnetic gradiometer). Two such groups of stations were installed in 1999 southwest and southeast of Tokyo at a distance of ～150 km from each other. It has been indicated that, several months before strong earthquakes (M > 5), the values of gradients and phase velocities start anomalously changing, and directions toward sources of impending strong earthquakes appear in the distribution of gradient vector directions. Directions from sources of impending earthquakes appear in the distribution of phase velocity vector directions. It is proposed to use gradients and phase velocities of ULF and ELF geomagnetic disturbances as one of the factors in a short-term prediction of strong earthquakes.     

Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Some new features of electron density irregularities over SHAR during strong spread F

An RH-560 rocket flight was conducted from Sriharikota rocket range (SHAR) (14°N, 80°E, dip latitude 5.5°N) to study electron density and electric field irregularities during spread F. The rocket was launched at 2130 local time (LT) and it attained an apogee of 348 km. Results of electron density fluctuations are presented here. Two extremely sharp layers of very high electron density were observed at 105 and 130 km. The electron density increase in these layers was by a factor of 50 in a vertical extent of 10 km. Large depletions in electron density were observed around 175 and 238 km. Both sharp layers as well as depletions were observed also during the descent. The presence of sharp layers and depletions during the ascent and the descent of the rocket as well as an order of magnitude less electron density, in 150/300 km region during the descent, indicate the presence of strong large-scale horizontal gradients in the electron density. Some of the valley region irregularities (165/178 km), in the intermediate scale size range, observed during this flight, show spectral peaks at 2 km and can be interpreted in terms of the image striation theory suggested by Vickrey et al. The irregularities at 176 km do not exhibit any peak at kilometer scales and appear to be of new type. The growth rate of intermediate scale size irregularities, produced through generalized Rayleigh Taylor instability, was calculated for the 200/330 km altitude, using observed values of electron density gradients and an assumed vertically downward wind of 20 ms^–1. These growth rate calculations suggest that the observed irregularities could be produced by the gradient drift instability.     

Geomagnetic solar flare effect on February 4 and 6，1986 at the china Antarctic Great wall Station

ＧｅｏｍａｇｎｅｔｉｃｓｏｌａｒｆｌａｒｅｅｆｆｅｃｔｏｎＦｅｂｒｕａｒｙ４ａｎｄ６，１９８６ａｔｔｈｅｃｈｉｎａＡｎｔａｒｃｔｉｃＧｒｅａｔｗａｌｌＳｔａｔｉｏｎＸｉａｏ－ＰｉｎｇＺＥＮＧａｎｄＹｕｎ－ＦａｎｇＬＩＮ（曾小苹，林云芳）（Ｉｎｓｔｉｔ...