Field-aligned conductance values estimated from Maxwellian and kappa distributions in quiet and disturbed events using Freja electron data

We study the question of what difference it makes for the derived field-aligned conductance (K) values if one uses Maxwellian or kappa distributions for the fitting of low-orbiting satellite electron flux spectra in the auroral region. This question has arisen because sometimes a high-energy tail is seen in the spectra. In principle, the kappa fits should always be better, because the kappa distribution is a generalization of the Maxwellian. However, the physical meaning of the parameters appearing in the Maxwellian is clearer. It therefore makes sense to study under which circumstances it is appropriate to use a Maxwellian. We use Freja electron data (TESP and MATE) from two events. One of the events represents quiet magnetospheric conditions (stable arc) and the other represents disturbed conditions (surge). In these Freja events, at least, using kappa rather than Maxwellian fitting gives a better fit to the observed distribution, but the difference in K values is not large (usually less than 20%). The difference can be of either sign. However, sometimes even the kappa distribution does not provide a good fit, and one needs a more complicated distribution such as two Maxwellians. We investigate the relative contributions of the two Maxwellians to the total field-aligned conductance value in these cases. We find that the contribution of the high-energy population is insignificant (usually much less than 20%). This is because K is proportional to n/E_c E_c, where n is the source plasma density and E_c is the characteristic energy.     

Effects of a kappa distribution function of electrons on incoherent scatter spectra

In usual incoherent scatter data analysis, the plasma distribution function is assumed to be Maxwellian. In space plasmas, however, distribution functions with a high energy tail which can be well modeled by a generalized Lorentzian distribution function with spectral index kappa (kappa distribution) have been observed. We have theoretically calculated incoherent scatter spectra for a plasma that consists of electrons with kappa distribution function and ions with Maxwellian neglecting the effects of the magnetic field and collisions. The ion line spectra have a double-humped shape similar to those from a Maxwellian plasma. The electron temperatures are underestimated, however, by up to 40% when interpreted assuming Maxwellian distribution. Ion temperatures and electron densities are affected little. Accordingly, actual electron temperatures might be underestimated when an energy input maintaining a high energy tail exists. We have also calculated plasma lines with the kappa distribution function. They are enhanced in total strength, and the peak frequencies appear to be slightly shifted to the transmitter frequency compared to the peak frequencies for a Maxwellian distribution. The damping rate depends on the electron temperature. For lower electron temperatures, plasma lines for electrons with a distribution function are more strongly damped than for a Maxwellian distribution. For higher electron temperatures, however, they have a relatively sharp peak.     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

Analytical One-Way Plane-Wave Solution in the 1-D Anisotropic “Simplified Twisted Crystal” Model

Analytical expressions for the exact 2 × 2 one-way propagator matrix of a plane S wave, propagating along the axis of spirality in the simple 1-D anisotropic simplified twisted crystal model, are presented. The analytical equations are useful in testing the applicability and accuracy of various approximate wavefield modelling methods, especially of the coupling ray theory and of its various quasi-isotropic approximations and various numerical implementations.In addition to the exact analytical solution of the elastodynamic equation in the simplified twisted crystal model, the analytical solutions of the equations of the four ray methods are given. The ray methods are (a) the coupling ray theory, (b) the coupling ray theory with the quasi-isotropic perturbation of travel times, (c) the anisotropic ray theory, (d) the isotropic ray theory. These four approximate solutions of the elastodynamic equation are roughly compared with the exact solution. Both the exact analytical solution and the analytical ray-theory solutions in the simplified twisted crystal model are also helpful in debugging computer codes for various approximate wavefield modelling methods, especially for the coupling ray theory.     

Approximate expression for the hilbert transform of a certain class of functions and their application to the ray theory of seismic waves

Summary Approximate expressions for the Hilbert transform of the functionf(t)=exp(- ₀ ² t ²/²) cos( ₀ t+v) are determined. This function, given a suitable choice of the three parameters ₀, and v, approximates a wide class of seismic signals very well. The approximate expressions for the Hilbert transform enable very simple formulae to be given for the elementary seismograms of the individual seismic body waves (in the zero approximation of the ray theory). This accelerates the computation of ray theoretical seismograms considerably.     

Analytical study of the energy rate balance equation for the magnetospheric storm-ring current

We present some results of the analytical integration of the energy rate balance equation, assuming that the input energy rate is proportional to the azimuthal interplanetary electric field, E_y, and can be described by simple rectangular or triangular functions, as approximations to the frequently observed shapes of E_y, especially during the passage of magnetic clouds. The input function is also parametrized by a reconnection-transfer efficiency factor (which is assumed to vary between 0.1 and 1). Our aim is to solve the balance equation and derive values for the decay parameter compatible with the observed Dst peak values. To facilitate the analytical integration we assume a constant value for through the main phase of the storm. The model is tested for two isolated and well-monitored intense storms. For these storms the analytical results are compared to those obtained by the numerical integration of the balance equation, based on the interplanetary data collected by the ISEE-3 satellite, with the values parametrized close to those obtained by the analytical study. From the best fit between this numerical integration and the observed Dst the most appropriate values of are then determined. Although we specifically focus on the main phase of the storms, this numerical integration has been also extended to the recovery phase by an independent adjust. The results of the best fit for the recovery phase show that the values of may differ drastically from those corresponding to the main phase. The values of the decay parameter for the main phase of each event, _m, are found to be very sensitive to the adopted efficiency factor, , decreasing as this factor increases. For the recovery phase, which is characterized by very low values of the power input, the response function becomes almost independent of the value of and the resulting values for the decay time parameter, _r, do not vary greatly as varies. As a consequence, the relative values of between the main and the recovery phase, _m/_r, can be greater or smaller than one as varies from 0.1 to 1.     

Single- and inter-station transfer functions for non-synchronous geomagnetic arrays-II. the results of the field data analysis

Summary The theory of methods of computing single- and inter-station transfer functions in both the spectral and time domains was developed in paper[1]. Both approaches are applied to the variation data recorded at field stations along two non-simultaneous profiles traversing the eastern margin of the Bohemian Massif, where a zone of anomalous induction seems to mark an important geological boundary of formations with different histories of development. The results of both analyses are found to coincide within reasonable bounds of 20–30% in the principal induction characteristics.

---

u m¶rt; ama ¶rt;-u -mau n¶rt;am u nma u am a ua am[1]. am nua m am a n¶rt;¶rt;a nu ¶rt; aaua ¶rt;aaum auau aumua a n mau ¶rt; u nu, nau m au aua,¶rt; aa a aa u¶rt;uu. a, u¶rt;u, mamau a¶rt;a ¶rt; ¶rt; ¶rt;uuau au umuu aumu. mam aau nma u am auam a 20–30% ¶rt; u¶rt;u naam.

     

Semiannual and annual variations in the height of the ionospheric F2-peak

Ionosonde data from sixteen stations are used to study the semiannual and annual variations in the height of the ionospheric F2-peak, hmF2. The semiannual variation, which peaks shortly after equinox, has an amplitude of about 8 km at an average level of solar activity (10.7 cm flux = 140 units), both at noon and midnight. The annual variation has an amplitude of about 11 km at northern midlatitudes, peaking in early summer; and is larger at southern stations, where it peaks in late summer. Both annual and semiannual amplitudes increase with increasing solar activity by day, but not at night. The semiannual variation in hmF2 is unrelated to the semiannual variation of the peak electron density NmF2, and is not reproduced by the CTIP and TIME-GCM computational models of the quiet-day thermosphere and ionosphere. The semiannual variation in hmF2 is approximately isobaric, in that its amplitude corresponds quite well to the semiannual variation in the height of fixed pressure-levels in the thermosphere, as represented by the MSIS empirical model. The annual variation is not isobaric. The annual mean of hmF2 increases with solar 10.7 cm flux, both by night and by day, on average by about 0.45 km/flux unit, rather smaller than the corresponding increase of height of constant pressure-levels in the MSIS model. The discrepancy may be due to solar-cycle variations of thermospheric winds. Although geomagnetic activity, which affects thermospheric density and temperature and therefore hmF2 also, is greatest at the equinoxes, this seems to account for less than half the semiannual variation of hmF2. The rest may be due to a semiannual variation of tidal and wave energy transmitted to the thermosphere from lower levels in the atmosphere.     

Approximate Conditions for the Off-Axis Triplication in Transversely Isotropic Media

The presented approximate formulas yield a critical value of anisotropy parameter , for which an incipient off-axis SV-wave triplication occurs in transversely isotropic media. The formulas are simple but approximate the exact solution with a high accuracy. The best results are obtained using the third-order approximation, which yields accuracy at least 30 times higher than the formulas presented by Thomsen and Dellinger (2003). The formula works safely for parameters  = a ₃₃/a ₄₄ > 2 and 0.2 >  = (a ₁₁ – a ₃₃)/2a ₃₃ > –0.2, and yields critical values of from 0.1 to 0.7. Outside this interval, it is recommended to use an exact solution.     

The Terrain Correction in a Moving Tangent Space

Conventionally the terrain/topographic reduction is based on the Bouguer Plate, which is flat and extends in the local tangent plane/horizontal plane to infinity. Here we aim at an error estimate of such a planar approximation of the Newton integral of the type of a disturbing potential and gravitational disturbance as linearized forms of the gravitational potential and the modulus of gravitational field intensity. To effect this quality control of the conventional terrain reduction, we first transform the spherical Newton functional from an equatorial frame of reference to an oblique meta-equatorial frame of reference with the evaluation point as a meta-North pole, and then by means of an oblique equiareal map projection of the azimuthal type to a tangent plane which moves at the evaluation point. The first term of these transformed Newton functionals is the planar approximation. The difference between the exact Newton kernels and their planar approximation are plotted and tabulated in Tables 1-3. Three configurations are studied in detail: for points at radius r = 10 km around the evaluation point the systematic error varies from 0.26% for a spherical height difference of the order of H – H* = 5 km, more than 0.80% for a spherical height difference of the order of H – H* = 1 km, and more than 1.60% for a spherical height difference of H – H* = 500 m. In contrast, the systematic error for spherical height difference H – H* = 1 km at a distance of r = 1000 km from the evaluation point increases to 44%. Indeed, the newly derived exact Newton kernels which are of the convolution type and are represented in the tangent space moving with the evaluation point can be preferably used with little extra computational effort.     

A statistical analysis for pattern recognition of small cloud particles sampled with a PMS-2DC probe

Although small particles (size between 25 m and 200 m) are frequently observed within ice and water clouds, they are not generally used properly for the calculation of structural, optical and microphysical quantities. Actually neither the exact shape nor the phase (ice or water) of these particles is well denned since the existing pattern recognition algorithms are only efficient for larger particle sizes. The present study describes a statistical analysis concerning small hexagonal columns and spherical particles sampled with a PMS-2DC probe, and the corresponding images are classified according to the occurrence probability of various pixels arrangements. This approach was first applied to synthetic data generated with a numerical model, including the effects of diffraction at a short distance, and then validated against actual data sets obtained from in-cloud flights during the pre-ICE89 campaign. Our method allows us to differentiate small hexagonal columns from spherical particles, thus making possible the characterization of the three dimensional shape (and consequently evaluation of the volume) of the particles, and finally to compute e.g., the liquid or the ice water content.     

Magnitudes of detectable and localizable shocks in nw bohemia

Summary The effectiveness of recording seismic phenomena in the Kruné hory (Mts.) region in NW Bohemia by selected stations in the CSR, GDR and Poland has been estimated. Magnitude isolines of the weakest earthquakes, which can be localized and detected with an 0.9 probability, were calculated on the basis of the level of seismic disturbances at the individual stations and of the empirical dependence of the attenuation of seismic waves with distance.

---

a a mum umauu uu u amu ¶rt; ana¶rt; uu uau mauu a mumuu , u a a uu n a m¶rt; mau u nuu auumu amau uu m amu u auma uuuu aum¶rt; a a mu, m mm 0.9 auuam u aum.

     

Elektronenproduktion,Rekombination, Elektronendichte und Absorption in der ionosphärischenD-Region und ihre jahreszeitlichen Variationen

Zusammenfassung Es wird die Gleichung für die Elektronenproduktionq(z) abgeleitet, die die meteorologischen Elemente der Mesosphäre berücksichtigt. Nach Angaben über die mit Satelliten und Raketen gemessene Röntgenstrahlung mit 8 Å wird das Differentialspektrum des ionisierenden Energieflusses für eine mittlere Sonnenaktivität konstruiert. Auf dieser Grundlage und nach der bekannten Intensität der Strahlung Ly- sowie nach Angaben über dieElektronenproduktion der kosmischen Strahlung werden die Profileq _Rö(z),q _Ly-(z) undq _CR(z) für mittlere geographische Breiten und Standardatmosphäre entwickelt. Nach eingehender Analyse der vollständigen Gleichung für den effektiven Rekombinationskoeffizienten wird für die Verhältnisse in der tiefen Ionosphäre der Beitrag jeder einzelnen Komponente der Gleichung bestimmt. ist eine recht veränderliche Grösse, die von den aeronomischen und meteorologischen Verhältnissen und der Sonnenzenitdistanz abhängt. Aus den fürq(z) und (z) erhaltenen Angaben werden zwei ElektronendichteprofileN(z) für =30° und 75° erhalten. Das ProfilN(z) bei =30° wird mit dem gemittelten Profil einer umfangreichen Gruppe experimentell gefundener VerteilungenN(z) verglichen; das Profil bei =75° wird durch Messung der deviativen und nondeviativen Absorption für eine längere Zeitperiode überprüft. In beiden Fällen hat sich die Richtigkeit der theoretisch erhaltenen Profile bestätigt. Die jahreszeitlichen Variationen der nondeviativen Absorption in derD-Region sind ausschliesslich durch die Variationen der meteorologischen Parameter im Bereich der Mesopause bei konstantem Energiefluss der ionisierenden Strahlung bedingt.

---

Summary An equation about the electron production is deduced in which the meteorological elements of the mesosphere are taken into account. The differential spectrum of the ionizing energy flux with 3 Å for average solar activity is constructed on evidence from rocket and satelitc measurements. The profilesq _Rö(z),q _Ly-(z) andq _CR(z) for mean geographical latitudes and standard atmosphere are plotted on that basis as well as on data fot the known intensity of the Ly- emission and the electron production of the cosmic rays. An exhaustive analysis is made of the full equation for the effective recombination coefficient and the contribution of all its components at lower ionosphere conditions is determined. is a rather variable quantity, dependent on the aeronomical and meteorological condition of the area under consideration, as well as on the solar zenith angle. Two profiles for the electron concentrationN(z) at =30° and 75° are drawn on the basis of data forq(z) and (z). The profileN(z) at =30° is compared with the averaged profile of a large group experimentally obtained distributionsN(z); the profile at =75° is checked by measurements of the deviative and nondeviative absorption taken for a lengthy period. Both checks are in good agreement with the theoretically obtained profiles. The seasonal variations of the nondeviative absorption in theD region could be completely explained with the variations of the meteorological parameters in the mesopause area at constant energy flux of the ionizing radiation.

     

The vibration magnetometer and the processing of its output signal strongly affected by noise

¶rt;aam ¶rt;a m¶rt;a amu uaa, u u , u auam u ma mum nuu ¶rt; uau aumma. u n u ama u ma n¶rt; ma m¶rt;a nm a umam u naa u ¶rt; nuu. mm m¶rt; n au auu u aau n¶rt;mam a au nm aum nu m unu a a.     

Deep seismically active fracture zones in ecuador and northern peru

Summary A system of 8 seismically active fracture zones was delineated on the basis of the distribution of earthquake foci in the continental lithosphere of Ecuador. The position and width of the outcrop, thickness, dip and maximum depth of the individual fracture zones were estimated and correlated with surface geological and tectonic phenomena, volcanism and hydrothermal manifestations. The existence and strike of the fracture zones was independently confirmed by the occurrence of historical disastrous earthquakes.

---

uma 8 uu amu a a aa a auu an¶rt;u a mu muma um a¶rt;a. u n¶rt; nu u uua a nmu, mua, u auaaua m¶rt; a . mau mu auu n¶rt;m¶rt;am a¶rt;u umuu aum mu u nmau nmu, mmu, au uu¶rt;mau nuu.

---

---

Visiting professors at Instituto Geofísico and Facultad de Geología, Minas y Petróleos, Escuela Politécnica Nacional, Quito (Ecuador).     

Spectral theory of electromagnetic induction in a radially and laterally inhomogeneous earth

Summary The partial differential equations of electromagnetic induction in a 3-D Earth of inhomogeneous conductivity are reduced to a system of ordinary differential equations of the 2nd order for the spectral coefficients of the field.

---

au am nu¶rt; ¶rt; maum u¶rt;uu u m ¶rt;¶rt; n n¶rt;umu n¶rt; um ¶rt;uua au m n¶rt;a ¶rt; nma uum n.

     

Sudden ionospheric disturbances and flare activity in the course of the solar cycle no 20

a¶rt;a ma nu an u u (SID) u nu n a (a n¶rt; au) mu 1965–1975 . u u SID na¶rt;am nau, mum u nm R.     

Effects of the imf sector structure in the midlatitude troposphere and stratosphere

Summary The data on geopotential heights and temperatures at 7 pressure levels between 1000-10 hPa above Berlin(52.5 °N, 13.4 °E) are analysed for the winters of 1963–1973. No demonstrable effect of the interplanetary magnetic field sector boundary crossing (IMF SBC) is found in the lower and middle stratosphere, but there is a demonstrable effect in the middle troposphere at the 500 hPa level. This effect is less important than the IMF SBC effect in the tropospheric vorticity area index and seems to be of a different type.

---

auum ¶rt;a nnmua m u mnam a 7 nm ¶rt;au ¶rt; 1000-10 a a¶rt; u(52,5 °.., 13,4 °.¶rt;.) ¶rt; u 1963–1973. ua ¶rt;aam m nu mau nam aum n( ) ¶rt;a amu u u ¶rt; mam, ma m a¶rt; ¶rt; mn a 500 a. mm m a, m u¶rt; na¶rt;u aumu am, u am m ¶rt; muna.

     

Motion of particles in turbulent flow: Hyperbolicity,strange attractors and dynamic stochasticity of a system

¶rt;naa, m ma um uu maunuu m am muu ¶rt;uauu um. nua a mau ammama a, n¶rt; mnu ma u u au uu u¶rt;mu.     

Scattering ofP andS waves by a spherically symmetric inclusion

Scattering of an arbitrary elastic wave incident upon a spherically symmetric inclusion is considered and solutions are developed in terms of the spherical vector system of Petrashen, which produces results in terms of displacements rather than displacement potentials and in a form suitable for accurate numerical computations. Analytical expressions for canonical scattering coefficients are obtained for both the cases of incidentP waves and incidentS waves. Calculations of energy flux in the scattered waves lead to elastic optical theorems for bothP andS waves, which relate the scattering cross sections to the amplitude of the scattered fields in the forward direction. The properties of the solutions for a homogeneous elastic sphere, a sphere filled by fluid, and a spherical cavity are illustrated with scattering cross sections that demonstrate important differences between these types of obstacles. A general result is that the frequency dependence of the scattering is defined by the wavelength of the scattered wave rather than the wavelength of the incident wave. This is consistent with the finding that the intensity of thePS scattering is generally much stronger than theSP scattering. When averaged over all scattering angles, the mean intensity of thePS converted waves is2V _p ² /V _s ⁴ times the mean intensity of theSP converted waves, and this ratio is independent of frequency. The exact solutions reduce to simple and easily used expressions in the case of the low frequency (Rayleigh) approximation and the low contrast (Rayleigh-Born) approximation. The case of energy absorbing inclusions can also be obtained by assigning complex values to the elastic parameters, which leads to the result that an increase in attenuation within the inclusion causes an increased scattering cross section with a marked preference for scatteredS waves. The complete generality of the results is demonstrated by showing waves scattered by the earth's core in the time domain, an example of high-frequency scattering that reveals a very complex relationship between geometrical arrivals and diffracted waves.