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.     

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.     

Correlations between soil gas and seismic activity in the Generalized Haiyuan Fault Zone,north-central China

Radon and mercury concentrations were measured in 10 fault gas profiles in Generalized Haiyuan Fault. This paper aims to predetermine the potential seismic risk in different segments of the fault zone from the perspective of geochemistry. The background value and anomaly threshold were adopted and synthesized using the maximum value method and average method to calculate concentration intensity values of radon and mercury. Fault soil gas mercury and radon concentrations show a decreasing gradient from NW to SE indicating evident segmentation. Higher values are mostly distributed in the Maomao Mountain–Tiger Mountain fault and Jingtai area. Combined with the seismotectonic background of historical and recent earthquakes and the spatial distribution characteristics of b-values, the fault soil gas concentration intensity shows a close correlation with earthquake activity within the fault zone. Concentrations of fault gas are higher and the b-value lower in areas of strong seismic activity, and regions with weak seismic activity correspond to lower fault gas concentrations and higher b-values. It is thus considered that the Jingtai area may be more dangerous than the other areas. This paper could provide vital background information for earthquake prediction in the Generalized Haiyuan Fault Zone.     

A laboratory study of anisotropic geomaterials incorporating recent micromechanical understanding

This paper presents an experimental investigation revisiting the anisotropic stress–strain–strength behaviour of geomaterials in drained monotonic shear using hollow cylinder apparatus. The test programme has been designed to cover the effect of material anisotropy, preshearing, material density and intermediate principal stress on the behaviour of Leighton Buzzard sand. Experiments have also been performed on glass beads to understand the effect of particle shape. This paper explains phenomenological observations based on recently acquired understanding in micromechanics, with attention focused on strength anisotropy and deformation non-coaxiality, i.e. non-coincidence between the principal stress direction and the principal strain rate direction. The test results demonstrate that the effects of initial anisotropy produced during sample preparation are significant. The stress–strain–strength behaviour of the specimen shows strong dependence on the principal stress direction. Preloading history, material density and particle shape are also found to be influential. In particular, it was found that non-coaxiality is more significant in presheared specimens. The observations on the strength anisotropy and deformation non-coaxiality were explained based on the stress–force–fabric relationship. It was observed that intermediate principal stress parameter b(b = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) has a significant effect on the non-coaxiality of sand. The lower the b-value, the higher the degree of non-coaxiality is induced. Visual inspection of shear band formed at the end of HCA testing has also been presented. The inclinations of the shear bands at different loading directions can be predicted well by taking account of the relative direction of the mobilized planes to the bedding plane.     

Energetic and spatial characterization of seismicity in the Algeria–Morocco region

We estimate the energetic and spatial characteristics of seismicity in the Algeria–Morocco region using a variety of seismic and statistical parameters, as a first step in a detailed investigation of regional seismic hazard. We divide the region into five seismotectonic regions, comprising the most important tectonic domains in the studied area: the Moroccan Meseta, the Rif, the Tell, the High Plateau, and the Atlas. Characteristic seismic hazard parameters, including the Gutenberg–Richter b-value, mean seismic activity rate, and maximum possible earthquake magnitude, were computed using an extension of the Aki–Utsu procedure for incomplete earthquake catalogs for each domain, based on recent earthquake catalogs compiled for northern Morocco and northern Algeria. Gutenberg–Richter b-values for each zone were initially estimated using the approach of Weichert (Bull Seismol Soc Am 70:1337–1346, 1980): the estimated b-values are 1.04 ± 0.04, 0.93 ± 0.10, 0.72 ± 0.03, 0.87 ± 0.02, and 0.77 ± 0.02 for the Atlas, Meseta, High Plateau, Rif, and Tell seismogenic zones, respectively. The fractal dimension D ₂ was also estimated for each zone. From the ratio D ₂/b, it appears that the Tell and Rif zones, with ratios of 2.09 and 2.12, respectively, have the highest potential earthquake hazard in the region. The Gutenberg–Richter relationship analysis allows us to derive that in the Tell and Rif, the number of earthquake with magnitude above Mw 4.0, since 1925 normalized to decade and to square cell with 100-km sides is equal to 2.6 and 1.91, respectively. This study provides the first detailed information about the potential seismicity of these large domains, including maximum regional magnitudes, characteristics of spatial clustering, and distribution of seismic energy release.     

Modeling of shear failure in fractured reservoirs with a porous matrix

A finite volume-based numerical modeling framework using a hierarchical fracture representation (HFR) has been developed to compute flow-induced shear failure. To accurately capture the mechanics near fracture manifolds, discontinuous basis functions are employed which ensure continuity of the displacement gradient across fractures. With these special basis functions, traction and compressive forces on the fracture segment can be calculated without any additional constraints, which is extremely useful for estimating the irreversible displacement along the fracture (slip) based on a constitutive friction law. The method is further extended to include slip-dependent hydraulic aperture change and grid convergent results are obtained. Further, the change in hydraulic aperture is modeled using an asymptotic representation which respects the experimentally observed behavior of pore volume dilation due to shear slip. The model allows the initial rapid increase in hydraulic aperture due to shear slip and asymptotically approaches a finite value after repeated shearing of a fracture segment. This aperture increase is the only feedback for mechanics into the fluid flow for a linear elastic mechanics problem. The same model is also extended to include poroelastic relations between flow and mechanics solver. The grid convergence result in the case of poroelastic flow-mechanics coupling for flow-induced shear failure is also obtained. This proves the robustness of the numerical and analytical modeling of fracture and friction in the extended finite volume method (XFVM) set-up. Finally, a grid convergent result for seismic moment magnitude for single fracture and fracture network with random initial hydraulic and friction properties is also obtained. The b-value, which represents the slope of seismic moment occurrence frequency decay vs seismic moment magnitude, which is approximately constant in a semi-logarithmic plot, is estimated. The numerical method leads to converged b-values for both single fracture and fracture network simulations, as grid and time resolutions are increased. For the resulting linear system, a sequential approach is used, that is, first, the flow and then the mechanics problems are solved. The new modeling framework is very useful to predict seismicity, permeability, and flow evolution in geological reservoirs. This is demonstrated with numerical simulations of enhancing a geothermal system.     

Influences of loading direction and intermediate principal stress ratio on the initiation of strain localization in cross-anisotropic sand

Since cross-anisotropic sand behaves differently when the loading direction or the stress state changes, the influences of the loading direction and the intermediate principal stress ratio (b = (σ ₂ ? σ ₃)/(σ ₁ ? σ ₃)) on the initiation of strain localization need study. According to the loading angle (angle between the major principal stress direction and the normal of bedding plane), a 3D non-coaxial non-associated elasto-plasticity hardening model was proposed by modifying Lode angle formulation of the Mohr–Coulomb yield function and the stress–dilatancy function. By using bifurcation analysis, the model was used to predict the initiation of strain localization under plane strain and true triaxial conditions. The predictions of the plane strain tests show that the major principal strain at the bifurcation points increases with the loading angle, while the stress ratio decreases with the loading angle. According to the loading angle and the intermediate principal stress ratio, the true triaxial tests were analyzed in three sectors. The stress–strain behavior and the volumetric strain in each sector can be well captured by the proposed model. Strain localization occurs in most b value conditions in all three sectors except for those which are close to triaxial compression condition (b = 0). The difference between the peak shear strength corresponding to the strain localization and the ultimate shear strength corresponding to plastic limit becomes obvious when the b value is near 0.4. The influence of bifurcation on the shear strength becomes weak when the loading direction changes from perpendicular to the bedding plane to parallel. The bifurcation analysis based on the proposed model gives out major principal strain and peak shear strength at the initiation of strain localization; the given results are consistent with experiments.     

China’s regional social vulnerability to geological disasters: evaluation and spatial characteristics analysis

The Andaman–Sumatra is one of the seismically active subduction zones and experienced three largest earthquakes in the recent past and rupturing more than 1,600-km-long portion of the plate boundary. The seismicity analysis of these large earthquakes source region (5°S–15°N latitude and 90°E–103°E longitude) has been carried out by several researchers and quantified the spatial and temporal variation of b-value which is a proxy to differential stress conditions and fractal dimension which is an indicator of material heterogeneity and strength. The results of all these studies clearly bring out the low b-value and low fractal dimension corresponding to locales were sizable magnitude earthquakes have occurred. Further locales of high stress regions are identified.     

A Criterion for Brittle Failure of Rocks Using the Theory of Critical Distances

This paper presents a new analytical criterion for brittle failure of rocks and heavily over-consolidated soils. Griffith’s model of a randomly oriented defect under a biaxial stress state is used to keep the criterion simple. The Griffith’s criterion is improved because the maximum tensile strength is not evaluated at the boundary of the defect but at a certain distance from the boundary, known as half of the critical distance. This fracture criterion is known as the point method, and is part of the theory of critical distances, which is utilised in fracture mechanics. The proposed failure criterion has two parameters: the inherent tensile strength, σ ₀, and the ratio of the half-length of the initial crack/flaw to the critical distance, a/L. These parameters are difficult to measure but they may be correlated with the uniaxial compressive and tensile strengths, σ _c and σ _t. The proposed criterion is able to reproduce the common range of strength ratios for rocks and heavily overconsolidated soils (σ _c/σ _t = 3–50) and the influence of several microstructural rock properties, such as texture and porosity. Good agreement with laboratory tests reported in the literature is found for tensile and low-confining stresses.     

Atmospheric elemental abundances for the components of the multiple system ADS 11061. 41 Draconis

We obtained speckle interferometric and spectroscopic observations of the system 41 Dra during its periastron passage in 2001. The components’ lines are resolved in the spectral interval 3700–9200 Å. The observed wavelength dependence of the brightness difference between the components is used to estimate the B-V indices separately for each of the components: B-V = 0.511 for component a and B-V = 0.502 for component b. We derived improved effective temperatures of the components from their B-V values and hydrogen-line profiles. The observations can be described with the parameters for the components T _eff ^a = 6370 K, log g^a = 4.05 and T _eff ^b = 6410 K, log g^b = 4.20. The iron, carbon, nitrogen, and oxygen abundances in the atmospheres of the components are log N(Fe)^a = 7.55, log N(Fe)^b = 7.60, log N(C)^a = 8.52, log N(C)^b = 8.58, log N(N)^a = 8.05, log N(N)^b = 7.99, log N(O)^a = 8.73, log N(O)^b = 8.76.     

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.     

Atmospheric chemical abundances of the components of the quadruple system ADS 11061. 40 Draconis

As part of our study of the components of the hierarchic quadruple system ADS 11061, we acquired spectroscopic observations of the binary 40 Dra. Echelle spectra showing the separation of the components’ lines were obtained in the spectral range 3700–9200 Å. Effective temperatures and surface gravities were derived for the components from BV photometry and the hydrogen-line profiles. The components of the 40 Dra system have parameters close to T _eff ^a = 6420 K, log g _a = 4.17, T _eff ^b = 6300 K, and log g _b = 4.20. We find the microturbulence velocity in the component atmospheres to be V _t = 2.6 km/s. The abundances of iron, carbon, nitrogen, and oxygen in the atmospheres of both components are estimated to be log N(Fe)^a = 7.50, log N(Fe)^b = 7.46, log N(C)^a = 8.39, log N(C)^b = 8.45, log N(N)^a = 8.12, log(N)^b = 8.15, log N(O)^a = 8.77, log N(O)^b = 8.74.     

Rheological properties of clayey soils originating from flow-like landslides

Flow-like landslides in clayey soils represent serious threats for populations and infrastructures and have been the subject of numerous studies in the past decade. However, despite the rising need for landslide mitigation with growing urbanization, the transient mechanisms involved in the solid-fluid transition are still poorly understood. One way of characterizing the solid-fluid transition is to carry out rheometrical tests on clayey soil samples to assess the evolution of viscosity with the shear stress. In this study, we carried out geotechnical and rheometrical tests on clayey samples collected from six flow-like landslides in order to assess if these clayey soils exhibit similar characteristics when they fluidize (solid-fluid transition). The results show that (1) all tested soils except one exhibit a yield-stress fluid behavior that can be associated with a bifurcation in viscosity (described by the critical shear rate \( \dot{\gamma_c} \)) and in shear modulus G; (2) the larger the amplitude of the viscosity bifurcation, the larger the associated drop in G; and (3) the water content (w) deviation from the Atterberg liquid limit (LL) seem a key parameter controlling a common mechanical behavior of these soils at the solid-fluid transition. We propose exponential laws describing the evolution of the critical shear stress τ_c, the critical shear rate \( \dot{\gamma_c} \), and the shear modulus G as a function of the deviation w-LL.     

Sequence signal processing of 2012 Ahar–Varzaghan earthquake (<Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> 6.4) of NW Iran

Before starting seismic cycle of Ahar–Varzaghan 2012 event, a partial gap in the form of a pre-seismic calm sequence (seismicity rate, r = 0.46 event/year, b = 1.4) with duration of 303 days spatially has dominated over the entire seismogenic area. From April 17, 2012, to May 31, 2012, r significantly increased to 2.16, indicating strong foreshock sequence, and b value changed to 1.9, remarkably. In the last two months before the mainshock, foreshocks have partially migrated toward the earthquake fault (with a decrease in size, b = 2.0). Significantly, high rate of seismicity and low V _P /V _S (1.64) in the foreshocks sequence and also very high seismicity rate (17.3) and high V _P /V _S (1.76) in the aftershocks sequence make substantial differences between the seismic cycle and the background seismicity. Moreover, a significant E–W migration of the microseismicity was confirmed in the study area.     

Dependence of the absorption-line profiles and radial-velocity curve of the optical star in an X-ray binary on the orbital inclination and component-mass ratio

Theoretical absorption-line profiles and radial-velocity curves for tidally deformed optical stars in X-ray binary systems are calculated assuming LTE. The variations in the profile shapes and radial-velocity curve of the optical star are analyzed as a function of the orbital inclination of the X-ray binary system. The dependence of the shape of the radial-velocity curve on the orbital inclination i increases with decreasing component-mass ratio q = m _x /m _v . The integrated line profiles and radial-velocity curves of the optical star are calculated for the Cyg X-1 binary, which are then used to estimate the orbital inclination and mass of the relativistic object: i < 43° andm _x = 8.2–12.8 M_⊙. These estimates are in good agreement with earlier results of fitting the radial-velocity curve of Cyg X-1 using a simpler model (i < 45°, m _x = 9.0–13.2 M_⊙).     

Simplified method for calculating inflow into a deep excavation with consideration of the effects of cutoff walls

Water inflow into deep excavations is a key parameter in the evaluation of environmental impact, and a simplified method is needed to calculate the inflow with adequate accuracy, especially for a deep excavation with a cutoff wall around it. This paper presents a model for calculating inflow based on the method of fragments, assuming two-dimensional sectional flow. The effects of cutoff wall thickness (w), wall penetration, and excavation width on the inflow are well considered in this model. Explicit formulas for the form factor are given. The accuracy and scope of this method are discussed. The method can be applied in both flooding and no flooding cases. Extensive verification shows that, for the flooding case and when one of four conditions is satisfied, the error of the proposed method will be within 10%. These conditions are b/d?≥?0.8 or w/d?≥?0.03 or s/d?≥?0.2 or s₁/d?≥?0.2, for which b is the half width of the excavation, d is the thickness of the soil layer underlying the wall, and s and s₁ are the wall penetration depths on the non-excavation side and excavation side, respectively. For the case of no flooding and when the ratio of the distance from the recharge boundary to the aquifer thickness (a/T) is?≥2, the error does not exceed 15%.     

Pile reinforcement mechanism of soil slopes

Stabilizing piles are widely used as an effective and economic reinforcement approach for slopes. Reasonable designs of pile reinforcement depend on the understanding of reinforcement mechanism of slopes. A series of centrifuge model tests were conducted on the pile-reinforced slopes and corresponding unreinforced slopes under self-weight and vertical loading conditions. The deformation of the slope was measured using image-based analysis and employed to investigate the pile reinforcement mechanism. The test results showed that the piles significantly reduced the deformation and changed the deformation distribution of the slope, and prevented the failure occurred in the unreinforced slope. The pile influence zone was determined according to the inflection points on the distribution curves of horizontal displacement, which comprehensively described the features of the pile–slope interaction and the characteristics of reinforced slopes. The concepts of anti-shear effect and compression effect were proposed to quantitatively describe the restriction features of the piles on the deformation of the slope, namely the reduction in the shear deformation and the increase in the compression deformation, respectively. The pile reinforcement effect mainly occurred in the pile influence zone and decreased with increasing distance from the piles. There was a dominated compression effect in the vicinities of the piles. The compression effect developed upwards in the slope with a transmission to the anti-shear effect. The anti-shear effect became significantly dominated near the slip surface and prevented the failure that occurred in the unreinforced slope.     

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.     

Accuracy assessment of near-shore bathymetry information retrieved from Landsat-8 imagery

The improvement in the capabilities of Landsat-8 imagery to retrieve bathymetric information in shallow coastal waters was examined. Landsat-8 images have an additional band named coastal/aerosol, Band 1: 435–451 nm in comparison with former generation of Landsat imagery. The selected Landsat-8 operational land image (OLI) was of Chabahar Bay, located in the southern part of Iran (acquired on February 22, 2014 in calm weather and relatively low turbidity). Accurate and high resolution bathymetric data from the study area, produced by field surveys using a single beam echo-sounder, were selected for calibrating the models and validating the results. Three methods, including traditional linear and ratio transform techniques, as well as a novel proposed integrated method, were used to determine depth values. All possible combinations of the three bands [coastal/aerosol (CB), blue (B), and green (G)] have been considered (11 options) using the traditional linear and ratio transform techniques, together with five model options for the integrated method. The accuracy of each model was assessed by comparing the determined bathymetric information with field measured values. The standard error of the estimates, correlation coefficients (R ² ) for both calibration and validation points, and root mean square errors (RMSE) were calculated for all cases. When compared with the ratio transform method, the method employing linear transformation with a combination of CB, B, and G bands yielded more accurate results (standard error = 1.712 m, R ² _calibration = 0.594, R ² _validation = 0.551, and RMSE =1.80 m). Adding the CB band to the ratio transform methodology also dramatically increased the accuracy of the estimated depths, whereas this increment was not statistically significant when using the linear transform methodology. The integrated transform method in form of Depth = b ₀  + b ₁ X _CB  + b ₂ X _B  + b ₅ ln(R _CB )/ln(R _G ) + b ₆ ln(R _B )/ln(R _G ) yielded the highest accuracy (standard error = 1.634 m, R ² _calibration = 0.634, R ² _validation = 0.595, and RMSE = 1.71 m), where R _i (i = CB, B, or G) refers to atmospherically corrected reflectance values in the i ^th band [X _i  = ln(R _i -R _{deep water})].     

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.