Interaction of interplanetary shocks with the bow shock

Fast forward interplanetary (IP) shocks have been identified as a source of large geomagnetic disturbances. However, the shocks can evolve in the solar wind, they are modified by interaction with the bow shock and during their propagation through the magnetosheath. A few previous papers refer the inclination and deceleration of the IP shock front in this region. Our contribution continues this effort and presents the study of an IP shock interaction with the bow shock. Since the bow shock is a reversed fast shock, the interaction of the IP shock and bow shock is a problem of interaction of two fast MHD shocks.

We compare profiles of magnetic field and plasma parameters observed by several spacecraft in the solar wind and magnetosheath with the profiles of the same parameters resulting from the MHD numerical model. The MHD model suggests that the interaction of an IP shock with the bow shock results in an inward bow shock displacement that is followed by its outward motion. Such motion will result in an indentation propagating along the bow shock surface. This scenario is confirmed by multipoint observations. Moreover, the model confirms also previous suggestions on the IP shock deceleration in the magnetosheath.     

Trajectories of ions specularly reflected from non-planar collisionless shocks

An important contribution to the thermalization of the solar wind ions at the Earth's bow shock for high Mach numbers comes from the reflection of a fraction of these ions from the shock. Previous studies have examined the trajectories of the reflected ions assuming the shock to be an infinite plane. In this paper a model is developed to describe the trajectories of particles after reflection for a variety of shock geometries. Of particular interest are the initial conditions which allow the particle to return to the shock with a greater normal velocity than at first encounter, or to return to the shock at all. The effects of the magnetic field direction and the curvature of the shock on particle trajectories are discussed for cylindrical and spherical shock geometries and compared to those for a planar shock.     

II型超新星瞬发爆炸模型的激波传播

本文专门讨论II型超新星瞬发爆炸模型中的激波形成点、激波波阵面的位置、激波速度和激波能量等一系列的问题 .研究发现在外星核区除存在首级激波外 ,还存在次级激波 ;激波波阵面的准确位置应在熵变化的最大点和外星核区中人为粘性压的最大点之间 ;并给出激波能量的不同定义量之间的大小关系     

Numerical simulation of slow shocks in the solar wind

The evolutionary state of slow forward shock waves is examined with the use of two MHD numerical codes. Our study is intended to be exploratory rather than a detailed parametric one. The first code is one-dimensional (with three components of velocity and magnetic field) which is used to follow a slow shock that propagates into a positive gradient of density versus distance. It is found that the slow shock evolves into an extraneous (intermediate) shock wave. The second code has a spherical, one-dimensional, planar geometry (with two velocity and magnetic field components) which is used to follow a spiral interplanetary magnetic field. It is found that a slow shock type perturbation can generate a forward slow shock; a fast forward shock is generated in the front of the slow shock; a contact discontinuity is formed behind the slow shock, and a compound nonlinear MHD wave is formed behind the contact discontinuity with a fast reverse shock formed further behind. Thus, we demonstrate that the evolution of a slow shock into (solely) a fast shock, as suggested by Whang (1987), is much more complicated.     

Experimental investigation of shock metamorphic effects in a metapelitic granulite: The importance of shock impedance contrast between components

Shock recovery experiments were performed at 12.5, 25, 34, 40, and 56 GPa at 25 °C, and at 18 and 25 GPa at 400 °C, on a high‐grade, migmatitic, garnet‐cordierite metapelite from the Etivé aureole, Scotland. Objectives for this study were to (1) characterize shock effects in a complex polymineralic rock with a significant proportion of hydrous ferromagnesian minerals, both as a function of variable shock pressure and preshock temperature, and (2) to explore the effects of shock impedance contrast between component minerals on the respective abundances and distribution of these features. At any shock pressure, the order of decreasing intensity of shock metamorphic effects in component phases is: cordierite (Crd)→biotite (Bt)→plagioclase (Pl)→K‐feldspar (Kfs)→quartz (Qtz)→garnet (Grt)→orthopyroxene (Opx). Samples shocked to pressures below 40 GPa (25 °C) were typically characterized by marked heterogeneous distribution of shock effects on both intragranular and intergranular scales. Shock heterogeneity is mainly attributed to shock impedance contrast between contiguous phases, and manifests as shock amplification locally where shock impedance contrast is greatest, and shock suppression where impedance contrast is least. The heterogeneous distribution of shock metamorphic effects in both experiments and natural rocks is a signature of extreme disequilibrium at the submillimeter scale. The heterogeneous distribution of shock metamorphic effects mitigates against the use of shock effects in minerals exclusively as regional shock pressure barometers, and ought to be augmented by additional constraints on shock pressure from numerical models.     

Shock re‐equilibration of fluid inclusions in crystalline basement rocks from the Ries crater,Germany

Abstract— This study examines the effects of shock metamorphism on fluid inclusions in crystalline basement target rocks from the Ries crater, Germany. The occurrence of two‐phase fluid inclusions decreases from shock stage 0 to shock stage 1, while single‐phase inclusions increase, likely as a result of re‐equilibration. In shock stages 2 and 3, both two‐phase and single‐phase inclusions decrease with increasing shock stage, indicating that fluid inclusion vesicles are destroyed due to plastic deformation and phase changes in the host minerals. However, quartz clasts entrained in shock stage 4 melts contain both single‐phase and two‐phase inclusions, demonstrating the rapid quenching of the melt and the heterogeneous nature of impact deformation. Inclusions in naturally shocked polycrystalline samples survive at higher shock pressures than those in single crystal shock experiments. However, fluid inclusions in both experimental and natural samples follow a similar trend in re‐equilibration at low to moderate shock pressures leading to destruction of inclusion vesicles in higher shock stages. This suggests that shock processing may lead to the destruction of fluid inclusions in many planetary materials and likely contributed to shock devolatilization of early planetesimals.     

Natural occurrence of reidite in the Xiuyan crater of China

The high‐pressure minerals of reidite and coesite have been identified in the moderately shock‐metamorphosed gneiss (shock stage II, 35–45 GPa) and the strongly shock‐metamorphosed gneiss (shock stage III, 45–55 GPa), respectively, from the polymict breccias of the Xiuyan crater, a simple impact structure 1.8 km in diameter in China. Reidite in the shock stage II gneiss displays lamellar textures developed in parental grains of zircon. The phase transformation of zircon to reidite likely corresponds to a martensitic mechanism. No coesite is found in the reidite‐bearing gneiss. The shock stage III gneiss contains abundant coesite, but no reidite is identified in the rock. Coesite occurs as acicular, dendritic, and spherulitic crystals characteristic of crystallization from shock‐produced silica melt. Zircon in the rock is mostly recrystallized. The postshock temperature in the shock stage III gneiss is too high for the preservation of reidite, whereas reidite survives in the shock stage II gneiss because of relatively low postshock temperature. Reidite does not occur together with coesite because of difference in shock‐induced temperature between the shock stage II gneiss and the shock stage III gneiss.     

The Injection, Acceleration, and Dynamical Influence of Interstellar Pickup Ions at the Solar Wind Termination Shock

The interaction of interstellar pickup ions with the solar wind termination shock is reviewed and assessed. The pickup ions mass and momentum load the wind and increase its pressure, effects which decrease the strength of the shock and its distance from the Sun. The pickup hydrogen may contribute substantially to the "reflected" ion population, which should provide most of the dissipation at the supercritical quasi-perpendicular shock. A fraction of the pickup ions impinging on the shock is "injected" into the process of diffusive shock acceleration to form the anomalous cosmic ray component. An injection mechanism which accounts for the apparent absence of solar wind ions in the anomalous component is "shock surfing", in which pickup ions which approach the shock slowly may be trapped between the upstream Lorentz force and the shock potential and accelerated in the motional electric field beyond the energy threshold for diffusive shock acceleration. However, the simplest interpretation of shock surfing would favor less massive pickup ion species, in contradiction with Voyager observations of anomalous component composition. A possible extension of the shock surfing mechanism is considered, as well as other injection mechanisms. Finally, the pressure of the anomalous component may modify the structure of the termination shock, which in turn may influence injection rates. This revised version was published online in July 2006 with corrections to the Cover Date.     

Shock surfing acceleration

Analytical and numerical analysis identify shock surfing acceleration as an ideal pre-energization mechanism for the slow pick-up ions at quasiperpendicular shocks. After gaining sufficient energy by shock surfing, pick-up ions undergo diffusive acceleration to reach their observed energies. Energetic ions upstream of the cometary bow shock, acceleration of solar energetic particles by magnetosonic waves in corona, ion enhancement in interplanetary shocks, generation of anomalous cosmic rays from interstellar pick-up ions at the termination shock are some of the cases where shock surfing acceleration apply. Inclusion of the lower-hybrid wave turbulence into the laminar model of shock surfing can explain the preferential acceleration of heavier particles as observed by Voyager at the termination shock. At relativistic energies, unlimited acceleration of ions is theoretically possible; because for sufficiently strong shocks main limitation of the mechanism, caused by the escape of accelerated particles downstream of the shock during acceleration no longer exists.     

A shock hodograph in collisionless plasma

In this paper, the shock hodograph (polar) for a collisionless transverse shock has been developed in the plane of the flow deflection angle and total pressure jump. The sonic point on the hodograph lies closer to the characteristic than to the point of maximum flow deflection for an attached shock. This hodograph is particularly useful in the analysis of three shock confluences and refraction of shock waves at gas interfaces. The first analysis is fully described in this paper. It is observed that the third shock wave is forward facing. It is pertinent to note that the limited region of supersonic flow also restricts the occurrences of three shock interactions.     

Statistical inference in the presence of an inclination effect in laboratory radiative shock experiments

A laser-driven experiment produces images of dense shocked material by x-ray transmission. The post-shock material is sufficiently dense that no significant signal passes through the dense layer, and therefore the shock compression cannot be directly measured by comparing transmitted intensities. One could try to determine the shock compression ratio by observing the ratio of the total distance travelled by the shock to the dense post-shock layer width, but small deviations of the angle of the shock with respect to the angle of imaging create large asymmetric errors in observation. A statistical approach to recovering shock compression by appropriately combining data from several experiments is developed, using fits to a simple model for the shock and shock tube geometry.     

Shock metamorphism of quartz in nature and experiment: II. Significance in geoscience*

Abstract— The occurrence of shock metamorphosed quartz is the most common petrographic criterion for the identification of terrestrial impact structures and lithologies. Its utility is due to its almost ubiquitous occurrence in terrestrial rocks, its overall stability and the fact that a variety of shock metamorphic effects, occurring over a range of shock pressures, have been well documented. These shock effects have been generally duplicated in shock recovery experiments and, thus, serve as shock pressure barometers. After reviewing the general character of shock effects in quartz, the differences between experimental and natural shock events and their potential effects on the shock metamorphism of quartz are explored. The short pulse lengths in experiments may account for the difficulty in synthesizing the high-pressure polymorphs, coesite and stishovite, compared to natural occurrences. In addition, post-shock thermal effects are possible in natural events, which can affect shock altered physical properties, such as refractive index, and cause annealing of shock damage and recrystallization. The orientations of planar microstructures, however, are unaffected by post-impact thermal events, except if quartz is recrystallized, and provide the best natural shock barometer in terms of utility and occurrence. The nature of planar microstructures, particularly planar deformation features (PDFs), is discussed in some detail and a scheme of variations in orientations with shock pressure is provided. The effect of post-impact events on PDFs is generally limited to annealing of the original glass lamellae to produce decorated PDFs, resulting from the exsolution of dissolved water during recrystallization. Basal (0001) PDFs differ from other PDF orientations in that they are multiple, mechanical Brazil twins, which are difficult to detect if not partially annealed and decorated. The occurrence and significance of shock metamorphosed quartz and its other phases (namely, coesite, stishovite, diaplectic glass and lechatelierite) are discussed for terrestrial impact structures in both crystalline (non-porous) and sedimentary (porous) targets. The bulk of past studies have dealt with crystalline targets, where variations in recorded shock pressure in quartz have been used to constrain aspects of the cratering process and to estimate crater dimensions at eroded structures. In sedimentary targets, the effect of pore space results in an inhomogeneous distribution in recorded shock pressure and temperature, which requires a different classification scheme for the variation of recorded shock compared to that in crystalline targets. This is discussed, along with examples of variations in the relative abundances of planar microstructures and their orientations, which are attributed to textural variations in sedimentary target rocks. Examples of the shock metamorphism of quartz in distal ejecta, such as at the K/T boundary, and from nuclear explosions are illustrated and are equivalent to that of known impact structures, except with respect to characteristics that are due to long-term, post-shock thermal effects. Finally, the differences between the deformation and phase transformation of quartz by shock and by endogenic, tectonic and volcanic processes are discussed. We confirm previous conclusions that they are completely dissimilar in character, due to the vastly different physical conditions and time scales typical for shock events, compared to tectonic and volcanic events. Well-characterized and documented shock effects in quartz are unequivocal indicators of impact in the natural environment.     

Laboratory Observation of Secondary Shock Formation Ahead of a Strongly Radiative Blast Wave

We have previously reported the experimental discovery of a second shock forming ahead of a radiative shock propagating in Xe. The initial shock is spherical, radiative, with a high Mach number, and it sends a supersonic radiative heat wave far ahead of itself. The heat wave rapidly slows to a transonic regime and when its Mach number drops to two with respect to the downstream plasma, the heat wave drives a second shock ahead of itself to satisfy mass and momentum conservation in the heat wave reference frame. We now show experimental data from a range of mixtures of Xe and N₂, gradually changing the properties of the initial shock and the environment into which the shock moves and radiates (the radiative conductivity and the heat capacity). We have successfully observed second shock formation over the entire range from 100% Xe mass fraction to 100% N₂. The formation radius of the second shock as a function of Xe mass fraction is consistent with an analytical estimate.     

Laboratory Simulations of Supernova Shockwave Propagation

Supernovae launch spherical shocks into the circumstellar medium (CSM). These shocks have high Mach numbers and may be radiative. We have created similar shocks in the laboratory by focusing laser pulses onto the tip of a solid pin surrounded by ambient gas; ablated material from the pin rapidly expands and launches a shock through the surrounding gas. Laser pulses were typically 5 ns in duration with ablative energies ranging from 1–150 J. Shocks in ambient gas pressures of ~1 kPa were observed at spatial scales of up to 5 cm using optical cameras with schlieren. Emission spectroscopy data were obtained to infer electron temperatures (< 10 eV). In this experiment we have observed a new phenomena; at the edge of the radiatively heated gas ahead of the shock, a second shock forms. The two expanding shocks are simultaneously visible for a time, until the original shock stalls from running into the heated gas. The second shock remains visible and continues to expand. A minimum condition for the formation of the second shock is that the original shock is super-critical, i.e., the temperature distribution ahead of the original shock has an inflexion point. In a non-radiative control experiment the second shock does not form. We hypothesize that a second shock could form in the astrophysical case, possibly in radiative supernova remnants such as SN1993J, or in shock-CSM interaction.     

Shock‐darkening in ordinary chondrites: Determination of the pressure‐temperature conditions by shock physics mesoscale modeling

We determined the shock‐darkening pressure range in ordinary chondrites using the iSALE shock physics code. We simulated planar shock waves on a mesoscale in a sample layer at different nominal pressures. Iron and troilite grains were resolved in a porous olivine matrix in the sample layer. We used equations of state (Tillotson EoS and ANEOS) and basic strength and thermal properties to describe the material phases. We used Lagrangian tracers to record the peak shock pressures in each material unit. The post‐shock temperatures (and the fractions of the tracers experiencing temperatures above the melting point) for each material were estimated after the passage of the shock wave and after the reflections of the shock at grain boundaries in the heterogeneous materials. The results showed that shock‐darkening, associated with troilite melt and the onset of olivine melt, happened between 40 and 50 GPa with 52 GPa being the pressure at which all tracers in the troilite material reach the melting point. We demonstrate the difficulties of shock heating in iron and also the importance of porosity. Material impedances, grain shapes, and the porosity models available in the iSALE code are discussed. We also discuss possible not‐shock‐related triggers for iron melt.     

Effect of azimuthal magnetic field on the propagation of cylindrical shock waves in self-gravitating gas

A propagation of diverging cylindrical shock in a self-gravitating gas, having an initial density and azimuthal magnetic field distributions variable, has been studied for the two cases (i) when the shock is weak and (ii) when it is strong. Analytical relations for shock velocity and shock strength have been obtained. Lastly, the expressions for the pressure, the density and the particle velocity immediately behind the shock have been also obtained for both cases.     

Quantitative in situ XRD measurement of shock metamorphism in Martian meteorites using lattice strain and strain‐related mosaicity in olivine

All Martian meteorites have experienced shock metamorphism to some degree. We quantitatively determined shock‐related strain in olivine crystals to measure shock level and peak shock pressure experienced by five Martian meteorites. Two independent methods employing nondestructive in situ micro X‐ray diffraction (μXRD) are applied, i.e., (1) the lattice strain method, in which the lattice strain value (ε) for each olivine grain is derived from a Williamson–Hall plot using its diffraction pattern (peak width variation with diffraction angle) with reference to a best fit calibration curve of ε values obtained from experimentally shocked olivine grains; (2) the strain‐related mosaicity method, allowing shock stage to be estimated by measuring the streaking along the Debye rings of olivine grain diffraction spots to define their strain‐related mosaic spread, which can then be compared with olivine mosaicity in ordinary chondrites of known shock stage. In this study, both the calculated peak shock pressures and the estimated shock stages for Dar al Gani 476 (45.6 ± 0.6 GPa), Sayh al Uhaymir 005/8 (46.1 ± 2.2 GPa), and Nakhla (18.0 ± 0.6 GPa) compare well with literature values. Formal shock assessments for North West Africa 1068/1110 (53.9 ± 2.1 GPa) and North West Africa 6234 (44.6 ± 3.1 GPa) have not been reported within the literature; however, their calculated peak shock pressures fall within the range of peak shock pressures defining their estimated shock stages. The availability of nondestructive and quantitative μXRD methods to determine shock stage and peak shock pressure from olivine crystals provides a key tool for shock metamorphism analysis.     

Energetic particle fluxes at heliospheric shocks: Evidences of superdiffusion and comparison between analytical and numerical modeling

Energetic particles fluxes measured by spacecraft in the heliosphere are frequently observed to peak during interplanetary shock crossings, suggesting the shock to be the source of acceleration. It has been shown that the shape of the energetic particle fluxes, upstream and downstream of the shock, is affected by energetic particle transport properties. In this study we make a comparison among a number of shock crossings observed by the ACE spacecraft, and the energetic particle fluxes derived by a test-particle numerical model in the vicinity of a planar shock. We find that observations are in good agreement with a particle density profile obtained in the simulation by assuming superdiffusive transport both upstream and downstream of the shock region.     

Analytical solution of magnetogasdynamic cylindrical shock waves in self-gravitating and rotating gas,II

Using the C.C.W. method, propagation of diverging cylindrical shock wave in a self-gravitating and rotating gas under the influence of a constant axial magnetic field has been studied for two cases of weak and strong shocks. Medium ahead of the shock is supposed to be homogeneous. Analytical relations for shock velocity and shock strength along with the expressions for the pressure, density, and particle velocity just behind the shock wave have been also obtained for both cases.     

Cosmic-ray diffusive acceleration at shock waves with finite upstream and downstream escape boundaries

In the present paper we discuss the modifications introduced into the first-order Fermi shock acceleration process due to a finite extent of diffusive regions near the shock or due to boundary conditions leading to an increased particle escape upstream and/or downstream of the shock. In the simple example of the planar shock wave considered we idealize the escape phenomenon by imposing a particle escape boundary at some distance from the shock. The presence of such a boundary (or boundaries) leads to coupled steepening of the accelerated particle spectrum and decreasing of the acceleration time scale. It allows for a semi-quantitative evaluation and, in some specific cases, also for modelling of the observed steep particle spectra as a result of the first-order Fermi shock acceleration. We also note that the particles close to the upper energy cut-off are younger than the estimate based on the respective acceleration time scale. In Appendix A we present a new time-dependent solution for infinite diffusive regions near the shock allowing for different constant diffusion coefficients upstream and downstream of the shock.