Spatial inhomogeneity of variations in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer critical frequencies under the conditions of winter 1995

The daily samples of the hourly measurements of the foF2 critical frequency, obtained on January 5–21, 1995, at the midlatitude and high-latitude automated ionospheric stations (geographic latitude higher than 60°), are considered. The {fo} sets are transformed into the {δfoF2} sets of relative variations, for which asymmetry (A) and excess (E) are calculated. The selected stations are grouped into 20 pairs of automated ionospheric stations (AISs) located at distances of 200–10 000 km from one another. Sign estimates of the cross-correlation coefficients between the sets of 16 A and E values for different pairs of stations are used. Two types of structures of the statistical invariant spatial distribution are established: the structures with a scale of about 300 km, invariant with respect to latitude, and with a scale of about 6000 km (for only high latitudes).     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

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.     

Estimation of the Maximum Earthquake Magnitude, <Emphasis Type="Italic">m</Emphasis><Subscript>max</Subscript>

This paper provides a generic equation for the evaluation of the maximum earthquake magnitude m_max for a given seismogenic zone or entire region. The equation is capable of generating solutions in different forms, depending on the assumptions of the statistical distribution model and/or the available information regarding past seismicity. It includes the cases (i) when earthquake magnitudes are distributed according to the doubly-truncated Gutenberg-Richter relation, (ii) when the empirical magnitude distribution deviates moderately from the Gutenberg-Richter relation, and (iii) when no specific type of magnitude distribution is assumed. Both synthetic, Monte-Carlo simulated seismic event catalogues, and actual data from Southern California, are used to demonstrate the procedures given for the evaluation of m_max.The three estimates of m_max for Southern California, obtained by the three procedures mentioned above, are respectively: 8.32 ± 0.43, 8.31 ± 0.42 and 8.34 ± 0.45. All three estimates are nearly identical, although higher than the value 7.99 obtained by Field et al. (1999). In general, since the third procedure is non-parametric and does not require specification of the functional form of the magnitude distribution, its estimate of the maximum earthquake magnitude m_max is considered more reliable than the other two which are based on the Gutenberg-Richter relation.     

Distribution and significance of C~_(40) alkanes in the extracts of Cretaceous source rocks from the Songliao Basin

Source rock extracts and crude oils from the Songliao Basin were analyzed by high-temperature gas chromatography (HTGC), gas chromatography-mass spectrometry (HTGC-MS) and gas chromatography-isotope ratio-mass spectrometry (GC-IRMS), for high molecular-weight alkanes. The distributions of n-alkanes in the Nenjiang Formation extracts are in the C14―C63 range; a bimodal distribution occurs in the C-21 and C21―40 regions. The C30―C37 n-alkanes are accompanied by C29―C35 hopanes, whereas the high molecular-weight C45―C47 n-alkanes co-occur with abundant isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The high δ 13C values of the n-alkanes and the microscopic maceral compositions indicate a highly diversified organic source input for the Nenjiang Formation source rocks, ranging from aquatic plants, blue alge-bacteria, to land plant material. In contrast, n-alkanes in the rock extracts of the Qingshankou Formation are characterized by a single modal distribution, with relatively low abundances of C29―C35 hopanes, but high molecular-weight isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The relatively low δ 13C values of C22―C44 n-alkanes and organic material compositions indicate that the source rocks in the Qingshankou Formation contain dominantly type I algal organic matter. The relative abundance of C 40 compounds in source rocks changes little at low maturity stage, but decreases drastically at higher maturity levels, with a concurrent reduction in the odd/even carbon predominance. In crude oils, in contrast, the relative abundance of C 40 compounds appears to relate closely with the oil source and oil viscosity.     

Estimation of velocity in the uppermost crust in a part of the western Gulf of Corinth,Greece, from the inversion of <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> arrival times using the neighbourhood algorithm

We determine the velocities in an upper crustal model, composed of three homogeneous layers, for one subregion of the western part of the Gulf of Corinth, NE of the town of Aigion, Greece. We have used local events that occurred there in the year 2001 and were recorded by the Corinth Rift Laboratory Network. Weighted P and S arrival time residuals are minimized using the Neighbourhood Algorithm of Sambridge (1999), combined with the grid search for source locations. The resolution of the inversion is tested by delete-one jackknifing. The model obtained is compared with some other models derived or applied to the subregion. A fast velocity increase between depths of 5 and 7 km is confirmed as the major structural element.     

On the interpretation of large gravimagnetic data by the modified method of <Emphasis Type="Italic">S</Emphasis>-approximations

The modified method of S-approximations applied for processing large and superlarge gravity and magnetic prospecting data is considered. The modified S-approximations of the elements of gravitational field are obtained due to the efficient block methods for solving the system of linear algebraic equations (SLAEs) to which the geophysically meaningful problem is reduced. The results of the mathematical experiment are presented.     

Experimental study on N<Subscript>2</Subscript>O and CH<Subscript>4</Subscript> fluxes from the dark coniferous forest zone soil of the Gongga Mountain,China

The static closed chamber technique is used in the study on the CH₄ and N₂O fluxes from the soils of primeval Abies fabri forest, the succession Abies fabri forest and the clear-cut areas of mid-aged Abies fabri forest in the Gongga Mountain from May 1998 to September 1999. The results indicate the following: (i) The forest soil serves as the source of atmospheric N₂O at the three measurement sites, while the fluxes of CH₄ are all negative, and soil is the sink of atmospheric CH₄. The comparative relations of N₂O emissions between the three sites are expressed as primeval Abies fabri forest < clear-cut areas < succession Abies fabri forest, and those of CH₄ consumption fluxes are primeval Abies fabri forest < succession Abies fabri forest < clear-cut areas, (ii) Significant seasonal variations of N₂O emission at various sites were observed, and two emission peaks of N₂O occurr during summer (July—August) and spring (February—March), whereas N₂O emission is relatively low in winter and spring (mid March—April). Seasonal variations of CH₄ consumption at each measurement site fluctuate drastically with unclear regularities. Generally, CH₄ consumption fluxes of succession Abies fabri forest and clear-cut areas are higher from mid May to late July but lower in the rest of sampling time, while the CH₄ flux keeps a relatively high value even up to September in primeval Abies fabri forest. In contrast to primeval Abies fabri forest, the CH₄ absorbabilities of succession Abies fabri forest and clear-cut areas of mid-aged Abies fabri forest are weaker. Particularly, the absorbability of the clear-cut areas is even weaker as compared with the other two sites, for the deforestation reduces the soil absorbability of atmospheric CH₄. (iii) Evident diurnal variation regularity exists in the N₂O emissions of primeval Abies fabri forest, and there is a statistic positive correlation between the fluxes of N₂O and air temperature (R=0.95, n=11, α· 0.01), and also the soil temperature of 5-cm layer (R=0.81, n=11, α> 0.01), whereas the CH₄ diurnal variation regularities are unclear and have no significant correlation with the soil temperature of 5-cm layer and air temperature.     

Updating the Probabilistic Seismic Hazard Values of Northern Algeria with the 21 May 2003 <Emphasis Type="Italic">M</Emphasis> 6.8 Algiers Earthquake Included

The occurrence of the Algiers earthquake (M 6.8) of May 21, 2003, has motivated the necessity to reassess the probabilistic seismic hazard of northern Algeria. The fact that this destructive earthquake took place in an area where there was no evidence of previous significant earthquakes, neither instrumental nor historical, strongly encourages us to review the seismic hazard map of this region. Recently, the probabilistic seismic hazard of northern Algeria was computed using the spatially smoothed seismicity methodology. The catalog used in the previous computation was updated for this review, and not only includes information until June 2003, but also considers a recent re-evaluation of several historical earthquakes. In this paper, the same methodology and seismicity models are utilized in an effort to compare this methodology against an improved and updated seismic catalog. The largest mean peak ground acceleration (PGA) values are obtained in northernmost Algeria, specifically in the central area of the Tell Atlas. These values are of the order of 0.48 g for a return period of 475 years. In the City of Algiers, the capital of Algeria, and approximately 50 km from the reported epicenter of this latest destructive earthquake, a new mean PGA value of 0.23 g is obtained for the same return period. This value is 0.07 g greater than that obtained in the previous computation. In general, we receive greater seismic hazard results in the surrounding area of Algiers, especially to the southwest. The main reason is not this recent earthquake by itself, but the significant increase in the m_max magnitude in the seismic source where the city and the epicenter are included.     

Parameters of the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer as a source of information on trends in thermospheric dynamics

The problem of determining trends in thermospheric dynamics parameters (horizontal winds) based on analysis of trends in various combinations of ionospheric F ₂-layer parameters is formulated. The previous attempts of the authors in this direction are briefly described. It is shown that all studied parameters lead to the same result: after the “boundary date” (approximately 1980) a systematic change in these parameters (a long-term trend) is observed, this fact manifesting changes in the dynamical regime of the thermosphere because of cooling and contraction of the entire middle and upper atmosphere. The results of a search for trends in the hmF2 height for the moment (T(ss) + 2 h) are described. These trends are found higher than the hmF2 trends obtained earlier by various authors analyzing the hmF2 behavior at fixed moments of local time.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

The instability of the <Emphasis Type="Italic">M</Emphasis><Subscript>max</Subscript> parameter and an alternative to its using

The work presents statistical methods for estimating the distribution parameters of rare, strong earthquakes. Using the two main theorems of extreme value theory (EVT), the distribution of T-maximum (the maximum magnitude over the time period T). Two methods for estimating the parameters of this distribution are proposed using the Generalized Pareto Distribution (GPD) and the General Extreme Value Distribution (GEV). In addition, the that allow the determination of the distribution of the T-maximum for an arbitrary value of T are proposed. The approach being used clarifies the nature of the instability of the widely accepted M max parameter. In the work, instead of unstable values of the M max parameter, the robust parameter Q _T (q), the q level quantile for the distribution of the T-maximum, is proposed to be used. The described method has been applied to the Harvard Catalogue of Seismic Moments of 1977–2006 and to the Magnitude Catalogue for Fennoscandia in 1900–2005. Moreover, the estimates of parameters of the corresponding GPD and GEV distributions, in particular, the most interesting shape parameter and the values of the M _max and Q _T (q) parameters are given.     

Improved method to determine the molar volume and compositions of the NaCl-H<Subscript>2</Subscript>O-CO<Subscript>2</Subscript> system inclusion

On the basis of Parry’s method (1986), an improved method was established to determine the molar volume (Vm) and compositions (X) of the NaCl-H2O-CO2 (NHC) system inclusion. To use this method, the determination of Vm-X only requires three microthermometric data of a NHC inclusion: partial homog-enization temperature (Th ,CO2), salinity (S) and total homogenization temperature (Th). Theoretically, four associated equations are needed containing four unknown parameters: X CO2, XNaCl, Vm and F (volume fraction of CO2 phase in total inclusion when occurring partial homogenization). When they are released, the Vm-X are determined. The former three equations, only correlated with Th ,CO2, S and F, have simplified expressions:XCO2=f1(Th,CO2,S,F),XNaCl=f2(Th,CO2,S,F),Vm=f3(Th,CO2,S,F). The last one is the thermodynamic relationship of X CO2, XNaCl, Vm and Th:f4(XCO2,XNaCl,Vm,Th)=0.Since the above four associated equations are complicated, it is necessary to adopt iterative technique to release them. The technique can be described by:(i) Freely input a F value (0≤F≤1),with Th ,CO2 and S, into the former three equations. As a result,X CO 2,XNaCl and the molar volume value recorded as Vm1 are derived. (ii) Input the X CO2 and XNaCl gotten in the step above into the last equation, and another molar volume value recorded as Vm2 is determined. (iii) If Vm1 is unequal to Vm2, the calculation will be restarted from “(i)”. The iteration is completed until Vm1 is equal to Vm2, which means that the four associated equations are released. Compared to Parry’s (1986) solution method, the improved method is more convenient to use, as well as more accurate to determine X CO 2. It is available for a NHC inlusion whose partial homogenization temperature is higher than clatherate melting temperature and there are no solid salt crystals in the inclusion at parital homogenization.     

Systematic examination of 'precursor waves' of strong earthquakes by using the wavelet method

Based on the borehole strainmeter recordings immediately before the strong earthquakes (&lt;i&lt;M&lt;/i&lt;&lt;sub&lt;S&lt;/sub&lt;≥?7.5) during the period from 2001 to 2005 observed at the Taian observation station, Shandong Province, the authors make a systematic and objective examination of the precursor waves of the quakes. The effects of Earth tides with periods larger than 128 min are eliminated through high-pass filtering; and atmospheric-pressure inferences are removed by linear regression. The 2?128 min signals are then separated into six frequency bands by employing the wavelet method. Results indicate that the wavelet method is capable of picking out information of weak variations in the signals. According to the characteristics of the 'precursor waves' obtained from wavelet transformation, the method of overrun ratio analysis is put forward for examination. All the detailed components of the wavelets have been analyzed. For the time series of the overrun rate in all these components, statisitical calculation has been made for the slopes of fit curves, and mean values and standard deviations were obtained and positive-negative slope ratios were analyzed. The three statistical data show that 'precursor waves' are not widely recorded by borehole strainmeter within 15 days before remote strong earthquakes.     

Model prediction of the most effective frequency for the large-scale modification of the midlatitude ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer by powerful HF radiowaves

The response of the midlatitude F ₂ layer to the effect of powerful HF radiowaves is studied using the numerical model of the ionosphere. The large-scale modification of the F ₂ layer over the Sura heating facility near Nizhni Novgorod is considered for autumnal conditions. The calculations are performed for various cases when the heating wave has different frequencies under the daytime and nighttime conditions. The calculation results show that large-scale changes in the electron temperature and density in the F ₂ layer caused by the artificial heating should substantially depend on the heating radiowave frequency. It is found that there should exist such, most effective, heating wave frequency at which a decrease in the electron density at the F ₂ layer maximum height over the heating facility should be maximal.     

Relation of long-period variations in the critical frequencies of the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer to geomagnetic activity

The relation of the long-period variations in the midnight and noon values of the critical frequency of the ionospheric F ₂ layer at three midlatitude stations (Irkutsk, Moscow, and Boulder) to the daily mean index of geomagnetic activity in years of different solar activity has been studied. It has been found that the correlation coefficients between the above parameters depend on time of day, season, and solar activity level. The correlation coefficients are higher at night than in the daytime, especially at low solar activity. The highest absolute values of the correlation coefficient most often appear during equinoxes: April–May and September–October. It has been shown that the variability of the critical frequencies of the midlatitude ionospheric F ₂ layer depends not only on geomagnetic activity but also (to a considerable degree) on the effect of the lower atmosphere.     

Geochemical evidence of the renewal of volcanic activity inferred from CO<Subscript>2</Subscript> soil and SO<Subscript>2</Subscript> plume fluxes: the 2007 Stromboli eruption (Italy)

On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO₂ and SO₂) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO₂ flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO₂ fluxes of the degassing plume and the CO₂ fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO₂ flux) and in soil degassing (CO₂ flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.     

The crystal structure of (Fe<Subscript>4</Subscript>Cr<Subscript>4</Subscript>Ni)<Subscript>9</Subscript>C<Subscript>4</Subscript>

(Fe4Cr4Ni)9C4 is a metal carbide mineral formed by combination of Fe, Cr and Ni with C. It occurs in a chromite deposit in the Luobusha ophiolite, Tibet. Based on the determination of its crystal structure, the empirical formula is (Fe4.12Cr3.84Ni0.96)8.92C3.70 and the simplified formula is (Fe4Cr4Ni)4C9. The mineral is hexagonal with a = 1.38392(2) nm, c = 0.44690(9) nm, pace group P63 m c, Z=6 and the calculated specific gravity Dx = 7.089 g/cm3. Fe, Cr and Ni occupy different crystallographic sites and their coordination numbers are approximately 12, forming an alternate stacking sequence of flat and puckered layers along the c axis. Some metallic atoms have a defect structure. The interatomic distances of Fe, Cr and Ni are 0.2525-0.2666 nm, and the distances between Fe, Cr, Ni and C are 0.1893-0.2169 nm. The coordination number of carbon is 6. It occurs in interstices of the metallic atoms Fe, Cr and Ni to form trigonalprismatically coordinated polyhedra. These coordination polyhedra are linked with each other via shared corners or shared edges into a new type of metal carbide structure.     

Vibrationally excited N<Subscript>2</Subscript> and O<Subscript>2</Subscript> in the upper atmosphere: A review

The present-day state of studies of vibrationally excited N₂ and O₂ in the upper atmosphere and the mechanisms by which these molecules affect the electron density and temperature in the ionosphere and neutral densities in the upper atmosphere are considered.