Quasilinear heating of electrons in the earth's plasmasphere

In the framework of a quasilinear theory we examine the interaction between thermal electrons and ion-cyclotron waves (ICW) in the outer plasmasphere of the earth. For this type of wave-particle interaction, a simplification of the quasilinear diffusion integral in a magneto-active plasma under plasmaspheric conditions is given. Under the assumption of a Maxwellian distribution of electrons we have calculated the collision frequency and the heating source as the electrons are scattered by ICW. The obtained values of intrinsic parameters of the outer plasmasphere may exceed greatly, accordingly, the Coulomb frequency of collisions and the heating source due to suprathermal electrons. ICW-heating causes the ionosphereward thermal flux to increase, and this must lead to an increase in electron temperature in lower-lying plasmaspheric regions and in the subauroral ionosphere. A quantitative estimation of the electron temperature for the hot zone, made in this paper, is consistent with available experimental data.     

A Local Least-Squares Modification of Stokes’ Formula

The combination of Stokes formula and an Earth Gravity Model (EGM) for geoid determination has become a standard procedure. However, the way of modifying Stokes formula vary from author to author, and numerous methods of modification exist. Most methods are deterministic, with the primary goal of reducing the truncation bias committed by limiting the area of Stokes integration around the computation point, but there are also some stochastic methods with the explicit goal to reduce the global mean square error of the geoid height estimator stemming from the truncation bias as well as the random errors of the EGM and the gravity data. The latter estimators are thus, at least from a theoretical point of view, optimal in a global mean sense, but in a local sense they may be far from optimality.Here we take advantage of the error variance-covariance matrices of the EGM and the terrestrial gravity data to derive the modification parameters of Stokes kernel in a local least-squares sense. The solution is given for the unbiased type of modification of Stokes formula of Sjöberg (1991).     

Critical remarks on the use of medieval eclipse records for the determination of long-term changes in the Earth's rotation

The use of some Arabic medieval solar and lunar eclipse records for the determination of secular changes in the Earth's rotation is critically reviewed. The published results derived from these data suggest a non-uniform decrease in the Earth's rotation rate over the last 27 cy. There is, however, up to this day no sound physical explanation for the deduced non-tidal oscillations, with an apparent period of about 1500 yr and a semi-amplitude of some 4 ms in the l.o.d., which overlayed to a constant secular tidal change in the Earth's rotation rate produce a net non-uniform deceleration of the Earth's rotation. In this paper we discuss a set of observations, which were executed by professional Arabic astronomers. We show by our analysis the way in which the non-uniform deceleration of the Earth's rotation was constructed. A correct reading of the Arabic medieval observations shows that they do not contradict a secular constant decrease in the Earth's rotation rate of nearly -4.6 10-22 rad s-2. This value is in accordance with other similar ones derived from ancient eclipse records and from satellite tracking data.     

Community preparedness for lava flows from Mauna Loa and Hualālai volcanoes,Kona, Hawai‘i

Lava flows from Mauna Loa and Huallai volcanoes are a major volcanic hazard that could impact the western portion of the island of Hawaii (e.g., Kona). The most recent eruptions of these two volcanoes to affect Kona occurred in a.d. 1950 and ca. 1800, respectively. In contrast, in eastern Hawaii, eruptions of neighboring Klauea volcano have occurred frequently since 1955, and therefore have been the focus for hazard mitigation. Official preparedness and response measures are therefore modeled on typical eruptions of Klauea.The combinations of short-lived precursory activity (e.g., volcanic tremor) at Mauna Loa, the potential for fast-moving lava flows, and the proximity of Kona communities to potential vents represent significant emergency management concerns in Kona. Less is known about past eruptions of Huallai, but similar concerns exist. Future lava flows present an increased threat to personal safety because of the short times that may be available for responding.Mitigation must address not only the specific characteristics of volcanic hazards in Kona, but also the manner in which the hazards relate to the communities likely to be affected. This paper describes the first steps in developing effective mitigation plans: measuring the current state of peoples knowledge of eruption parameters and the implications for their safety. We present results of a questionnaire survey administered to 462 high school students and adults in Kona. The rationale for this study was the long lapsed time since the last Kona eruption, and the high population growth and expansion of infrastructure over this time interval. Anticipated future growth in social and economic infrastructure in this area provides additional justification for this work.The residents of Kona have received little or no specific information about how to react to future volcanic eruptions or warnings, and short-term preparedness levels are low. Respondents appear uncertain about how to respond to threatening lava flows and overestimate the minimum time available to react, suggesting that personal risk levels are unnecessarily high. A successful volcanic warning plan in Kona must be tailored to meet the unique situation there.     

Channel overflows of the Pōhue Bay flow,Mauna Loa,Hawai'i: examples of the contrast between surface and interior lava

A number of overflows from a large lava channel and tube system on the southwest rift zone of Mauna Loa were studied. Initial overflows were very low viscosity gas-rich phoehoe evidenced by flow-unit aspect ratios and vesicle sizes and contents. Calculated volumetric flow-rates in the channel range between 80 and 890 m³/s, and those of the overflows between 35 and 110 m³/s. After traveling tens to hundreds of meters the tops of these sheet-like overflows were disrupted into a surface composed of clinker and phoehoe fragments. After these 'a' overflows came to rest, lava from the interiors was able to break out on to the surface as phoehoe. The surface structure of a lava flow records the interaction between the differential shear rate (usually correlated with the volumetric flow-rate) and viscosity-induced resistance to flow. However, the interior of a flow, being better insulated, may react differently or record a later set of emplacement conditions. Clefts of toothpaste lava occurring within fields of clinker on proximal-type 'a' flows also record different shear rates during different times of flow emplacement. The interplay between viscosity and shear rate determines the final morphological lava type, and although no specific portion of lava ever makes a transition from 'a' back to phoehoe, parts of a flow can appear to do so.     

Origin of the fumarolic fluids of Vulcano Island,Italy and implications for volcanic surveillance

Variations in D and ¹⁸O values with H₂O contents and outlet temperatures indicate that the fumaroles of La Fossa crater have discharged mixtures of magmatic water and marine hydrothermal water, since 1979. The contribution of meteoric water was low in the period 1979–1982 and very low afterwards. The ¹⁸O values of the marine-hydrothermal component of +5 to +7.2 are due to isotopic exchange with the ¹⁸O-rich silicates of the rocks under high-temperature and low-permeability conditions. The ¹⁸O value of the magmatic end-member is generally +3.5 to +4.3, although values as high as +5.5 to +6.5 were reached in the summer of 1988, when magma degassing appears to have extended into the core of the magma body. The D values of the end-member were close to -20, typical of andesitic waters. Both the isotopic values and chemical data strongly support a dry model, consisting of a central magmatic gas column and a surrounding hydrothermal envelope, in which marine hydrothermal brines move along limited fracture zones to undergo total evaporation on approaching the conduits of magmatic fluids. The vents at the eastern and western boundaries of the fumarolic field are fed by fluids whose pressure is governed by the coexistence of vapor, liquid and halite, giving rise to a high risk of phreato magmatic explosions, should magma penetrate into these wet environments. Most La Fossa eruptions were triggered by an initial hydrothermal blast and continued with a series of phreatomagmatic explosions. The fluids discharged by the Forgia Vecchia fumaroles are mixed with meteoric water, which is largely evaporated, although subordinate loss of condensed steam may be responsible for scrubbing most of the acidic gas species. The temperatures and pressures, and the risk of a sudden pressure increase, are low. A boiling hydrothermal aquifer at 230° C is present underneath the Baia di Levante beach. This area has a minor risk of hydrothermal explosions.     

Mean Field Dynamos with Algebraic and Dynamic α-Quenchings

Calculations for mean field dynamo models (in both full spheres and spherical shells), with both algebraic and dynamic -quenchings, show qualitative as well as quantitative differences and similarities in the dynamical behaviour of these models. We summarise and enhance recent results with extra examples. Overall, the effect of using a dynamic appears to be complicated and is affected by the region of parameter space examined.     

Cenozoic Intraplate Alkaline Volcanism of Western Bohemia

Three independent volcanic suites have been recognised in W Bohemia: (i) the old unimodal alkaline ol. nephelinite-tephrite (29-19 Ma) in the Ohe Rift, (ii) two contemporaneous weakly (trachybasalt/trachyandesite-trachyte/rhyolite; 13-11 Ma) and strongly (ol. nephelinite-tephrite/basanite; 12-8 Ma) alkaline series in the flank of the Cheb-Domalice Graben formed by the Teplá Highland and (iii) the young unimodal ol. melilitite/ol. nephelinite alkaline suite (2.0-0.12 Ma) at the intersection of the above mentioned structures in the Cheb Basin. The magmas of all the suites are mantle-derived and, in the case of the Cheb-Domalice Graben series, associated with the AFC process. Two main fault systems: (i) ENE-WSW and (ii) NNW-SSE are developed in W Bohemia, corresponding to the directions of the two prominent taphrogenic structures. The southwesterly continuation of the Ohe Rift across the Mariánské Lázn Fault is marked by volcanics only.     

Integral Equations for Kinematic Dynamo Models

A new technique for the treatment of the kinematic dynamo problem is presented. The method is applicable when the dynamo is surrounded by a medium of finite conductivity and is based on a reformulation of the induction equation and boundary conditions at infinity into an integral equation. We show that the integral operator involved here is compact in the case of homogeneous conductivity, which is important for both mathematical and numerical treatment. A lower bound for the norm of then yields a necessary condition for the generation of magnetic fields by kinematic dynamos. Numerical results are presented for some simple ²-dynamo models. The far-field asymptotics for stationary and time-dependent field modes are discussed.     

Comparison between Pollutant Concentration in the Samples of Fog and Rime Water Collected at Mt. Milešovka

The article aims at showing the differences in concentration of pollutants that are contained in the samples of fog and/or low cloud water in comparison with the water from rime. The results follow from the fog and rime measurements made at Mt. Mileovka (eské Stedohoí Mountains). They are compared with the results of other studies that also report the differences in fog and rime chemistry.     

Tilt Observations in the Normal Mode Frequency Band at the Geodynamic Observatory Cueva de los Verdes,Lanzarote

We have conducted observations with the aid of a seismo-tiltmeter station, which is based on the Ostrovsky pendulum and installed at the Geodynamic Observatory Cueva de los Verdes at Lanzarote Island since 1995. In this station the signal is separated into two frequency bands – tidal tilts (from 0 to 5 mHz) and ground oscillations in the frequency range of free Earths normal modes (from 0.2 to 5 mHz). The later band, called accelerometer channel, has additional amplification. We analyzed the background records in the frequency range of Earths free oscillations from August 2000 to September 2001, as well as, Earths normal modes after strong earthquakes. We found several distinctive persistent peaks in the spectra of background oscillations. Both amplitudes of distinguished peaks and noises have seasonal variations. We found that spectra of background oscillations are different in the frequency interval between 1.4 and 2.5 mHz for North- South and East-West components.     

Lagrangian Solution of Dynamo Equations and Path Integrals

We discuss an explicit solution of the Cauchy problem for induction equation and suggest its generalization for equations of ²-dynamo. These solutions are based on concepts of multiplicative, Wiener path, and stochastic integrals. Obtained explicit solution can be useful as a tool in investigations of a dynamo with fluctuating helicity.     

A method of searching for radial velocity anomalies in the lower mantle using interfering signals

Radial velocity anomalies in the lower mantle that give rise to triplications in the travel-time curve for short-periodP waves will produce arrivals havingdT/d values that differe by roughly 0.2–0.5 s/deg. The first two arrivals associated with such triplications will be separated by less than one second over a distance range of 4°–10° they may not, therefore, be separable visually on single seismograms, so that their presence can only be inferred from some measurable property that depends on their mutual interference. If there are lateral variations in the regions of anomalous velocity gradients, the interfering signals will also have different azimuths of arrival. Using two synthetic wavelets we have investigated the effect of interference on bothdT/d and azimuth measurements at the Yellowknife Array. We found that if the interfering pulses have a dominant frequencyv, there is a range of time separations (0.30/v0.55/v) over which the measureddT/d and azimuth values may fluctuate by much more than the differences indT/d and azimuth between the interfering signals. We have evaluated the following empirically defined functions for three different primary signals, and for three different relative amplitudes of the interfering signals:f (t), the drift function, which expresses how the measured slownesses,p, and azimuths, , differ from the slownesses and azimuths of the primary wavelets; f(), the range function, which describes the behaviour of the upper and lower bounds ofp and as a function of the difference in arrival times of the signals, andf , studied the properties of these functions, and have outlined how these properties provide criteria based on the numerical and statistical characteristics of the arrival vectors, and on the waveform of the signal that will enable small radial velocity anomalies to be more clearly delineated.Contribution No. 863 from the Earth Physics Branch.     

A note on the general concept of wave breaking for Rossby and gravity waves

A recently proposed general definition of wave breaking is further discussed, in order to deal with some points on which misunderstanding appears to have arisen. As with surface and internal gravity waves, the classification of Rossby waves into breaking and not breaking is a generic classification based on dynamical considerations, and not a statement about any unique signature or automatically recognizable shape. Nor is it a statement about passive tracers uncorrelated with potential vorticity on isentropic surfaces. A strong motivation for the definition is that proofs of the nonacceleration theorem of wave, mean-flow interaction theory rely, explicitly or implicitly, on a hypothesis that the waves do not break in the sense envisaged.The general definition refers to the qualitative behaviour of a certain set of material contours, namely those, and only those, which would undulate reversibly, with small slopes, under the influence of the waves' restoring mechanism, in those circumstances for which linearized, nondissipative wave theory is a self-consistent approximation to nonlinear reality. The waves' restoring mechanism depends upon the basic-state vertical potential density gradient in the case of gravity waves, and upon the basic-state isentropic gradient of potential vorticity in the case of Rossby waves. In the usual linearized theory of planetary scale Rossby waves on a zonal shear flow, the relevant material contours lie along latitude circles when undisturbed.     

A statistical study of the upstream intermediate ion boundary in the Earth’s foreshock

A statistical investigation of the location of onset of intermediate and gyrating ion populations in the Earths foreshock is presented based on Fixed Voltage Analyzer data from ISEE 1. This study reveals the existence of a spatial boundary for intermediate and gyrating ion populations that coincides with the reported ULF wave boundary. This boundary position in the Earths foreshock depends strongly upon the magnetic cone angle _BX and appears well defined for relatively large cone angles, though not for small cone angles. As reported in a previous study of the ULF wave boundary, the position of the intermediate-gyrating ion boundary is not compatible with a fixed growth rate of the waves resulting from the interaction between a uniform beam and the ambient plasma. The present work examines the momentum associated with protons which travel along this boundary, and we show that the variation of the boundary position (or equivalently, the associated particle momentum) with the cone angle is related to classical acceleration mechanisms at the bow shock surface. The same functional behavior as a function of the cone angle is obtained for the momentum predicted by an acceleration model and for the particle momentum associated with the boundary. However, the model predicts systematically larger values of the momentum than the observation related values by a constant amount; we suggest that this difference may be due to some momentum exchange between the incident solar-wind population and the backstreaming particles through a wave-particle interaction resulting from a beam plasma instability.     

Estimation of the Source Time Function Based on Blind Deconvolution with Gaussian Mixtures

It is demonstrated that the blind deconvolution method is fully capable of recovering the unknown Greens function and of estimating the source time functions from observed seismic data of small earthquakes. Based on the assumption of the Gaussian-mixture model of the Greens function, the newly-formulated algorithm is evaluated using synthetic seismic data along with those of the May 8, 1996 Mexico earthquake (M_c = 4.6). Since the estimated results closely match the theoretical input very well, the method is then employed to analyze the source time functions of the July 7, 1995 Pu-Li, Taiwan earthquake (M_L = 5.3). The stations triggered by this event were azimuthally well covered. Using the estimated source time functions, information pertaining to the directivity effect is readily obtained, and the actual fault plane of this event is identified, thus clearly indicating that this method provides a most efficient way to estimate the source time function of a small earthquake.Acknowledgment The authors would like to express their thanks to two anonymous reviewers for their valuable suggestions and Dr. I. Santamarias courteous assistance. They also appreciate the efforts of Drs. H.C. Chiu and R.J. Rau, who provided the seismic data and the fault plane solutions. The National Science Council, Taiwan, has supported this research (NSC 91-2119-M-194-011).     

A daily rainfall occurrence process

A model for the periodic (annual cycle), discrete rainfall occurrence process is presented. Using this model the probabilistic properties of the process in -day intervals can be investigated. In such an interval the rainfall occurrence process is approximated by some stationary processa _t ,tIN. The processa _t ,tIN is described by the distributions of the lengths of wet and dry sequences. It is assumed that the lengths of successive wet and dry sequences are independent. For this process the distribution of the number of wet days in -day intervals is calculated. The model is fitted to 50-year rainfall data from Wroclaw, Poland. Rainfall amounts of 0.1, 1.0 and 2.0 mm are considered as thresholds defining a wet day. To estimate the distribution of the length of wet and dry sequences the family of Pascal distribution is chosen.     

Large-scale structure of magnetospheric plasma

Recent investigations of magnetospheric plasma structure are summarized under the broad categories of empirical models, transport across boundaries, formation, and dynamics of the plasma sheet. This report reviews work in these areas during the period 1991 to 1993. Fully three-dimensional empirical models and simulations have become important contributors to our understanding of the magnetospheric system. Some new structural concepts have appeared in the literature: the entry boundary and geopause, the plasma sheet region 1 vortices, the low-energy layer, the adiabaticity boundary or wall region, and a region in the tail to which we refer as the injection port. Traditional structural concepts have also been the subject of recent study, notably the plasmapause, the magnetopause, and the plasma sheet. Significant progress has been made in understanding the nature of plasma sheet formation and dynamics, but the acceleration of electrons to high energy remains somewhat mysterious.     

Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.