极区Hall电流经度差异特征的统计学研究:CHAMP卫星观测

本文利用2001至2010年间CHAMP(CHallenging Minisatellite Payload)卫星标量磁力仪(Overhauser Magnetometer)观测的磁场数据,反演得到电离层霍尔(Hall)电流,并且对极区电离层Hall电流的特征进行了统计学研究.研究主要关注平静期,即重联电场小于2 mV·m~(-1)条件下,在磁纬60°至90°范围内的Hall电流在不同太阳活动、季节、磁经度、磁地方时等条件下的变化特征.研究发现:Hall电流具有明显的经度差异,在南北半球呈现显著一波经度结构,而且南北半球反相,即北半球电流密度呈现峰-谷-峰结构,而南半球呈现谷-峰-谷结构.Hall电流密度的经度结构与太阳活动紧密相关,太阳活动高年经度差异最大,太阳活动中年经度差异次之,太阳活动低年经度差异最小.研究发现,在磁地方时为10-14MLT的白天,影响Hall电流的因素主要是太阳辐射;而磁地方时为21-03MLT的夜晚,除了电导率的影响之外,可能存在其他的重要的物理过程影响着Hall电流的经度分布.本文还研究了与电流相关的焦耳热的经度分布情况,发现其在南北半球分别呈现单峰、单谷结构,经度差异亦十分明显.     

The main features of solar wind plasma correlations of importance to space weather strategy

In this work solar wind measurements from several spacecraft were used to investigate the correlations of solar wind plasma parameters. These results provide a test of the concept of predicting space weather by monitoring the condition of the solar wind at a large distance (up to 230R_e, the L1 point) upstream from the Earth.We compared the ion flux and bulk velocity time behavior measured by widely-separated spacecraft: the spacecraft pairs INTERBALL-1 and IMP 8 (separations up to 30R_e), INTERBALL-1 and WIND, and IMP 8 and WIND (both with separations up to 250R_e). The average value of the ion flux correlation coefficient is about 0.73. But in some cases the plasma parameters from two spacecraft are very different in both behavior and value, so correlations are very poor.The technique of multifactorial analysis was used to obtain the physical dependences of the correlations on the spacecraft separation and on different plasma and magnetic field parameters. We found that the correlation values have a weak but significant dependence on the separation perpendicular to the Sun–Earth line (YZ_se-separations up to 90R_e).The most important factors influencing the correlation level are density (or ion flux) variability, the direction of the IMF vector to the Sun–Earth line (cone angle), and the solar wind bulk velocity.     

Water retention effects on elastic properties of Opalinus shale

Shales play an important role in many engineering applications such as nuclear waste, CO₂ storage and oil or gas production. Shales are often utilized as an impermeable seal or an unconventional reservoir. For both situations, shales are often studied using seismic waves. Elastic properties of shales strongly depend on their hydration, which can lead to substantial structural changes. Thus, in order to explore shaly formations with seismic methods, it is necessary to understand the dependency of shale elastic properties on variations in hydration. In this work, we investigate structural changes in Opalinus shale at different hydration states using laboratory measurements and X-ray micro-computed tomography. We show that the shale swells with hydration and shrinks with drying with no visible damage. The pore space of the shale deforms, exhibiting a reduction in the total porosity with drying and an increase in the total porosity with hydration. We study the elastic properties of the shale at different hydration states using ultrasonic velocities measurements. The elastic moduli of the shale show substantial changes with variations in hydration, which cannot be explained with a single driving mechanism. We suggest that changes of the elastic moduli with variations in hydration are driven by multiple competing factors: (1) variations in total porosity, (2) substitution of pore-filling fluid, (3) change in stiffness of contacts between clay particles and (4) chemical hardening/softening of clay particles. We qualitatively and quantitatively analyse and discuss the influence of each of these factors on the elastic moduli. We conclude that depending on the microstructure and composition of a particular shale, some of the factors dominate over the others, resulting in different dependencies of the elastic moduli on hydration.     

Solar Irradiance Models and Measurements: A Comparison in the 220–240?nm wavelength band

Solar irradiance models that assume solar irradiance variations to be due to changes in the solar surface magnetic flux have been successfully used to reconstruct total solar irradiance on rotational as well as cyclical and secular time scales. Modelling spectral solar irradiance is not yet as advanced, and also suffers from a lack of comparison data, in particular on solar cycle time scales. Here, we compare solar irradiance in the 220–240?nm band as modelled with SATIRE-S and measured by different instruments on the UARS and SORCE satellites. We find good agreement between the model and measurements on rotational time scales. The long-term trends, however, show significant differences. Both SORCE instruments, in particular, show a much steeper gradient over the decaying part of cycle 23 than the modelled irradiance or that measured by UARS/SUSIM.     

Phenological Responses to ENSO in the Global Oceans

Phenology relates to the study of timing of periodic events in the life cycle of plants or animals as influenced by environmental conditions and climatic forcing. Phenological metrics provide information essential to quantify variations in the life cycle of these organisms. The metrics also allow us to estimate the speed at which living organisms respond to environmental changes. At the surface of the oceans, microscopic plant cells, so-called phytoplankton, grow and sometimes form blooms, with concentrations reaching up to 100 million cells per litre and extending over many square kilometres. These blooms can have a huge collective impact on ocean colour, because they contain chlorophyll and other auxiliary pigments, making them visible from space. Phytoplankton populations have a high turnover rate and can respond within hours to days to environmental perturbations. This makes them ideal indicators to study the first-level biological response to environmental changes. In the Earth’s climate system, the El Niño–Southern Oscillation (ENSO) dominates large-scale inter-annual variations in environmental conditions. It serves as a natural experiment to study and understand how phytoplankton in the ocean (and hence the organisms at higher trophic levels) respond to climate variability. Here, the ENSO influence on phytoplankton is estimated through variations in chlorophyll concentration, primary production and timings of initiation, peak, termination and duration of the growing period. The phenological variabilities are used to characterise phytoplankton responses to changes in some physical variables: sea surface temperature, sea surface height and wind. It is reported that in oceanic regions experiencing high annual variations in the solar cycle, such as in high latitudes, the influence of ENSO may be readily measured using annual mean anomalies of physical variables. In contrast, in oceanic regions where ENSO modulates a climate system characterised by a seasonal reversal of the wind forcing, such as the monsoon system in the Indian Ocean, phenology-based mean anomalies of physical variables help refine evaluation of the mechanisms driving the biological responses and provide a more comprehensive understanding of the integrated processes.     

Cyclic variations in the Earth’s flattening and questions of seismotectonics

For more than a decade, the global network of GPS stations whose measurements are part of the International GPS Service (IGS) have been recording cyclic variations in the radius vector of the geodetic ellipsoid with a period of one year and amplitude of ~10 mm. The analysis of the figure of the Earth carried out by us shows that the observed variations in the vertical component of the Earth’s surface displacements can induce small changes in the flattening of the Earth’s figure which are, in turn, caused by the instability of the Earth’s rotation. The variations in the angular velocity and flattening of the Earth change the kinetic energy of the Earth’s rotation. The additional energy is ~10²¹ J. The emerging variations in the flattening of the Earth’s ellipsoid lead to changes in the surface area of the Earth’s figure, cause the development of deformations in rocks, accumulation of damage, activation of seismotectonic processes, and preparation of earthquakes. It is shown that earthquakes can be caused by the instability of the Earth’s rotation which induces pulsations in the shape of the Earth and leads to the development of alternating-sign deformations in the Earth’s solid shell.     

The Grad-Shafranov reconstruction in twenty years: 1996–2016

We review and summarize the applications of the Grad-Shafranov (GS) reconstruction technique to space plasma structures in the Earth’s magnetosphere and in the interplanetary space. We organize our presentations following the branches of the “academic family tree” rooted on Prof. Bengt U. Ö. Sonnerup, the inventor of the GS method. Special attentions are paid to validations of the GS reconstruction results via (1) the direct validation by co-spatial in-situ measurements among multiple spacecraft, and (2) indirect validation by implications and interpretations of the physical connection between the structures reconstructed and other related processes. For the latter, the inter-comparison and interconnection between the large-scale magnetic flux ropes (i.e., Magnetic Clouds) in the solar wind and their solar source properties are presented. In addition, we also summarize various GS-type (or -like) reconstruction and an extension of the GS technique to toroidal geometry. In particular, we point to a possible advancement with added complexity of “helical symmetry” and mixed helicity, in the hope of stimulating interest in future development. We close by offering some thoughts on appreciating the scientific merit of GS reconstruction in general.     

Consistent parameter constraints for soil hydraulic functions

Parameters of functions to describe soil hydraulic properties are derived from measurements by means of parameter estimation. Of crucial importance here is the choice of correct constraints in the parameter space. Often, the parameters are mere shape parameters without physical meaning, giving flexibility to the model. A fundamental requirement is that the hydraulic functions are monotonic: the retention function and the conductivity function can only decrease as the capillary suction increases. A stricter physical requirement for the conductivity function is that its decrease with respect to saturation is at least linear. This linear decrease would occur if all pores of a capillary bundle had an equal radius. In the first part of this contribution, we derive constraints for the so-called tortuosity parameter of the Mualem conductivity model, which allow highest possible flexibility on one hand and guarantee physical consistency on the other hand. In combination with the retention functions of Brooks and Corey, van Genuchten, or Durner, such a constraint can be expressed as a function of the pore-size distribution parameters. In the second part, we show that a common modification of retention models, which is applied to reach zero water content at finite suction, can lead to the physically unrealistic case of increasing water content with increasing suction. We propose a solution for this problem by slightly modifying these models and introducing a correct parameter constraint.     

Total solar irradiance variations

Total solar irradiance has been monitored from space for nearly two decades. These space-borne observations have established conclusively that total solar irradiance changes over a wide range of periodicities—from minutes to the 11-year solar cycle. Since the total energy flux of the Sun is the principal driver for all Earths atmospheric phenomena, the accurate knowledge of the solar radiation received by the Earth and its variations is an extremely important issue. In this paper we review the long-term variations of total solar irradiance during solar cycles 21 and 22. We conclude that, within the current accuracy and precision of the measurements, the minimum level of total solar irradiance is about the same for both solar cycles 21 and 22.     

New approach to the construction of the Earth’s upper atmosphere model using the satellite local data on the atmosphere

The present-day models of the Earth’s upper atmosphere make it possible to construct the spatial-temporal pattern of variations in the atmospheric parameters on the planetary scale in essence in the averaged form. The set of data on the satellite deceleration in the atmosphere, probe measurements aboard geophysical rockets, and radiowave incoherent scatter measurements in the Earth’s atmosphere are used to construct these standard models. The current level of the space studies makes it possible to use a new method to study the Earth’s upper atmosphere: to study the upper atmosphere by measuring the absorption of the solar XUV radiation by the Earth’s atmosphere during the solar disk observations.     

Correlation Properties of the Samplings of Variations in the <Emphasis Type="Italic">F</Emphasis>2-Layer Critical Frequency in Various Solar and Geophysical Conditions

The temporal and spatial correlation of variations in the F2-layer critical frequency in various solar and geophysical conditions are considered based on a large array of experimental data. The possible extrapolation of ionospheric observational data is studied based on the obtained results on the stability of the correlation in time and space. It is shown that it is possible to use the obtained results for rapid short-term ionospheric forecast. Moreover, the data on the temporal and spatial correlation radius (at the level of a few hours and thousand kilometers) are used to study the irregular structure of the ionosphere and to determine the characteristic dimension of large-scale irregularities and apparent velocities of their motion.     

The influence of geomagnetic and solar variabilities on lower thermosphere density

Atmospheric density measurements near 200 km from the Satellite Electrostatic Triaxial Accelerometer (SETA) experiment are analyzed for geomagnetic and solar flux variability effects. Data from the SETA experiment, onboard two satellites, are available for the periods of May to November 1982, and July 1983 to March 1984. The data utilized the span ±79.5° latitude, and are available for both day (1030 LT) and night (2230 LT). Annual and semiannual density variations are removed and regression analyses are performed on the residuals using a series of lagged 3 h K_p indices to determine and remove geomagnetic fluctuations. Densities are found to increase by as much as 134% in response to an increase in the K_p index from 1 to 6. Monthly curves are generated for the K_p regression coefficients to delineate seasonal-latitudinal and day/night dependences, which reflect the effects of mean meridional advection of disturbances from high to low latitudes. Further analyses are performed comparing measured densities with MSISE-90 predictions. Results show that the model is able to capture many of the prominent features, but does not fully predict the level of variability for the individual disturbance periods analyzed. After the geomagnetic effects are removed, the residual densities are interpreted in terms of solar flux variability. The daily-averaged SETA density residuals are strongly correlated with long-term solar flux variability, and exhibit a much greater dependence on the 27-day solar rotation period than MSISE-90 predictions. Variations in residual density of the order of 10–20% occur in association with day-to-day and 27-day solar flux variations. The MSIS model does not accurately predict the magnitude of these short-term density variations in response to solar activity.     

Transition of solar cycle length in association with the occurrence of grand solar minima indicated by radiocarbon content in tree-rings

In this paper, we review the variation of the 11-year solar cycle since the 15th century revealed by the measurement of radiocarbon content in single-year tree-rings of Japanese cedar trees. Measurements of radiocarbon content in absolutely dated tree-rings provide a calibration curve for accurate dating of archaeological matters, but at the same time, enable us to examine the variations of solar magnetic activity in the pre-historical period. The Sun holds several long-term quasi-cyclic variations in addition to the fundamental 11-year sunspot activity cycle and the 22-year polarity reversal cycle, and it is speculated that the property of the 11-year and the 22-year solar cycle varies in association with such long-term quasi-cycles. It is essential to reveal the details of solar variations around the transition time of solar dynamo for illuminating the mechanisms of the long-term solar variations. We therefore have investigated the property of the 11-year and 22-year cycles around the two grand solar minima; the Maunder Minimum (1645–1715 AD) and the Spoerer Minimum (1415–1534 AD), the periods of prolonged sunspot minima. As a result, slight stretching of the “11-year” and the “22-year” solar cycles was found during these two grand solar activity minima; continuously during the Maunder Minimum and only intermittently during the Spoerer Minimum. On the contrary, normal or slightly shortened 11-year cycles were detected during the interval period of these two minima. It suggests the inverse correlation between the solar cycle length and solar magnetic activity level, and also the change of meridional flow during the grand solar activity minima. Further measurements for the beginning of the grand solar minima will provide a clue to the occurrence of such prolonged sunspot disappearance. We also discuss the effect of solar variations to radiocarbon dating.     

Polar cap response to the solar wind density jump under constant southward IMF

Sharp changes of the solar wind parameters determining the dynamic pressure jump lead to strong magnetosphere-ionosphere disturbances. Here the effect on the Earth’s ionospheric high latitudes of the solar wind dynamic pressure pulse caused only by the increase of the interplanetary plasma density under southward constant IMF is considered. We investigate reaction of the cross-polar cap potential on the increase of AL index and/or jump of the solar wind density. It is found that for the case of 10 January 1997 the main contribution to the polar cap potential drop increase gave the growth of AL index relative to the input of the solar wind density jump. We also study the influence of the solar wind density increase on the crosspolar cap potential for the quiet magnetospheric conditions. It occurred that the polar cap potential difference decreases with the great increase of the interplanetary plasma density. For the disturbed magnetosphere the main role in the polar cap potential drop increase plays increase of AL. Thus, we found the change of the cross-polar cap potential due to the AL index variations and/or the solar wind density drop even in a case when the interplanetary electric field is constant.     

A solar pattern in the longest temperature series from three stations in Europe

We analyze the longest temperature series from Prague, Bologna and Uccle. We partition daily minimum and maximum temperatures and their differences in two subsets as a function of high vs low solar activity, using the superimposed epochs method. Differences display patterns with significant amplitudes and time constants ～3 months. These are recognized in all stations and are stable against a change in the analyzed period. Amplitude of variations is ～1 °C. Differences between average annual values corresponding to high vs low activity periods are also ～1 °C. Solar activity may account for these long-term temperature variations. These variations also present local characteristics, which may render identification of a global correlation delicate. We discuss possible physical mechanisms by which solar variation could force climate changes (e.g. through solar activity itself, the EUV part of the solar flux, cosmic rays, the downward ionosphere-earth current density, etc.).     

Solar flare effect on the geomagnetic field and ionosphere

This paper studies the ionospheric and geomagnetic response to an X6.2 solar flare recorded at 14:30 UT on December 13, 2001, in quiet geomagnetic conditions which allow the variations in the geomagnetic field and ionosphere measurements to be easily related to the solar flare radiation.By using measurements from the global positioning system (GPS) and geomagnetic observatories, the temporal evolution of ionospheric total electron content variation, vTEC_V, and geomagnetic field variations, δB, as well as their rates of variation, were obtained around the subsolar point at different solar zenith angles. The enhancement of both parameters was recorded one to three minutes later than the Geostationary Operational Environmental Satellite (GOES) programme recording; such delay tends to depend on the latitude, longitude, and solar zenith angle of the observatory's observations.The vTEC_V is related to the local time and the δB to the intensity and position of the ionospheric currents.The vTEC_V′s maximum value is always recorded later than the maximum values reached by δB and the X-ray intensity. The maximum δB is larger in the local morning than in the afternoon.The rates of vTEC_V and δB have two maximum values at the same time as the maximum values recorded by Hα (for each ribbon).This work shows the quantitative and qualitative relations between a solar flare and the ionospheric and geomagnetic variations that it produces.     

Variations of daytime and nighttime electron temperature and heat flux in the upper ionosphere,topside ionosphere and lower plasmasphere for low and high solar activity

A database of the electron temperature (T_e) comprising of most of the available LEO satellite measurements is used for studying the solar activity variations of T_e. The T_e data are grouped for two levels of solar activity (low LSA and high HSA), five altitude ranges between 350 and 2000 km, and day and night. By fitting a theoretical expression to the T_e values we obtain variation of T_e along magnetic field lines and heat flux for LSA and HSA. We have found that T_e increases with increase in solar activity at low and mid-latitudes during nighttime at all altitudes studied. During daytime the T_e response to solar activity depends on latitude, altitude, and season. This analysis shows existence of anti-correlation between T_e and solar activity at mid-latitudes below 700 km during the equinox and winter day hours. Heat fluxes show small latitudinal dependence for daytime but substantial for nighttime.     

Growth in the atmospheric aerosol concentration as a climate forcing agent

The level of scattered radiation is analyzed using data of diurnal solar coronograph observations at the Mountain Astronomic Station in the period 1957?C2010 around coronal spectral lines at 5303 ? and 6374 ?. The observations were performed near the solar limb and were normalized by the intensity of the solar disk center. The measurements revealed variations on different timescales: seasonal variations, local maxima on timescales of a few years, and long-term trends. The local changes in the level of scattered radiation were found to be probably due to volcanic eruptions. An analysis revealed a tendency towards an increased in scattered radiation by approximately 40% during the last 50 years. The variations in the level of scattered radiation are compared with the concentration of atmospheric aerosols. The long-term growth in scattered radiation compares well with changes in the Earth??s near-surface temperature and is possibly associated with global climate change.     

GPU-accelerated computing of three-dimensional solar wind background

High-performance computational models are required to make the real-time or faster than real-time numerical prediction of adverse space weather events and their influence on the geospace environment. The main objective in this article is to explore the application of programmable graphic processing units (GPUs) to the numerical space weather modeling for the study of solar wind background that is a crucial part in the numerical space weather modeling. GPU programming is realized for our Solar-Interplanetary-CESE MHD model (SIP-CESE MHD model) by numerically studying the solar corona/interplanetary solar wind. The global solar wind structures are obtained by the established GPU model with the magnetic field synoptic data as input. Meanwhile, the time-dependent solar surface boundary conditions derived from the method of characteristics and the mass flux limit are incorporated to couple the observation and the three-dimensional (3D) MHD model. The simulated evolution of the global structures for two Carrington rotations 2058 and 2062 is compared with solar observations and solar wind measurements from spacecraft near the Earth. The MHD model is also validated by comparison with the standard potential field source surface (PFSS) model. Comparisons show that the MHD results are in good overall agreement with coronal and interplanetary structures, including the size and distribution of coronal holes, the position and shape of the streamer belts, and the transition of the solar wind speeds and magnetic field polarities.