Parameter estimation of SWAT and quantification of consequent confidence bands of model simulations

Soil and Water Assessment Tool (SWAT) is a river basin scale model widely used to study the impact of land management practices in large, complex watersheds. Even though model output uncertainties are generally recognized to affect watershed management decisions, those uncertainties are largely ignored in model applications. The uncertainties of SWAT simulations are quantified using various methods, but simultaneous attempt to calibrate a model so as to reduce the uncertainty are seldom done. This study aims to use an uncertainty reduction procedure that helps calibrate the SWAT model. The shuffled complex evolutionary metropolis algorithm for uncertainty analysis is employed for this purpose, and is demonstrated using the data from the St. Joseph River basin, USA. The values of the performance indices, the r² and the Nash–Sutcliffe efficiency (NSE) for the simulations during calibration period was found to be 0.81 (same for r² and NSE) and 0.79 for validation period indicating a good simulation by the model. The results also indicate that the algorithm helps reduce the uncertainty (percentage of coverage?=?62% and average width?=?19.2 m³/s), and also identifies the plausible range of parameters that simulate the processes with less uncertainty. The confidence bands of simulations are obtained that can be employed in making uncertainty-based decisions on watershed management practices.     

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

Evaluation of the SWAT model for water balance study of a mountainous snowfed river basin of Nepal

In this study, a semi-distributed hydrologic model Soil and Water Assessment Tool (SWAT) has been employed for the Karnali River basin, Nepal to test its applicability for hydrological simulation. Further, model was evaluated to carry out the water balance study of the basin and to determine the snowmelt contribution in the river flow. Snowmelt Runoff Model (SRM) was also used to compare the snowmelt runoff simulated from the SWAT model. The statistical results show that performance of the SWAT model in the Karnali River basin is quite good (p-factor = 0.88 and 0.88, for daily calibration and validation, respectively; r-factor = 0.76 and 0.71, for daily calibration and validation, respectively). Baseflow alpha factor (ALPHA_BF) was found most sensitive parameter for the flow simulation. The study revealed that the average annual runoff volume available at the basin outlet is about 47.16 billion cubic metre out of which about 12% of runoff volume is contributed by the snowmelt runoff. About 25% of annual precipitation seems to be lost as evapotranspiration. The results revealed that both the models, SWAT and SRM, can be efficiently applied in the mountainous river basins of Nepal for planning and management of water resources.     

Satellite-based rainfall estimation and discharge measurement of Middle Indus River,Pakistan

The impacts of floods and droughts are intensified by climate change, lack of preparedness, and coordination. The average rainfall in study area is ranging from 200 to 400 mm per year. Rain gauge generally provides very accurate measurement of point rain rates and the amounts of rainfall but due to scarcity of the gauge locations provides very general information of the area on regional scale. Recognizing these practical limitations, it is essential to use remote sensing techniques for measuring the quantity of rainfall in the Middle Indus. In this research, Tropical Rainfall Measuring Mission (TRMM) estimation can be used as a proxy for the magnitude of rainfall estimates from classical methods (rain gauge), quantity, and its spatial distribution for Middle Indus river basin. In order to use TRMM satellite data for discharge measurement, its accuracy is determined by statistically comparing it with in situ gauged data on daily and monthly bases. The daily R ² value (0.42) is significantly lower than monthly R ² value (0.82), probably due to the time of summation of TRMM 3-hourly precipitation data into daily estimates. Daily TRMM data from 2003 to 2012 was used as input forcing in Soil and Water Assessment Tool (SWAT) hydrological model along with other input parameters. The calibration and validation results of SWAT model give R ² = 0.72 and 0.73 and Nash-Sutcliffe coefficient of efficiency = 0.69 and 0.65, respectively. Daily and monthly comparison graphs are generated on the basis of model discharge output and observed data.     

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

SWAT manual calibration and parameters sensitivity analysis in a semi-arid watershed in North-western Morocco

Being a laborious approach, manual calibration of hydrologic model in a semi-arid context requires in-depth knowledge of the watershed and as much as possible field input data to obtain reliable simulations. In this study, manual calibration and relative sensitivity analysis approaches of the SWAT model (Soil and Water Assessment Tool) were applied for water balance in a 1993 km² watershed (on the R’dom river) located in North-western Morocco. The watershed is located in a semi-arid area dominated by agro-forestry activities. The objectives of this study were (i) to perform a local sensitivity analysis of the SWAT model taking into consideration the watershed characteristics and (ii) to implement a detailed methodology of manual calibration and validation of the model in a semi-arid context. Sensitivity analysis has been carried out on 12 different SWAT input parameters, and has revealed that 4 input parameters only were the most influential ones on flow components of the R’dom watershed. Model manual calibration was conducted along 2006 and 2007 by comparing measured and predicted monthly and daily discharges and taking Nash-Sutcliffe coefficient (NSE), determination coefficient (R ²), and percent bias (PBIAS) as goodness-of-fit indicators. Validation has been performed by the same approach through 2008 and 2009 period. All final NSE values were above 0.5, R ² values exceeded 0.7, and PBIAS lower than 25% demonstrating satisfactory model performances over the study watershed conditions. The SWAT model set-up with measured input data, manually calibrated and validated, reflects well the real hydrologic processes occurring in the R’dom watershed and can be used to assess current and future conditions and to evaluate alternative management practices.     

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

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.     

Is NGC 7006 a globular cluster with a primarily red horizontal branch?

A CCD BV R photometric study of the central region (15″ ≤ r ≤ 100″) of the globular cluster NGC 7006 based on color-magnitude diagrams is presented. We find for the main parameters of the cluster [Fe/H] = ?1.62, Y = 0.21, E _B?V = 0.15 ^m , V _HB = 18.84 ^m , M _V ^HB =+0.56 ^m , R _⊙ = 37.1 kpc). Two previously unknown RR Lyr variables were discovered in the central region of the cluster. The morphological index of the horizontal branch for the entire region studied indicates that the red stellar population dominates, consistent with previous studies: HB _mi = ?0.13. Such anomalously negative morphological indices are possessed by a whole group of Ool clusters with intermediate metallicities, which also display a characteristic distribution of stars along the horizontal branch. There is a radial dependence for the horizontal-branch morphology, with the color becoming primarily blue with approach toward the cluster center. One possible origin for this behavior could be the effect of inner dynamical processes on the spatial distribution of hot stars.     

Early Eocene planktonic and benthic foraminifera from the Khangiran formation (northeast of Iran)

The Kopet-Dagh sedimentary basin is located in northeast of Iran and southeast of Turkmenistan in the Middle East. The Khangiran formation represents the last marine deposition in the Kopet-Dagh sedimentary basin. The early Eocene planktonic and benthic foraminifera from the Khangiran formation which belongs to the lower 376 m thickness of this formation are identified and biostratigraphically evaluated. Due to rarity of Morozovella species, planktonic foraminiferal zonation was difficult to determine in this formation. The determination of upper part of the late Paleocene for the lower 124 m of this formation is according to the lowest occurrence of Acarinina sibaiyaensis species. From E₅ to near middle of E₇ biozone, increasing trend of Acarinina frequency and a peak in Morozovella species and decreasing trend of frequency of Subbotina and Pseudohastigerina species indicate the warm and oligotrophic condition of the seawater during sedimentation of the studied interval. Toward the Ypresian-Lutetian boundary, the increasing trend of Pseudohastigerina and Subbotina species and decreasing in frequency of Acarinina species suggest the low-oxygen level, eutrophic and intermediate condition of the seawater. In this formation, the high abundance of the epifaunal taxa such as Anomalinoides spp., Cibicidoides spp., Gyroidinoides spp., and Lenticulina spp. from the base (late Paleocene sediments) up to E₆ biozone reflects oligotrophic and oxic shallow water conditions. The occurrence of several peaks in abundance of Bulimina and Uvigerina species at the Ypresian/Lutetian boundary suggests eutrophic condition. These paleoecological conditions could be correlated with fluctuations in the numbers of the planktonic foraminifera.     

Modeling and mapping of solar radiation using geostatistical analysis methods in Iran

Because of economic and technical limitations, measuring solar energy received at ground level (R _s ) isn’t possible in all parts of the country, and in only 12% of synoptic stations is this parameter measured and recorded. Thus, it should be estimated and modeled spatially based on other climatic variables using mathematical methods. In this research, many attempts have been made to introduce an air temperature-based model for Rs estimation, and then, based on the output of the mentioned models, several geostatistical methods have been tested, and finally an elegant spatial model is proposed for (Rs) zoning in Iran. In this regard, the relationships between the measured amounts of monthly solar radiation and other climatic parameters, such as a monthly average, maximum and minimum temperature, precipitation, relative humidity, and the number of sunny hours during the period 1970–2010, are examined and modeled. It was revealed that based on the linear relationship between the monthly average air temperatures and solar radiation values recorded in each of the stations, that the best-fit linear model, with R ²  = 0.822, MAE = 1.81, RMSE = 2.51%, and MAPE = 10.08, can be introduced for Rs estimation. Then, using the outputs of the proposed model, the amounts of (R _s ) are estimated in another 171 meteorological stations (a total of 192 stations), and eight geostatistical methods (IDW, GPI, RBF, LPI, OK, SK, UK, and EBK) were investigated for zoning. Comparing the resulting variograms showed that in addition to proof of spatial correlation between solar radiation data, they can be applied for modeling changes in various directions. Analyzing the ratio of the nugget effect on the roof of the variograms showed that the Gaussian model with the lowest ratio (Co/Co + C = 0.883) and (R ²  = 0.972), could model the highest correlation between the data and, therefore, it was used for data interpolation. To select the best geostatistical model, R2, MAE, and RMSE were used. On this basis, it was found that the RBF method with R ²  = 0.904, MAE = 3.02, RMSE = 0.39% is the most effective. Also, the IDW method with R ²  = 0.90, MAE = 3.08, RMSE = 0.391%, compared to other methods is the most effective. In addition, for data validation, correlations between observed and estimated values of solar radiation were studied and found R ²  = 0.86.     

Prediction of strength parameters of sedimentary rocks using artificial neural networks and regression analysis

Accurate laboratory measurement of geo-engineering properties of intact rock including uniaxial compressive strength (UCS) and modulus of elasticity (E) involves high costs and a substantial amount of time. For this reason, it is of great necessity to develop some relationships and models for estimating these parameters in rock engineering. The present study was conducted to forecast UCS and E in the sedimentary rocks using artificial neural networks (ANNs) and multivariable regression analysis (MLR). For this purpose, a total of 196 rock samples from four rock types (i.e., sandstone, conglomerate, limestone, and marl) were cored and subjected to comprehensive laboratory tests. To develop the predictive models, physical properties of studied rocks such as P wave velocity (V_p), dry density (γ_d), porosity, and water absorption (A_b) were considered as model inputs, while UCS and E were the output parameters. We evaluated the performance of MLR and ANN models by calculating correlation coefficient (R), mean absolute error (MAE), and root-mean-square error (RMSE) indices. The comparison of the obtained results revealed that ANN outperforms MLR when predicting the UCS and E.     

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

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

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.     

Analytical methods for measuring the parameters of interstellar gas using methanol observations

The excitation of methanol in the absence of external radiation is analyzed, and LTE methods for probing interstellar gas considered. It is shown that rotation diagrams correctly estimate the gas kinetic temperature only if they are constructed using lines whose upper levels are located in the same K-ladders, such as the J₀?J_?1E lines at 157 GHz, the J₁?J₀E lines at 165 GHz, and the J₂?J₁E lines at 25 GHz. The gas density must be no less than 10⁷ cm^?3. Rotation diagrams constructed from lines with different K values for their upper levels (e.g., 2_K?1_K at 96 GHz, 3_K?2_K at 145 GHz, 5_K?4_K at 241 GHz) significantly underestimate the temperature, but enable estimation of the density. In addition, diagrams based on the 2_K?1_K lines can be used to estimate the methanol column density within a factor of about two to five. It is suggested that rotation diagrams should be used in the following manner. First, two rotation diagrams should be constructed, one from the lines at 96, 145, or 241 GHz, and another from the lines at 157, 165, or 25 GHz. The former diagram is used to estimate the gas density. If the density is about 10⁷ cm^?3 or higher, the latter diagram reproduces the temperature fairly well. If the density is around 10⁶ cm^?3, the temperature obtained from the latter diagram should be multiplied by a factor of 1.5–2. If the density is about 10⁵ cm^?3 or lower, then the latter diagram yields a temperature that is lower than the kinetic temperature by a factor of three or more, and should be used only as a lower limit for the kinetic temperature. The errors in the methanol column density determined from the integrated intensity of a single line can be more than an order of magnitude, even when the gas temperature is well known. However, if the J₀?(J ? 1)₀E lines, as well as the J₁?(J ? 1)₁A⁺ or A^? lines are used, the relative error in the column density is no more than a factor of a few.     

Vertical stability and the Brunt-Väisäla frequency of deep natural waters by the example of Lake Baikal,Lake Tanganyika,and the World Ocean

Theoretical analysis, calculations, and comparison with the results of observations in Lake Baikal, Lake Tanganyika, and the World Ocean are performed for the vertical stability E and the Brunt-Väisäla frequency N in the form of N ² with regard to all components (at the constant temperature T and the salinity S, the common adiabatic form at T, S Const). The adiabatic stability E _ad and the Väisäla frequency N in the form of N _ad ² are always positive; at a change from the inverse to the direct temperature stratification, they have deep minimums reaching 10^?16 m^?1 and 10^?15 s^?2 and less; the minimums have the form of a special point, a reversal point of the first kind called a “cusp.” The reality of these reversal points is confirmed by the analysis of the investigation procedure, comparison with the results of previous theoretical (Sherstyankin, et al., 2007), and experimental (observations in Baikal, Shimaraev et al., 1994) works. The features of vertical profiles of E _ad , E andN _ad ² , N ², as well as the layers where the Brunt-Väisäla frequency is less than the inertial frequency, are studied. The analysis with regard to all components of the stability E _ad and the Brunt-Väisäla frequency N makes a great contribution to understanding of mixing processes in theoretical and experimental investigations; it is valid in all reservoirs of the Earth with inverse and direct temperature stratification, including Lake Baikal, Lake Tanganyika, and the World Ocean.     

Seyfert nuclei as sources of ultrahigh-energy cosmic rays

Particles can be accelerated to ultrahigh energies E≈10²¹ eV in moderate Seyfert nuclei. This acceleration occurs in shock fronts in relativistic jets. The maximum energy and chemical composition of the accelerated particles depend on the magnetic field in the jet, which is not well known; fields in the range ～5–1000 G are considered in the model. The highest energies of E≈10²¹ eV are acquired by Fe nuclei when the field in the jet is B≈16 G. When B～(5–40) G, nuclei with Z<10 are accelerated to E≤10²⁰ eV, while nuclei with Z≥10 acquire energies E≥2×10²⁰ eV. Only particles with Z≥23 acquire energies E≤10²⁰ eV when B～1000 G. Protons are accelerated to E<4×10¹⁹ eV, and do not fall into the range of energies of interest for any magnetic field B. The particles lose a negligible amount of their energy in interactions with infrared photons in the accretion disk; losses in the thick gas-dust torus are also negligible if the luminosity of the galaxy is L≤10⁴⁶ erg/s and the angle between the normal to the galactic plane and the line of sight is sufficiently small, i.e., if the axial ratio of the galactic disk is comparatively high. The particles do not lose energy to curvature radiation if their deviations from the jet axis do not exceed 0.03–0.04 pc at distances from the center of R≈40–50 pc. Synchrotron losses are small, since the magnetic field frozen in the galactic wind at R≤40–50 pc is directed (as in the jet) primarily in the direction of motion. If the model considered is valid, the detected cosmic-ray protons could be either fragments of Seyfert nuclei or be accelerated in other sources. The jet magnetic fields can be estimated both from direct astronomical observations and from the energy spectrum and chemical composition of cosmic rays.     

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.     

The Influence of Benthic Macrofauna on the Erodibility of Intertidal Sediments with Varying mud Content in Three New Zealand Estuaries

Fine sediment inputs can alter estuarine ecosystem structure and function. However, natural variations in the processes that regulate sediment transport make it difficult to predict their fate. In this study, sediments were sampled at different times (2011–2012) from 45 points across intertidal sandflat transects in three New Zealand estuaries (Whitford, Whangamata, and Kawhia) encompassing a wide range in mud (≤63 μm) content (0–56 %) and macrofaunal community structure. Using a core-based erosion measurement device (EROMES), we calculated three distinct measures of sediment erosion potential: erosion threshold (? _c ; N m^?2), erosion rate (ER; g m^?2 s^?1), and change in erosion rate with increasing bed shear stress (m _e ; g N^?1 s^?1). Collectively, these measures characterized surface (? _c and ER) and sub-surface (m _e ) erosion. Benthic macrofauna were grouped by functional traits (size and motility) and data pooled across estuaries to determine relationships between abiotic (mud content, mean grain size) and biotic (benthic macrofauna, microbial biomass) variables and erosion measures. Results indicated that small bioturbating macrofauna (predominantly freely motile species <5 mm in size) destabilized surface sediments, explaining 23 % of the variation in ? _c (p ≤ 0.01) and 59 % of the variation in ER (p ≤ 0.01). Alternatively, mud content and mean grain size cumulatively explained 61 % of the variation in m _e (p ≤ 0.01), where increasing mud and grain size stabilized sub-surface sediments. These results highlight that the importance of biotic and abiotic predictors vary with erosion stage and that functional group classifications are a useful way to determine the impact of benthic macrofauna on sediment erodibility across communities with different species composition.