Strmömgren zones of type II supernovae

The emission measures EM in the directions of supernova remnants and pulsars are considered as functions of their ages t. The resulting plot has a well-defined lower boundary, which can be approximated by the expression EM_min∝1/t. The quantity EM_min increases with decreasing age t and does not level off or reach a maximum until t?500 yr. It is concluded that the bulk of the radiative energy that goes into ionizing and heating the interstellar gas is released at early stages of the supernova remnant’s evolution. We suggest that most of the kinetic energy of the supernova shell is converted into thermal energy and radiated at remnant ages t<100 yr, when the supernova shell, which is expanding at an enormous speed (about 10⁴ km/s), overtakes the shell produced by the presupernova in the supergiant stage. We have estimated the ionization energy E?10⁵¹ erg, diameter L?60 pc, and electron density N_e?7 cm^?3 of the HII regions around the supernovae (the supernova Strömgren zones). A list of objects that can be reliably identified as Strömgren zones of type II supernovae is presented. The plot of pulsar pulse broadening τ as a function of the pulsar age t also has a well-defined lower boundary, for which τ∝t^?2 when t≥1000 yr. This suggests that turbulence develops during the first thousand years after the supernova outburst. It is also concluded that turbulence plays an important role in the formation and evolution of the Strömgren zones of type II supernovae.     

Application of cubic spline in soil erosion modeling from Narmada Watersheds,India

Soil erosion by water is ubiquitous, exhibits spatio-temporal variability, and is fundamental to determining sediment yield which is key to proper watershed management. In this study, we propose a relationship between the curve number and sediment yield index (SYI) using cubic splines. Using field data from four watersheds, the relation between observed and computed SYI is found to have a coefficient of determination (R²) value from 0.63 to 0.88 suggesting that such a relation can be used to determine SYI from the available CN value. It is found that cubic splines perform satisfactorily with Nash-Sutcliff efficiency ranging from 60.18 to 64.01%, absolute prediction error from 1.35 to 5.56%, integral square error from 1.21 to 5.82%, coefficient of correlation from 79.32 to 93.78%, and degree of agreement from 0.87 to 0.99%.     

Prograde garnet growth along complex <Emphasis Type="Italic">P–T–t</Emphasis> paths: results from numerical experiments on polyphase garnet from the Wölz Complex (Austroalpine basement)

Garnet in metapelites from the Wölz Complex of the Austroalpine crystalline basement east of the Tauern Window characteristically consists of two growth phases, which preserve a comprehensive record of the geothermal history during polymetamorphism. From numerical modelling of garnet formation, detailed information on the pressure–temperature–time (P–T–t) evolution during prograde metamorphism is obtained. In that respect, the combined influences of chemical fractionation associated with garnet growth, modification of the original growth zoning through intragranular diffusion and the nucleation history on the chemical zoning of garnet as P and T change during growth are considered. The concentric chemical zoning observed in garnet and the homogenous rock matrix, which is devoid of chemical segregation, render the simulation of garnet growth through successive equilibrium states reliable. Whereas the first growth phase of garnet was formed at isobaric conditions of ～3.8 kbar at low heating/cooling rates, the second growth phase grew along a Barrovian P–T path marked with a thermal peak of ～625°C at ～10 kbar and a maximum in P of ～10.4 kbar at ～610°C. For the heating rate during the growth of the second phase of garnet, average rates faster than 50°C Ma^?1 are obtained. From geochronological investigations the first growth phase of garnet from the Wölz Complex pertains to the Permian metamorphic event. The second growth phase grew in the course of Eo-Alpine metamorphism during the Cretaceous.     

Evaluation of Curve Fitting Techniques in Deriving p–y Curves for Laterally Loaded Piles

The p–y method is one of the most popular methods for the analysis and design of laterally loaded piles. The mathematical relationship it provides between the bending moment, which can be easily measured at strain gauges along the pile, and the soil resistance and lateral pile displacement, facilitates the construction of p–y curves. Numerical techniques are required to fit smooth continuous curves to the discrete bending moment data in order to improve the accuracy of subsequent differentiation and integration operations. Due to the lack of guidance on the optimum positioning of strain gauges and the reliability and accuracy of curve fitting methods, a unifying study, inclusive of small (0.61 m) and large (3.8 and 7.5 m) diameter piles in clay, was carried out using 18 strain gauge layouts and cubic spline, cubic to quintic B-spline and 3rd to 10th degree global polynomial techniques. Bending moment data was obtained using 3D finite element analysis. Through a comprehensive evaluation, the cubic and cubic B-spline methods were found to be consistently accurate in deriving p–y curves for both the small and large diameter piles.     

Coda <Emphasis Type="Italic">Q</Emphasis> in Qeshm Island,south of Iran,using aftershocks of the Qeshm earthquake of November 27, 2005

Coda wave attenuation is estimated for Qeshm Island which is located in the southeastern part of Zagros. For this purpose, the aftershocks of Qeshm earthquake in November 27, 2005, recorded within an epicentral distance less than 100 km, have been used. More than 829 earthquakes were recorded by a local temporary network consisting of 16 short period stations installed after a week after the main shock for ～10 weeks. The coda quality factor, Q _c, was estimated using the single-backscattering model in frequency bands of 0.5–24 Hz. In this research, lateral and vertical variations of coda Q in Qeshm Island are explored. In Qeshm Island, absence of significant lateral variation of coda Q is observed. To investigate the attenuation variation with depth, the coda Q value was calculated for coda time windows with different lengths (5, 10, 15, 20, 25, and 30 s). It is observed that coda Q increases with depth. However, in our study area, the rate of increase of coda Q with depth is not uniform. Beneath Qeshm Island, the rate of increase of coda Q is greater at depths less than ～40 km compared with those of larger depths. This is indicating the existence of a low attenuation anomalous structure under the ～40-km depth which may be correlated with the Moho depth in this region. The average frequency relation for this region is Q _c = 36 ± 1.2f ^{0.94 ± 0.039} at a 5 s-lapse time window length and Q _c = 110 ± 1.8f ^{0.88 ± 0.09} at a 30-s lapse time window length.     

Detailed comparison of nine intact rock failure criteria using polyaxial intact coal strength data obtained through PFC<Superscript>3D</Superscript> simulations

Study of intact rock failure criteria is an important topic in rock mechanics. In this study, applicability of nine different intact rock failure criteria is investigated for intact coal strength data. PFC^3D modeling was used to simulate the laboratory polyaxial tests for cubic intact coal blocks of side dimension 110 mm under different confining stress combinations. A modified grid search procedure is proposed and used to find the best-fitting parameter values and to calculate the coefficient of determination (R ²) values for each criterion. Detailed comparisons of the nine criteria are made using the following aspects: R ² values, σ ₁ ? σ ₂ plots for different σ ₃, shapes on the deviatoric plane, linearity or nonlinearity on the meridian planes. Through the comparisons of R ² values, σ ₁ ? σ ₂ plots and meridian lines, the modified Wiebols–Cook and modified Lade criteria were found to fit the intact coal strength data best. The nine failure criteria are categorized into three types based on the appearances on the deviatoric plane.     

Evaluation of ferrous and ferric Mössbauer fractions. Part II

Mössbauer fractions f are reported for various ferrous- and/or ferric-containing oxides, hydroxides, silicates, and phosphates to extend the list previously reported by De Grave and Van Alboom (1991). The f fractions were evaluated from the experimental temperature dependencies of their center shifts, assuming the Debye model for the lattice vibrations. For most Fe²⁺ sites the characteristic Mössbauer or lattice temperatures Θ_M are in the range 300–400 K, while those for Fe³⁺ sites are close to or exceed 500 K, implying significantly higher f fractions for Fe³⁺ than for Fe²⁺, in particular at room temperature. A correlation between Θ_M and the coordination type, or, for a given valence state and coordination type, between Θ_M and the mineral type is, however, not obvious.     

Experimental Investigation into Hydraulic Fracture Network Propagation in Gas Shales Using CT Scanning Technology

Multistage fracturing of the horizontal well is recognized as the main stimulation technology for shale gas development. The hydraulic fracture geometry and stimulated reservoir volume (SRV) is interpreted by using the microseismic mapping technology. In this paper, we used a computerized tomography (CT) scanning technique to reveal the fracture geometry created in natural bedding-developed shale (cubic block of 30 cm × 30 cm × 30 cm) by laboratory fracturing. Experimental results show that partially opened bedding planes are helpful in increasing fracture complexity in shale. However, they tend to dominate fracture patterns for vertical stress difference Δσ _v  ≤ 6 MPa, which decreases the vertical fracture number, resulting in the minimum SRV. A uniformly distributed complex fracture network requires the induced hydraulic fractures that can connect the pre-existing fractures as well as pulverize the continuum rock mass. In typical shale with a narrow (<0.05 mm) and closed natural fracture system, it is likely to create complex fracture for horizontal stress difference Δσ _h  ≤ 6 MPa and simple transverse fracture for Δσ _h  ≥ 9 MPa. However, high naturally fractured shale with a wide open natural fracture system (>0.1 mm) does not agree with the rule that low Δσ _h is favorable for uniformly creating a complex fracture network in zone. In such case, a moderate Δσ _h from 3 to 6 MPa is favorable for both the growth of new hydraulic fractures and the activation of a natural fracture system. Shale bedding, natural fracture, and geostress are objective formation conditions that we cannot change; we can only maximize the fracture complexity by controlling the engineering design for fluid viscosity, flow rate, and well completion type. Variable flow rate fracturing with low-viscosity slickwater fluid of 2.5 mPa s was proved to be an effective treatment to improve the connectivity of induced hydraulic fracture with pre-existing fractures. Moreover, the simultaneous fracturing can effectively reduce the stress difference and increase the fracture number, making it possible to generate a large-scale complex fracture network, even for high Δσ _h from 6 MPa to 12 MPa.     

The transformation features of impurity defects in natural diamonds of various habits under high <Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">T</Emphasis> conditions

The results of the investigations of the transformation of impurity defects in natural diamonds of various habits at the stage of high-temperature annealing at P = 6 GPa and T = 2200°C are presented. The studies conducted allowed us to ascertain that the transformations of Aand B-defects in diamonds of octahedral and cubic habits follow general regularities. This fact shows that most of the diamonds of cubic habit with low degree of aggregation of nitrogen centers were not really annealed over a long-term interval. Unlike octahedral diamonds, those of cubic habit are characterized by a pronounced increase in the peak of H-containing defects (3107 cm^–1) after annealing.     

The stress regime in a Rotliegend reservoir of the Northeast German Basin

In-situ stresses have significant impact, either positive or negative, on the short and long term behaviour of fractured reservoirs. The knowledge of the stress conditions are therefore important for planning and utilization of man-made geothermal reservoirs. The geothermal field Groß Schönebeck (40 km north of Berlin/Germany) belongs to the key sites in the northeastern German Basin. We present a stress state determination for this Lower Permian (Rotliegend) reservoir by an integrated approach of 3D structural modelling, 3D fault mapping, stress ratio definition based on frictional constraints, and slip-tendency analysis. The results indicate stress ratios of the minimum horizontal stress S _hmin being equal or increasing 0.55 times the amount of the vertical stress S _V (S _hmin ≥ 0.55S _V ) and of the maximum horizontal stress S _Hmax ≤ 0.78–1.00S _V in stress regimes from normal to strike slip faulting. Thus, acting stresses in the 4,100-m deep reservoir are S _V  = 100 MPa, S _hmin = 55 MPa and S _Hmax = 78?100 MPa. Values from hydraulic fracturing support these results. Various fault sets of the reservoir are characterized in terms of their potential to conduct geothermal fluids based on their slip and dilatation tendency. This combined approach can be adopted to any other geothermal site investigation.     

Spectroscopic monitoring of SS 433: A search for long-term variations of kinematic model parameters

Between 1994 and 2006, we obtained uniform spectroscopic observations of SS 433 in the region of Hα. We determined Doppler shifts of the moving emission lines, Hα ₊ and Hα _?, and studied various irregularities in the profiles for the moving emission lines. The total number of Doppler shifts measured in these 13 years is 488 for Hα _? and 389 for Hα ₊. We have also used published data to study possible long-term variations of the SS 433 system, based on 755 Doppler shifts for Hα _? and 630 for Hα ₊ obtained over 28 years. We have derived improved kinematic model parameters for the precessing relativistic jets of S S 433 using five-and eight-parameter models. On average, the precession period was stable during the 28 years of observations (60 precession cycles), at 162.250^d ± 0.003^d. Phase jumps of the precession period and random variations of its length with amplitudes of ≈6% and ≈1%, respectively, were observed, but no secular changes in the precession period were detected. The nutation period, P _nut = 6.2876^d ± 0.00035^d, and its phase were stable during 28 years (more than 1600 nutation cycles). We find no secular variations of the nutation cycle. The ejection speed of the relativistic jets, v, was, on average, constant during the 28 years, β = v/c = 0.2561 ± 0.0157. No secular variation of β is detected. In general, S S 433 demonstrates remarkably stable long-term characteristics of its precession and nutation, as well as of the central “engine” near the relativistic object that collimates the plasma in the jets and accelerates it to v = 0.2561c. Our results support a model with a “slaved” accretion disk in S S 433, which follows the precession of the optical star’s rotation axis.     

Physical and optical characteristics of atmospheric aerosols during ICARB at Manora Peak,Nainital: A sparsely inhabited,high-altitude location in the Himalayas

Collocated measurements of the optical and physical properties of columnar and near-surface aerosols were carried out from Manora Peak, Nainital (a sparsely inhabited, high altitude location, ～2 km above mean sea level, in the Himalayas), during the Integrated Campaign for Aerosols, gases and Radiation Budget (ICARB) under the Geosphere Biosphere Programme of the Indian Space Research Organization (ISRO-GBP). Under this, observational data of spectral aerosol optical depths (AOD), mass concentration of aerosol black carbon (M _B ), mass concentration (M _T ) and number concentration (N _t ) of composite (total) aerosols near the surface and meteorological parameters were collected during the period February 15 to April 30, 2006. Though very low (<0.1 at 500 nm) AODs were observed during clear days, as much as a four-fold increase was seen on hazy days. The Ångström exponent (α), deduced from the spectral AODs, revealed high values during clear days, while on hazy days α was low; with an overall mean value of 0.69 ± 0.06 for the campaign period. BC mass concentration varied between 0.36 and 2.87 μg m^?3 and contributed in the range 0.7 to 1.8% to the total aerosol mass. Total aerosol number concentration and BC mass concentration showed diurnal variation with a midnight and early morning minimum and a late afternoon maximum; a pattern quite opposite to that seen in low altitude stations. These are attributed to the dynamics of the atmospheric boundary layer.     

Influence of oxygen fugacity and temperature on the redox state of iron in natural silicic aluminosilicate melts

The influence of oxygen fugacity (fO₂) and temperature on the valence and structural state of iron was experimentally studied in glasses quenched from natural aluminosilicate melts of granite and pantellerite compositions exposed to various T-fO₂ conditions (1100–1420°C and 10^?12–10^?0.68 bar) at a total pressure of 1 atm. The quenched glasses were investigated by Mössbauer spectroscopy. It was shown that the effect of oxygen fugacity on the redox state of iron at 1320–1420°C can be described by the equation log(Fe³⁺/Fe²⁺) = k log(fO₂) + q, where k and q are constants depending on melt composition and temperature. The Fe³⁺/Fe²⁺ ratio decreases with decreasing fO₂ (T = const) and increasing temperature (fO₂ = const). The structural state of Fe³⁺ depends on the degree of iron oxidation. With increasing Fe³⁺/Fe²⁺ ≥ 1, the dominant coordination of Fe³⁺ changes from octahedral to tetrahedral. Ferrous iron ions occur in octahedral (and/or five-coordinated) sites independent of Fe³⁺/Fe²⁺.     

H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].     

The space distribution of quasars

The space distribution of quasars from the 2dF and SDSS DR5 catalogs in the redshift interval 0.3 < z < 1.9 is analyzed. The distributions of quasars in both catalogs are found to have the following common features: (1) when the distance between the nearest objects exceeds 35h ^?1 Mpc (where h = H ₀/100 km/s Mpc is the dimensionless Hubble constant), the distribution of quasars virtually coincides with a uniform three-dimensional distribution; (2) on scale lengths of (5–35)h ^?1 Mpc, the fractal dimension of the quasar distribution is 2.3; (3) the amplitude of quasar clustering and the average distance between neighboring quasars increase slowly with z (at a significance level of about 1.5σ). Twenty large groups of quasars with sizes of (50–150)h ^?1 Mpc can be identified in the 2dF catalog at the 4σ significance level. These groups are incipient superclusters (two earlier known groups are confirmed). The space density of these groups is of the order of 7h ³ Gpc^?3.     

A many-year photometric and spectroscopic variability study of planetary nebulae (1968–2008). The young planetary nebula IC 4997

We present the results of our many-year (1968–2008) photoelectric photometry and absolute spectrophotometry of the young variable planetary nebula IC 4997, performed in uniform systems. Integrated line energy fluxes in the range 3727–5007 Å are tabulated, along with the integrated (star + nebula) continuum flux at 4500 Å, and time variations investigated. We study the time behavior of fluxes in the hydrogen, HeI, [OIII], [OII], [NeIII] emission lines. Variations of the line intensity ratio R = F([OIII]4363 Å)/F(Hγ) between 1938 and 2008 are presented. We estimate the nebula’s electron density N _e and electron temperature T _e for 1972–1992: the mean nebula N _e increased by the factor of five, from ～4 × 10⁵ to ～2 × 10⁶ cm^?3, while the mean T _e increased from 12 000 to 14 000 K. The color temperature of the exciting star, measured from the 4000–4900 Å continuum, increased from 37 000–40 000 K to 47 000 K during the time interval covered by our observations, as is confirmed by the growth of ionization of HeI, NeIII, and OIII (from the 4363 Å line). V-band variations of the object’s integrated (nebula + star) light with an amplitude about 0.5 ^m can essentially be fully explained as being due to changing contributions from variable fluxes in the [OIII] 5007 and 4959 Å lines. The V magnitude in 2009 (after 40 years) happened to be the same as at the start of our observations in 1968. A period of the order of 50 years can also be noted in log R(t). This may provide evidence for binarity of the central star and be related to its orbital period.     

Photometric elements and apsidal motion of the eclipsing binary ES Lac

We have obtained high-accuracy photoelectric measurements of ES Lac, an eclipsing binary with an elliptical orbit (B9III + B9III; P = 4.459^d, e = 0.198) in 1985–2004 at the Sternberg Astronomical Institute’s Tien Shan High-Altitude Observatory. Our detailed analysis of the 19-year uniform series of measurements has yielded the first photometric elements for this system, as well as a self-consistent set of physical and geometrical parameters for the binary. The virtually identical components (M ₁ = M ₂ = 3.0 M _⊙; R ₁ = R ₂ = 4.12 R _⊙) are appreciably separated from the main sequence, and are located on the giant branch: their age is t = (3.5 ± 0.2) × 10⁸ yrs. An analysis of our observations together with previously published times of minima has enabled a considerable refinement of the period of the apsidal motion, U = 355 ± 20 years, and a first determination of the apsidal parameter reflecting the radial density distributions for the components stars: k ₂ ^obs = 0.00213(18). This value is in a good agreement with the value expected theoretically for current evolutionary models of such stars: k ₂ ^th = 0.00257(15).     

Factors of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O fractionation in corundum estimated from the calculation of isotopic shifts on vibration frequencies

The β-factors of corundum were estimated on the basis of DFT calculations of vibrational frequency changes due to ¹⁶O–¹⁸O isotope substitution in a harmonic approximation using an all-electron Gaussian-type basis set and the B3LYP hybrid functional (the CRYSTAL program). Calculations were performed accounting for eight phonon wave vectors within the first Brillouin zone. The results are approximated by the relation 1000ln β _crn = 9.19874x–0.12326x ² + 0.00213x ³ (x = 10⁶/T(K)², 400 < T(K) < 1300), which can be used in isotope geochemical studies in combination with the known temperature effects on the β-factors of other phases.     

Homogenization of carbonate-bearing microinclusions in diamond at <Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">T</Emphasis> parameters of the upper mantle

The staged high-pressure annealing of natural cubic diamonds with numerous melt microinclusions from the Internatsional’naya kimberlite pipe was studied experimentally. The results mainly show that the carbonate phases, the daughter phases in partially crystallized microinclusions in diamonds, may undergo phase transformations under the mantle P–T conditions. Most likely, partial melting and further dissolution of dolomite in the carbonate–silicate melt (homogenization of inclusions) occur in inclusions. The experimental data on the staged high-pressure annealing of diamonds with melt microinclusions allow us to estimate the temperature of their homogenization as 1400–1500°C. Thus, cubic diamonds from the Internatsional’naya pipe could have been formed under quite high temperatures corresponding to the lithosphere/asthenosphere boundary. However, it should be noted that the effect of selective capture of inclusions with partial loss of volatiles in relation to the composition of the crystallization medium is not excluded during the growth. This may increase the temperature of their homogenization significantly between 1400 and 1500°C.