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.     

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_⊙.     

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.     

Magnetic field models for HD 116458 and HD 126515

We have modeled the magnetic fields of the slowly rotating stars HD 116458 and HD 126515 using the “magnetic charge” technique. HD 116458 has a small angle between its rotation axis and dipole axis (β = 12°), whereas this angle is large for HD 126515 (β = 86°). Both stars can be described with a decentered-dipole model, with the respective displacements being r = 0.07 and r = 0.24 in units of the stellar radius. The decentered-dipole model is able to satisfactorily explain the phase relations for the effective field, B_e(P), and the mean surface field, B_s(P), for both stars, along with the fact that the B_e(P) phase relation for HD 126515 is anharmonic. We discuss the role of systematic measurement errors possibly resulting from instrumental or methodical effects in one or both of the phase relations. The displacement of the dipole probably reflects real asymmetry of the stellar field structure, and is not due to measurement errors. Using both phase relations, B_e(P) and B_s(P), in the modeling considerably reduces the influence of the nonuniform distribution of chemical elements on the stellar surface.     

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

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

Physicochemical model for the crystallization of rocks of the calc-alkaline series

The Al-rich region of the CaO-MgO-Al₂O₃-SiO₂ system was experimentally studied at pressures of 1.0–2.8 GPa and temperatures of 1300–1535°C. The slopes of the lines of the monovariant reactions An + Sp = Cpx + Cor + (Ga) and L = Cpx + Ga + Cor + Sp and the compositions of the phases involved in these reactions are determined. The results are utilized in the topological analysis of the aluminous region of the CaO-MgO-Al₂O₃-SiO₂ system. On this basis, the principal structure of the phase diagram is analyzed, and a phase diagram is constructed for the junction region of the quaternary system and the CaO-Al₂O₃-SiO₂ ternary system. A continuous series of the monovariant eutectic: L = Cpx + Opx + Fo + An, L = Cpx + Opx + An + Sp, L = Cpx + (Ga) + An + Sp, L = Cpx + Cor + (Ga) + An, L = An + Ga + Cpx + Ky and L = Ga + Cpx + Ky + Qz is examined within the pressure range from atmospheric to 3.0 GPa and higher. Analogous “telescoped” eutectic series are of fundamental character for interpreting the evolution of magmatic melts. A physicochemical model is suggested for the evolution of magmatic melts that produce rocks of the calc-alkaline series, with this model underlain by the fact that a change in the composition of magmatic melt at a pressure decrease should correspond to the minimum melting temperatures, i.e., to melts in the fundamental series of eutectic reactions. The comparison of our physicochemical model and rocks of the calc-alkaline series shows that the compositions of rocks of the calc-alkaline series are close to the compositions determined for the eutectic equilibria, and the mineralogical composition of xenoliths and megacrysts in volcanic and dike varieties of the rocks are similar to the subsolidus phases of the established fundamental eutectic system.     

Photosynthetic capability and Fe,Mn, Cu,and Zn contents in two Moraceae species under different phosphorus levels

The strong adaptability of Broussonetia papyrifera (L.) Vent. to low phosphorus (P) conditions can be attributed to the large amount of root-exuded organic acids and the high efficiency of P extraction. However, microelement contents are influenced by low-P stress, and their effects on the photosynthetic capability of B. papyrifera remain unknown. In this study, we investigated the effects of low-P treatment on net photosynthetic rate (P _N); chlorophyll a fluorescence (ChlF) characteristics; and Fe, Mn, Cu, and Zn contents of B. papyrifera and Morus alba L. seedlings. Results show that B. papyrifera exhibited better photosynthetic capability under moderate P deficiency (0.125, 0.063, and 0.031 mmol/L P treatments), whereas the photosynthetic capability of M. alba decreased under moderate and severe P deficiency (0.016 and 0 mmol/L P treatments). Under moderate P deficiency, the decrease in Cu and Zn contents in B. papyrifera was lower than that in M. alba. Under severe P deficiency, a considerable decrease of photosynthetic capability in B. papyrifera and M. alba was associated with low Cu and Zn contents. The P _N of the two Moraceae species exhibited a better correlation with Cu and Zn contents than with Fe or Mn content. P deficiency could not only decrease cyclic photophorylation and photosynthetic efficiency, but could also affect the stability of thylakoid membrane structure and electron transport efficiency by influencing the contents of Cu or Zn, thereby affecting photosynthesis.     

Field testing of one-way and two-way cyclic lateral responses of single and jet-grouting reinforced piles in soft clay

Piles supporting transmission towers, offshore structures (such as wind turbines), or infrastructures in seismic areas are frequently subjected to either one-way or two-way cyclic lateral loadings. Relatively little attention, however, has been paid to compare and understand the effects of different loading regimes (one-way or two-way cycling) on lateral responses of piles in soft clay. For this reason, a series of field tests in soft clay are carried out to compare one-way and two-way cyclic responses of single piles and of jet-grouting reinforced piles. The field tests reveal that the single pile subjected to two-way cycling experiences much more rapid degradation in lateral stiffness and capacity, but accumulates much smaller residual pile deflection (δ _p), than the single pile under one-way cycling. This is because the reverse part of the two-way cycling also generates plastic strain, causing additional softening and strength reduction in the soil surrounding the pile. After each cycling, non-zero bending moment (i.e. locked in moment, or M _L) is retained in the single piles, and the M _L increases with the δ _p. The one-way cycling leads to two times larger M _L than the two-way cycling, as it causes greater δ _p. The maximum M _L in the pile after one-way cycling can be up to 40% of the maximum bending moment induced during the previous cyclic loading stage. After application of jet-grouting surrounding the upper part of the single pile, it greatly reduces degradation of lateral pile stiffness, accumulation of δ _p and therefore development of M _L. Compared to the field measurements, the API (API RP 2A-WSD, recommended practice for planning, designing, and constructing fixed offshore platform-working stress design, 21st edn. API, Washington, 2000) code underestimates the lateral stiffness of the pile under one-way cycling, while overestimates that of the pile under two-way cycling, leading to a non-conservative prediction of bending moment in the latter pile.     

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.     

New dwarf galaxy candidates in the Leo-I Group

We have carried out a search for low-surface-brightness dwarf galaxies in the region of the Leo-I Group (M96) in images of the second Palomar Sky Survey. We found a total of 36 likely dwarf members of the group with typical magnitudes B _t ～18^m–19^m in an area of sky covering 120 square degrees. Half of these galaxies are absent from known catalogs and lists of galaxies. The radial-velocity dispersion calculated for 19 galaxies is 130 km/s. The Leo-I Group has a mean distance from the Sun of 10.4 Mpc, a mean projected radius of 352 kpc, an integrated luminosity of 6.7 × 10¹⁰L_⊙, a virial mass-to-luminosity ratio of 107 M_⊙/L_⊙, and a mean crossing time of 2.7 Gyr. The group shows evidence for a radial segregation of the galaxies according to morphological type and luminosity, suggesting that the group is in a state of dynamical relaxation. The subsystem of bright galaxies in the Leo-I Group is smaller in size (250 kpc) and has a lower velocity dispersion (92 km/s), resulting in a lower virial mass-to-luminosity ratio (34 M_⊙/L_⊙), as is typical of the Local Group and other nearby groups of galaxies.     

Earthquake hazard and risk assessment based on unified scaling law for earthquakes: Altai–Sayan Region

We apply the general concept of seismic risk analysis based on morphostructural analysis of the territory, pattern recognition of earthquake-prone nodes, and the Unified Scaling Law for Earthquakes, USLE, in another seismic region of Russia to the west from Lake Baikal, i.e., Altai–Sayan Region. The USLE generalizes the empirical Gutenberg–Richter relationship making use of apparently fractal distribution of earthquake sources of different size: \( \log_{10} N\left( {M,L} \right)\, = \,A\, + \,B \cdot \left( {5\, - \,M} \right)\, + \,C \cdot \log_{10} L, \) where N (M, L) is the expected annual number of earthquakes of a certain magnitude M within an seismically prone area of linear dimension L. The local estimates of A, B, and C allow determination of the expected maximum credible magnitude in a given time interval and the associated spread around ground shaking parameters (e.g., peak ground acceleration, PGA, or macroseismic intensity, I₀). Compilation of the corresponding seismic hazard map of Altai–Sayan Region and its rigorous testing against the available seismic evidences in the past is used to model regional maps of specific earthquake risks for population, cities, and infrastructures.     

Prediction of longitudinal dispersion coefficient using multivariate adaptive regression splines

In this paper, multivariate adaptive regression splines (MARS) was developed as a novel soft-computing technique for predicting longitudinal dispersion coefficient (D_L) in rivers. As mentioned in the literature, experimental dataset related to D_L was collected and used for preparing MARS model. Results of MARS model were compared with multi-layer neural network model and empirical formulas. To define the most effective parameters on D_L, the Gamma test was used. Performance of MARS model was assessed by calculation of standard error indices. Error indices showed that MARS model has suitable performance and is more accurate compared to multi-layer neural network model and empirical formulas. Results of the Gamma test and MARS model showed that flow depth (H) and ratio of the mean velocity to shear velocity (u/u^?) were the most effective parameters on the D_L.     

Properties of galactic disks at optical and near-IR wavelengths

We have analyzed the radial scales, central surface brightnesses, and colors of 400 disks of various types of galaxies. For nine galaxies, the brightness decrease and the central disk brightness were obtained via a two-dimensional decomposition of the U BV RI J H K photometric images into bulge and disk components. We used published disk parameters for 392 of the galaxies. The central surface brightness μ _0,i ⁰ and linear (disk) scale length h vary smoothly along the Hubble sequence of galaxies within a rather narrow interval. The disks of relatively early-type galaxies display higher central K surface brightnesses, higher central surface densities, higher central mass-to-luminosity ratios M/L(B), smaller sizes (relative to the diameter of the galaxy D ₂₅), redder integrated colors, and redder central colors. The color gradient normalized to the radius of the galaxy and the “blue” central surface brightness of the disk, μ _0,i/0(B), are both independent of the galaxy type. The radial disk scales in different photometric bands differ less in early-type than in late-type galaxies. A correlation between the central disk surface brightness and the total luminosity of the galaxy is observed. We also consider the influence of dust on the photometric parameters of the disks.     

Supermassive black holes: Relation to dark halos

Estimates of the masses of supermassive black holes (M _bh ) in the nuclei of disk galaxies with known rotation curves are compared with estimates of the rotational velocities V _m and the “indicative” masses of the galaxies M _i . Although there is a correlation between M _bh and V _m or M _i , it is appreciably weaker than the correlation with the central velocity dispersion. The values of M _bh for early-type galaxies (S0-Sab), which have more massive bulges, are, on average, higher than the values for late-type galaxies with the same rotational velocities. We conclude that the black-hole masses are determined primarily by the properties of the bulge and not the rotational velocity or the mass of the galaxy.     

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

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.     

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.     

Accounting for mudflow genesis in preliminary assessment of the maximum volume of solid mudflow sediments in the North Caucasus

The area of investigation was subject to multifactor analysis of the relationship between the maximum amount of mudflow solid sediments (W) and parameters such as the mudflow basin area (S), average channel slope (α), and mudflow channel length (L). They were used to obtain analytical expressions in order to approximate the W(S, L, α) relation based on the mudflow genesis and source height. Statistical data on mudflow manifestations in different basins in the North Caucasus covering more than fifty years were used to obtain the analytical expressions in order to assess the maximum volume of mudflow solid sediments.     

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 morphological characteristics of gully systems and watersheds in Dry-Hot Valley,SW China

Gully systems and watersheds are geomorphic units with clear boundaries that are relatively independent of basin landscapes and play an important role in natural geography. In order to explore the morphological characteristics of gully systems and watersheds in the Dry-Hot Valley [South West (SW) China], gullies are interpreted from online Google images with high resolution and watersheds are extracted from digital elevation model at a scale of 1:50,000. The results show that: (1) There are 17,382 gullies (with a total area of 1141.66 km²) and 42 watersheds in the study area. (2) The average gully density of the study area (D) is 4.29 km/km², gully frequency (F) is 14.39 gullies/km², the branching ratio (B) is 5.13, the length ratio (L) is 3.12, and the coefficient of the main and tributary gullies (M) is 0.06. The degree of gully erosion is strong to extremely strong, the main development intensity of gully erosion ranges from intense to moderate, and the type of gully system is tributary. (3) The watershed areas (A) are between 0.39 and 96.43 km², the relief ratio (R) is from 0.10 to 0.19, the circularity ratio (C) is from 0.30 to 0.83, the texture ratio (T) is from 0.82 to 39.35, and the dominant geomorphological texture type is fine. (4) There is a quantitative relationship between F and D:F?=?0.624D² (R?=0.84) and T is closely related to D, F, M (R²?>?0.7). A, R and C are related to M (R²?>?0.5). The development of gully systems is the result of coupling effects between multiple factors. In this area, the degree of erosion and the condition of the main and tributary gullies can be controlled by the degree of topographic breakage in the watershed, which provides some theoretical basis for the evaluation of gully erosion by the latter. In addition, the scale, relief, and shape have a significant impact on the locations of the main and tributary gullies. For tributary gullies, attention should be paid to the interception and control of runoff and sediment in the small confluence branches in order to prevent gully expansion and head advance. These features can inform the development of targeted measures for the control of soil erosion.