Properties of the aftershocks sequences of the 2000 and 2002 earthquakes in Aksehir-Afyon graben,west-central Anatolia,Turkey

The b-value of the Gutenberg–Richter’s frequency–magnitude relation and the p-value of the modified Omori law, which describes the decay rate of aftershock activity, were investigated for more than 500 aftershocks in the Aksehir-Afyon graben (AAG) following the 15 December 2000 Sultandagi–Aksehir and the 3 February 2002 Çay–Eber and Çobanlar earthquakes. We used the Kandilli Observatory’s catalog, which contains records of aftershocks with magnitudes ≥2.5. For the Çobanlar earthquake, the estimated b-values for three aftershock sequences are in the range 0.34 ≤  b ≤ 2.85, with the exception of the one that occurred during the first hour (4.77), while the obtained p-values are in the range 0.44 ≤ p ≤ 1.77. The aftershocks of the Sultandagi earthquake have a high p-value, indicating fast decay of the aftershock activity. A regular increase of b can be observed, with b < 1.0 after 0.208 days for the Çay–Eber earthquake. A systematic and similar increase and decrease pattern exists for the b- and p-values of the Çobanlar earthquakes during the first 5 days.     

Galactic parameters derived from open-cluster data

The components U₀ and V₀ of the solar motion and the Oort constant A₀ are determined using the data of a homogeneous open-cluster catalog with corrected distance moduli. The results are based on a sample of 146 open clusters with known radial velocities located in the Galactic plane (b<7°) within 4 kpc of the Sun. The solar Galactocentric distance R₀ is determined using two kinematic methods. The following results are obtained: A₀=17.0±0.9 km/s kpc, U₀=10.5±1.0 km/s, V₀=11.5±1.1 km/s, R₀=8.3±0.3 pc.     

Surface irrigation performance of date palms under water scarcity in arid irrigated lands

In this paper, a study on the performance of surface irrigation of date palms in a Tunisian arid area (Douz oasis) is presented. The study is conducted in 16 plots with various sizes and soil textures over a 4-year period (2012–2015). In the first step, an assessment of total water requirements of the date palms is carried out. Then, the surface irrigation performance is analyzed using three indicators, i.e., the relative water supply (RWS) indicator, the uniformity index of water distribution (D _U ), and the water application efficiency (E _a ). Finally, the irrigation management problems are identified. The results indicate that in the arid Tunisian Saharan oases, the soil texture, plot size, and farmers’ practices (especially irrigation duration) have significant effects on surface irrigation performance. The average annual net irrigation requirements of date palms are about 2400 mm. The RWS increases from 1.8 in the smaller plots (0.5 ha) to 3.6 in the largest plots (2.5 ha), implying that the increase in the plot size requires an excessive water supply. D _U decreases from 80.7 in the 0.5 ha plots to 65.4 in the 2.5 ha plots; however, no significant difference in the E _a is observed. The results show that the soil texture has no influence on the RWS and D _U , but the E _a is significantly higher in the loamy-sand soils (46.7%) compared to the sandy soils (36.3%). Overall, RWS indicator is higher than 1 (RWS?=?2.6) implying excessive irrigation supply to the system. Although D _U is relatively uniform (>?60%), E _a is relatively low (<?50%) indicating that the current irrigation management is inefficient. These findings have a paramount importance for improving irrigation water management in the Tunisian Saharan oases.     

The dust envelope of R CrB in June 2001: Detection of a compact IR source

We present the results of a joint analysis of JHK interferometric and UBVJHKLM photometric observations of RCrB acquired in June 2001. The baseline for the IOTA interferometer was 21.18 m. During the observations, the star was in its bright state in the V band and near its maximum brightness in the L band. Our analysis reveals an IR source that is considerably smaller than the extended dust envelope discovered earlier. We identify this compact IR source with the emission from a group of dust clouds. The linear scale (diameter) of the IR source was d _in,c ≈ 13.5D_* (its angular diameter is θ_in,c≈6.4 mas). About 7% of the star’s radiation was obscured by this group of clouds, which contributed ～14% of the total IR excess of R CrB and ～22% of the K-band flux. The color temperature of the compact source was only ～300 K higher than the color temperature of the extended dust envelope. The inner boundary of the extended dust envelope had a diameter of d _in,e ≈ 90D_* (θ_in,e≈43 mas).     

Mineralogy,Mössbauer spectra and electrical conductivity of triphylite Li(Fe<Superscript>2+</Superscript>,Mn<Superscript>2+</Superscript>) PO<Subscript>4</Subscript>

The electrical conductivity of monocrystalline triphylite, Li(Fe²⁺,Mn²⁺)PO₄, with the orthorhombic olivine-type structure was measured parallel (∥) to the [010] direction and ∥ [001] (space group Pnma), between ～400 and ～700 K. Electrical measurements on triphylite are of technological interest because LiFePO₄ is a promising electrode material for rechargeable Li batteries. Triphylite was examined by electron microprobe, ICP atomic emission spectroscopy, X-ray diffraction, Mössbauer spectroscopy and microscopic analysis. The DC conductivity σ_DC was determined from AC impedance data (20 Hz–1 MHz) extrapolating to zero frequency. Triphylite shows σ_DC with activated behavior measured ∥ [010] between ～500 and ～700 K during the first heating up, with activation energy of E _A = 1.52 eV; on cooling E _A = 0.61 eV was found down to ～400 K and extrapolated σ_DC (295 K) ～10^?9 Ω^?1cm^?1; ∥ [001] E _A = 0.65 eV and extrapolated σ_DC(295 K) ～10^?9 to 10^?10 Ω^?1cm^?1, measured during the second heating cycle. The enhanced AC conductivity relative to σ_DC at lower temperatures indicates a hopping-type charge transport between localized levels. Conduction during the first heating up is ascribed to ionic Li⁺ hopping. DC polarization experiments showed conduction after the first heating up to be electronic related to lowered activation energy. Electronic conduction appears to be coupled with the presence of Li⁺ vacancies and Fe³⁺, formed by triphylite alteration. For comparison, σ_DC was measured on the synthetic compound LiMgPO₄ with olivine-type structure, where also an activated behavior of σ_DC with E _A ～1.45 eV was observed during heating and cooling due to ionic Li⁺ conduction; here no oxidation can occur associated with formation of trivalent cations.     

Shear heating by translational brittle reverse faulting along a single,sharp and straight fault plane

Shear heating by reverse faulting on a sharp straight fault plane is modelled. Increase in temperature (T _i ) of faulted hangingwall and footwall blocks by frictional/shear heating for planar rough reverse faults is proportional to the coefficient of friction (μ), density and thickness of the hangingwall block (ρ). T _i increases as movement progresses with time. Thermal conductivity (K _i ) and thermal diffusivity (\(k_{\mathrm {i}}^{\prime }\)) of faulted blocks govern T _i but they do not bear simple relation. T _i is significant only near the fault plane. If the lithology is dry and faulting brings adjacent hangingwall and footwall blocks of the same lithology in contact, those blocks undergo the same rate of increase in shear heating per unit area per unit time.     

Infrared photometry of the unique object FG Sge in 1985–2001

Our analysis of many years of infrared photometry of the unique object FG Sge indicates that the dust envelope formed around the supergiant in August 1992 is spherically symmetrical and contains compact, dense dust clouds. The emission from the spherically symmetrical dust envelope is consistent with the observed radiation from the star at 3.5–5 µm, and the presence of the dust clouds can explain the radiation observed at 1.25–2.2 µm. The mean integrated flux from the dust envelope in 1992–2001 was ～(1.0±0.2)×10^?8 erg s^?1cm^?2. The variations of its optical depth in 1992–2001 were within 0.5–1.0. The maximum density of the dust envelope was recorded in the second half of 1993 and corresponded to mean optical depths as high as unity. Several times in the interval from 1992 to 2001, the dusty material of the envelope partially dissipated and was then replenished. For example, the optical depth of the dust cloud at λ=1.25 µm during the last brigthness minimum in the J band was τ_1.25≈4.3, which is much higher than the optical depth of the dust envelope of FG Sge. During maxima of the J brightness, the mean spectral energy distribution at 0.36–5 µm can be represented as a combination of radiation from a G0 supergiant that is attenuated by a dust envelope with a mean optical depth of 0.65±0.15 and emission from the spherically symmetrical dust envelope itself, with the temperature of the graphite grains being 750±150 K. At minima of the J brightness, only radiation from the dust envelope is observed at 1.65–5 µm, with the radiation from the supergiant barely detectable at 1.25 µm. As a result, the integrated flux during J minima is almost half that during J maxima. The mean mass of the spherically symmetrical dust envelope of FG Sge in 1992–2001 was (3 ± 1) × 10^?7M_⊙. This envelope’s mass varied by nearly a factor of two during 1992–2001, in the range (2 – 4) × 10^?7M_⊙. In Autumn 1992, the mass-loss rate from the supergiant exceeded 2 × 10^?7M_⊙/yr. The average rate at which matter was injected into the envelope during 1993–2001 was 10^?8M_⊙/yr. The mean rate of dissipation of the dust envelope was about 1 × 10^?8M_⊙/yr. During 1992–2001, the supergiant lost about 8.7 × 10^?7M_⊙. The parameters of the dust envelope were relatively constant from 1999 until the middle of 2001.     

A numerical method for the solution of the nonlinear turbulent one-dimensional free surface flow equations

Free surface flow of an incompressible fluid over a shallow plane/undulating horizontal bed is characteristically turbulent due to disturbances generated by the bed resistance and other causes. The governing equations of such flows in one dimension, for finite amplitude of surface elevation over the bed, are the Continuity Equation and a highly nonlinear Momentum Equation of order three. The method developed in this paper introduces the “discharge” variable q = η U, where η = elevation of the free surface above the bed level, and U = average stream-wise forward velocity. By this substitution, the continuity equation becomes a linear first-order PDE and the momentum equation is transformed after introduction of a small approximation in the fifth term. Next, it is shown by an invertibility argument that q can be a function of η: q = F(η), rendering the momentum equation as a first order, second degree ODE for F(η), that can be be integrated by the Runge-Kutta method. The continuity equation then takes the form of a first order evolutionary PDE that can be integrated by a Lax-Wendroff type of scheme for the temporal evolution of the surface elevation η. The method is implemented for two particular cases: when the initial elevation is triangular with vertical angle of 120 ^° and when it has a sinusoidal form. The computations exhibit the physically interesting feature that the frontal portion of the propagating wave undergoes a sharp jump followed by tumbling over as a breaker. Compared to other discretization methods, the application of the Runge-Kutta and an extended version of the Lax-Wendroff scheme is much easier.     

New models to predict rock failure by linking the volume dilatant point with the peak stress point

The rock mass failure process can be divided into several distinct deformation stages: the compaction stage, elastic stage, stable failure stage, accelerated failure stage, and post-peak stage. Although each stage has been well studied, the relationship among the stages has not been established. Here, we establish two models which are the Strain model Q and Energy density model S by using the renormalization group theory and investigate the mechanical relationship between the volume dilatant point and peak stress point on the rock stress-strain curve. Our models show that the strain ratio (ε _f /ε _c ) and energy ratio (E _f /E _c ) at the volume dilatant point and peak stress point are solely functions of the shape parameter m. To verify our models, we further studied the failure process of rock specimens through several uniaxial compression experiments and found that the relationship between ε _f /ε _c or E _f /E _c and m shares a notably similar pattern to that from our theoretical model. However, the ε _f /ε _c and E _f /E _c values in our experiments are slightly smaller than those predicted by the models. In brief, we demonstrate that our models can be used to predict the failure process of the laboratory-scale hard brittle rock samples.     

IR photometry and models for the dust envelopes of two carbon stars

The results of JHKLM photometry of two carbon stars are presented: the irregular variable NQ Cas and the Mira star BD Vul. Data on the mean fluxes supplemented with mid-IR observations with the IRAS, AKARI, andWISE satellites are used to compute spherically symmetrical model dust envelopes for the stars, consisting of particles of amorphous carbon and silicon carbide. The optical depth in the visible for the comparatively cool dust envelope of BD Vul, with a dust temperature at its inner boundary T₁ = 610 K, is fairly low: τ_V = 0.13. The dust envelope of NQ Cas is appreciably hotter (T₁ = 1550 K), and has τ_V = 0.32. The estimated mass-loss rates are 1.5 × 10^?7M_⊙/yr for NQ Cas and 5.9 × 10^?7M_⊙/yr for BD Vul.     

Mixing of metals during stripping of galactic gaseous halos

The mixing of metals in the intergalactic gas when a galaxy with a metal-rich envelope moves through the intergalactic medium is analyzed. Two simple models for the initial distribution of metals are considered. In the first case, the metals are concentrated in a fairly thin envelope with thickness ΔR _s =1 kpc, outer radius R _s =31 kpc, and metallicity Z=10^?3. In the second case, material with the same metallicity uniformly fills an entire spherical region of radius R _s . After 2.85 Gyr, the metals are distributed over a fairly extended volume with a typical size of ?200 kpc in the direction of the motion of the intergalactic gas, with a mean metallicity of ?4.6×10^?4 in metal-enriched regions. However, the distribution of metals remains extremely nonuniform, so that the main contribution to the overall metallicity is provided by metal-rich islands Z?6×10^?4 that occupy only ～10% of the total mixing volume. Moreover, metal-free regions remain in this volume.     

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.     

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.     

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.     

IR photometry of the symbiotic star BF Cyg: Detection of the red giant’s ellipsoidal brightness variability

We present the results of our IR photometric observations of the classical symbiotic star BF Cyg acquired in 1978–2003. The variability range in the J and K bands was ～0.2^m. A periodic component in the cool star’s brightness variations is clearly visible, its period being half the orbital one and its J amplitude being ～0.15^m. This component is associated with the ellipsoidal shape of the red giant, which model calculations show fills its Roche lobe. This is required in order to reproduce ellipsoidal brightness variability with such a large amplitude: the calculated amplitude for a red giant filling 90% of its Roche lobe is half the observed value. At the same time, it was not possible to confidently chose the optimum component-mass ratio, q = M _giant /M_hot, and orbital inclination, i, from possible values in the ranges q = 2–4, i = 70°–90°. Including the contribution from the hot radiation sources (the hot component and ionized envelope), which vary with a period equal to the orbital period, has a considerable influence on the estimated parameters associated with the red giant’s ellipsoidal brightness variations, and this contribution cannot be neglected. The deviations of the observed from the calculated light curve are irregular, with the rms deviation being σ(O-C) ≈ 0.04^m.     

Paleoweathering and Paleoclimatic Reconstructions using Neoproterozoic Paleosol in the Eastern Margin of the Pranhita-Godavari Basin,India

Paleogene sediments of the inner fold belt, Naga hills, have very well preserved ichnofossils. 16 ichnospecies have been documented among 13 ichnogenera such as Arenicolites isp., Chondrites targionii, Cylindrichnus isp., Diplocraterion parallelum, Gyrochorte isp., Ophiomorpha annulata, O. irregulaire, O. nodosa, O. rudis, Palaeophycus tubularis, Planolites beverleyensis, Scolicia palaeobullia, Skolithos linearis, Trypinites weisei, Thalassinoides horizontalis and Zoophycos isp. The ichnofossil assemblages comprise mostly domichnia and fodinichnia benthos of the Skolithos and Cruziana ichnofacies. A shallow marine nearshore to offshore marine environment with fluctuating energy condition has been envisaged.     

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

Bearing strength surface for bridge caisson foundations in frictional soil under combined loading

The problem of estimating the bearing capacity of massive caisson foundations in frictional soil under combined vertical (N), horizontal (Q) and moment (M) loading is examined numerically by means of three-dimensional finite element analyses. The analysis is performed with due consideration to the foundation’s depth-to-width ratio (D/B), the magnitude of the vertical load and the caisson-soil contact interface conditions. The constitutive law for soil behavior is appropriately validated against experimental results from 1-g small-scale tests, available in the literature. The ultimate limit states are presented in the form of a bearing strength surface in dimensionless and normalized form, while detailed discussion is provided on the physical and geometrical interpretation of the kinematic mechanisms that accompany failure. A generalized closed-form expression for the failure envelope in M–Q–N space is then fitted to the numerical results with use of an appropriately trained artificial neural network. An upper-bound limit equilibrium solution for a certain failure mechanism (designated as the “sliding” mechanism) associated with maximum horizontal bearing capacity is also developed for verification purposes. One of the originalities of the paper lies with respect to the post-failure response of the caissons, where it is shown that the incremental displacement vector is accurately reproduced by assuming normality on the bearing strength surface irrespective of the considered plastic flow rule (associative or non-associative) at the microscale (soil element).     

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.