Precursory seismicity rate changes prior to the large and major earthquakes along the Sagaing fault zone,Central Myanmar

In this study, the seismicity rate changes that can represent an earthquake precursor were investigated along the Sagaing Fault Zone (SFZ), Central Myanmar, using the Z value technique. After statistical improvement of the existing seismicity data (the instrumental earthquake records) by removal of the foreshocks and aftershocks and man-made seismicity changes and standardization of the reported magnitude scales, 3574 earthquake events with a M _w ≥ 4.2 reported during 1977–2015 were found to directly represent the seismotectonic activities of the SFZ. To find the characteristic parameters specifically suitable for the SFZ, seven known events of M _w ≥ 6.0 earthquakes were recognized and used for retrospective tests. As a result, utilizing the conditions of 25 fixed earthquake events considered (N) and a 2-year time window (T _w), a significantly high Z value was found to precede most of the M _w ≥ 6.0 earthquakes. Therefore, to evaluate the prospective areas of upcoming earthquakes, these conditions (N = 25 and T _w = 2) were applied with the most up-to-date seismicity data of 2010–2015. The results illustrate that the vicinity of Myitkyina and Naypyidaw (Z = 4.2–5.1) cities might be subject to strong or major earthquakes in the future.     

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.     

Climate extremes: an observation and projection of its impacts on food production in ASEAN

Multifractal behaviour of interevent time sequences is investigated for the earthquake events in the NW Himalaya, which is one of the most seismically active zones of India and experienced moderate to large damaging earthquakes in the past. In the present study, the multifractal detrended fluctuation analysis (MF-DFA) is used to understand the multifractal behaviour of the earthquake data. For this purpose, a complete and homogeneous earthquake catalogue of the period 1965–2013 with a magnitude of completeness M _w 4.3 is used. The analysis revealed the presence of multifractal behaviour and sharp changes near the occurrence of three earthquakes of magnitude (M _w ) greater than 6.6 including the October 2005, Muzaffarabad–Kashmir earthquake. The multifractal spectrum and related parameters are explored to understand the time dynamics and clustering of the events.     

Radial pulsations of helium stars and a model for BX CIR

The results of hydrodynamical calculations of radially pulsating helium stars with masses 0.5M_⊙≤M≤0.9M_⊙, bolometric luminosities 600L_⊙≤5×10³L_⊙, and effective temperatures 1.5×10⁴ K≤T_eff≤3.5×10⁴ K are presented. The pulsation instability of these stars is due to the effects of ionization of iron-group elements in layers with temperatures T～2×10⁵ K. The calculations were carried out using opacities for the relative mass abundances of hydrogen and heavy elements X=0 and Z=0.01, 0.015, and 0.02. Approximate formulas for the pulsation constant Q over the entire range of pulsation instability of the hot helium stars in terms of the mass M, radius R, effective temperature T_eff, and heavy-element abundance Z are derived. The instability of BX Cir to radial pulsations with the observed period Π=0.1066 d occurs only for a mass M≥0.55M_⊙, effective temperature T_eff≥23000 K, and heavy-element abundance Z≥0.015. The allowed mass of BX Cir is in the range 0.55M_⊙≤M≤0.8M_⊙, which corresponds to luminosities 800L_⊙≤M≤1400L_⊙ and mean radii 1.7R_⊙?R?2.1R_⊙.     

Site response of the Ganges basin inferred from re-evaluated macroseismic observations from the 1897 Shillong, 1905 Kangra,and 1934 Nepal earthquakes

We analyze previously published geodetic data and intensity values for the M _s = 8.1 Shillong (1897), M _s = 7.8 Kangra (1905), and M _s = 8.2 Nepal/Bihar (1934) earthquakes to investigate the rupture zones of these earthquakes as well as the amplification of ground motions throughout the Punjab, Ganges and Brahmaputra valleys. For each earthquake we subtract the observed MSK intensities from a synthetic intensity derived from an inferred planar rupture model of the earthquake, combined with an attenuation function derived from instrumentally recorded earthquakes. The resulting residuals are contoured to identify regions of anomalous intensity caused primarily by local site effects. Observations indicative of liquefaction are treated separately from other indications of shaking severity lest they inflate inferred residual shaking estimates. Despite this precaution we find that intensites are 1–3 units higher near the major rivers, as well as at the edges of the Ganges basin. We find evidence for a post-critical Moho reflection from the 1897 and 1905 earthquakes that raises intensities 1–2 units at distances of the order of 150 km from the rupture zone, and we find that the 1905 earthquake triggered a substantial subsequent earthquake at Dehra Dun, at a distance of approximately 150 km. Four or more M = 8 earthquakes are apparently overdue in the region based on seismic moment summation in the past 500 years. Results from the current study permit anticipated intensities in these future earthquakes to be refined to incorporate site effects derived from dense macroseismic data.     

Probabilistic seismic hazard assessment in the Constantine region,Northeast of Algeria

This paper presents a seismic hazard evaluation and develops an earthquake catalogue for the Constantine region over the period from 1357 to 2014. The study contributes to the improvement of seismic risk management by evaluating the seismic hazards in Northeast Algeria. A regional seismicity analysis was conducted based on reliable earthquake data obtained from various agencies (CRAAG, IGN, USGS and ISC). All magnitudes (M _l, m _b) and intensities (I ₀, I _MM, I _MSK and I _EMS) were converted to M _s magnitudes using the appropriate relationships. Earthquake hazard maps were created for the Constantine region. These maps were estimated in terms of spectral acceleration (SA) at periods of 0.1, 0.2, 0.5, 0.7, 0.9, 1.0, 1.5 and 2.0 s. Five seismogenic zones are proposed. This new method differs from the conventional method because it incorporates earthquake magnitude uncertainty and mixed datasets containing large historical events and recent data. The method can be used to estimate the b value of the Gutenberg-Richter relationship, annual activity rate λ(M) of an event and maximum possible magnitude M _max using incomplete and heterogeneous data files. In addition, an earthquake is considered a Poisson with an annual activity rate λ and with a doubly truncated exponential earthquake magnitude distribution. Map of seismic hazard and an earthquake catalogue, graphs and maps were created using geographic information systems (GIS), the Z-map code version 6 and Crisis software 2012.     

Seismic anisotropy of the mantle beneath Eastern Asia based on <Emphasis Type="Italic">ScS</Emphasis> and <Emphasis Type="Italic">S</Emphasis> waves from deep-focus earthquakes

The seismic anisotropy of the mantle is studied based on the data of S and ScS waves from earthquakes occurred in the mantle transition zone over the period of 2007–2013 and recorded by seismic stations in the continental margin of Asia, on Sakhalin Island, and in the southern part of the Kamchatka Peninsula. The measurements of the azimuths of polarization of the fast S and ScS waves in the continental margin of Asia show that they are predominantly oriented in the E–SE directions. Based on the distribution of the shear wave splitting parameters, the symmetry of the medium can be described in terms of a transversely isotropic model with a horizontal symmetry axis and may correspond to horizontal flow in the upper mantle beneath the Amur Plate. The fast azimuths of polarization of ScS wave, which were determined to be of N–NE directions in the northern area of Sakhalin Island and in the continental part of Asia, may correspond to an inclined flow under the conditions of oblique subduction and complex geometry of the downgoing Pacific Plate. In the south of the Kamchatka Peninsula, the S- and ScS-wave azimuths of polarization from the M 8.4 Sea of Okhotsk earthquake are determined to be oriented along the direction of the Pacific Plate motion. The fast-S-wave azimuths of polarization from the aftershocks of the Sea of Okhotsk earthquake and from other large events of 2008–2009 are determined to be nearly parallel to the motion trend of the Pacific Plate, but orthogonal to it for the events of 2008–2009. On the basis of the distribution of azimuths of polarization of the fast S waves, the symmetry of the medium can be described in terms of a transversely isotropic model with the symmetry axis inclined orthogonally to the plane of downgoing plate and oriented westward orthogonally to the trench strike.     

Absorption in the H(93–95)<Emphasis Type="Italic">α</Emphasis> and H(78–79)<Emphasis Type="Italic">α</Emphasis> recombination lines in the HII region of the quasar 3C 345

Absorption of the synchrotron emission of the quasar 3C 345 in the continuum and H(93–95)α and H(78–79)α radio recombination lines is studied. The upper limit for absorption in the H(93–95)α lines is T_al/T_ac < 0.7%; absorption in the H(78–79)α lines with antenna temperature T_al = 25 mK, linewidth Δf = 5.3 ± 0.08 MHz, and T_al/T_ac ≥ 0.3% has been detected. A correction to the redshift Δz = 0.00135 ± 0.00008 (z = 0.59365) has been determined.     

Masses of the visual components and black holes in X-ray novae: Effects of proximity of the components

It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V _rot sin i and underestimation of the component mass ratio q = M _x /M _v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M _x and visual components M _v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M _x /M _v , reducing M _v , while M _x changes only slightly. The resulting distribution of M _v attains its maximum near M _v ? 0.35M _⊙, in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M _x also differs strongly from the mass distribution for massive stars in the Galaxy.     

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.     

Structural systematics of Ge substitution in primitive pyroxenes

Pyroxenes of general stoichiometry Mg(Ge _x Si_1?x )O₃ were encountered in attempts to synthesise Ge-substituted talcs at 0.2 GPa, 650–700 °C. Orthopyroxenes (Pbca) of compositions x = 0.21, 0.30, and 0.34 were identified, and also a P2₁/c clinopyroxene of composition x = 0.63, and C2/c clinopyroxenes of compositions x = 0.91 and 1. End-member clinoenstatite MgSiO₃-P2₁/c synthesised at 16 GPa, 1300 °C and transformed from C2/c was also included in the study. Crystal structure refinements using single-crystal XRD data showed that unit-cell parameters vary linearly with Si–Ge for the Pbca and P2₁/c pyroxenes, both of which have two symmetrically non-equivalent tetrahedral chains. Refinement of Si–Ge occupancies at tetrahedral sites showed that the two chains of all primitive pyroxenes have very different compositions, with X_Ge(TB) ? X_Ge(TA). This difference arises from the greater flexibility of the B-chain to rotate in response to tetrahedral expansion due to increasing Ge content. The TA-M2 shared polyhedral edge imposes significant constraints on the flexibility of the A-chain, which can accommodate much less Ge than the B-chain. Linear trends of cell parameters, site occupancies, and structural parameters for the primitive pyroxenes, when extrapolated to published data for MgGeO₃–Pbca, extend across the entire Si–Ge join.     

Nonlinear radial pulsations of Wolf-Rayet stars

The evolution of Population I stars with initial masses 60 M _⊙ ≤ M _ZAMS ≤ 120 M _⊙ is computed up to the Wolf-Rayet stage, when the central helium abundance decreases to Y _c ≈ 0.05. Several models from evolutionary sequences in the core helium-burning stage were used as initial conditions when solving the equations of radiative hydrodynamics for self-exciting stellar radial pulsations. The low-density envelope surrounding the compact core during the core helium burning is unstable against radial oscillations in a wide range of effective temperatures extending to T _eff ～ 10⁵ K. The e-folding time of the amplitude growth is comparable to the dynamical time scale of the star, and, when the instability ceases growing, the radial displacement of the outer layers is comparable to the stellar radius. Evolutionary changes of the stellar radius and luminosity are accompanied by a decrease in the amplitude of radial pulsations, but, at the effective temperature T _eff ≈ 10⁵ K, the stellar oscillations are still nonlinear, with a maximum expansion velocity of the outer layers of about one-third the local escape velocity. The period of the radial oscillations decreases from 9 hr to 4 min as stellar mass decreases from M = 28 M _⊙ to M = 6 M _⊙ in the course of evolution. The nonlinear oscillations lead to a substantial increase of the radii of the Lagrangian mass zones compared to their equilibrium radii throughout the instability region. The instability of Wolf-Rayet stars against radial oscillations is due to the action of the κ mechanism in the iron-group ionization zone, which has a temperature of T ～ 2 × 10⁵ K.     

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.     

Application of time- and magnitude-predictable model for long-term earthquake prediction in Iran

The regional time- and magnitude-predictable model has been applied successfully in diverse regions of the world to describe the occurrence of main shocks. In the current study, the model has been calibrated against the historical and instrumental catalog of Iranian earthquakes. The Iranian plateau is divided into 15 seismogenic provinces; then, the interevent times for strong main shocks have been determined for each one. The empirical relations reported by Papazachos et al. (Tectonophysics 271:295–323, 1997a) for the Alpine–Himalayan belt (including Iran) were adopted except for the constant terms that were calculated separately for every seismotectonic area. By using the calibrated equations developed for the study area and taking into account the occurrence time and magnitude of the last main shocks in each seismogenic source, the time-dependent conditional probabilities of occurrence P(?t) of the next main shocks during next 10, 20, 30, 40 and 50 years as well as the magnitude of the expected main shocks (M _f) have been estimated. The immediate probability (within next 10 years) of a large main shock is estimated to be high and moderate (>35 %) in all regions except zones 9 (M _f = 5.8) and 15 (M _f = 6.1). However, it should be noted that the probabilities have been estimated for different M _f values in 15 regions. Comparing the model predictions with the actual earthquake occurrence rates shows the good performance of the model for Iranian plateau.     

Abundances of carbon,nitrogen, oxygen,and other elements in the atmospheres of the giants 25 Mon and HR 7389

An analysis of high-resolution CCD spectra of the giant 25 Mon, which shows signs of metallicity, and the normal giant HR 7389 is presented. The derived effective temperatures, gravitational accelerations, and microturbulence velocities are T_eff = 6700 K, log g = 3.24, and ξ _t = 3.1 km/s for 25 Mon and T_eff = 6630 K, log g = 3.71, and ξ _t = 2.6 km/s for HR 7389. The abundances (log ε) of nine elements are determined: carbon, nitrogen, oxygen, sodium, silicon, calcium, iron, nickel, and barium. The derived excess carbon abundances are 0.23 dex for 25 Mon and 0.16 dex for HR 7389. 25 Mon displays a modest (0.08 dex) oxygen excess, with the oxygen excess for HR 7389 being somewhat higher (0.15 dex). The nitrogen abundance is probably no lower than the solar value for both stars. The abundances of iron, sodium, calcium (for HR 7389), barium, and nickel exceed the solar values by 0.22–0.40 dex for both stars. The highest excess (0.62 dex) is exhibited by the calcium abundance for 25 Mon. Silicon displays a nearly solar abundance in both stars—small deficits of ?0.03 dex and ?0.07 dex for 25 Mon and HR 7389, respectively. No fundamental differences in the elemental abundances were found in the atmospheres of 25 Mon and HR 7389. Based on their T_eff and log g values, as well as theoretical calculations, A. Claret estimated the masses, radii, luminosities, and ages of 25 Mon (M/M_⊙ = 2.45, log(R/R_⊙) = 0.79, log(L/L_⊙) = 1.85, t = 5.3 × 10⁸ yr) and HR 7389 (M/M_⊙ = 2.36, log(R/R_⊙) = 0.50, log(L/L_⊙) = 1.24, t = 4.6 × 10⁸ yr), and also of the stars 20 Peg (M/M_⊙ = 2.36, log(R/R_⊙) = 0.73, log(L/L_⊙) = 1.79, t = 4.9 × 10⁸ yr) and 30 LMi (M/M_⊙ = 2.47, log(R/R_⊙) = 0.73, log(L/L_⊙) = 1.88, t = 4.8 × 10⁸ yr) studied by the author earlier.     

Seismic source characteristics in Kachchh and Saurashtra regions of Western India: <Emphasis Type="Italic">b</Emphasis>-value and fractal dimension mapping of aftershock sequences

Seismic source characteristics in the Kachchh rift basin and Saurashtra horst tectonic blocks in the stable continental region (SCR) of western peninsular India are studied using the earthquake catalog data for the period 2006–2011 recorded by a 52-station broadband seismic network known as Gujarat State Network (GSNet) running by Institute of Seismological Research (ISR), Gujarat. These data are mainly the aftershock sequences of three mainshocks, the 2001 Bhuj earthquake (M _w 7.7) in the Kachchh rift basin, and the 2007 and 2011 Talala earthquakes (M _w ≥ 5.0) in the Saurashtra horst. Two important seismological parameters, the frequency–magnitude relation (b-value) and the fractal correlation dimension (D _c) of the hypocenters, are estimated. The b-value and the D _c maps indicate a difference in seismic characteristics of these two tectonic regions. The average b-value in Kachchh region is 1.2 ± 0.05 and that in the Saurashtra region 0.7 ± 0.04. The average D _c in Kachchh is 2.64 ± 0.01 and in Saurashtra 2.46 ± 0.01. The hypocenters in Kachchh rift basin cluster at a depth range 20–35 km and that in Saurashtra at 5–10 km. The b-value and D _c cross sections image the seismogenic structures that shed new light on seismotectonics of these two tectonic regions. The mainshock sources at depth are identified as lower b-value or stressed zones at the fault end. Crustal heterogeneities are well reflected in the maps as well as in the cross sections. We also find a positive correlation between b- and D _c-values in both the tectonic regions.     

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.     

The formation and evolution of the most massive stars

We consider the formation of massive stars under the assumption that a young star accretes material from the protostellar cloud through its accretion disk while losing gas in the polar directions via its stellar wind. The mass of the star reaches its maximum when the intensity of the gradually strengthening stellar wind of the young star becomes equal to the accretion rate. We show that the maximum mass of the forming stars increases with the temperature of gas in the protostellar cloud T ₀, since the rate at which the protostellar matter is accreted increases with T ₀. Numerical modeling indicates that the maximum mass of the forming stars increases to ～900 M _⊙ for T ₀ ～ 300 K. Such high temperatures of the protostellar gas can be reached either in dense star-formation regions or in the vicinity of bright active galactic nuclei. It is also shown that, the lower the abundance of heavy elements in the initial stellar material Z, the larger the maximum mass of the star, since the mass-loss rate due to the stellar wind decreases with decreasing Z. This suggests that supermassive stars with masses up to 10⁶ M _⊙ could be formed at early stages in the evolution of the Universe, in young galaxies that are almost devoid of heavy elements. Under the current conditions, for T ₀ = (30–100) K, the maximum mass of a star can reach ～100M _⊙, as is confirmed by observations. Another opportunity for the most massive stars to increase their masses emerges in connection with the formation and early stages of evolution of the most massive close binary systems: the most massive stars can be produced either by coalescence of the binary components or via mass transfer in such systems.     

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.     

The probable quadruple systems FK Aql and FZ Del

Orbital-period variations of the eclipsing binaries FK Aql and FZ Del are analyzed. For each of the systems, a superposition of two cyclic variations of their orbital periods is found. FK Aql may be a quadruple system that contains two more bodies, besides the eclipsing binary, with masses M ₃ ? 1.75M _⊙ and M ₄ ? 1.47M _⊙, and the corresponding periods 15 and 82 yrs. This could also be a triple system with a third body of mass M ₃ ? 1.75M _⊙ and a period of the long-period orbit P ₃ = 15 yrs, or with a third body of mass M ₃ ? 1.30M _⊙ and a period of the long-period orbit P ₃ = 82 yrs. FZ Del may be a quadruple system with the additional componentmasses M ₃ ? 0.2M _⊙ and M ₄ ? 0.3M _⊙, with the periods 10.2 and 53.7 yrs. This could also be a triple system with a third-body mass M ₃ ? 0.2M _⊙ and a period of the long-period orbit P ₃ = 10.2 yrs. In both systems, the residual period variations could be due to magnetic cycles of the secondary. The period variations of the eclipsing binary FZ Del could also be due to apsidal motion, together with the influence of a third body or the effects of magnetic activity.