Carriers of magnetism and vertical inhomogeneities in ultramafic-mafic rocks of the Kondyor Massif (Aldan shield)

Dunite samples from a borehole drilled in the platiniferous concentrically-zoned Kondyor Massif are studied by electron spin resonance (ESR). The spectrum profiles, relative intensities I, and volume magnetic susceptibilities κ are analyzed. These values experience strong irregular variations, sometimes by an order of magnitude, in the upper and medium parts of the column, at depths from 100 to 400 m; and the variations decay at greater depths. The magnetic properties of the samples are determined by iron (II) ions in the olivine lattice and by iron (III) ions in the magnetite and pyrrhotite microphases and in the products of breakdown of the solid solution: chromiferous magnetite, chromoferrite, etc. The I and κ values are directly related: κ_max = 27.8 × 10⁻³ SI units, κ_min=2.63 × 10⁻³ SI units, and κ_mean = 12.7 × 10⁻³ SI units. The maximum κ values are found in the zones with elevated contents of magnetite and pyrrhotite particles, and the minimum ones, in zones with few medium and small clusters with Fe³⁺ ions. The uneven distributions of solid solutions and magnetic phases over depths are suggested to be related to the disturbances in the conditions of crystallization.     

Dependence of geomagnetic disturbances on extreme values of the solar wind <Emphasis Type="Italic">E</Emphasis><Subscript><Emphasis Type="Italic">y</Emphasis></Subscript> component

The dependence of the zonal geomagnetic indices (AE, Ap, Kp, Kn, and Dst) on the solar wind parameters (the electric field E _y component, dynamic pressure P _d and IMF irregularity σB) has been studied for two types of events: magnetic clouds and high-speed streams. Based on the empirical relationships, it has been established that the AE, Ap, Kp, and Kn indices are directly proportional to the E _y value at E _y < 12 mV m^?1 and are inversely proportional to this value at E _y > 12 mV m^?1 for the first-type events. On the contrary, the dependence of Dst on E _y is monotonous nonlinear. A linear dependence of all geomagnetic indices on E _y is typical of the second-type events. It has been indicated that the specific features of geoeffectiveness of magnetic clouds and high-speed solar wind streams are caused by the dependence of the electric field potential across the polar cap on the electric field, solar wind dynamic pressure, and IMF fluctuations.     

Sorption of Eu<Superscript>3+</Superscript>onto nano-size silica-water interfaces

The sorption of Eu species onto nano-size silica-water interfaces is investigated at pH range of 1―8.5 and the initial Eu concentrations (C_Eu) of 2×10⁻⁵, 2×10⁻⁴ and 2×10⁻³ M using fluorescence spectroscopy. The sorption rate of Eu is initially low, but significantly increases at pH > 4. For the initial C_Eu of 2×10⁻⁵, 2×10⁻⁴ and 2×10⁻³ M, the dissolved Eu species are completely sorbed onto silica-water interfaces at pH = 4.75, −5.8 and 6.6, respectively, with the respective sorption densities of −1.58×10⁻⁸, 1.58×10⁻⁷ and 1.58×10⁻⁶ mol/m². The sorbed Eu species at pH < 6 is aquo Eu³⁺, which is sorbed onto silica-water interfaces as an outer-sphere complex at pH < 5, but may be sorbed as an inner-sphere bidentate complex at 5 < pH < 6, due to the decrease of the N_H2O to −6 at pH = 6. At pH = 6 – 8, Eu(OH)²⁺, Eu(CO₃)⁺and Eu(CO₃)₂ ⁻ form in the solutions, and Eu(CO₃)⁺is dominant at pH = −7.5. These ions may be sorbed onto silica-water interfaces as inner-sphere bidentate complexes or multi-nuclear pre-cipitates.     

Region-specific <Emphasis Type="Italic">S</Emphasis>-wave attenuation for earthquakes in northwestern Iran

In this study, continuous wavelet transform is applied to estimate the frequency-dependent quality factor of shear waves, Q _S , in northwestern Iran. The dataset used in this study includes velocigrams of more than 50 events with magnitudes between 4.0 and 6.5, which have occurred in the study area. The CWT-based method shows a high-resolution technique for the estimation of S-wave frequency-dependent attenuation. The quality factor values are determined in the form of a power law as Q _S (f)?=?(147?±?16)f ^0.71?±?0.02 and (126?±?12)f ^0.73?±?0.02 for vertical and horizontal components, respectively, where f is between 0.9 and 12 Hz. Furthermore, in order to verify the reliability of the suggested Q _S estimator method, an additional test is performed by using accelerograms of Ahar-Varzaghan dual earthquakes on August 11, 2012, of moment magnitudes 6.4 and 6.3 and their aftershocks. Results indicate that the estimated Q _S values from CWT-based method are not very sensitive to the numbers and types of waveforms used (velocity or acceleration).     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

Anomalous variations in the structure of the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region at geomagnetic midlatitudes of the Southern and Northern hemispheres in going from summer to winter conditions at high solar activity

The structure and dynamics of the ionosphere and plasmasphere at high solar activity under quiet geomagnetic conditions of June 2–3, 1979, and January 5–6, 1980, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The plasma drift velocity, determined by comparing the calculated and measured heights of the F ₂ layer maximum (hmF2), and the correction of [N₂] and [O₂], found in the NRLMSISE-00 model, make it possible to coordinate the electron densities (NmF2) calculated at the hmF2 height and the measured anomalous variations in NmF2 over the Argentine Islands ionosonde as well as the calculated and measured NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that, if the interference of the diffusion velocities of O⁺(⁴S) and H⁺ ions is taken into account, the additional heating of plasmaspheric electrons leads to an increase in the flux of O⁺(⁴S) ions from the topside ionosphere to lower F ₂ layer altitudes, as a result of which an anomalous nighttime increase in NmF2 6, observed on January 6, 1980, over Millstone Hill station, is mainly produced. The second component of the formation of anomalous night-time NmF2 is the plasma drift along the magnetic field caused by the neutral wind, which shifts O⁺(⁴S) ions to higher altitudes where the recombination rate of O⁺(⁴S) with N₂ and O₂ is lower and slows down a decrease in NmF2 in the course of time. It has been shown that the influence of electronically excited O⁺ ions and vibrationally excited N₂ and O₂ molecules on electron density (N _e ) considerably differs under winter and summer conditions. This difference forms significant part of the winter anomaly in N _e at heights of the F ₂ region and topside ionosphere over Millstone Hill station.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

CO<Subscript>2</Subscript> flux evaluation over the evergreen coniferous and broad-leaved mixed forest in Dinghushan,China

The Dinghushan flux observation site, as one of the four forest sites of ChinaFLUX, aims to acquire long-term measurements of CO₂ flux over a typical southern subtropical evergreen coniferous and broad-leaved mixed forest ecosystem using the open path eddy covariance method. Based on two years of data from 2003 to 2004, the characteristics of temporal variation in CO₂ flux and its response to environmental factors in the forest ecosystem are analyzed. Provided two-dimensional coordinate rotation, WPL correction and quality control, poor energy-balance and underestimation of ecosystem respiration during nighttime implied that there could be a CO₂ leak during the nighttime at the site. Using daytime (PAR > 1.0 μmol^?1·m^?2·s^?1) flux data during windy conditions (u* > 0.2 m·s^?1), monthly ecosystem respiration (Reco) was derived through the Michaelis-Menten equation modeling the relationship between net ecosystem C0₂ exchange (NEE) and photosynthetically active radiation (PAR). Exponential function was employed to describe the relationship between Reco and soil temperature at 5 cm depth (Ts05), then Reco of both daytime and nighttime was calculated respectively by the function. The major results are: (i) Derived from the Michaelis-Menten equation, the apparent quantum yield (α) was 0.0027±0.0011 mgCO₂·μmol^?1 photons, and the maximum photosynthetic assimilation rate (Amax) was 1.102±0.288 mgCO₂·m^?2·s^?1. Indistinctive seasonal variation of α or Amax was consistent with weak seasonal dynamics of leaf area index (LAf) in such a lower subtropical evergreen mixed forest, (ii) Monthly accumulated Reco was estimated as 95.3±21.1 gC·m^?2mon^?1, accounting for about 68% of the gross primary product (GPP). Monthly accumulated WEE was estimated as ?43.2±29.6 gC·m^?2·mon^?1. The forest ecosystem acted as carbon sink all year round without any seasonal carbon efflux period. Annual NEE of 2003 and 2004 was estimated as ?563.0 and ?441.2 gC·m^?2·a^?1 respectively, accounting for about 32% of GPP.     

Anomalous variations in the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript>-layer structure at geomagnetic midlatitudes of the Southern and Northern hemispheres at the transition from summer to winter conditions under low solar activity

The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N₂ and O₂ molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O⁺ on N _e ions under winter and summer conditions forms not more than 11% of the N _e winter anomaly event in the F ₂ layer and topside ionosphere. The model without electronically excited O⁺ ions reduces the duration of the N _e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.     

Gradual late stage deepening of Gega ice-dammed lake,Tsangpo gorge,southeastern Tibet,indicated by preliminary sedimentary rock magnetic properties

Gega lake, in southeastern Tibet, was formed by the blocking of Yarlung Tsangpo gorge by a glacier and is a well-known example of geomorphological damming. However, the evolution of the damming process at the site is still not understood in detail. Here, we use measurements of multiple magnetic parameters of the sediments from the Yusong (YS) 3 section, which is well-dated by optically stimulated luminescence, to provide a detailed history of the late stage of Gega dammed lake since 17.0 ka. Low-frequency field magnetic susceptibility (χ_lf) increases upwards gradually from 25 to 79?×?10^?8 m³/kg above 5.5 m, but other magnetic properties, such as frequency-dependent magnetic susceptibility (χ_fd % and χ_fd), susceptibility of anhysteretic remanent magnetization (χ_ARM), and saturation isothermal remanent magnetization (SIRM) did not show a similar degree of enhancement. The magnetic grain size indicators of χ_ARM/χ_lf, χ_ARM/SIRM, and χ_lf/SIRM all indicate a trend of increasing magnetic grain size from the 5.5 m to the top of section; however, the bulk sediment grain size decreases gradually within the same interval. The total organic carbon is very low (0.2–0.7%), and thus, it is unlikely that the sedimentary environment is sufficiently strongly reducing to lead to the dissolution of magnetic minerals. Therefore, we infer that the coarsening-upwards of the magnetic grain size, and the increasing magnetic susceptibility from 5.5 m, reflect the gradual preferential preservation of magnetic minerals, caused by the deepening of the lake. Thus, we interpret the magnetic record of the section as reflecting the gradual deepening of the late stage of Gega ice-dammed lake. Last, the abrupt disappearance of the lake may have resulted in an outburst flood.     

Generation of superDreicer electric fields in the solar chromosphere

The electric field generation at the front of the current pulse, which originates in a coronal magnetic loop owing to the development of the Rayleigh–Taylor magnetic instability at loop footpoints, has been considered. During the τ_A ≈ l/V _A ≈ 5?25 s time (where l is the plasma plume height entering a magnetic loop as a result of the Rayleigh–Taylor instability), a disturbance related to the magnetic field tension B _?(r,t), “escapes” the instability region with the Alfvén velocity in this case. As a result, an electric current pulse I_z(z ? V _A t), at the front of which an induction magnetic field E _z, which is directed along the magnetic tube axis and can therefore accelerate particles, starts propagating along a magnetic loop with a characteristic scale of Δξ ≈ l. In the case of sufficiently large currents, when B _? ²/8π > p, an electric current pulse propagates nonlinearly, and a relatively large longitudinal electric field originates E _z ≈ 2I _z ³ V _A/c ⁴a²B_z ²l, which can be larger than the Dreicer field, depending on the electric current value.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

Seasonal and annual variation of CO<Subscript>2</Subscript> flux above a broad-leaved Korean pine mixed forest

Long-term measurement of carbon metabolism of old-growth forests is critical to predict their behaviors and to reduce the uncertainties of carbon accounting under changing climate. Eddy covariance technology was applied to investigate the long-term carbon exchange over a 200 year-old Chinese broad-leaved Korean pine mixed forest in the Changbai Mountains (128°28′E and 42°24′N, Jilin Province, P. R. China) since August 2002. On the data obtained with open-path eddy covariance system and CO₂ profile measurement system from Jan. 2003 to Dec. 2004, this paper reports (i) annual and seasonal variation of F _NEE, F _GPP and R _E; (ii) regulation of environmental factors on phase and amplitude of ecosystem CO₂ uptake and release Corrections due to storage and friction velocity were applied to the eddy carbon flux.L_AI and soil temperature determined the seasonal and annual dynamics of F_GPP and R_E separately. V_PD and air temperature regulated ecosystem photosynthesis at finer scales in growing seasons. Water condition at the root zone exerted a significant influence on ecosystem maintenance carbon metabolism of this forest in winter.The forest was a net sink of atmospheric CO₂ and sequestered ?449 g C·m^?2 during the study period; ?278 and ?171 gC·m^?2 for 2003 and 2004 respectively. F _GPP and F _RE over 2003 and 2004 were ?1332, ?1294 g C·m^?2. and 1054, 1124 g C·m^?2 respectively. This study shows that old-growth forest can be a strong net carbon sink of atmospheric CO₂.There was significant seasonal and annual variation in carbon metabolism. In winter, there was weak photosynthesis while the ecosystem emitted CO₂. Carbon exchanges were active in spring and fall but contributed little to carbon sequestration on an annual scale. The summer is the most significant season as far as ecosystem carbon balance is concerned. The 90 days of summer contributed 66.9, 68.9% of F _GPP, and 60.4, 62.1% of R _E of the entire year.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Seasonal dynamics of soil CO<Subscript>2</Subscript> effluxes with responses to environmental factors in lower subtropical forests of China

Seasonal metrics and environmental responses to forestry soil surface CO₂ emission effluxes among three types of lower subtropical forests were consistently monitored over two years with static chamber-gas chromatograph techniques among three types of lower subtropical forests. Results showed that annual CO₂ effluxes (S+L) reached 3942.20, 3422.36 and 2163.02 CO₂ g·m^?2·a^?1, respectively in the monsoon evergreen broadleaf forest, mixed broadleaf-coniferous forest and coniferous forest. All the three types of forests revealed the same characteristics of seasonal changes with the CO₂ effluxes peaking throughout June to August. During this peaking period, the effluxes were 35.9%, 38.1% and 40.2% of the total annual effluxes, respectively. The CO₂ emission process responding to the environmental factors displayed significantly different patterns in forestry soils of the three types of forests. The coniferous forest (CF) was more sensitive to temperature than the other two types. The Q ₁₀ values were higher, along with greater seasonal variations of the CO₂ efflux, indicating that the structurally unique forestry ecosystem has disadvantage against interferences. All the three types of forestry CO₂ effluxes showed significant correlation with the soil temperature (T _s), soil water content (M _s) and air pressure (P _a). However, stepwise regression analysis indicated no significant correlation between air pressure and the soil CO₂ efflux. With an empirical model to measure soil temperature and water content in 5 cm beneath the soil surface, the CO₂ effluxes accounting for 75.7%, 77.8% and 86.5% of the efflux variability respectively in soils of BF, MF and PF were calculated. This model can be better used to evaluate the CO₂ emission of soils under water stress and arid or semi-arid conditions.     

Estimation of storm-time level of day-side wave geomagnetic activity using a new <Emphasis Type="Italic">ULF</Emphasis> index

The level of wave geomagnetic activity in the morning and daytime sectors of auroral latitudes during strong magnetic storms with Dst _min varying from ?100 to ?150 nT in 1995–2002 have been studied using a new ULF index of wave activity proposed in [Kozyreva et al., 2007]. It has been detected that daytime Pc5 pulsations (2–6 mHz) are most intense during the main phase of a magnetic storm rather than during the recovery phase as was considered previously. It has been indicated that morning geomagnetic pulsations during the substorm recovery phase mainly contribute to daytime wave activity. The appearance of individual intervals with the southward IMF B _z component during the magnetic storm recovery phase results in increases in the ULF index values.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

Seasonal variation in carbon exchange and its ecological analysis over <Emphasis Type="Italic">Leymus chinensis</Emphasis> steppe in Inner Mongolia

Eddy covariance technique was used to measure carbon flux during two growing seasons in 2003 and 2004 over typical steppe in the Inner Mongolia Plateau, China. The results showed that there were two different CO₂ flux diurnal patterns at the grassland ecosystem. One had a dual peak in diurnal course of CO₂ fluxes with a depression of CO₂ flux after noon, and the other had a single peak. In 2003, the maximum diurnal uptake and emitting value of CO₂ were ?7.4 and 5.4 g·m^?2·d^?1 respectively and both occurred in July. While in 2004, the maximum diurnal uptake and release of CO₂ were ?12.8 and 5.8 g·m^?2·d^?1 and occurred both in August. The grassland fixed 294.66 and 467.46 g CO₂·m^?2 in 2003 and 2004, and released 333.14 and 437.17 g CO₂·m^?2 in 2003 and 2004, respectively from May to September. Water availability and photosynthetic active radiation (PAR) are two important factors of controlling CO₂ flux. Consecutive precipitation can cause reduction in the ability of ecosystem carbon exchange. Under favorable soil water conditions, daytime CO₂ flux is dependent on PAR. CO₂ flux, under soil water stress conditions, is obviously less than those under favorable soil water conditions, and there is a light saturation phenomena at PAR=1200 μmol·m^?2·s^?1. Soil respiration was temperature dependent when there was no soil water stress; otherwise, this response became accumulatively decoupled from soil temperature.