The crystal structure of (Fe_4Cr_4Ni)_9C_4

The Luobusha podiform chromite deposit occursin the Luobusha ophiolite exposed at the In-dus-Yarlung Zangbo suture[1]. The wall rocks of chro-mitites are harzburgites. Many mantle minerals werefound in chromitites and harzburgites. Among of them,there are 120 grains of diamond[2], many elementalmetals (mainly Au, Cu, Fe, Ni, Cr, Al, W, Zn, Pb, Sn,Os, Ir, Ru, C, etc.) and their intermetallic compounds.Besides diamond and graphite, the element C isalso combined with Fe, Cr, Ni, Ti and …     

Mineralogy and crystal structure determination of Mg-fillowite

Accurate crystal structure determination of complex phosphate Mg-fillowite.—Na₂Ca(Mn₄Mg₂Fe)₇[PO₄]₆ has been finished. Mg-fillowite occurs in muscovite-pegmatite in the Altay area, Xinjiang, China. It is of the trigonal lattice, with space group R-3, unit cell a=1.5143(3)nm, c=4.3191(2) nm, V=8.5736 nm³, Z = 18. The R-factor of the determination accuracy is R (I<2 σ(I)) = 0.0776. The cationic polyhedra consist of two kinds of structural units: compound column and screw column. In three-dimensional space, they are linked with each other through corner sharing or edge sharing to build the entire structure. There are 45 atoms in an asymmetric unit and 720 atoms in a unit cell. In terms of types and connection manners of the coordination polyhedra, the degree of complicity is rarely seen in all mineral crystal structures. Study of such a structure is significant for the classification of phosphate and the building rule of coordination polyhedra in a complex crystal structure.     

Variations in carbon isotope ratios of C<Subscript>3</Subscript> plants and distribution of C<Subscript>4</Subscript> plants along an altitudinal transect on the eastern slope of Mount Gongga

Variations in carbon isotopic ratios (δ ¹³C) of C₃ plants and distribution of C₄ plants were investigated along an altitudinal transect on the eastern slope of Mount Gongga, and the environmental effects on them were discussed. It is shown that plants with C₄ photosynthetic pathway mainly occur at altitudes below 2100 m a.s.l., suggesting that the low summer temperature is responsible for the distributional pattern. In addition, δ ¹³C of C₃ plants increases with elevation at the region above 2000 m a.s.l. with the characteristics of humid climate, and the increase rate in δ ¹³C for C₃ plants is about 1.3% per kilometer. Temperature determines the altitudinal trend of δ ¹³C.     

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.     

The 20-second regional magnitude <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">S</Emphasis></Subscript>(20R) for the Russian Far East

A new modified magnitude scale M _S (20R) is elaborated. It permits us to extend the teleseismic magnitude scale M _S (20) to the regional epicenter distances. The data set used in this study contains digital records at 12 seismic stations of 392 earthquakes that occured in the northwest Pacific Ocean in the period of 1993–2008. The new scale is based on amplitudes of surface waves of a narrow range of the periods (16–25 s) close to the period of 20 s, for distances of 80–3000 km. The digital Butterworth filter is used for processing. On the basis of the found regional features concerning distance dependence for seismic wave attenuation, all the stations of the region have been subdivided into two groups, namely, “continental” and “island-arc.” For each group of stations, its own calibration function is proposed. Individual station corrections are used to compensate for the local features.     

<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.     

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.     

Dispersion Measurements of <Emphasis Type="Italic">P</Emphasis> Waves and their Implications for Mantle <Emphasis Type="Italic">Q</Emphasis><Subscript><Emphasis Type="Italic">p</Emphasis></Subscript>

We analyze the anelasticity of the earth using group delays of P-body waves of deep (>200 km) events in the period range 4–32 s for epicentral distances of 5–85 degrees. We show that Time Frequency Analysis (TFA), which is usually applied to very dispersive surface waves, can be applied to the much less dispersive P-body waves to measure frequency-dependent group delays with respect to arrival times predicted from the CMT centroid location and PREM reference model. We find that the measured dispersion is due to: (1) anelasticity (described by the P-wave quality factor Q _p ), (2) ambient noise, which results in randomly distributed noise in the dispersion measurements, (3) interference with other phases (triplications, crustal reverberations, conversions at deep mantle boundaries), for which the total dispersion depends on the amplitude and time separation between the different phases, and (4) the source time function, which is dispersive when the wavelet is asymmetrical or contains subevents. These mechanisms yield dispersion ranging in the order of one to 10 seconds with anelasticity responsible for the more modest dispersion. We select 150 seismograms which all have small coda amplitudes extending to ten percent of the main arrival, minimizing the effect of interference. The main P waves have short durations, minimizing effects of the source. We construct a two-layer model of Q _p with an interface at 660 km depth and take Q _p constant with period. Our data set is too small to solve for a possible frequency dependence of Q _p . The upper mantle Q ₁ is 476 [299–1176] and the lower mantle Q ₂ is 794 [633–1064] (the bracketed numbers indicate the 68 percent confidence range of Q _p ^–1). These values are in-between the AK135 model (Kennett et al., 1995) and the PREM model (Dziewonski and Anderson, 1981) for the lower mantle and confirm results of Warren and Shearer (2000) that the upper mantle is less attenuating than PREM and AK135.     

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.     

Structure refinement of astrophyllite

The crystal structure of astrophyllite K₂Na(Fe, Mn, Mg,□)₇[Ti₂(Si₄O₁₂)₂|O₃](OH, F)₄ has been refined. The dimensions of the triclinic unit cell are: a = 0.5359(2) nm,b = 1.1614(4) nm, c = 1.1861(4) nm, α= 113.16(2)°, β= 103.04(2)°,γ= 94.56(2)°,V = 0.6495(5) nm³, Z= 1, space group P1, R=0.057 for 5308 reflections |F_o|>3σ|F_o|. According to structural and compositional differences the monoclinic astrophyllite K₂NaNa(Fe, Mn)₄Mg₂Ti2[Si₄O₁₂]₂(OH)₄(OH, F)₂ and astrophyllite should be considered as two different mineral species. Astrophyllite, monoclinic astrophyllite, bafertisite and lamprophyllite contain heteropolyhedral sheets which topologically are related with Si, O sheets of mica where one or several SiO₄ tetrahedra are replaced by TiO _n polyhedra. Therefore this heterophyllotitanosilicate series represents a kind of functional substitution in inorganic crystals.     

Response of dayside auroras to abrupt increases in the solar wind dynamic pressure at positive and negative polarity of the IMF <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> component

The optical observations on Heiss Island (Φ′ = 75.0°) have been used to study the characteristics of auroras in the near-noon MLT sector after abrupt increases in the solar wind dynamic pressure at negative and positive polarity of the IMF B _z component. It has been found out that the 427.8 and 557.7 nm emission intensities considerably increased at B _z < 0 both equatorward of the dayside red luminosity band and within this band. The value of the emission intensities at a red luminosity maximum (I ₆₃₀₀/I ₅₅₇₇ ∼ 0.5) indicates that energetic electron precipitation is of the magnetospheric origin. At B _z > 0, fluxes of harder (E > 1 keV) precipitating electrons were superimposed on the soft spectrum of precipitating particles in the equatorial part of the red luminosity band. This red band part was hypothetically caused by the low-latitude boundary layer (LLBL) on closed lines of the geomagnetic field, the estimated thickness of which is ∼3 R _e . The 557.7 nm emission intensity increased during 3–5 min after SC/SI and was accompanied by the displacement of the red band equatorward boundary toward lower latitudes. The displacement value was ∼150–200 km when the dynamic pressure abruptly increased by a factor of 3–5. After SC/SI, the 630.0 nm emission intensity continued increasing during 16–18 min. It is assumed that the time of an increase in the red line intensity corresponds to the time of saturation of the magnetospheric boundary layers with magnetosheath particles after an abrupt increase in their density.     

Experimental study on supercritical phenomena of WO<Subscript>3</Subscript> solubility in NaCl-H<Subscript>2</Subscript>O system

The solubility measurements of WO₃ in NaCl-H₂O system with 4.0% by weight salt were carried out to study the supercritical phenomena of solubility. These experiments were carried out using rapid-quench pressure vessels, with quantitative Na₂WO₄-2H₂O, HCl, NaCl and H₂O in Pt capsule to determine the solubility of WO₃ from supersaturated solution. The pressure was kept at 34 MPa, which was near the critical pressure 31.4 MPa and the temperatures varied from 250°C to 550°C. The experimental results indicate that the solubility of WO₃ has the features of supercritical phenomenon in critical region and is sensitive to the change of system temperature and solvent density. The temperature and pressure in this experiment are close to those of ore-forming fluids in many tungsten deposits, so the features of supercritical geofluids can be helpful to investigating the ore-forming mechanism of tungsten deposit.     

Elements in size-fractionated bottom sediments of the Elbe River in its Czech part

The concentrations of 45 elements, and the content of organic substances, amorphous inorganic substances and minerals were determined in 26 samples of Elbe River bottom sediments to determine the correlation between the element composition and the grain size in the ￡ \le 4, 4--8, 8--16, 16--32, 32--63 7m and bulk samples ￡ \le 63 7m fractions. The purpose of this study was to obtain information about the role of particle size and their mineral composition as a control of element concentrations in sediments, the potential remobilization of elements under changing environmental conditions, and the chemical speciation of the large excess concentrations of elements in this system.¶About 22--56% of the estimated element yield was associated with sediments with a grain size of less than 4 7m and 82--97% of the element yield was found in the fraction smaller than 32 7m. Sequential extraction carried out for 20 potentially reactive elements in the ￡ \le 63 7m bulk samples and the ￡ \le 4 7m fraction showed the increased importance of exchangeable ions in the order Mn = Ni < Mg < Sr < Cu < Zn < Ca < K < Cd < Na < Mo, reducible fractions Zn < Sr < Mo < Co < Cu < Ni < Pb < Mn < Cr < Be < V = Fe < As < Ag, oxidizable chemical phases As = Ca < Cd < K < Sr < Fe = Mn < Zn = Mg < Mo < Co = Ag < Pb < Li < Cu = V < Ni < Cr = Na < Be with small concentrations from lithogenic phases Mn < As < Pb < Sr < Zn = Cu = Na < Be < V < Ni < Cr < Fe < K < Mg < Li. The major accumulative phases were identified for the individual elements.     

Spatial distribution and bonding forms of heavy metals in sediments along the middle course of the River Elbe (km 287ċ390)

Sediments of eight groyne fields along the middle course of the River Elbe (river km 287?390) were geochemically studied. The 78 sediment samples were analysed for pH and grain size distribution. The grain size fraction < 2 μm was used for mineralogical and chemical analysis: semiquantitative clay mineral analysis; total element content (Al, Ca, Fe, K, Mn, Cd, Cr, Cu, Ni, Pb, Zn), total inorganic and organic C, and bonding form fractionation with a six‐step sequential chemical extraction. The latter was performed on selected samples (n = 32). The sediments along the Elbe's course are still contaminated with heavy metals far above the local geogenic background level. An enrichment factor of more than 15 was calculated for zinc. Cr and Ni are the elements with the lowest enrichment. The bonding form analysis of selected heavy metals shows a dominance of relatively immobile bonding forms, e. g. the moderately reducible and the residual fraction, which implies a relative low mobility potential. Only Zn poses a higher potential threat to the environment, since it has a higher percentage of the first three extracted phases: adsorbed, carbonate, and easily reducible fraction.     

The moment magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> and the energy magnitude <Emphasis Type="Italic">M</Emphasis><Subscript>e</Subscript>: common roots and differences

Starting from the classical empirical magnitude-energy relationships, in this article, the derivation of the modern scales for moment magnitude M _w and energy magnitude M _e is outlined and critically discussed. The formulas for M _w and M _e calculation are presented in a way that reveals, besides the contributions of the physically defined measurement parameters seismic moment M ₀ and radiated seismic energy E _S, the role of the constants in the classical Gutenberg–Richter magnitude–energy relationship. Further, it is shown that M _w and M _e are linked via the parameter Θ = log(E _S/M ₀), and the formula for M _e can be written as M _e = M _w + (Θ + 4.7)/1.5. This relationship directly links M _e with M _w via their common scaling to classical magnitudes and, at the same time, highlights the reason why M _w and M _e can significantly differ. In fact, Θ is assumed to be constant when calculating M _w. However, variations over three to four orders of magnitude in stress drop Δσ (as well as related variations in rupture velocity V _R and seismic wave radiation efficiency η _R) are responsible for the large variability of actual Θ values of earthquakes. As a result, for the same earthquake, M _e may sometimes differ by more than one magnitude unit from M _w. Such a difference is highly relevant when assessing the actual damage potential associated with a given earthquake, because it expresses rather different static and dynamic source properties. While M _w is most appropriate for estimating the earthquake size (i.e., the product of rupture area times average displacement) and thus the potential tsunami hazard posed by strong and great earthquakes in marine environs, M _e is more suitable than M _w for assessing the potential hazard of damage due to strong ground shaking, i.e., the earthquake strength. Therefore, whenever possible, these two magnitudes should be both independently determined and jointly considered. Usually, only M _w is taken as a unified magnitude in many seismological applications (ShakeMap, seismic hazard studies, etc.) since procedures to calculate it are well developed and accepted to be stable with small uncertainty. For many reasons, procedures for E _S and M _e calculation are affected by a larger uncertainty and are currently not yet available for all global earthquakes. Thus, despite the physical importance of E _S in characterizing the seismic source, the use of M _e has been limited so far to the detriment of quicker and more complete rough estimates of both earthquake size and strength and their causal relationships. Further studies are needed to improve E _S estimations in order to allow M _e to be extensively used as an important complement to M _w in common seismological practice and its applications.     

Dayside auroras during unusually large positive values of the IMF <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> component

The characteristics of dayside auroras during the large (16–24 nT) positive values of the IMF B _z component, observed on January 14, 1988, during the interaction between the Earth’s magnetosphere and the body of the interplanetary magnetic cloud, have been studied based on the optical observations on Heiss Island. A wide band of diffuse red luminosity with an intensity of 1–2 kilorayleigh (kR) was observed during 6 h in the interval 1030–1630 MLT at latitudes higher than 75° CGL. Rayed auroral arcs, the brightness of which in the 557.7 nm emission sharply increased to 3–7 kR in the postnoon sector immediately after the polarity reversal of the IMF B _y component from positive to negative, were continuously registered within the band. Bright auroral arcs were observed at the equatorward edge of red luminosity. It has been found out that the red auroral intensity increases and the band equatorward boundary shifts to lower latitudes with increasing solar wind dynamic pressure. However, a direct proportional dependence of the variations in the auroral features on the dynamic pressure variations has not been found. It has been concluded that the source of bright discrete auroras is located in the region of the low-latitude boundary layer (LLBL) on closed geomagnetic field lines. The estimated LLBL thickness is ∼3 R _e . It has been concluded that the intensity of the dayside red band depends on the solar wind plasma density, whereas the position of the position equatorward boundary depends on the dynamic pressure value and its variations.     

Soil organic carbon storage and soil CO<Subscript>2</Subscript> flux in the alpine meadow ecosystem

High-resolution sampling, measurements of organic carbon contents and ¹⁴C signatures of selected four soil profiles in the Haibei Station situated on the northeast Tibetan Plateau, and application of ¹⁴C tracing technology were conducted in an attempt to investigate the turnover times of soil organic carbon and the soil-CO₂ flux in the alpine meadow ecosystem. The results show that the organic carbon stored in the soils varies from 22.12×10⁴ kg C hm⁻² to 30.75×10⁴ kg C hm⁻² in the alpine meadow ecosystems, with an average of 26.86×10⁴ kg C hm⁻². Turnover times of organic carbon pools increase with depth from 45 a to 73 a in the surface soil horizon to hundreds of years or millennia or even longer at the deep soil horizons in the alpine meadow ecosystems. The soil-CO₂ flux ranges from 103.24 g C m⁻² a⁻¹ to 254.93 gC m⁻² a⁻¹, with an average of 191.23 g C m⁻² a⁻¹. The CO₂ efflux produced from microbial decomposition of organic matter varies from 73.3 g C m⁻² a⁻¹ to 181 g C m⁻² a⁻¹. More than 30% of total soil organic carbon resides in the active carbon pool and 72.8%281.23% of total CO₂ emitted from organic matter decomposition results from the topsoil horizon (from 0 cm to 10 cm) for the Kobresia meadow. Responding to global warming, the storage, volume of flow and fate of the soil organic carbon in the alpine meadow ecosystem of the Tibetan Plateau will be changed, which needs further research. Supported by the National Natural Science Foundation of China (Grant Nos. 40231015, 40471120 and 40473002) and the Guangdong Provincial Natural Science Foundation of China (Grant No. 06300102)     

Comparison of SO<Subscript>2</Subscript> column retrievals from BRD and DOAS algorithms

Atmospheric SO₂ has a significant impact on the urban environment and global climate. Band Residual Difference Algorithm (BRD) and Differential Optical Absorption Spectroscopy (DOAS) were used respectively by NASA and ESA science team to derive SO₂ columns from satellite observations, but there are few studies on the comparison and validation of BRD and DOAS SO₂ retrievals under the same observation conditions. In this study, the radiative transfer model SCIATRAN was firstly used to validate the accuracies of BRD and DOAS SO₂ retrievals, and analyse the uncertainty of SO₂ retrieval caused by band selection, O3 absorption, aerosol, surface reflectance, solar and viewing zenith angle. Finally, BRD and DOAS algorithms were applied to the same radiances from satellite observations, and comparisons of BRD and DOAS SO₂ retrievals were conducted over volcanic eruption and North China. Results show that, for the case with low SO₂ columns, BRD SO₂ retrievals have higher retrieval accuracy than DOAS, but typical seasonal variation with high SO₂ column in winter and low in summer can be more clearly discernible in DOAS SO₂ retrievals than BRD from satellite observations. For the case with high SO₂ columns, the differences between BRD (310.8–314.4 nm) and DOAS (315–327 nm) retrievals are large, and the value and accuracy of BRD (310.8–314.4 nm) SO₂ retrievals are lower than those of DOAS (315–327 nm) retrievals. Compared with the SO₂ inputs in forward model, both BRD (310.8–314.4 nm) and DOAS (315–327 nm) SO₂ retrievals are underestimated for the case with high SO₂ columns. The selection of wavelength range can significantly affect the accuracy of SO₂ retrieval. The error of BRD SO₂ retrieval from 310.8–314.4 nm is lower than other bands in the ultraviolet spectral region (306–327 nm). The increase of wavelength in the ultraviolet spectral region 306–330 nm can reduce the underestimation of DOAS SO₂ retrievals in the case of high SO₂ column, but slight overestimation of SO₂ retrieval is found from the 315–327 nm range in the case of low SO₂ column. The values of BRD and DOAS SO₂ retrieval decrease with atmospheric O₃ column and aerosol optical depth increasing, but increase with surface reflectance increasing. Large solar zenith angle and viewing zenith angle can introduce more errors to the BRD and DOAS SO₂ retrievals. This study is important for the improvement of retrieval algorithm and the application of SO₂ products from satellite observations.     

Calibration of Local Magnitude M<Subscript>L</Subscript> in the Azores Archipelago Based on Recent Digital Recordings

—The development of the digital seismic network in the Azores Archipelago during recent years made it possible to obtain the amplitudes (waveform) of recorded motion in a large set of stations. With this new data, maximum amplitudes of the Wood Anderson seismograph are computed, for each station/component, which, together with epicentral distances, allows for the estimation of local magnitude M _L . We used data recorded in 8 digital permanent three-component stations, with inter-stations distances up to 300 km, in the period June 1998 – June 2000, corresponding to a set of 1315 events with magnitude (M _L or M _D ) 2<M<5.8 and epicenters located in the Azores region, to estimate the coefficients of the equation to compute M _L , as well as to determine the corrections to be applied to each station. The new set of parameters, formed by attenuation coefficients and station corrections, were introduced in the calculations of the M _L , leading to smaller dispersions in the analyzed dataset. We also conclude that the attenuation in the first 150 km is similar to the California values, although higher for longer distances.     

First <Superscript>13</Superscript>C/<Superscript>12</Superscript>C isotopic characterisation of volcanic plume CO<Subscript>2</Subscript>

We describe analytical details and uncertainty evaluation of a simple technique for the measurement of the carbon isotopic composition of CO₂ in volcanic plumes. Data collected at Solfatara and Vulcano, where plumes are fed by fumaroles which are accessible for direct sampling, were first used to validate the technique. For both volcanoes, the plume-derived carbon isotopic compositions are in good agreement with the fumarolic compositions, thus providing confidence on the method, and allowing its application at volcanoes where the volcanic component is inaccessible to direct sampling. As a notable example, we applied the same method to Mount Etna where we derived a δ¹³C of volcanic CO₂ between −0.9 ± 0.27‰ and −1.41 ± 0.27‰ (Bocca Nuova and Voragine craters). The comparison of our measurements to data reported in previous work highlights a temporal trend of systematic increase of δ¹³C values of Etna CO₂ from ~ −4‰, in the 1970’s and the 1980’s, to ~ −1‰ at the present time (2009). This shift toward more positive δ¹³C values matches a concurrent change in magma composition and an increase in the eruption frequency and energy. We discuss such variations in terms of two possible processes: magma carbonate assimilation and carbon isotopic fractionation due to magma degassing along the Etna plumbing system. Finally, our results highlight potential of systematic measurements of the carbon isotopic composition of the CO₂ emitted by volcanic plumes for a better understanding of volcanic processes and for improved surveillance of volcanic activity.