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.     

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.     

Temporal shoreline change and infrastructure influences along the southern Red Sea coast of Saudi Arabia

Jizan is one of the Saudi Arabian coastal cities endowed with diverse natural settings, which includes Ash Shuqayq in the north, Turfah in the centre and Jizan in the South. This work analysed specific environmental characteristics, such as spits, sabkhas and wadis. Assessments used Landsat imagery to examine coastal change between 1973 and 2011. The cumulative temporal change identified regression trends given by coefficients of determination that explained a significant percentage of data variation for Jizan (R ² = 69%) and Turfah (R ² = 72%), while Ash Shuqayq was insignificant (R ² = 14%). Inter-survey results predicted future change, although trends were not significant, i.e. Jizan (R ² = 22%), Turfah (R ² = 14%) and Ash Shuqayq (R ² = 3 and 61% with outlying value removed). Aerial photos showed regional coastal changes, which included a maximum accretion of 36.4 m and maximum erosion of 12.9 m. These are scientifically effective techniques to monitor regional coastal change, i.e. erosion and accretion and identified rates of 0.59, 1.80 and 3.53 myr^?1 for Ash Shuqayq, Turfah and Jizan. Changes were linked to infrastructure developments, e.g. tourism, port development and natural causes, e.g. spit formations and wadi outfalls.     

Determination of atrazine isotherms on agricultural soils

For determination of atrazine isotherms in agricultural soils of Fars Province, composite soil samples from 0 to 5 cm depth with textures of silty clay loam, clay loam and loam were collected. In order to form the atrazine isotherms, 10, 50 and 100 µg atrazine g^?1 soil was added to the soil samples. Soluble atrazine in water:soil ratios of 10:1, 50:1 and 200:1 was measured after 3-h shaking. Finally, for each cases of applied atrazine, water extractable atrazine was determined and quantified using gas chromatography instrument. The results indicated that there was a linear relationship between the logarithms of water extractable atrazine and added atrazine for different water:soil ratios. A general equation of WEA = K(WS) ^α (AA) ^β is obtained experimentally between water extractable atrazine, µg g^?1(WEA), and added atrazine, µg g^?1 (AA), where K, α and β are absorption constants; WS is the water:soil ratio, g g^?1. For the loam, silty clay loam and clay loam soil textures, the α were 0.49, 0.23 and 0.13, respectively, the β were 0.55, 0.806 and 0.21, respectively, and the K were 1.44, 0.78 and 25.38, respectively.     

Detailed comparison of nine intact rock failure criteria using polyaxial intact coal strength data obtained through PFC<Superscript>3D</Superscript> simulations

Study of intact rock failure criteria is an important topic in rock mechanics. In this study, applicability of nine different intact rock failure criteria is investigated for intact coal strength data. PFC^3D modeling was used to simulate the laboratory polyaxial tests for cubic intact coal blocks of side dimension 110 mm under different confining stress combinations. A modified grid search procedure is proposed and used to find the best-fitting parameter values and to calculate the coefficient of determination (R ²) values for each criterion. Detailed comparisons of the nine criteria are made using the following aspects: R ² values, σ ₁ ? σ ₂ plots for different σ ₃, shapes on the deviatoric plane, linearity or nonlinearity on the meridian planes. Through the comparisons of R ² values, σ ₁ ? σ ₂ plots and meridian lines, the modified Wiebols–Cook and modified Lade criteria were found to fit the intact coal strength data best. The nine failure criteria are categorized into three types based on the appearances on the deviatoric plane.     

Comparative study of the absorptive potential of raw and activated carbon <Emphasis Type="Italic">Acacia nilotica</Emphasis> for Reactive Black 5 dye

Acacia nilotica was used for the adsorption of Reactive Black 5 (RB5) dye from an aqueous solution. Both the raw and activated (with H₃PO₄) carbon forms of Acacia nilotica (RAN and ANAC, respectively) were used for comparison. Various parameters (including dye concentration, contact time, temperature, and pH) were optimized to obtain the maximum adsorption capacity. RAN and ANAC were characterized using Fourier transform infrared spectroscopy, scanning electron microscopy, and energy-dispersive X-ray spectroscopy. The maximum experimental adsorption capacities for RAN and ANAC were 34.79 and 41.01 mg g^?1, respectively, which agreed with the maximum adsorption capacities predicted by the Langmuir, Freundlich, and Dubinin–Radushkevich equilibrium isotherm models. The adsorption data of ANAC showed a good fit to the isotherm models based on the coefficient of determination (R ²): Langmuir type II (R ² = 0.99) > Freundlich (R ² = 0.9853) > Dubinin–Radushkevich (R ² = 0.9659). This result suggested monolayer adsorption of RB5 dye. The adsorption of RB5 dye followed pseudo-second-order kinetics. The RAN adsorbent reflected an exothermic reaction (enthalpy change, ΔH = ?0.006 kJ mol^?1) and increased randomness (standard entropy change, ΔS = 0.038 kJ mol^?1) at the solid–solution interface. In contrast, ANAC reflected both exothermic [?0.011 kJ mol^?1 (303–313 K)] and endothermic [0.003 kJ mol^?1 (313–323 K)] reactions. However, the ΔS value of ANAC was lower when the RB5 adsorption increased from 313 to 323 K. The negative values for the Gibbs free energy change at all temperatures indicated that the adsorption of RB5 dye onto RAN and ANAC was spontaneous in the forward direction.     

The effects of petroleum-contaminated soil on photosynthesis of <Emphasis Type="Italic">Amorpha fruticosa</Emphasis> seedlings

A pot experiment was conducted to monitor the dynamic response of photosynthesis of Amorpha fruticosa seedlings to different concentrations of petroleum-contaminated soils from April to September. The results showed that the photosynthetic rates, stomatal conductance and transpiration rate of seedlings significantly decreased in 5–20 g kg^?1 petroleum-contaminated soil during the three given sampling period of July 31 (early), August 30 (mid-term) and September 29 (late). However, the intercellular CO₂ concentration significantly increased in 10 g kg^?1 contaminated soil, while declined in 20 g kg^?1 contaminated soil during the early sampling period as well as in 20 g kg^?1 contaminated soil during the late sampling period. The leaf relative water content of seedlings significantly increased in 20 g kg^?1 contaminated soil during the early sampling period, while it dropped dramatically in 15–20 g kg^?1 contaminated soil during the late sampling period. The contents of chlorophyll a, chlorophyll b and the total chlorophyll of seedlings showed a sharp decline during the three sampling periods in contaminated soil. Comprehensively, considering the negative effects of petroleum on the photosynthesis, growth performance and remediation effect on petroleum of A. fruticosa seedlings, this plant was tolerant of petroleum-contaminated soil and was potentially useful for the phytoremediation of petroleum-contaminated sites in northern Shaanxi, China.     

Geospatial approach in mapping soil erodibility using CartoDEM – A case study in hilly watershed of Lower Himalayan Range

Soil erodibility is one of the most important factors used in spatial soil erosion risk assessment. Soil information derived from soil map is used to generate soil erodibility factor map. Soil maps are not available at appropriate scale. In general, soil maps at small scale are used in deriving soil erodibility map that largely generalized spatial variability and it largely ignores the spatial variability since soil map units are discrete polygons. The present study was attempted to generate soil erodibilty map using terrain indices derived from DTM and surface soil sample data. Soil variability in the hilly landscape is largely controlled by topography represented by DTM. The CartoDEM (30 m) was used to derive terrain indices such as terrain wetness index (TWI), stream power index (SPI), sediment transport index (STI) and slope parameters. A total of 95 surface soil samples were collected to compute soil erodibility factor (K) values. The K values ranged from 0.23 to 0.81 t ha^?1R^?1 in the watershed. Correlation analysis among K-factor and terrain parameters showed highest correlation of soil erodibilty with TWI (r ²= 0.561) followed by slope (r ²= 0.33). A multiple linear regression model was developed to derive soil erodibilty using terrain parameters. A set of 20 soil sample points were used to assess the accuracy of the model. The coefficient of determination (r ²) and RMSE were computed to be 0.76 and 0.07 t ha^?1R^?1 respectively. The proposed methodology is quite useful in generating soil erodibilty factor map using digital elevation model (DEM) for any hilly terrain areas. The equation/model need to be established for the particular hilly terrain under the study. The developed model was used to generate spatial soil erodibility factor (K) map of the watershed in the lower Himalayan range.     

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

The human impact on the transformation of juniper forest landscape in the western part of the Pamir-Alay range (Tajikistan)

Detailed analyses were conducted of human impact on juniper forest landscapes occurring within the Zarafshan Range (Pamir-Alay). Juniperus seravschanica and J. semiglobosa belong to forest-forming species in Central Asia. At present, juniper forests all over Tajikistan are seriously threatened as a result of excessive logging and cattle grazing. The aim of this paper is to present juniper forest transformation as a result of human activities as well as the diversity of soil properties in the organic and humus horizons in the altitudinal system of soil zonation. Three groups of phytocoenoses were distinguished: those with a dominant share of Juniperus seravschanica; those with a dominant share of J. semiglobosa; and mixed. Associations with Juniperus seravschanica and J. semiglobosa feature several variants of phytocoenoses with dominant species: Artemisia lehmanniana, A. dracunculus, Eremurus olgae, Festuca sulcata, Ligularia thomsonii, Stipa turkestanica, Thymus seravschanicus, and Ziziphora pamiroalaica. The collected soil samples differ in their granulometric composition. Gravelly cobble fractions >2 mm are dominant; the share of sandy particles <2 mm is much lower (about 10–20%). Fraction 0.5–0.05 attains 35% on average. The Corg content of the soil varied from 0.26 to 11.40% in the humus horizon (A) and from 4.3 to 25% in the organic (O). Similar relationships were reported in the case of Ntot concentration. A clear relationship can be observed between concentrations of Corg and Ntot. Soil pH varied, ranging from very low acidic (pH 5.5) to neutral (pH 8.5). The content of available P varied; high concentrations were noted in organic (O) (40.46–211 mg kg^?1) and mixed horizons (OA) (2.61–119 mg kg^?1). Maximum accumulations of Pavail (1739.6 mg kg^?1) and Ptot (9696 mg kg^?1) were observed at a site heavily affected by intense grazing. Concentrations of Mgavail varied from 116 to 964 mg kg^?1. Most of the analysed soil profiles lacked an organic horizon; only thin humus occurred.     

Dynamics of carbon fluxes with responses to vegetation,meteorological and terrain factors in the south-eastern Tibetan Plateau

Tibetan Plateau (TP) is the highest and most extensive plateau in the world and has been known as the roof of the world, and it is sensitive to climate change. The researches of CO₂ fluxes (F _C) in the TP region play a significant role in understanding regional and global carbon balance and climate change. Eddy covariance flux measurements were conducted at three sites of south-eastern TP comprising Dali (DL, cropland ecosystem), LinZhi (LZ, alpine meadow ecosystem) and Wenjiang (WJ, cropland ecosystem); amongst those DL and LZ are located in plateau region, while WJ is in plain region. Dynamics of F _C and influences of vegetation, meteorological (air temperature, photosynthetically active radiation, soil temperature and soil water content) and terrain factors (altitude) were analysed on the basis of data taken during 2008. The results showed that, in the cool sub-season (March, April, October and December), carbon sink appeared even in December with fluxes of (?0.021 to ?0.05) mg CO₂ m^?2 s^?1 and carbon source only in October (0.03 ± 0.0048) mg CO₂ m^?2 s^?1 in DL and WJ site. In LZ site, carbon sink was observed in April: (?0.036 ± 0.0023) mg CO₂ m^?2 s^?1 and carbon sources in December and March (0.008–0.010 mg CO₂ m^?2 s^?1). In the hot sub-season (May–August), carbon source was observed only in May with (0.011 ± 0.0022), (0.104 ± 0.0029) and (0.036 ± 0.0017) fluxes in LZ, DL and WJ site, respectively, while carbon sinks with (?0.021 ± 0.0041), (?0.213 ± 0.0007) and (?0.110 ± 0.0015) mg CO₂ m^?2 s^?1 fluxes in LZ, DL, and WJ, respectively. Comparing with plain region (WJ), carbon sinks in plateau region (DL and LZ) lasted for a longer time, and the absorption sum was large and up to (–357.718 ± 0.0054) and (?371.111 ± 0.0039) g C m^?2 year^?1, respectively. The LZ site had the weakest carbon sink with (?178.547 ± 0.0070) g C m^?2 year^?1. Multivariate analysis of covariance showed that altitude (AL) as an independent factor explained 39.5 % of F _C (P < 0.026). F _C had a quadratic relationship with Normalized difference vegetation index (NDVI) (R ² ranges from 0.485 to 0.640 for three sites), an exponential relationship with soil temperature at 5-cm depth (ST ₅) at night time and a quadratic relationship with air temperature (T _a) at day time. Path analysis indicated that photosynthetically active radiation (PAR), sensible heat fluxes (H) and other factors all had direct or indirect effects on F _C in all of the three tested sites around the south-eastern TP.     

A modified normalized model for predicting effective soil thermal conductivity

Effective soil thermal conductivity (λ _eff) describes the ability of a multiphase soil to transmit heat by conduction under unit temperature gradient. It is a critical parameter for environmental science, earth and planetary science, and engineering applications. Numerous models are available in the literature, but their applicability is generally restricted to certain soil types or water contents (θ). The objective of this study was to develop a new model in the similar form of the Johansen 1975 model to simulate the λ _eff(θ) relationship of soils of various soil textures and water contents. An exponential type model with two parameters is developed and a new function for calculating dry soil thermal conductivity is presented. Performance of the new model and six other normalized models were evaluated with published datasets. The results show that the new model is able to well mimic λ _eff(θ) relationship of soils from sand to silt loam and from oven dry to full saturation. In addition, it has the best performance among the seven models under test (with root-mean-square error of 0.059 W m^?1 °C^?1, average deviations of 0.0009 W m^?1 °C^?1, and Nash–Sutcliffe efficiency of 0.994). The new model has potential to improve the reliability of soil thermal conductivity estimation and be incorporated into numerical modeling for environmental, earth and engineering studies.     

Impact of rainfall intensity on the hydrological performance of erosion control geotextiles

The impact of erosion control geotextiles on the surface runoff from slopes is quite variable and depends strongly on site-specific conditions (soil characteristics, slope morphology, climate, etc.), as has been shown in several earlier studies. In addition, little is known about the proportion of runoff reduction that is caused by the geotextile and the proportion that is caused by soil characteristics. To shed more light on this issue, an experiment was carried out to test the impact of 500 g m^?2 jute nets (J500) and 400 g m^?2; 700 g m^?2 coir nets (C400, C700) on the surface runoff from simulated rainfall of four different intensities (I ₁ = 18.7; I ₂ = 27.2; I ₃ = 53.6; I ₄ = 90.5 mm h^?1). Data on runoff volume, peak discharge and time to peak discharge were collected from 40 simulated rainfall events. An impermeable “no-soil” subgrade was used to examine the impact of the geotextile on runoff without any influence of soil. All tested geotextiles significantly reduced runoff (volume, peak discharge) at all rainfall intensities, with the exception of C400 and C700 during simulated rainfall intensity I ₄. J500 seemed to have the most effective runoff reduction performance at all rainfall intensities. In general, as the rainfall intensity increased, the effectiveness of the geotextiles decreased. Interesting behaviour was observed for J500 under simulated rainfall intensity I ₄—the effectiveness of the geotextile increased with the duration of the rainfall.     

Behavior of 10-Å phase at low temperatures

The thermal evolution of 10-Å phase Mg₃Si₄O₁₀(OH)₂·H₂O, a phyllosilicate which may have an important role in the storage/release of water in subducting slabs, was studied by X-ray single-crystal diffraction in the temperature range 116–293 K. The lattice parameters were measured at several intervals both on cooling and heating. The structural model was refined with intensity data collected at 116 K and compared to the model refined at room temperature. As expected for a layer silicate on cooling in this temperature range, the a and b lattice parameters undergo a small linear decrease, α _a  = 1.7(4) 10^?6 K^?1 and α _b  = 1.9(4) 10^?6 K^?1, where α is the linear thermal expansion coefficient. The greater variation is along the c axis and can be modeled with the second order polynomial c _T  = c ₂₉₃(1 + 6.7(4)10^?5 K^?1ΔT + 9.5(2.5)10^?8 K^?2(ΔT)²) where ΔT = T ? 293 K; the monoclinic angle β slightly increased. The cell volume thermal expansion can be modeled with the polynomial V _T  = V ₂₉₃ (1 + 8.0 10^?5 K^?1 ΔT + 1.4 10^?7 K^?2 (ΔT)²) where ΔT = T ? 293 is in K and V in ų. These variations were similar to those expected for a pressure increase, indicating that T and P effects are approximately inverse. The least-squares refinement with intensity data measured at 116 K shows that the volume of the SiO₄ tetrahedra does not change significantly, whereas the volume of the Mg octahedra slightly decreases. To adjust for the increased misfit between the tetrahedral and octahedral sheets, the tetrahedral rotation angle α changes from 0.58° to 1.38°, increasing the ditrigonalization of the silicate sheet. This deformation has implications on the H-bonds between the water molecule and the basal oxygen atoms. Furthermore, the highly anisotropic thermal ellipsoid of the H₂O oxygen indicates positional disorder, similar to the disorder observed at room temperature. The low-temperature results support the hypothesis that the disorder is static. It can be modeled with a splitting of the interlayer oxygen site with a statistical distribution of the H₂O molecules into two positions, 0.6 Å apart. The resulting shortest O_bas–O_W distances are 2.97 Å, with a significant shortening with respect to the value at room temperature. The low-temperature behavior of the H-bond system is consistent with that hypothesized at high pressure on the basis of the Raman spectra evolution with P.     

Thermal equation of state of natural tourmaline at high pressure and temperature

Synchrotron-based in situ angle-dispersive X-ray diffraction experiments were conducted on a natural uvite-dominated tourmaline sample by using an external-heating diamond anvil cell at simultaneously high pressures and temperatures up to 18 GPa and 723 K, respectively. The angle-dispersive X-ray diffraction data reveal no indication of a structural phase transition over the P–T range of the current experiment in this study. The pressure–volume–temperature data were fitted by the high-temperature Birch–Murnaghan equation of state. Isothermal bulk modulus of K ₀ = 96.6 (9) GPa, pressure derivative of the bulk modulus of \(K_{0}^{\prime } = 12.5 \;(4)\), thermal expansion coefficient of α ₀ = 4.39 (27) × 10^?5 K^?1 and temperature derivative of the bulk modulus (?K/?T) _P  = ?0.009 (6) GPa K^?1 were obtained. The axial thermoelastic properties were also obtained with K _a0 = 139 (2) GPa, \(K_{a0}^{\prime }\) = 11.5 (7) and α _a0 = 1.00 (11) × 10^?5 K^?1 for the a-axis, and K _c0 = 59 (1) GPa, \(K_{c0}^{\prime }\) = 11.4 (5) and α _c0 = 2.41 (24) × 10^?5 K^?1 for the c-axis. Both of axial compression and thermal expansion exhibit large anisotropic behavior. Thermoelastic parameters of tourmaline in this study were also compared with that of the other two ring silicates of beryl and cordierite.     

Adsorption study on municipal solid waste leachate using <Emphasis Type="Italic">Moringa oleifera</Emphasis> seed

Effects of initial concentrations of Moringa oleifera seed coagulant for removing Chemical Oxygen Demand and Total Dissolved Solids from municipal solid waste leachate have been evaluated at an optimum pH of 7 and temperature of 318 K. The kinetic data obtained from the experiments were fitted to the pseudo first-order, pseudo second-order, Elovich and intraparticle diffusion models. Based on a regression coefficient (R ²), the equilibrium (kinetic) data were best fitted with the Elovich model (R ² = 0.993 for Chemical Oxygen Demand and R ² = 0.996 for Total Dissolved Solids) than that of other models. The results of the kinetic models study indicated that the adsorption capacity of M. oleifera seed as a coagulant for removing Chemical Oxygen Demand and Total Dissolved Solids in a leachate increased up to 100 mg L^?1, beyond which the adsorption capacity got reduced. Finally, the present study concluded that M. oleifera seed coagulant could be employed effectively for the removal of Chemical Oxygen Demand and Total Dissolved Solids in a municipal solid waste leachate.     

Short- and Long-Term Vegetative Propagation of Two <Emphasis Type="Italic">Spartina</Emphasis> Species on a Salt Marsh in Southern Brazil

Spartina alterniflora and Spartina densiflora are native salt marsh plants from the Atlantic coast; their habitats in Patos Lagoon estuary (southern Brazil) are characterized by a microtidal regime (<0.5 m) and, during El Niño events, high estuarine water levels and prolonged flooding due to elevated freshwater discharge from a 200,000-km² watershed. During and between El Niño events, the vegetative propagation of these two Spartina species in the largest estuary of southern Brazil (Patos Lagoon) was evaluated by monitoring transplanted plants for 10 years (short-term study) and interpreting aerial photos of natural stands for 56 years (long-term study). During the short-term study, S. alterniflora quickly occupied mud flats (up to 208 cm year^?1) by elongation of rhizomes, whereas S. densiflora showed a modest lateral spread (up to 13 cm year^?1) and generated dense circular-shaped stands. However, moderate and strong El Niño events can promote excessive flooding and positive anomalies in the estuarine water level that reduce the lateral spread and competitive ability of S. densiflora. During the long-term study, natural stands of S. alterniflora and S. densiflora had steady lateral spread rates of 152 and 5.2 cm year^?1, respectively, over mud flats. In the microtidal marshes of the southwest Atlantic, the continuous long-term lateral expansion of both Spartina species embodies periods of intense flooding stress (moderate and strong El Niños), when there is a decrease of vegetative propagation and less stressful low water periods of fast spread over mud flats (non-El Niño periods and weak intensity El Niños).     

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.     

Spatial information of soil hydraulic conductivity and performance of cokriging over kriging in a semi-arid basin scale

Unlike the studies in small parcels by systematic measurements, the spatial variability of soil properties is expected to increase in those over relatively large areas or scales. Spatial variability of soil hydraulic conductivity (K _h) is of significance for the environmental processes, such as soil erosion, plant growth, transport of the plant nutrients in a soil profile and ground water levels. However, its variability is not much and sufficiently known at basin scale. A study of testing the performance of cokriging of K _h compared with that of kriging was conducted in the catchment area of Sarayköy II Irrigation Dam in Cank?r?, Turkey. A total of 300 soil surface samples (0–10 cm) were collected from the catchment with irregular intervals. Of the selected soil properties, because the water-stable aggregates (WSA) indicated the highest relationship with the hydraulic conductivity by the Pearson correlation analysis, it is used as an auxiliary variable to predict K _h by the cokriging procedure. In addition, the sampling density was reduced randomly to n = 175, n = 150, n = 75 and n = 50 for K _h to determine if the superiority of cokriging over kriging would exist. Statistically, the results showed that all reduced K _h was as good as the complete K _h when its auxiliary relations with WSA were used in cokriging. Particularly, the results of the “Relative Reduction in MSE” (RMSE) revealed that the reduced data set of n = 75 produced the most accurate map than the others. In this basin-scaled study, there was a clear superiority of the cokriging procedure by the reduction in data although a very undulating topography and topographically different aspects, two different land uses with non-uniform vegetation density, different parent materials and soil textures were present in the area. Hence, using the statistically significant auxiliary relationship between K _h and WSA might bring about a very useful data set for watershed hydrological researches.