Lateral lithospheric conductivity patterns,from European synchronous observatory data

Much information on the regional lithospheric structure may come from MTS data acquired by synchronous 2D arrays and processed with regard to the nonlocal response of a laterally inhomogeneous subsurface. We suggest to invert the nonlocal MT responses applying correlation of all surface horizontal and vertical components of the geomagnetic field recorded simultaneously at all stations. The inversion algorithm has been applied to 2004–2005 European observatory data of diurnal S_q variations for first five harmonics and yielded lateral conductivity patterns for different periods. The maps show spatial correlation between conductivity maxima and lithospheric thickness minima and, specifically, highlight the contours of the Pannonian basin, where lithosphere is as thin as ～50 km, from seismic data.     

Upper mantle lateral inhomogeneity estimated using global MVS

We processed synchronous observatory MVS data at periods from 1.2 to 6 days to estimate lateral conductivity variations in the upper mantle. The inverse problem for depth-dependent conductivity distribution was solved using correlation of the measured surface components of the geomagnetic field. We applied spherical harmonic expansion of the complex amplitudes of time harmonics and multiple averaging of the amplitudes obtained for different SH sets. The chosen approximation approach was tested by comparing the EM field simulated for a smoothly inhomogeneous Earth with the exact forward solution. The SHA-analyzed observatory data for fifteen magnetic storms were used to map global apparent conductivity variations.     

Variations of the solar-wind stream structure in the region of subsonic flow during the 11-year solar cycle

The large-scale stream structure of the solar wind near the Sun and its evolution during the 11-year solar activity cycle are investigated. The study is based on observations of scattering of the radiation from compact natural radio sources at radial distances R≤14R _S (R _S is the solar radius). Regular observations were conducted in 1981–1998 on the RT-22 and DKR-1000 radio telescopes of the Russian Academy of Sciences at Pushchino, at λ=1.35 cm and 2.7 m, respectively. The radial dependences of the interplanetary scintillations m(R) and the scattering angle 2?(R) are considered together with the structure of large-scale magnetic fields in the solar corona at R=2.5R _S. The entire range of variations in the level of scattering and the associated heliolatitude flow structures in the subsonic solar wind forms over the 11-year solar cycle, as a direct result of the large-scale structure of the evolving magnetic fields at the source of the solar-wind streamlines.     

Hydrochemical and PCO2 variations of a cave stream in a subtropical karst area,Chongqing, SW China: piston effects,dilution effects,soil CO2 and buffer effects

The role of CO₂ in karst has been of interest for decades, and emphasized by IGCP 379, International Geoscience Programme, UNESCO started in 1995. There are still open questions about the dynamics of carbon in karst systems, particularly the flux of carbon between the surface and subsurface and between different components in the karst subsurface. This research report focuses on the variations of hydrochemistry and PCO₂ (partial pressures of carbon dioxide) in subtropical karst groundwater, using high-resolution auto-monitoring hydrochemical data (15-min intervals). The aim of this study was to understand how hydrochemistry and PCO₂ in karst systems respond to recharge over different time scales and what the controlling factors are. An auto-monitoring hydrochemistry station was installed about 300 m upstream from the exit in the active stream channel of Xueyu Cave, a typical subtropical karst cave. Four years of high-resolution continuous pH, specific conductivity (Spc), temperature and water-level data were collected. A thermodynamic model was used to link the continuous data to monthly water quality data, allowing the calculation of CO₂ partial pressures and calcite saturation (SIc) levels on a continuous basis. Seasonal, diurnal and storm-scale variations were captured in the hydrochemistry and calculated PCO₂ records, indicating that the cave stream is a dynamic and variable system. Seasonal features (higher specific conductivity and lower pH in summer; lower specific conductivity and higher pH in winter) tend to covary with temperature which influences the production of CO₂ in soils, thus being the driving force for the variations (the soil CO₂ effect). Due to the buffer effect of a thick vadose zone and large void cave space, diurnal variations are not obvious compared with epikarst springs in SW China. Storm-scale fluctuations due to storm events occur during the summer rainy season. Piston flow effects, dilution and soil CO₂ effects determine the variations in different storm events. At the beginning of the rains, the piston effect drives the variations, characterized by increase in Spc, SIc and pH in the cave stream and decrease in PCO₂. With heavy rainfall, decrease in Spc shows control by the dilution effect, while decrease in SIc and pH and increase in PCO₂ indicates the greater influence of soil CO₂. These results imply that the soil and cave voids are important factors influencing the hydrochemical evolution of karst groundwater. Future works need to use such high-resolution technology widely for tracing the PCO₂ and hydrochemical variations in different karst aquifers.     

Magnetotelluric and audiomagnetotelluric measurements in Finland

Joint Finnish—Hungarian MT (magnetotelluric) and AMT (audiomagnetotelluric) measurements were carried out in Finland in the framework of the international ELAS project. The conditions for MT measurements are favorable at these latitudes. Five MT and 150 AMT stations gave information on the electrical conductivity distribution in the area: AMT results guided the choice of MT sites with minimal near-surface distortion effects and helped the interpretation of the MT soundings; the MT measurements indicate the presence of large conductivity anomalies and can be best interpreted as lateral induction effects of near-surface dyke structures. This result is confirmed by a certain correspondence between the directions of the maximum impedances and of the tectonic zones of the area.Any information about the upper mantle would require the use of S_q harmonics because of the crustal conductivity anomalies detected by the MT measurements.     

The large-scale solar magnetic field and 11-year activity cycles

Magnetic Hα synoptic maps of the Sun for 1915–1999 are analyzed and the intensities of spherical harmonics of the large-scale solar magnetic field computed. The possibility of using these Hα maps as a database for investigations of long-term variations of solar activity is demonstrated. As an example, the magnetic-field polarity distribution for the Hα maps and the analogous polarity distribution for the magnetographic maps of the Stanford observatory for 1975–1999 are compared. An activity index A(t) is introduced for the large-scale magnetic field, which is the sum of the magnetic-moment intensities for the dipole and octupole components. The 11-year cycle of the large-scale solar magnetic field leads the 11-year sunspot cycle by, on average, 5.5 years. It is concluded that the observed weak large-scale solar magnetic field is not the product of the decay of strong active-region fields. Based on the new data, the level of the current (23rd) solar-activity cycle and some aspects of solar-cycle theory are discussed.     

Analysis of pumping test data for determining unconfined-aquifer parameters: Composite analysis or not?

Recently, composite analysis (CA), which simultaneously analyzes all drawdown data from multiple observation wells, has been applied to determine the hydraulic parameters of an unconfined aquifer. Moench (1994) claimed that the value of specific yield (S _y) determined from non-composite analysis (nonCA) is sometimes unrealistically low as compared with that obtained by water-balance calculation, and results from CA are better representative of aquifer properties than those from nonCA. To examine the validity of this assertion, the drawdown data from a pumping test conducted at Cape Cod, Massachusetts, USA, were analyzed using both nonCA and CA methods. The results show that the mean estimates of hydraulic conductivity and S _y determined from CA are close to those determined from nonCA. In some cases the analysis based on CA also results in low estimates of S _y as compared with those determined based on nonCA. A hypothetical case study is presented, which examines the effect of measurement errors on the estimated parameters. The results indicate that the CA method also gives poorer estimates of S _y than the nonCA method if the pumping test data contain measurement errors. Moench AF (1994) Specific yield as determined by type-curve analysis of aquifer-test data. Ground Water, 32(6):949–957.     

Horizontal magnetic fields in the solar photosphere

Two-dimensional simulations of time-dependent solar magnetogranulation are used to analyze the horizontal magnetic fields and the response of the synthesized Stokes profiles of the IR FeI λ1564.85 nm line to the magnetic fields. The 1.5-h series of MHD models used for the analyses reproduces a region of the magnetic network in the photosphere with an unsigned magnetic flux density of 192 G at the solar surface. According to the magnetic-field distribution obtained, the most probable absolute strength of the horizontal magnetic field at an optical depth of τ ₅ = 1(τ ₅ denotes τ at λ = 500 nm) is 50 G, while the mean value is 244 G. On average, the horizontal magnetic fields are stronger than the vertical fields to heights of about 400 km in the photosphere due to their higher density and the larger area they occupy. The maximum factor by which the horizontal fields are greater is 1.5. Strong horizontal magnetic flux tubes emerge at the surface as spots with field strengths of more than 500 G. These are smaller than granules in size, and have lifetimes of 3–6 min. They form in the photosphere due to the expulsion of magnetic fields by convective flows coming from deep subphotospheric layers. The data obtained qualitatively agree with observations with the Hinode space observatory.     

Radio and optical spectra of objects from three complete samples of radio sources

We present classifications, optical identifications, and radio spectra for 19 radio sources from three complete samples, with declinations 4°–6° (B1950, S _{3.9 GHz} > 200 mJy), 10°–12°30′ (J2000, S _{4.85 GHz} > 200 mJy), and 74°–75° (J2000, S _{4.85 GHz} > 100 mJy). We also present corresponding information for the radio source J0527+0331. The right ascensions are 0–24^h and the Galactic latitudes |b| > 15° for all the samples. Our observations were obtained with the 6 m telescope from the Special Astrophysical Observatory in the range 4000–9000 Å or 4000–7500 Å and the RATAN-600 radio telescope at frequencies in the range 0.97–21.7 GHz. We obtained flux densities for the radio sources and optical spectra for their optical counterparts. Nine objects were classified as quasars with redshifts from z = 1.029 to 3.212; nine objects are emission-line galaxies with redshifts from 0.172 to 0.546, and one is a galaxy with burstlike star formation at z = 0.156, and one is a BL Lac object with z = 0.509. The spectra of five radio sources were decomposed into extended and compact components. The radio source J0527+0331, identified with a BL Lac object, displays significant variations of time scales from several days to several years. Data on flux variations are presented for 11 radio sources, as well as their spectra at several epochs.     

Restricted utility of δC of bulk organic matter as a record of paleovegetation in some loess-paleosol sequences in the Chinese Loess Plateau

Molecular stratigraphic analyses using gas chromatograph-mass spectrometry have been performed in the upper section (S₀, L₁, S₁) of the Yuanbo loess-paleosol sequences in northwest China, with a record extending from the last interglaciation through the present interglaciation. The CPI (Carbon Preference Index) values of both n-alkanols and n-alkan-2-ones display variations between loess deposits and paleosols, showing a correlation with the magnetic susceptibility record, an indicator of the East Asian summer monsoon. The observed variations in the indexes in relation to changes in lithology/paleoclimate are proposed to result from microbial degradation of higher plant lipids in the paleosols. The CPI values of n-alkanes, n-alkanols, and n-alkan-2-ones are negatively correlated with δ¹³C of bulk organic matter. The correlations suggest that the observed glacial-interglacial variations of δ¹³C data in the loess stratigraphy reflect the relative importance of the contribution of paleovegetation compared with microorganisms (including both the degradation and the addition of organic matter) and allochthonous loess/soil parent materials. It is thus necessary to evaluate the contributions of the latter two before the paleovegetation can be reconstructed based on the δ¹³C analysis of bulk organic matter in some loess-paleosol sequences of the Chinese Loess Plateau.     

A new generalized soil thermal conductivity model for sand–kaolin clay mixtures using thermo-time domain reflectometry probe test

As the demand of exploitation and utilization of geothermal energy increases, more geothermal-related earth structures occur recently. The design of the structures depends upon an accurate prediction of soil thermal conductivity. The existing soil thermal conductivity models were mostly developed by empirical fits to datasets of soil thermal conductivity measurements. Due to the gaps in measured thermal conductivities between any two tested natural soils, the models may not provide accurate prediction for other soils, and the predicted thermal conductivity might not be continuous over the entire range of soil type. In this research, a generalized soil thermal conductivity model was proposed based on a series of laboratory experiments on sand, kaolin clay and sand–kaolin clay mixtures using a newly designed thermo-time domain reflectometry probe. The model was then validated with respect to k _dry–n (thermal conductivity of dry soils and porosity) and k _r–S _r (normalized thermal conductivity and degree of saturation) relationships by comparing with previous experimental studies. The predicted thermal conductivities were found to be in a good agreement with the experimental data collected from both this study and the other literatures with at least 85% confidence interval. It is concluded that the proposed model accounts for the effects of both environmental factors (i.e., moisture content and dry density) and compositional factors (i.e., quartz content and soil type) on soil thermal conductivity, and it has a great potential in predicting soil thermal conductivity more accurately for geothermal applications.     

Paleostress reconstructions of Jabal Hafit structures,Southeast of Al Ain City,United Arab Emirates (UAE)

This paper presents the first paleostress results obtained from displacement and fracture systems within the Lower Eocene sediments at Jabal Hafit, Abu Dhabi Emirate, UAE. Detailed investigation of Paleogene structures at Jabal Hafit reveal the existence of both extensional structures (normal faults) and compressional structures (strike-slip and reverse faults). Structural analysis and paleostress reconstructions show that the Paleogene kinematic history is characterized by the succession of four paleostress stages. Orientation of principal stresses was found from fault-slip data using an improved right-dihedra method, followed by rotational optimisation (TENSOR program).The paleostress results confirm four transtensional tectonic stages (T1–T4) which affected the study area. The first tectonic stage (T1) is characterized by S_Hmax NW–SE σ₂-orientation. This stage produced NW–SE striking joints (tension veins) and E–W to ENE–WSW striking dextral strike-slip faults. The proposed age of this stage is Early Eocene. The second stage (T2) had S_Hmax N–S σ₂-orientation. N–S striking joints and NNE–SSW striking sinistral strike-slip faults, E–W striking reverse faults and N–S striking normal faults were created during this stage. The T2 stage is interpreted to be post-Early Eocene in age. The third stage (T3) is characterized by S_Hmax E–W σ₂-orientation. This stage reactivated the E–W reverse faults as sinistral strike-slip faults and created E–W striking joints and NE–SW reverse faults. The proposed age for T3 is post-Middle Eocene. During the T3 (S_Hmax E–W σ₂-orientation) stage the NNW-plunging Hafit anticline was formed. The last tectonic stage that affected the study area (T4) is characterized by S_Hmax NE–SW σ₂-orientation. During this stage, the ENE–WSW faults were reactivated as sinistral strike-slip and reverse faults. NE–SW oriented joints were also created during the T4 (S_Hmax NE–SW σ₂-orientation) stage. The interpreted age of this stage is post-Middle Miocene time but younger than T3 (S_Hmax E–W σ₂-orientation) stage.     

Temperature-driven meltwater production and hydrochemical variations at a glaciated alpine karst aquifer: implication for the atmospheric CO<Subscript>2</Subscript> sink under global warming

About two hydrological years of continuous data of discharge, temperature, electrical conductivity and pH have been recorded at the Glarey spring in the Tsanfleuron glaciated karst area in the Swiss Alps, to understand how glaciated karst aquifer systems respond hydrochemically to diurnal and seasonal recharge variations, and how calcite dissolution by glacial meltwater contributes to the atmospheric CO₂ sink. A thermodynamic model was used to link the continuous data to monthly water quality data allowing the calculation of CO₂ partial pressures and calcite saturation indexes. The results show diurnal and seasonal hydrochemical variations controlled chiefly by air temperature, the latter influencing karst aquifer recharge by ice and snowmelt. Karst process-related atmospheric CO₂ sinks were more than four times higher in the melting season than those in the freezing season. This finding has implication for understanding the atmospheric CO₂ sink in glaciated carbonate rock terrains: the carbon sink will increase with increasing runoff caused by global warming, i.e., carbonate weathering provides a negative feedback for anthropogenic CO₂ release. However, this is a transient regulation effect that is most efficient when glacial meltwater production is highest, which in turn depends on the future climatic evolution.     

Geomagnetic induction responses of anisotropic conducting mantle

Phase change of dielectric magnesiowustite in the lower mantle may leave signatures in geomagnetic records of the globally distributed array of observatories. The related features appear in EM induction responses of lower mantle, which are studied theoretically. The surface EM field corresponding to a response of the earth with conductivity anisotropy in a mantle spherical layer is presented as the sum of the magnetic and electric modes. Equations for the fields of both modes and their relationship in a weakly anisotropic earth are obtained by the perturbation method. The two field modes are analyzed jointly and separately to characterize the conductivity tensor of the anisotropic lower mantle. The tensor elements corresponding to the tangential components of the field can be estimated from the magnetic mode alone recorded currently by the global network of geomagnetic observatories. For the tensor data to be complete, observatory data on lateral variations of the electric field are required in addition to three-component geomagnetic records.     

Variability of the gravitational lens UM 673 in 2003–2005

Results of VRI photometry of two components of the gravitationally lensed quasar UM 673 in 2003–2005 are presented. The observational data were obtained on the 1.5-m telescope of the Maidanak observatory. During the monitoring of the system, considerable brightness variations in components A and B with amplitudes of about 0.2 ^m were recorded, demonstrating variability of the quasar. The increase in the brightness of the components was accompanied by reddening of their V-R color indices. Analysis of the brightness and color variations of the components shows no variations related to microlensing.     

Sources of the global magnetic field of the Sun

Data on the global magnetic field (GMF) of the Sun as a star for 1968–1999 are used to determine the correlation of the GMF with the radial component of the interplanetary magnetic field (IMF) |B _r|; all data were averaged over a half year. The time variations in the GMF |H| are better correlated with variations in |B _r|; than the results of extrapolating the field from the “source surface” to the Earth’s orbit in a potential model based on magnetic synoptic maps of the photosphere. Possible origins for the higher correlation between the GMF and IMF are discussed. For both the GMF and IMF, the source surface actually corresponds to the quiet photosphere—i.e., background fields and coronal holes—rather than to a spherical surface artificially placed ≈2.5 R _⊙ from the center of the Sun, as assumed in potential models (R _⊙ is the solar radius). The mean effective strength of the photospheric field is about 1.9 G. There is a nearly linear dependence between |H| and |B _r|. The strong correlation between variations in |H| and |B _r| casts doubt on the validity of correcting solar magnetic fields using the so-called “saturation” factor δ^?1 (for magnetograph measurements in the λ 525.0 nm FeI line).     

Large-scale stratification of the ionosphere during earthquake preparation

Modulation of ionospheric plasma density by acoustic waves (AW) is investigated. A rigorous electrodynamic model of the interaction of acoustic waves, generated by the lithospheric displacements during the seismic event preparation, with the cold dense ionospheric plasma in the presence of a magnetic field is constructed in conjunction with the actual height variations of the plasma particles' mobility and conductivity tensors in the ionosphere. A mechanism of creation of ionospheric inhomogeneities with different scales and with various degree of disturbance ΔN/N₀=10⁻⁴–10⁻⁵ (N₀ is the concentration of plasma electrons (ions) in background ionospheric plasma) is presented.     

Excitation of electromagnetic lower hybrid instability in the region of bouncing electron beams

The region of bouncing electron beams in the earth’s magnetosphere can be unstable against a non-resonant electromagnetic lower hybrid instability. The instability is purely growing in the rest frame of the plasma, and can be excited either by the temperature anisotropy or the drift velocity of the bouncing electron beams. The growth rates of the instability decrease with the increase of cold electron density. Consequently the growth rate is maximum at the equator where the cold electron density is minimum. The intense turbulence generated by this instability could broaden the bouncing electron beams thereby explaining the observed wider cone width of the beams at the equator. The instability could generate magnetic pulsations in the frequency range of orderPc ₁?Pc₃ with typical wavelength ≈ (3–10) km in the, magnetosphere during magnetic storms or substorms.     

On the correlation between spectra of solar microwave bursts and proton fluxes near the Earth

Studies of the extreme solar proton event of January 20, 2005 intensified the contest over of a long-standing problem: are solar cosmic rays arriving at the Earth accelerated by solar flares or by shocks preceding rapidly moving coronal mass ejections? Among the most important questions is the relationship between the energy spectra of the solar cosmic rays and the frequency spectra of flare microwave bursts. Some studies of previous solar-activity cycles have shown that such a relationship does exist, in particular, for protons with energies of tens of MeV. The present work analyzes this relation using data for 1987–2008. For flare events observed in the western half of the disk, there is a significant correlation between the index δ, which is equivalent to the power-law index of the integrated energy spectrum of 10–100 MeV protons detected near the Earth’s orbit, and radio burst parameters such as a ratio of peak fluxes S at two frequencies (for example, at 9 and 15 GHz) and a microwave peak frequency f _m . Proton fluxes with hard (flat) energy spectra (δ ≤ 1.5) correspond to hard microwave frequency spectra (S ₉/S ₁₅ ≤ 1 and f _m ≥ 15 GHz), while flares with soft radio spectra (S ₉/S ₁₅ ≥ 1.5 and f _m ≤ 5 GHz) result in proton fluxes with soft (steep) energy spectra (δ ≥ 1.5–2). It is also shown that powerful high-frequency bursts with the hardest radio spectra (f _m ≈ 30 GHz) can point at acceleration of significant proton fluxes in flares occurring in strong magnetic fields. These results argue that solar cosmic rays (or at least their initial impulses) are mainly accelerated in flares associated with impulsive and post-eruptive energy release, rather than in shocks driven by coronal mass ejections.     

Estimating hydraulic parameters of a heterogeneous aquitard using long-term multi-extensometer and groundwater level data

The classical aquitard-drainage model COMPAC has been modified to simulate the compaction process of a heterogeneous aquitard consisting of multiple sub-units (Multi-COMPAC). By coupling Multi-COMPAC with the parameter estimation code PEST++, the vertical hydraulic conductivity (K _v) and elastic (S _ske) and inelastic (S _skp) skeletal specific-storage values of each sub-unit can be estimated using observed long-term multi-extensometer and groundwater level data. The approach was first tested through a synthetic case with known parameters. Results of the synthetic case revealed that it was possible to accurately estimate the three parameters for each sub-unit. Next, the methodology was applied to a field site located in Changzhou city, China. Based on the detailed stratigraphic information and extensometer data, the aquitard of interest was subdivided into three sub-units. Parameters K _v, S _ske and S _skp of each sub-unit were estimated simultaneously and then were compared with laboratory results and with bulk values and geologic data from previous studies, demonstrating the reliability of parameter estimates. Estimated S _skp values ranged within the magnitude of 10^?4 m^?1, while K _v ranged over 10^?10–10^?8 m/s, suggesting moderately high heterogeneity of the aquitard. However, the elastic deformation of the third sub-unit, consisting of soft plastic silty clay, is masked by delayed drainage, and the inverse procedure leads to large uncertainty in the S _ske estimate for this sub-unit.