Surface latent heat flux anomalies quasi-synchronous with ionospheric disturbances before the 2007 Pu’er earthquake in China

One of various mechanisms of pre-earthquake lithosphere–atmosphere–ionosphere coupling as possible explanation of the seismo-ionospheric effects is connected with the release of latent heat. Abnormal variations of ionospheric electromagnetic parameters possibly related to the 2007 Ms 6.4 Pu’er earthquake in China were reported. This paper attempts to examine whether there were abnormal changes of surface latent heat flux (SLHF) linked with this pre-earthquake ionospheric disturbances. The spatio-temporal statistical analyzes of multi-years SLHF data from USA NCEP/NCAR Reanalysis Project reveal that local SLHF enhancements appeared 11, 10 and 7 days before the Pu’er earthquake, respectively. As contrasted to the formerly reported local ionospheric Ne enhancement 9 days before the shock observed by DEMETER satellite, it is discovered that the SLHF and Ne anomalies are quasi-synchronous and have good spatial correspondence with the epicentre and the local active faults. This is valuable for understanding the seismogenic coupling processes and for recognizing earthquake anomaly with multiple parameters from integrated Earth observation system.     

Model for the VLF/LF radio signal anomalies formation associated with earthquakes

The influence of quasi-static electric field of seismic origin on the characteristics of the internal gravity waves (IGWs) in the Earth’s ionosphere is considered. The electric field in the ionosphere arises due to the injection of charged aerosols into the atmosphere, formation of an EMF in the near Earth atmosphere and perturbation of the conductive electric current in the global electric circuit. Amplification of the electric current in seismic zone is accompanied by the formation of perturbation of the lower ionosphere that affects the amplitude and phase of VLF/LF signals. The action of the electric field on the IGWs is connected with the appearance of the Ampere’s force in the ionosphere. In the spectral range of these waves the latter acts on the neutral component of the ionosphere plasma. As the result of this interaction the ionosphere starts to support the discrete spectrum of oscillations. Periods of their maximums increase as numbers of natural sequence. The existence of such peculiarities of the waves in the ionosphere is confirmed by observations.     

Long-term (50 years) measurements of cosmic ray fluxes in the atmosphere

Since the middle of 1957 till present time the group of researchers of P.N. Lebedev Physical Institute of the Russian Academy of Sciences has carried out the regular balloon borne measurements of charged particle fluxes in the atmosphere. The measurements are performed at polar (northern and southern) and middle latitudes and cover the interval of heights from the ground level up to 30–35 km. Standard detectors of particles (gas-discharged counters) have been used. More than 80,000 measurements of cosmic ray fluxes in the atmosphere have been performed to the present time. In the data analysis the geomagnetic field and the Earth’s atmosphere are used as cosmic ray spectrometers.     

Three-dimensional magnetic reconnection regimes: A review

The magnetic field in many astrophysical plasmas – such as the solar corona and Earth’s magnetosphere – has been shown to have a highly complex, three-dimensional structure. Recent advances in theory and computational simulations have shown that reconnection in these fields also has a three-dimensional nature, in contrast to the widely used two-dimensional (or 2.5-dimensional) models. Here we discuss the underlying theory of three-dimensional magnetic reconnection. We also review a selection of new models that illustrate the current state of the art, as well as highlighting the complexity of energy release processes mediated by reconnection in complicated three-dimensional magnetic fields.     

Remote sensing using GNSS signals: Current status and future directions

The refracted, reflected and scattered signals of global navigation satellite systems (GNSS) have been successfully used to remotely sense the Earth’s surface and atmosphere. It has demonstrated its potential to sense the atmosphere and ionosphere, ocean, land surfaces (including soil moisture) and the cryosphere. These new measurements, although in need of refinement and further validation in many cases, can be used to complement existing techniques and sensors, e.g., radiosonde, ionosonde, radar altimetry and synthetic aperture radar (SAR). This paper presents the current status and new developments of remote sensing using GNSS signals as well as its future directions and applications. Some notable emerging applications include monitoring sea ice, dangerous sea states, ocean eddy and storm surges. With the further improvement of the next generation multi-frequency GNSS systems and receivers and new space-based instruments utilizing GNSS reflections and refractions, new scientific applications of GNSS are expected in various environment remote sensing fields in the near future.     

Pion and electromagnetic contribution to dose: Comparisons of HZETRN to Monte Carlo results and ISS data

Recent work has indicated that pion production and the associated electromagnetic (EM) cascade may be an important contribution to the total astronaut exposure in space. Recent extensions to the deterministic space radiation transport code, HZETRN, allow the production and transport of pions, muons, electrons, positrons, and photons. In this paper, the extended code is compared to the Monte Carlo codes, Geant4, PHITS, and FLUKA, in slab geometries exposed to galactic cosmic ray (GCR) boundary conditions. While improvements in the HZETRN transport formalism for the new particles are needed, it is shown that reasonable agreement on dose is found at larger shielding thicknesses commonly found on the International Space Station (ISS). Finally, the extended code is compared to ISS data on a minute-by-minute basis over a seven day period in 2001. The impact of pion/EM production on exposure estimates and validation results is clearly shown. The Badhwar–O’Neill (BO) 2004 and 2010 models are used to generate the GCR boundary condition at each time-step allowing the impact of environmental model improvements on validation results to be quantified as well. It is found that the updated BO2010 model noticeably reduces overall exposure estimates from the BO2004 model, and the additional production mechanisms in HZETRN provide some compensation. It is shown that the overestimates provided by the BO2004 GCR model in previous validation studies led to deflated uncertainty estimates for environmental, physics, and transport models, and allowed an important physical interaction (π/EM) to be overlooked in model development. Despite the additional π/EM production mechanisms in HZETRN, a systematic under-prediction of total dose is observed in comparison to Monte Carlo results and measured data.     

Planet temperatures with surface cooling parameterized

A semigray (shortwave and longwave) surface temperature model is developed from conditions on Venus, Earth and Mars, where the greenhouse effect is mostly due to carbon dioxide and water vapor. In addition to estimating longwave optical depths, parameterizations are developed for surface cooling due to shortwave absorption in the atmosphere, and for convective (sensible and latent) heat transfer. An approximation to the Clausius–Clapeyron relation provides water–vapor feedback. The resulting iterative algorithm is applied to three “super-Earths” in the Gliese 581 system, including the “Goldilocks” planet g (Vogt et al., 2010). Surprisingly, none of the three appear habitable. One cannot accurately locate a star’s habitable zone without data or assumptions about a planet’s atmosphere.     

Ion-acoustic waves in unmagnetized collisionless weakly relativistic plasma of warm-ion and isothermal-electron using time-fractional KdV equation

Collisionless unmagnetized plasma consisting of a mixture of warm ion-fluid and isothermal-electron is considered, assuming that the ion flow velocity has a weak relativistic effect. The reductive perturbation method has been employed to derive the Korteweg–de Vries (KdV) equation for small – but finite-amplitude electrostatic ion-acoustic waves in this plasma. The semi-inverse method and Agrawal’s method lead to the Euler–Lagrange equation that leads to the time fractional KdV equation. The variational-iteration method given by He is used to solve the derived time fractional KdV equation. The calculations show that the fractional order may play the same rule of higher order dissipation in KdV equation to modulate the soliton wave amplitude in the plasma system. The results of the present investigation may be applicable to some plasma environments, such as space-plasmas, laser-plasma interaction, plasma sheet boundary layer of the earth’s magnetosphere, solar atmosphere and interplanetary space.     

EUV observation from the Earth-orbiting satellite,EXCEED

An Earth-orbiting small satellite “EXtreme ultraviolet spectrosCope for ExosphEric Dynamics” (EXCEED) which will be launched in 2012 is under development. The mission will carry out spectroscopic and imaging observation of EUV (Extreme Ultraviolet: 60–145 nm) emissions from tenuous plasmas around the planets (Venus, Mars, Mercury, and Jupiter). It is essential for EUV observation to put on an observing site outside the Earth’s atmosphere to avoid the absorption. It is also essential that the detection efficiency must be very high in order to catch the faint signals from those targets. In this mission, we employ cesium iodide coated microchannel plate as a 2 dimensional photon counting devise which shows 1.5–50 times higher quantum detection efficiency comparing with the bared one. We coat the surface of the grating and entrance mirror with silicon carbides by the chemical vapor deposition method in order to archive the high diffraction efficiency and reflectivity. The whole spectrometer is shielded by the 2 mm thick stainless steel to prevent the contamination caused by the high energy electrons from the inner radiation belt. In this paper, we will introduce the mission overview, its instrument, and their performance.     

Consequences of the application of the streamer fluid model to the study of the sprite inception mechanism

We present the results of a streamer-fluid model used to investigate the electrodynamical coupling between the troposphere and upper atmosphere due to the penetration of lightning electric fields into the mesosphere and the lower ionosphere, generating sprites. The model solves the continuity equation for electrons and ions coupled to Poisson equation. The dominant physical response of the atmosphere is the formation of a screening-ionization wave. The wave shields the atmosphere above it from the action of the lightning field and, together with the conductivity reduction below it due to attachment, the wave amplifies the total field below it, allowing for the penetration of intense electric fields in the mesosphere as it propagates downwards into regions of higher density that compress the wave. This is the key physical mechanism for sprite inception. We evaluated the effects of the thundercloud charge geometry, lightning current waveshape, atmospheric conductivity, via different electron density profiles, and the effect of ionization, attachment and electron mobility coefficients in the electrical breakdown process, related to halo production, and sprite streamer initiation. The results showed that electrons with higher mobility are more efficient in shielding the lightning electric field before breakdown, causing delay, and they contribute to the formation of the streamer seed after breakdown, anticipating the sprite streamer inception. Similarly, a higher effective ionization rate, produced by modifications in the attachment and ionization coefficients, anticipates sprite inception. The simulations with 6 different electron density profiles, and therefore conductivities, spanning 4 orders of magnitude, showed that the altitude of breakdown and sprite initiation, as well as their time delays from the lightning discharge are directly related to atmospheric conductivity: higher conductivities produce halo and sprite inception at lower altitudes with longer delays and may hinder sprite formation. We document that variations of 30 times in the lightning current leads to sprite initiation altitudes in the range 66.0–73.5 km, with delays between 1.550 and 34.500 ms, while variations of 4 orders of magnitude in the conductivity profile lead to initiation altitudes 61.0–70.6 km, with delays in the range 3.825–9.825 ms. Consequently, we suggest that lightning characteristics dominate over atmospheric parameters in determining sprites’ initiation altitude and delay. The simulation of a −CG, with a constant current of 30 kA, did not produce a sprite seed, confirming an asymmetry in the response of the atmosphere to positive and negative lightning. This is due to the free electron drift direction that is away from the screening ionization wave, preventing the formation of the streamer seed for the great majority of −CGs. The same does not apply to halos, which depend on the occurrence of breakdown and can be produced by discharges of both polarities.     

Vertical motions in Thailand after the 2004 Sumatra–Andaman Earthquake from GPS observations and its geophysical modelling

Following previous findings from ongoing GPS research in Thailand since 2004 we continue to exploit the GPS technique to monitor and model land motions induced by the Sumatra–Andaman Earthquake. Our latest results show that up to the end of 2010, Thailand has been co-seismically displaced and is subsequently undergoing a post-seismic horizontal deformation with total displacements (co-seismic plus post-seismic) ranging from 10.5 to 74.7 cm. We observed the largest horizontal displacements in the southern part of Thailand and moderate and small displacements in the central and northern parts. In addition to horizontal displacements throughout Thailand, continuous GPS measurements show that large parts of Thailand are subsiding at rates up to 1 cm/yr. It is the first time that such vertical post-seismic deformations at large distances (650–1500 km away from the Earthquake’s epicentre) have been recorded. We have investigated the physical processes leading to the observed subsidence. While after-slip on the subduction interface induces negligible or even slightly positive vertical motions, relaxation in the asthenosphere is associated with a sizable subsidence. Predictions from a 3D finite element model feature an asthenosphere with an effective viscosity of the order of 3 * 10¹⁸ Pas, fit the horizontal post-seismic data and the observed subsidence well. This model is then used to predict the subsidence over the whole seismic cycle. The subsidence should go on with a diminishing rate through the next two decades and its final magnitude should not exceed 10 cm in the Bangkok area.     

A numerical model approximating extreme energetic electron events involved in the physical and chemical processes of the middle atmosphere

Relativistic electrons (with energies >150 keV) which originate in the outer radiation belt and detected by the Russian ‘Meteor’ series of satellites have been correlated with the atmospheric total ozone data compiled by almost 90 stations located around the world within the latitude zone 40°–70°N. In more than 60% of the stations examined we have detected a clear decrease of the ozone 3–5 days after the electron flux excess. A numerical model has been applied to approximate this effect based on relativistic electron initiated nitric oxides creation in the upper mesosphere with subsequent atmospheric transport (both vertical and horizontal) towards the upper stratosphere. A first attempt of local and temporal prediction of ozone depletion because of energetic electrons impact in the middle atmosphere has been illustrated.     

Global tectonics and metallogeny: Deep roots of some ore-controlling fracture zones. A possible relation to small-scale convective cells at the base of the lithosphere?

Criteria suggest that some of the ore-controlling fracture zones have deep roots, some extending into the upper mantle. The question is analyzed of whether a pattern of deep-seated fracture zones, extending to the base of a continental lithosphere, may be related to, or interact with, the small-scale convective cells (rolls) which originate [1,2] in the mantle beneath a moving lithospheric plate.The analysis follows a previous application of the same concept by others [3,4] who explained the origin of a chain of volcanoes by the movement of an oceanic lithospheric plate over a pattern of small-scale convective cells of the underlying mantle.The Canadian Shield is used to discuss the above question. A correlation has been found between: (1) the direction of two sets of long wavelength gravity anomalies detected by Stephenson and Beaumont, [8] in a study of the isostatic response of the Canadian Shield, and (2) the pattern of N-S and E-W trending trajectories of ore-controlling fracture zones, postulated by the author [5,6] in the southern part of the Canadian Shield. The Hudson Bay Paleolineament [7] correlates very well with the N-S set of gravity anomalies.The above correlation suggests that the ore-controlling lineaments and the long wavelength anomalies of the Canadian Shield are related. If the gravity anomalies reflect, as Stephenson and Beaumont [8] tentatively suggest, a pattern of convective cells beneath the base of the lithosphere, then the volcanic activity and metallogenesis may be related to the boundaries and corners of the cells. The author suggests that the convective cells could, in this case, originate by melting of the mantle material proceeding preferentially along the intersecting deep-seated fracture zones.The pattern of deep-seated fracture zones, compiled for the western United States [9] also shows a relationship of major ore deposits and ore clusters to the corners of rectangular blocks, defined by mutual intersection of the E-W and N-S fracture zones. In this case, the size of the blocks, measured in an E-W direction, is about 530 km, and in a N-S direction about 600 km. These figures are significantly close to the distance from the seismic discontinuity at a depth of 650 km, which is considered by others [2,4] as the lower boundary of the small-scale convection.     

Ionospheric perturbations in plasma parameters before global strong earthquakes

Based on the electron density (Ne) and temperature (Te) data from DEMETER, the ionospheric perturbations before 82 Ms ? 7.0 earthquakes (EQs) during 2005–2010 were studied, using moving median and space difference methods within 10 days before and 2 days after these events in local nighttime. It was found that the plasma parameters disturbances appeared before 49 EQs, in which more disturbances were detected before shallow-focus earthquakes than deep ones, and there was little difference between continental and oceanic ones, both exceeding 1/2 percentage. For the disturbed time, more perturbations were seen in 1, 3, 5, 6, 8 days before EQs and 1 day after EQs. For the spatial distribution, the anomalies before EQs were not just above the epicenters, but shifted equatorward with several degrees to almost twenty degrees. Most of the abnormities were positive ones, which demonstrate that Ne increases before EQs at the altitude of 670 km of DEMETER. Perturbations of Ne were more than that of Te, which illustrates that Ne is much more sensitive to seismic activity than Te.     

Model of energetic populations at Ganymede,implications for an orbiter

A model is developed to study the energetic particle populations in Ganymede’s magnetosphere. The main objective is to estimate to what extent the moon could protect an orbiter from radiations. Using Liouville’s theorem, the phase space density of particles coming from Jupiter’s magnetosphere is calculated at any point of Ganymede’s environment. Up to energies of ∼50–100 keV for ions and ∼10–20 MeV for electrons, Ganymede’s magnetic field appears to be able to form distinctive populations as loss-cones over the polar caps and radiation belts. At larger energies, these features are blurred by Larmor radius effects; the moon absorption simply creates a quasi-isotropic layer of ∼500 km thickness where the flux is reduced by ∼40–50%. The predictions are compared to Galileo measurements. In particular, we demonstrate the importance of the moon sweeping in reducing the flux over the polar caps. Interestingly, this can be accounted for by assuming that the particles bouncing between Jupiter and Ganymede are ideally scattered in pitch angle and permanently re-fill the loss-cone, which increases the precipitation on Ganymede’s polar cap. In overall, it is estimated that the radiation dose received by an orbiter of Ganymede will be reduced by more than 50–60% compared to the expected dose at Jupiter/Ganymede distance. This should have a positive impact on the design of a future orbiter of Ganymede.     

Modeling of the outer heliosphere with the realistic solar cycle

Time-dependent kinetic-continuum model of the solar wind interaction with the two-component local interstellar cloud (LIC) has been developed recently [Izmodenov, V., Malama, Y.G., Ruderman, M.S. Solar cycle influence on the interaction of the solar wind with local interstellar cloud. Astron. Astrophys. 429, 1069–1080, 2005a.]. Here, we adopted this model to the realistic solar cycle, when the solar wind parameters at the Earth’s orbit are taken from space data. This paper focuses on the results related to the termination shock (TS) excursion with the solar cycle that may help to understand Voyager 1 data obtained at and after the crossing of the termination shock and to predict the time of the TS crossing by Voyager 2.     

MHD modeling of the outer heliosphere: Achievements and challenges

An overview is presented of magnetic-field-related effects in the solar wind (SW) interaction with the local interstellar medium (LISM) and the different theoretical approaches used in their investigation. We discuss the possibility that the interstellar magnetic field (ISMF) introduces north–south and east–west asymmetries of the heliosphere, which might explain observational data obtained by the Voyager 1 and Voyager 2 spacecraft. The SW–LISM interaction parameters that are responsible for the deflection of the interstellar neutral hydrogen flow from the direction of propagation of neutral helium in the inner heliosheath are outlined. The possibility of a strong ISMF, which increases the heliospheric asymmetry and the H–He flow deflection, is discussed. The effect of the combination of a slow-fast solar wind during solar minimum over the Sun’s 11-year activity cycle is illustrated. The consequences of a tilt between the Sun’s magnetic and rotational axes are analyzed. Band-like areas of an increased magnetic field distribution in the outer heliosheath are sought in order to discover regions of possible 2–3 kHz radio emission.     

A facility for simulating Titan’s environment

As a result of measurements acquired by the Cassini–Huygens mission of Titan’s near surface atmospheric composition and temperature, Titan conditions can now be simulated in the laboratory and samples can subsequently be subjected to those conditions. Titan demonstrates an active hydrological-like cycle with its thick atmosphere, dynamic clouds, polar lakes of methane and ethane, moist regolith, and extensive fluvial erosive features. Unlike Earth, Titan’s hydrological-like cycle likely involves several constituents, primarily methane and ethane. Here the properties of a new Titan simulation facility are presented, including conceptual methodology, design, implementation, and performance results. The chamber maintains Titan’s surface temperature and pressure, and the sample cryogenic liquids undergoing experimentation are condensed within the chamber itself. During the experiments, the evaporation rates of the sample liquids are directly determined by continually measuring mass. Constituents are analyzed utilizing a Fourier Transform Infrared Spectroscopy (FTIR), and vapor concentrations are determined using a gas chromatograph fitted with a Flame Ionization Detector (FID). All pertinent data is logged via computer. Under laboratory conditions, the direct measurements of the evaporation rates of methane, ethane, and mixtures thereof can be achieved. Among the processes to be studied are the effects of regolith on transport from the subsurface to the atmosphere, the freezing point depression effects of dissolved nitrogen, and the solubility of various relevant organic compounds.     

Measurements of the radiation quality factor Q at aviation altitudes during solar minimum (2006–2008)

In radiation protection, the Q-factor has been defined to describe the biological effectiveness of the energy deposition or absorbed dose to humans in the mixed radiation fields at aviation altitudes. This particular radiation field is generated by the interactions of primary cosmic particles with the atoms of the constituents of the Earth’s atmosphere. Thus the intensity, characterized by the ambient dose equivalent rate H∗(10), depends on the flight altitude and the energy spectra of the particles, mainly protons and alpha particles, impinging on the atmosphere. These charged cosmic projectiles are deflected both by the interplanetary and the Earth’s magnetic field such that the corresponding energy spectra are modulated by these fields. The solar minimum is a time period of particular interest since the interplanetary magnetic field is weakest within the 11-year solar cycle and the dose rates at aviation altitudes reach their maximum due to the reduced shielding of galactic cosmic radiation. For this reason, the German Aerospace Center (DLR) performed repeated dosimetric on-board measurements in cooperation with several German airlines during the past solar minimum from March 2006 to August 2008. The Q-factors measured with a TEPC range from 1.98 at the equator to 2.60 in the polar region.     

An ion analyzer for the lunar surface with E-parallel–B

We present a novel instrument concept to measure the energy and mass spectra of ions incident on the lunar surface, based on the E-parallel–B or Thomson-parabola device used extensively as a diagnostic in the plasma fusion community. The Apollo-era Suprathermal Ion Detector Experiment (SIDE) was the first instrument package to perform in-situ measurements of ions incident on the lunar surface. The ions can originate from a variety of sources, including the solar wind, the Earth’s magnetotail, and photoionization of the thin lunar atmosphere. The species and energy distribution of ions arriving at the lunar surface depend in a complicated and poorly-understood fashion on the phase of the lunar day, the position of the Moon with respect to the Earth, and on the local plasma environment.