The Kelvin-Helmholtz instability on the magnetopause

Conditions for the development of Kelvin-Helmholtz (K-H) waves on the magnetopause have been known for more than 15 years; more recently, spacecraft observations have stimulated further examination of the properties of K-H waves. For amagnetopause with no boundary layer, two different modes of surface waves have been identified and their properties have been investigated for various assumed orientations of magnetic field and flow velocity vectors. The power radiated into the magnetosphere from the velocity shear at the boundary has been estimated. Other calculations have focused on the consequences of finite thickness boundary layers, both uniform and non-uniform. The boundary layer is found to modify the wave modes present at the magnetopause and to yield a criterion for the wavelength of the fastest growing surface waves. The paper concludes by questioning the extent to which the inferences from boundary layer models are model dependent and identifies areas where further work is needed or anticipated.     

Boltysh,another end‐Cretaceous impact

Abstract— The Chixculub impact occurred at the Cretaceous/Tertiary (K/T) boundary, and although several other Late Cretaceous and Paleogene impact craters have, at times, been linked with the K/T boundary, isotope geochronology has demonstrated that all have significantly different ages. The currently accepted age of the 24 km diameter Boltysh crater, a K‐Ar whole‐rock age, places it in the Coniacian at 88 ± 3 Ma. However, comprehensive Ar‐Ar dating of a range of melt samples yields a mean age of 65.17 ± 0.64 Ma, within errors of the K/T boundary. Several of the fresh samples exhibit signs of excess argon but this seems to be concentrated in rapidly crystallized glass‐rich samples. The Ar‐Ar age confirms an earlier fission track measurement and thus two dating techniques have yielded an age within errors of the K/T boundary for this crater. Crucially, although the ages of Boltysh and Chixculub are within errors, they may not have formed synchronously. Craters of 24 km diameter occur much more commonly than impacts of Chixculub dimensions, but their proximity does raise the important question of how many impacts there might have been close to the K/T boundary.     

Calculating and Testing Nonlinear Force-Free Fields

Improvements to an existing method for calculating nonlinear force-free magnetic fields (Wheatland, Solar Phys. 238, 29, 2006) are described. In particular a solution of the 3-D Poisson equation using 2-D Fourier transforms is presented. The improved nonlinear force-free method is demonstrated in application to linear force-free test cases with localized nonzero values of the normal component of the field in the boundary. These fields provide suitable test cases for nonlinear force-free calculations because the boundary conditions involve localized nonzero values of the normal components of the field and of the current density, and because (being linear force-free fields) they have more direct numerical solutions. Despite their simplicity, fields of this kind have not been recognized as test cases for nonlinear methods before. The examples illustrate the treatment of the boundary conditions on current in the nonlinear force-free method, and in particular the limitations imposed by field lines that connect outside of the boundary region.     

3D Solar Null Point Reconnection MHD Simulations

Numerical MHD simulations of 3D reconnection events in the solar corona have improved enormously over the last few years, not only in resolution, but also in their complexity, enabling more and more realistic modeling. Various ways to obtain the initial magnetic field, different forms of solar atmospheric models as well as diverse driving speeds and patterns have been employed. This study considers differences between simulations with stratified and non-stratified solar atmospheres, addresses the influence of the driving speed on the plasma flow and energetics, and provides quantitative formulas for mapping electric fields and dissipation levels obtained in numerical simulations to the corresponding solar quantities. The simulations start out from a potential magnetic field containing a null-point, obtained from a Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) magnetogram magnetogram extrapolation approximately 8?hours before a C-class flare was observed. The magnetic field is stressed with a boundary motion pattern similar to?–?although simpler than?–?horizontal motions observed by SOHO during the period preceding the flare. The general behavior is nearly independent of the driving speed, and is also very similar in stratified and non-stratified models, provided only that the boundary motions are slow enough. The boundary motions cause a build-up of current sheets, mainly in the fan-plane of the magnetic null-point, but do not result in a flare-like energy release. The additional free energy required for the flare could have been partly present in non-potential form at the initial state, with subsequent additions from magnetic flux emergence or from components of the boundary motion that were not represented by the idealized driving pattern.     

Boundary layer ion composition at Jupiter during the inbound pass of the Ulysses flyby

During the inbound segment of the Ulysses flyby of Jupiter, there were multiple incursions into the dawnside low-latitude boundary layer, as identified by Bame et al. (Science257, 1539–1542, 1992) using plasma electron data. In the present study, ion composition and spectral measurements provide independent collaborative evidence for the existence of distinct boundary layer regions. Measurements are taken in the energy-per-charge range of 0.6–60 keV/e and involve mass as well as mass-per-charge identification by the Ulysses/SWICS experiment. Ion species of Jovian magnetospheric origin (including O⁺, O²⁺, S²⁺, S³⁺) and sheath origin (including He²⁺ and high charge state CNO) have been directly identified for the first time in the Jovian magnetospheric boundary layer. Protons of probably mixed origin and He⁺ of possibly sheath (ultimately interstellar pickup) origin were also observed in the boundary layer. Sheath-like ions are observed throughout the boundary layer; however, the Jovian ions are depleted or absent for portions of two boundary layer cases studied. Ions of solar wind origin are observed within the outer magnetosphere. and ions of magnetospheric origin are found within the sheath, indicating that transport across the magnetopause boundary can work both ways, at least under some conditions. Although their source cannot be uniquely identified, the proton energy spectrum in the boundary layer suggests a sheath origin for the lower energy protons.     

High latitude,outer zone boundary observations of electrons and protons

The diurnal variation of the high latitude outer zone boundary at 1400km has been determined for electrons ?140keV electrons, and for protons in two energy intervals: 0.56?E?1.1 MeV, 1.1?E?3.2 MeV, from detectors aboard the NOAA-2 satellite. The dependence of the 140 keV electron boundary on D_st has been examined as well?. A wel?l defined correlation of boundary position with D_st is found to exist during the main phase of disturbances, together with an evident local time dependence. All the boundaries are found to be consistent witn the supposition of adiabatic drift and demonstrate the stability of the boundary position over approximately ten years of comparable observation. No statistically significant hemispheric differences in boundary location were observed to occur.     

Unrestricted second-order tensor virial equations for linear oscillations of magnetic configurations

This paper deals with the second-order tensor virial equations for the linear oscillations of a gaseous mass in the presence of a magnetic field. It is shown that the commonly used linearized versions of the tensor virial equations are restricted integral equations that incorporate the linearized equation of motion but not the boundary condition. These restricted equations only allow trial functions that fulfil the boundary condition and are of limited practical value.The unrestricted variational principle for the linear oscillations of a magnetic configuration is used to derive a more general formulation of the second-order tensor virial equations so that the linear trial function _i =X _ij x _j can be used to study the oscillations of a configuration with a magnetic field that extends in the exterior vacuum. The unrestricted virial equations have been applied to Ferraro's model and approximate results for the eigenfrequencies and eigenfunctions have been obtained for nine oscillation modes.     

Relative Velocities and Linewidths in a Coronal Hole and Outside

Relative Doppler velocities and spectral linewidths in a coronal hole and in the quiet Sun region outside have been obtained from Solar and Heliospheric Observatory (SOHO)/Coronal Diagnostic Spectrometer (CDS) observations. Five strong emission lines in the CDS wavelength range (namely, O? iii 599 Å, O?v 630 Å, Ne?vi 562.8 Å, He?ii 304 Å, and Mg?ix 368 Å), whose formation temperatures represent different heights in the solar atmosphere from the lower transition region to the inner corona, have been used in the study. As reported earlier, relative velocities in the coronal hole are generally blueshifted with respect to the quiet Sun, and the magnitude of the blueshifts increases with height. It has been found that the polar coronal hole has larger relative velocities than the equatorial extension in the inner corona. Several localized velocity contours have been found mainly on network brightenings and in the vicinity of the coronal hole boundary. The presence of velocity contours on the network may represent network outflows whereas the latter could be due to localized jets probably arising from magnetic reconnection at the boundary. All spectral lines have larger widths in the coronal hole than in the quiet Sun. In O?v 630 Å an extended low-linewidth region is seen in the coronal hole?–?quiet Sun boundary, which may indicate fresh mass transfer across the boundary. Also polar coronal holes have larger linewidths in comparison with the equatorial extension. Together with larger relative velocities, this suggests that the solar wind emanating from polar hole regions is faster than that from equatorial hole regions.     

On the cause of geomagnetic activity at Hale sector boundary

The solar wind velocity and density have been studied around the interplanetary magnetic field sector boundaries (+ to ? and ? to +) during 1965–1974, separating the data into autumnal and vernal periods. It is noticed that the solar wind velocity shows a sharp increase around the Hale type of sector boundary in both northern and southern heliosphere indicating a more favourable condition for the high speed stream after Hale type of sector boundary crossing than non-Hale boundary.     

Study on Two Methods for Nonlinear Force-Free Extrapolation Based on Semi-Analytical Field

In this paper, two semi-analytical solutions of force-free fields (Low and Lou, Astrophys. J. 352, 343, 1990) have been used to test two nonlinear force-free extrapolation methods. One is the boundary integral equation (BIE) method developed by Yan and Sakurai (Solar Phys. 195, 89, 2000), and the other is the approximate vertical integration (AVI) method developed by Song et al. (Astrophys. J. 649, 1084, 2006). Some improvements have been made to the AVI method to avoid the singular points in the process of calculation. It is found that the correlation coefficients between the first semi-analytical field and extrapolated field using the BIE method, and also that obtained by the improved AVI method, are greater than 90% below a height 10 of the 64×64 lower boundary. For the second semi-analytical field, these correlation coefficients are greater than 80% below the same relative height. Although differences between the semi-analytical solutions and the extrapolated fields exist for both the BIE and AVI methods, these two methods can give reliable results for heights of about 15% of the extent of the lower boundary.     

The formation of boundary clinopyroxenes and associated glass veins in type B1 CAIs

Abstract— We used focused ion beam thin section preparation and scanning transmission electron microscopy (FIB/STEM) to examine the interfacial region between spinel and host melilite for spinel grains in type B1 inclusions from the Allende and Leoville carbonaceous chondrites. Boundary clinopyroxenes decorating spinel surfaces have compositions similar to those of coarser clinopyroxenes from the same region of the inclusion, suggesting little movement after formation. Host melilite displays no anomalous compositions near the interface and late‐stage minerals are not observed, suggesting that boundary pyroxenes did not form by crystallization of residual liquid. Allende spinels display either direct spinel‐melilite contact or an intervening boundary clinopyroxene between the two phases. Spinel‐melilite interfacial regions in a Leoville B1 are more complex, with boundary clinopyroxene, as observed in Allende, but also variable amounts of glass, secondary calcite, perovskite, and an Mg‐, Al‐, OH‐rich and Ca‐, Si‐poor crystalline phase that may be a layered double hydrate. One possible scenario of formation for the glass veins is that open system alteration of melilite produced a porous, hydrated aggregate of Mg‐carpholite or sudoite + aluminous diopside that was shock melted and quenched to a glass. The hydrated crystalline phase we observed may have been a shocked remnant of the precursor phase assemblage, but is more likely to have formed later by alteration of the glass. In the mantle, boundary clinopyroxenes may have been crystallized from Ti‐rich liquids formed by the direct dissolution of perovskite and an associated Sc‐Zr‐rich phase or as a reaction product between dissolving perovskite and liquid. In the core, any perovskite and associated Ti‐enriched liquids that may have originally been present disappeared before the growth of boundary clinopyroxene, and the observed boundary clinopyroxene may have nucleated and grown from the liquid, along with the larger core clinopyroxene.     

A re-determination of the boundaries of the cusp

Ariel 4 data (Goodallet al., 1973) of February 4 to May 5 1972 have been used to investigate the cusp boundaries. Data with the same local time (LT) and magnetic local time (MLT) around noontime and at high invariant latitudes (INL) have been used to represent the cusp area. (1) For noontime hours (MLT=1200 and LT=1200) high electron density values are observed. (2) the boundary towards the equator of the high density region starts at 76°–77° INL and the poleward boundary is at about 82° INL. Therefore, high electron densities are observed along an INL interval of 5°–6° at the Ariel 4 heights of 500–600 km, agreeing with other investigations.     

Coronal heating by magnetic reconnection

The theory of magnetic reconnection has advanced substantially over the past few years. There now exists a new generation of fast two-dimensional models known as almost-uniform reconnection and nonuniform reconnection, depending on the boundary conditions. Also, we are beginning to explore the uncharted region of three-dimensional reconnection, where regimes of “spine reconnection” and “fan reconnection” have been discovered. Furthermore, part of the coronal heating problem appears to have been solved with recent observational support for the Converging Flux Model in which heating is produced by coronal reconnection driven by footpoint motions.     

Radiation transfer through atmospheric aerosol media with anisotropic scattering

Radiation transfer in atmospheric aerosol media with general boundary conditions has been studied for anisotropic scattering. The considered aerosol medium assumed to have specular and diffused reflecting boundary surfaces and in the presence of internal source. The radiation transfer scattering parameters as single scattering albedo, asymmetry factor, scattering, absorption, extinction efficiencies and anisotropic scattering coefficient have been calculated using the Mie theory. The problem with general boundary conditions is solved in terms of the solution of source-free problem with simply boundary conditions. Pomraning-Eddington approximation is used to solve the source-free problem. For the sake of comparison, a weight function is introduced and used in two special forms. The calculated partial heat fluxes with the two methods are compared and showed good agreement. Some of our results are found in a good agreement with published data.     

A search for shocked quartz grains in the Allerød‐Younger Dryas boundary layer

The Younger Dryas impact hypothesis suggests that multiple airbursts or extraterrestrial impacts occurring at the end of the Allerød interstadial resulted in the Younger Dryas cold period. So far, no reproducible, diagnostic evidence has, however, been reported. Quartz grains containing planar deformation features (known as shocked quartz grains), are considered a reliable indicator for the occurrence of an extraterrestrial impact when found in a geological setting. Although alleged shocked quartz grains have been reported at a possible Allerød‐Younger Dryas boundary layer in Venezuela, the identification of shocked quartz in this layer is ambiguous. To test whether shocked quartz is indeed present in the proposed impact layer, we investigated the quartz fraction of multiple Allerød‐Younger Dryas boundary layers from Europe and North America, where proposed impact markers have been reported. Grains were analyzed using a combination of light and electron microscopy techniques. All samples contained a variable amount of quartz grains with (sub)planar microstructures, often tectonic deformation lamellae. A total of one quartz grain containing planar deformation features was found in our samples. This shocked quartz grain comes from the Usselo palaeosol at Geldrop Aalsterhut, the Netherlands. Scanning electron microscopy cathodoluminescence imaging and transmission electron microscopy imaging, however, show that the planar deformation features in this grain are healed and thus likely to be older than the Allerød‐Younger Dryas boundary. We suggest that this grain was possibly eroded from an older crater or distal ejecta layer and later redeposited in the European sandbelt. The single shocked quartz grain at this moment thus cannot be used to support the Younger Dryas impact hypothesis.     

Boundary layer plasmas as a source for high-latitude, early afternoon, auroral arcs

Simultaneous measurements of hot boundary layer plasma from PROGNOZ-7 and particle precipitation from the TIROS/NOAA satellite in nearly magnetically conjugate regions have been used to study the dynamo process responsible for the formation of high latitude, early afternoon, auroral arcs.

Characteristic for the PROGNOZ-7 observations in the dayside boundary layer at high latitudes is the frequent occurrence of regions with injected magnetosheath plasma embedded in a “halo” of antisunward flowing magnetosphere plasma. The injected magnetosheath plasma have several features which indicate that it also acts as a local source of EMF in the boundary layer. The process resembles that of a local MHD dynamo driven by the excess drift velocity of the injected magnetosheath plasma relative to the background magnetospheric plasma.

The dynamo region is capable of driving field-aligned currents that couple to the ionosphere, where the upward current is associated with the high latitude auroral arcs.

We demonstrate that the large-scale morphology as well as the detailed data intercomparison between PROGNOZ-7 and TIROS-N both agree well with a local injection of magnetosheath plasma into the dayside boundary layer as the main dynamo process powering the high-latitude, early afternoon auroral arcs.     

Dynamical problems of thermoelastic solids

The dynamical problem for thermal stresses in an infinite isotropic elastic cylinder of radius a with its axis along the z-axis, subject to fixed boundary conditions is studied. The Fourier heat conduction equation has been solved applying the Fourier transform and the theory of complex variable. The thermoelastic equation of motion has been separated into two wave equations which can be solved separately. The temperature, the displacement and the stress components have been obtained in analytical form as series involving Bessel function of first kind and of order zero.     

Planetary dynamics in the system α Centauri: The stability diagrams

The stability of the motion of a hypothetical planet in the binary system ?? Cen A?CB has been investigated. The analysis has been performed within the framework of a planar (restricted and full) three-body problem for the case of prograde orbits. Based on a representative set of initial data, we have obtained the Lyapunov spectra of the motion of a triple system with a single planet. Chaotic domains have been identified in the pericenter distance-eccentricity plane of initial conditions for the planet through a statistical analysis of the data obtained. We have studied the correspondence of these chaotic domains to the domains of initial conditions that lead to the planet??s encounter with one of the binary??s stars or to the escape of the planet from the system. We show that the stability criterion based on the maximum Lyapunov exponent gives a more clear-cut boundary of the instability domains than does the encounterescape criterion at the same integration time. The typical Lyapunov time of chaotic motion is ??500 yr for unstable outer orbits and ??60 yr for unstable inner ones. The domain of chaos expands significantly as the initial orbital eccentricity of the planet increases. The chaos-order boundary has a fractal structure due to the presence of orbital resonances.     

Lineated valley fill (LVF) and lobate debris aprons (LDA) in the Deuteronilus Mensae northern dichotomy boundary region, Mars: Constraints on the extent, age and episodicity of Amazonian glacial events

In order to assess the nature, degradational processes and history of the dichotomy boundary on Mars, we conducted a detailed morphological analysis of a 70,000 km² region of its northern portion (north-central Deuteronilus Mensae, south of Lyot, in the vicinity of Sinton Crater). This region is characterized by the distinctive sinuous ∼2 km-high plateau scarp boundary, outlying massifs to the north, and extensive fretted valleys dissecting the plateau to the south. These features represent the first-order modification and retreat of the dichotomy boundary, and are further modified by processes that form lineated valley fill (LVF) in the fretted valleys, and lobate debris aprons (LDA) along the dichotomy scarp and surrounding the outlying massifs. We use new high-resolution image and topography data to examine the nature and origin of LVF and LDA and to investigate the climatic and accompanying degradational history of the escarpment. On the basis of our analysis, we conclude that: (1) LVF and LDA deposits within the study region are comprised of the same material, show integrated flow patterns, and originate as debris-covered valley glaciers; a significant amount of ice (hundreds of meters) is likely to remain today beneath a thin cover of sublimation till. (2) There is depositional evidence to suggest glacial highstands at least 800 m above the present level, implying previous conditions in which the distribution of ice was much more widespread; this is supported by similar deposits within many other areas across the dichotomy boundary. (3) The timing of the most recent large-scale activity of the LDA/LVF in this area is about 100-500 million years ago, similar to ages reported elsewhere along the dichotomy boundary. (4) There is evidence for a secondary, but significantly limited phase of glaciation; the deposits of which are limited to the vicinity of the alcoves; similar later phases have also been reported elsewhere along the dichotomy boundary. (5) Modification of the fretted valleys of the dichotomy boundary has been substantial locally, but we find no evidence that the Amazonian glacial epochs caused retreat of the dichotomy boundary of the scale of tens to hundreds of kilometers. Our findings support the results of an analysis just to the east of the study region and of studies carried out elsewhere along the dichotomy boundary that find further evidence for the remnants of debris-covered glaciers and extensive valley glacial land systems.     

Propagation of solar-flare cosmic rays along corotating interplanetary flux-tubes

The propagation of solar-flare cosmic rays along a corotating spiral mean interplanetary magnetic flux-tube is investigated under a variety of initial, inner and outer boundary conditions, by solving the transport equations numerically. The inner boundary condition is shown to have negligible effects on the redistribution of the cosmic rays in interplanetary space, apart from a proportionality constant. The diffusion coefficients have been classified into type I and type II: for type I the peakalong-spiral stays at a finite heliocentric distance, while for type II it travels to infinity. The time dependence of the anisotropy vector is characteristically different in each type and in the case of type I is a function of the observation point. The general propagation characteristics with a free-escape outer boundary are like those of type I. It is shown that closely exponential decay can be produced without a free-escape boundary. The concept of an ‘equilibrium shape’ at late time is introduced. With a source that decays slowly, a maximum and minimum pair appears in the spatial distribution; the maximum subsequently propagates outward and the minimum inward. The impulsive models given do not reproduce satisfactorily all the observed characteristics of solar burst events and a solar source extended in time appears to be indicated.