Development of a technique for IR spectroscopic determination of nitrogen content and aggregation degree in diamond crystals

A technique for IR spectroscopic determination of the total nitrogen content N _S in the form of A-and B ₁-defects is suggested. It provides for the computer processing and decomposition of IR spectra into constituent bands, calculation of the total absorption band area S _N and individual areas S _A and S _B1 and their normalization with respect to the total area of the diamond intrinsic absorption S ₀, with the normalization coefficients K _S , K _A , and K _B1 being calculated. Based on the analysis of the IR spectra of 60 octahedral diamond crystals from the Mir and Yubileinaya pipes (Sakha-Yakutiya), the empirical functions N _S = 911.85 K _S ^0.9919 ppm (R ² = 0.9859), N _A = 1185.6 K _A ^1.1511 ppm (R ² = 0.8703), and N _B1 = 911.85 K _S ^0.9919 ? 1185.6 K _A ^1.1511 ppm have been defined.     

Nickelpicromerite,K<Subscript>2</Subscript>Ni(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>•6H<Subscript>2</Subscript>O,a new picromerite-group mineral from Slyudorudnik,South Urals,Russia

A new picromerite-group mineral, nickelpicromerite, K₂Ni(SO₄)₂?·?6H₂O (IMA 2012–053), was found at the Vein #169 of the Ufaley quartz deposit, near the town of Slyudorudnik, Kyshtym District, Chelyabinsk area, South Urals, Russia. It is a supergene mineral that occurs, with gypsum and goethite, in the fractures of slightly weathered actinolite-talc schist containing partially vermiculitized biotite and partially altered sulfides: pyrrhotite, pentlandite, millerite, pyrite and marcasite. Nickelpicromerite forms equant to short prismatic or tabular crystals up to 0.07 mm in size and anhedral grains up to 0.5 mm across, their clusters or crusts up to 1 mm. Nickelpicromerite is light greenish blue. Lustre is vitreous. Mohs hardness is 2–2½. Cleavage is distinct, parallel to {10–2}. D _meas is 2.20(2), D _calc is 2.22 g cm^?3. Nickelpicromerite is optically biaxial (+), α?=?1.486(2), β?=?1.489(2), γ?=?1.494(2), 2V_meas =75(10)°, 2V_calc =76°. The chemical composition (wt.%, electron-microprobe data) is: K₂O 20.93, MgO 0.38, FeO 0.07, NiO 16.76, SO₃ 37.20, H₂O (calc.) 24.66, total 100.00. The empirical formula, calculated based on 14 O, is: K_1.93Mg_0.04Ni_0.98S_2.02O_8.05(H₂O)_5.95. Nickelpicromerite is monoclinic, P2₁/c, a?=?6.1310(7), b?=?12.1863(14), c?=?9.0076(10) Å, β?=?105.045(2)°, V?=?649.9(1) ų, Z?=?2. Eight strongest reflections of the powder XRD pattern are [d,Å-I(hkl)]: 5.386–34(110); 4.312–46(002); 4.240–33(120); 4.085–100(012, 10–2); 3.685–85(031), 3.041–45(040, 112), 2.808–31(013, 20–2, 122), 2.368–34(13–3, 21–3, 033). Nickelpicromerite (single-crystal X-ray data, R?=?0.028) is isostructural to other picromerite-group minerals and synthetic Tutton’s salts. Its crystal structure consists of [Ni(H₂O)₆]²⁺ octahedra linked to (SO₄)^2? tetrahedra via hydrogen bonds. K⁺ cations are coordinated by eight anions. Nickelpicromerite is the product of alteration of primary sulfide minerals and the reaction of the acid Ni-sulfate solutions with biotite.     

Simulation of suspended sediment based on gamma test,heuristic, and regression-based techniques

In the present study, four different heuristic techniques viz. multi-layer perceptron (MLP), radial basis function (RBF), self-organizing maps (SOM), and co-active neuro-fuzzy inference system (CANFIS) with hyperbolic tangent and sigmoid transfer functions and two regression-based techniques, i.e., multiple linear regression (MLR) and sediment-rating curve (SRC), were used for suspended sediment modeling. Gamma test (GT), correlation function (CF), M test, and trail–error procedure were applied for estimation of appropriate input variables as well as training data length. The results of the GT and CF suggested the five input variables (Q_t, Q_t?1, Q_t?2, S_t?1, and S_t?2, where Q_t?1 and S_t?1 indicate the discharge and sediment values of one previous day) as the best combination. The optimal training data length (75% of total data) was estimated by M test and trail–error procedure for development of the applied models. The MLP with sigmoid transfer function (M-2) performed better than the all other models. The results of sensitivity analysis indicated that the present-day discharge (Q_t), 1-day lag discharge (Q_t?1) and 1-day lag suspended sediment (S_t?1) are the most influenced parameters in modeling current day suspended sediment (S_t).     

Analysis of magnetic activity of the rapidly rotating stars He 373 and AP 225

Spectroscopic and photometric data for the two rapidly rotating members of the α Persei cluster He 373 and AP 225 are analyzed. Improved estimates have been obtained for the projected equatorial rotation velocities: v sin i = 164 km/s for He 323 and v sin i = 129 km/s for AP 225. Multi-band photometric mapping is used to map the spot distributions on the surfaces of the two stars. The fractional spotted areas S and mean temperature difference ΔT between the unspotted photosphere and the spots are estimated (S = 7% and ΔT = 1000 K for He 373; S = 9% and ΔT = 800 K for AP 225). The H _α line profiles of both stars have variable emission components whose widths are used to deduce the presence of extended regions of emission reaching the corotation radius.     

Evaluating of the earthquake hazard parameters with Bayesian method for the different seismic source regions of the North Anatolian Fault Zone

The earthquake hazard parameters and earthquake occurrence probabilities are computed for the different regions of the North Anatolia Fault Zone (NAFZ) using Bayesian method. A homogenous earthquake catalog for M _S magnitude which is equal or larger than 4.0 is used for a time period between 1900 and 2015. Only two historical earthquakes (1766, M _S = 7. 3 and 1897, M _S = 7. 0) are included in Region 2 (Marmara Region) where a large earthquake is expected in the near future since no large earthquake has been observed for the instrumental period. In order to evaluate earthquake hazard parameters for next 5, 10, 20, 50, 100 years, M _max (maximum regional magnitude), β value, λ (seismic activity or density) are computed for the different regions of NAFZ. The computed M _max values are changed between 7.11 and 7.89. While the highest magnitude value is calculated in the Region 9 related to Tokat-Erzincan, the lowest value in the Region 10 including the eastern of Erzincan. The “quantiles” of “apparent” and “true” magnitudes of future time intervals of 5, 10, 20, 50, and 100 years are calculated for confidence limits of probability levels of 50, 70 and 90 % of the 10 different seismic source regions. The region between Tokat and Erzincan has earthquake hazard level according to the determined parameters. In this region the expected maximum earthquake size is 7.8 with 90 % occurrence probability in next 100 years. While the regional M _max value of Marmara Region is computed as 7.61, expected maximum earthquake size is 7.37 with 90 % occurrence probability in next 100 years.     

Bounding relationships between parameters of earthquakes and seismogenic clastic dikes (based on worldwide data)

This paper presents the first boundary equations describing the relationship between earthquake parameters (magnitude M _S and macroseismic intensity I _P at the observation point on the MSK-64 scale) and clastic dikes (having maximal thickness m _cd , visible height h _cd , and the index of manifestation intensity of dikes in the cross section I _cd ). As was expected, the maximal size of dikes grows with an increase in the earthquake magnitude and macroseismic intensity. Analysis of the dependences showed that it is better to use all three parameters for estimation of the minimal threshold M _S or I _P from clastic dikes, and, in the absence of data on seismogenic rupture, the maximal calculated value should be used. Some limitations in application and the advantages of the equations obtained are discussed with respect to characterizing earthquakes of the pre-instrumental period.     

Seismic hazard scenario and attenuation model of the Garhwal Himalaya using near-field synthesis from weak motion seismometry

In this paper, we present a seismic hazard scenario for the Garhwal region of the north-western Himalayan range, in terms of the horizontal Peak Ground Acceleration. The scenario earthquake of moment magnitude M _w 8.5 has a 10% exceedance probability over the next 50 years. These estimates, the first for the region, were calculated through a stepwise process based on:

• An estimation of the Maximum Credible Earthquake from the seismicity of the region and Global Seismic Hazard Assessment Program considerations, and
• four seismotectonic parameters abstracted from near field weak-motion data recorded at five stations installed in Chamoli District of the Garhwal region in the aftermath of the 1999 Chamoli earthquake. The latter include
• The frequency dependent power law for the shear wave quality factor, Q _S
• the site amplification at each station using horizontal-to-vertical-spectral ratio and generalized inversion technique
• source parameters of various events recorded by the array and application of the resulting relations between the scalar seismic moment M ₀ (dyne-cm) and moment magnitude M _w and the corner frequency, ? _c (Hz) and moment magnitude M _w to simulate spectral acceleration due to higher magnitude events corresponding to the estimated Maximum Credible Earthquake, and
• regional and site specific local spectral attenuation relations at different geometrically central frequencies in the low, moderate and high frequency bands.

     

Time constant of hydraulic-head response in aquifers subjected to sudden recharge change: application to large basins

Analytical formulae are proposed to describe the first-order temporal evolution of the head in large groundwater systems (such as those found in North Africa or eastern Australia) that are subjected to drastic modifications of their recharge conditions (such as those in Pleistocene and Holocene times). The mathematical model is based on the hydrodynamics of a mixed-aquifer system composed of a confined aquifer connected to an unconfined one with a large storage capacity. The transient behaviour of the head following a sudden change of recharge conditions is computed with Laplace transforms for linear one-dimensional and cylindrical geometries. This transient evolution closely follows an exponential trend exp(?t/τ). The time constant τ is expressed analytically as a function of the various parameters characterizing the system. In many commonly occurring situations, τ depends on only four parameters: the width a _c of the main confined aquifer, its transmissivity T _c, the integrated storage situated upstream in the unconfined aquifer M?=?S _u a _u, and a curvature parameter accounting for convergence/divergence effects. This model is applied to the natural decay of large aquifer basins of the Sahara and Australia following the end of the mid-Holocene humid period. The observed persistence of the resource is discussed on the basis of the time constant estimated with the system parameters. This comparison confirms the role of the upstream water reserve, which is modelled as an unconfined aquifer, and highlights the significant increase of the time constant in case of converging flow.     

The character of scaling earthquake source spectra for Kamchatka in the 3.5–6.5 magnitude range

The properties of the source spectra of local shallow-focus earthquakes on Kamchatka in the range of magnitudes M _w = 3.5–6.5 are studied using 460 records of S-waves obtained at the PET station. The family of average source spectra is constructed; the spectra are used to study the relationship between M _w and the key quasi-dimensionless source parameters: stress drop Δσ and apparent stress σ_a. It is found that the parameter Δσ is almost stable, while σ_a grows steadily as the magnitude M _w increases, indicating that the similarity is violated. It is known that at sufficiently large M _w the similarity hypothesis is approximately valid: both parameters Δσ and σ_a do not show any noticeable magnitude dependence. It has been established that M _w ≈ 5.7 is the threshold value of the magnitude when the change in regimes described occurs for the conditions on Kamchatka.     

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.     

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

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

Spots and the activity of M dwarfs from observations with the Kepler Space Telescope

A method for estimating the spottedness parameter S (the spotted area as a fraction of the surface of an active star) proposed earlier is applied to an analysis of activity in 1570 M dwarf stars. The analysis is based on observational material obtained with the Kepler Space Telescope, as well as data on the fluxes of the studied objects in the near and far ultraviolet (NUV and FUV) based on data from the GALEX space telescope. The variations of S with the ages of the stars are studied (four groups with different ages are distinguished), as well as variations of S with their rotational periods. A diagram characterizing the relationship between S and the Rossby number Ro resembles the classical dependence of the X-ray luminosities of active stars on Ro, and a saturation regime is attained at the same value, Ro = 0.13. Moreover, objects with ages of more than 100 million years do not form a single sequence (and stars older than 900 million years possess surface spottednesses of order 1%). The S?Ro dependence obtained could expand possibilities for analyzing the dependence of the X-ray luminosities of active stars on their Rossby numbers, and could also be applied to refine parameters characterizing the action of dynamo mechanisms, such as the dynamo number N _D . A comparison of the GALEX NUV and FUV brightness estimates with the activity parameters of the stars suggests that younger, more rapidly rotating active stars are brighter in the NUV, and that the FUV flux grows and the difference of the FUV and NUV brightnesses decreases with increasing spottedness S.     

Comparison of the magnetic properties of leading and following spots and the overlying ultraviolet emission

SDO/HMI and SDO/AIA data for the 24th solar-activity cycle are analyzed using a quicker and more accurate method for resolving π ambiguities in the transverse component of the photospheric magnetic field, yielding new results and confirming some earlier results on the magnetic properties of leading and following magnetically connected spots and single spots. The minimum inclination of the field lines to the positive normal to the solar surface α _min within umbrae is smaller in leading than in following spots in 78% of the spot pairs considered; the same trend is found for the mean angle 〈α〉 in 83% of the spot pairs. Positive correlations between the α _min values and the 〈α〉 values in leading and following spots are also found. On average, in umbrae, the mean values of 〈B〉, the umbra area S, and the angles α _min and 〈α〉 decrease with growth in the maximum magnetic field B _max in both leading and following spots. The presence of a positive correlation between B _max and S is confirmed, and a positive correlation between 〈B〉 and S in leading and following spots has been found. Themagnetic properties of the umbrae of magnetically connected pairs of spots are compared with the contrast of the He II 304 emission above the umbrae, C ₃₀₄. Spots satisfying certain conditions display a positive correlation between C _304?L and 〈α _L 〉 for the leading (L) spots, and between C _304?L /C _304?F and l _L /l _F , where l _L (l _F ) are the lengths of the field lines connecting leading (L) or following (F) spots from the corresponding spot umbrae to the apex of the field line.     

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

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

Development of Analytical Models to Estimate Pile Setup in Cohesive Soils Based on FE Numerical Analyses

This paper presents the numerical simulation of pile installation and the subsequent increase in the pile capacity over time (or setup) after installation that was performed using the finite element software Abaqus. In the first part, pile installation and the following load tests were simulated numerically using the volumetric cavity expansion concept. The anisotropic modified Cam-Clay and Dracker–Prager models were adopted in the FE model to describe the behavior of the clayey and sandy soils, respectively. The proposed FE model proposed was successfully validated through simulating two full-scale instrumented driven pile case studies. In the second part, over 100 different actual properties of individual soil layers distracted from literature were used in the finite element analysis to conduct parametric study and to evaluate the effect of different soil properties on the pile setup behavior. The setup factor A was targeted here to describe the pile setup as a function of time after the end of driving. The selected soil properties in this study to evaluate the setup factor A include: soil plasticity index (PI), undrained shear strength (S _u ), vertical coefficient of consolidation (C _v ), sensitivity ratio (S _r ), and over-consolidation ratio (OCR). The predicted setup factor showed direct proportion with the PI and S _r and inverse relation with S _u , C _v and OCR. These soil properties were selected as independent variables, and nonlinear multivariable regression analysis was performed using Gauss–Newton algorithm to develop appropriate regression models for A. Best models were selected among all based on level of errors of prediction, which were validated with additional nineteen different site information available in the literature. The results indicated that the developed model is able to predict the setup behavior for individual cohesive soil layers, especially for values of setup factor greater than 0.10, which is the most expectable case in nature.     

Mineral CuFe<Subscript>2</Subscript>S<Subscript>4</Subscript> from sulfide Copper–Nickel ores of the Lovnoozero deposit,Kola Peninsula

The unnamed mineral CuFe₂S₄ has been found from sulfide Cu–Ni ores of the Lovnoozero deposit in the Kola Peninsula, Russia. It occurs in norite composed of orthopyroxene (bronzite), Ca-rich plagioclase (66% An), pargasite, and phlogopite. The last two minerals are replaced by talc, chlorite and carbonates. Monoclinic pyrrhotite, pentlandite, chalcopyrite, and pyrite are associated ore minerals. Phase CuFe₂S₄ is enclosed predominantly in chalcopyrite, probably replacing it, and occurs in later carbonate veinlets together with redeposited sulfides. It is light yellow with a brownish tint and metallic luster. The Mohs hardness is 5–5.5; VHN 654 ± 86 kgs/mm². Density (calc.) = 4.524 g/cm³. The mineral is anisotropic, internal reflections are absent. Reflectance values (λ, nm R′ _g and R′ _p %) are: 440 30.3 29.5, 500 43.7 42.8, 560 50.9 49.6, 620 52.4 51.2, 640 52.6 51.4, 680 52.8 51.6, 700 52.7 51.4. CuFe₂S₄ is monoclinic, a = 6.260(4), b = 5.39(1), c = 13.19(1) Å, β = 94.88(7)°, V = 443(1) ų, Z = 4. The strongest reflections in the powder diffraction pattern are [d, Å (I) (hkl)]: 4.150 (10) (012), 3.559 (4) (\(11\bar 2\)), 3.020 (4) (\(10\bar 4\)), 2.560 (3) (\(21\bar 2\)), 2.500 (3) (\(10\bar 5\)), 2.340 (3) (\(12\bar 2\)), 1.817 (3) (215), 1.489 (3) (402). The chemical composition is as follows, wt %: 20.44 Cu, 35.85 Fe, 0.65 Ni, 0.14 Co, 43.15 S, total is 100.23. The empirical formula calculated on the basis of 7 atoms is Cu_0.969(Fe_1.934Ni_0.034Co_0.007)_1.975S_4.056. According to its mode of occurrence, the mineral was formed as a result of low temperature processes involving metamorphic hydrothermal solutions.     

Towards a better understanding of the evolution of the flood risk in Mediterranean urban areas: the case of Barcelona

Dam failure constitutes a grave threat to human life. However, there is still a lack of systematic and comprehensive research on the loss of life (L) caused by dam break in China. From the perspective of protecting human life, a new calculation method for L occurred in dam break floods is put forward. Fourteen dam failure cases in China are selected as the basic data by three-dimensional stratified sampling, balancing spatial, vertical elevation and temporal representations, as well as considering various conditions of the dam collapse. The method includes three progressive steps: Firstly, some impact factors of loss of life (IFL) are selected by literature survey, i.e., severity of dam break flood (S _F), population at risk (P _R), understanding of dam break (U _B), warning time (T _W) and evacuation condition (E _C). And the other IFL of weather during dam break (W _B), dam break mode (M _B), water storage (S _W), building vulnerability (V _B), dam break time (T _B) and average distance from affected area to dam (D _D) are also taken into account to get a more comprehensive consideration. According to disaster system and disaster risk, these eleven IFL are divided into four categories. Through the improved entropy method, eight key IFL are further selected out of the eleven. Secondly, four L modules are built based on four categories, which are L-causing factor module (M ₁), L-prone environment module (M ₂), affected body module (M ₃) and rescue condition module (M ₄). Eventually, by using two methods of multivariate nonlinear regression and leave-one-out cross-validation in combination with coupled four modules, the calculation method for L is established. Compared with the results of Graham method and D&M method, the result of the proposed one is much closer to the actual value and performs better in fitting effect and regional applicability. In the application, L calculation and consequence assessment are carried out in the example of Hengjiang reservoir that has already broken down. At the same time, L calculation and risk prediction are used in the analysis of Yunshan reservoir, which is under planning. The proposed method can not only be applied to estimate L and its rate (f _L ) under various types of dam break conditions in China, but also provide a reliable consequence assessment and prediction approach to reduce the risk of L.     

Fractal dimension of pore space in carbonate samples from Tushka area (Egypt)

To investigate inhomogeneous and porous structures in nature, the concept of fractal dimension was established. This paper briefly introduces the definition and measurement methods of fractal dimension. Three different methods including mercury injection capillary pressure (MICP), nuclear magnetic resonance (NMR), and nitrogen adsorption (BET) were applied to determine the fractal dimensions of the pore space of eight carbonate rock samples taken from West Tushka area, Egypt. In the case of fractal behavior, the capillary pressure P _c and cumulative fraction V _c resulting from MICP are linearly related with a slope of D-3 in a double logarithmic plot with D being the value of fractal dimension. For NMR, the cumulative intensity fraction V _c and relaxation time T ₂ show a linear relation with a slope of 3-D in a double logarithmic plot. Fractal dimension can also be determined by the specific surface area S _por derived from nitrogen adsorption measurements and the effective hydraulic radius. The fractal dimension D shows a linear relation with the logarithm of S _por . The fractal dimension is also used in models of permeability prediction. To consider a more comprehensive data set, another 34 carbonate samples taken from the same study area were integrated in the discussion on BET method and permeability prediction. Most of the 42 rock samples show a good agreement between measured permeability and predicted permeability if the mean surface fractal dimension for each facies is used.     

The impact of reflection effects on the parameters of the old pre-cataclysmic variables MS Peg and LM Com

We have determined the main parameters of the old precataclysmic variable stars MS Peg and LM Com. The radial velocities of the components, reflection effects in the spectra, and light curves of the systems are studied based on model stellar atmospheres subject to external irradiation. Forty-seven moderate-resolution spectra for MS Peg and 57 for LM Com obtained with the 6-m telescope of the Special Astrophysical Observatory are used to derive the refined orbital periods of 0.1736660 days and 0.2586873 days, respectively; the orbital eccentricities do not exceed e=0.04. The mass (M _w =0.49e_⊙) and radius (e _w =0.015R_⊙) of the MS Peg primary calculated using the gravitational redshift correspond to those for a cooling carbon white dwarf with a thin hydrogen envelope. The parameters of the red dwarf (M _r =0.19M_⊙, T_eff=3560 K, R _r =0.18R_⊙) are close to those derived from evolutionary tracks for main-sequence M stars with solar chemical composition. The radius (R _r =0.22R_⊙) and temperature (T_eff=3650 K) of the LM Com secondary exceed theoretical estimates for main-sequence stars with masses of M _r =0.17M_⊙. The luminosity excess of the red dwarf in LM Com can be explained by a prolonged (T>5×10⁶ yrs) relaxation of the M star to its normal state after the binary leaves the common-envelope stage. For both systems, theoretical U, B, V, and R light curves and spectra calculated using the adopted sets of parameters are generally consistent with the observations. This confirms the radiative origin of the hot spots, the unimportance of horizontal radiative transport, and the absence of large-scale velocity fields with high values (V_trans>50 km/s) at the surfaces of the secondaries. Most of the emission lines in the spectra of these objects are formed under conditions close to thermalization, enabling modeling of their pro files in an LTE approximation. A strong λ3905 Å emission line has been identified as the 3s²3p4s 1P⁰-3s²3p² 1S SiI λ3905.52 Å line formed in the atmosphere of the hot spot. The observed intensity can be explained by non-LTE “superionization” of SiI atoms by soft UV radiation from the white dwarf. We suggest a technique for identifying binaries whose cool components are subject to UV irradiation based on observations of λ3905 Å emission in their spectra.     

The transfer of oxygen isotopic signals from precipitation to drip water and modern calcite on the seasonal time scale in Yongxing Cave,central China

Stable isotope data of precipitation (δ¹⁸O_p and deuterium excess), drip water (δ¹⁸O_d), and modern calcite precipitates (δ¹⁸O_c and δ¹³C_c) from Yongxing Cave, central China, are presented, with monthly sampling intervals from June 2013 to September 2016. Moderate correlations between the monthly variation of δ¹⁸O_p values (from ??11.5 to ??0.7‰) and precipitation amount (r = ??0.59, n?=?34, p?<?0.01) and deuterium excess (r?=?0.39, n?=?31, p?<?0.01) imply a combined effect of changes in precipitation amount and atmospheric circulation. At five drip sites, the δ¹⁸O_d values have a much smaller variability (from ??9.1 to ??7.5‰), without seasonal signals, probably a consequence of the mixing in the karst reservoir with a deep aquifer. The mean δ¹⁸O_d value (??8.4‰) for all drip waters is significantly more negative than the mean δ¹⁸O_p value (??6.9‰) weighted by precipitation amount, but close to the wet season (May to September) mean value (??8.3‰), suggesting that a threshold of precipitation amount must be exceeded to provide recharge. Calculation based on the equilibrium fractionation factor indicates that the δ¹⁸O_c values are not in isotopic equilibrium with their corresponding drip waters, with a range of disequilibrium effects from 0.4 to 1.4‰. The δ¹⁸O_c and δ¹³C_c values generally increase progressively away from the locus of precipitation on glass plates. The disequilibrium effects in the cave are likely caused by progressive calcite precipitation and CO₂ degassing related to a high gradient of CO₂ concentration between drip waters and cave air. Our study provides an important reference to interpret δ¹⁸O_c records from the monsoon region of China.