In this study, daily averages of air quality parameters were measured in two stations (S1 and S2) of the organized industrial district in Samsun. The meteorological variables were measured at only one station (S1), such as temperature, relative humidity, wind speed, solar radiation, and ambient pressure in 2007, and the daily promised limit for nitrogen dioxide has been especially exceeded at 206 times for 1^st station. However, exceeds of the limit value in 2006 for 1^st station was reduced by approximately 3.5 times. The daily nitrogen dioxide concentration did not exceed the daily limit of WHO^[1] as for 2^nd station. The results obtained showed that under the influence of dominant wind direction, the second station measurement results are higher than that of the first station. To determine all of the possible environmental effects, the measurements should be analyzed from a multi-point perspective.

Richard’s equation was approximated by finite-difference numerical scheme to model water infiltration profile in variably unsaturated soil^[1]. The published data of Philip’s semi-analytical solution was used to validate the simulated results from the numerical scheme. A discrepancy was found between the simulated and the published semi-analytical results. Morris method as a global sensitivity tool was used as an alternative to local sensitivity analysis to assess the results discrepancy. Morris method with different sampling strategies were tested, of which Manhattan distance method has resulted a better sensitivity measures and also a better scan of input space than Euclidean method. Moreover, Morris method at p = 2 , r = 2 and Manhattan distance sampling strategy, with only 2 extra simulation runs than local sensitivity analysis, was able to produce reliable sensitivity measures (μ^*, σ). The sensitivity analysis results were cross-validated by Sobol’ variance-based method with 150,000 simulation runs. The global sensitivity tool has identified three important parameters, of which spatial discretization size was the sole reason of the discrepancy observed. In addition, a high proportion of total output variance contributed by parameters β and θ_s is suggesting a greater significant digits to reduce its input uncertainty range.

We studied Zeta potentials of nanoparticles titanium dioxides （nTiO₂） in different concentration of NaNO₃ and phosphate (P) solutions. In addition, the effect of flow rate on the transport of nTiO₂ in P was investigated at pH＝6.5. Experimental results show that the Zeta potential of nTiO₂ is compressed with the increasing ion concentration (IC) of NaNO₃ at pH＝6.5. The negative charge increases with the augment of P. Therefore, the high P and low NaNO3 induce the stabilization of nTiO₂ aggregates. The transport experiments suggest that the rapid flow rate is favorable for the transportability of nTiO₂ and soluble phosphate. The breakthrough transport curves (BTCs) of nTiO₂ in sand columns can be fitted well with two-site kinetic attachment model. The modeling results suggest that the values of first-order attachment rate coefficients (k₂) and detachment rate coefficients (k_2d) on site 2 and first-order attachment rate coefficients (k₁) on site 1 are responsible to the attaching efficiency of nTiO₂ on sands and their transportability.

Every year, hundreds of fires occur in the forests and rangelands across the world and damage thousands hectare of trees, shrubs, and plants which cause environmental and economic damages. This study aims to establish a real time forest fire alert system for better forest management and monitoring in Golestan Province. In this study, in order to prepare fire hazard maps, the required layers were produced based on fire data in Golestan forests and MODIS sensor data. At first, the natural fire data was divided into two categories of training and test samples randomly. Then, the vegetation moisture stresses and greenness were considered using six indexes of NDVI, MSI, WDVI, OSAVI, GVMI and NDWI in natural fire area of training category on the day before fire occurrence and a long period of 15 years, and the risk threshold of the parameters was considered in addition to selecting the best spectral index of vegetation. Finally, the model output was validated for fire occurrences of the test category. The results showed the possibility of prediction of fire site before occurrence of fire with more than 80 percent accuracy.

In this study, the development of rinnenkarren systems is analyzed. During the field studies, 36 rinnenkarren systems were investigated. The width and depth were measured at every 10 cm on the main channels and then shape was calculated to these places (the quotient of channel width and depth). Water flow was performed on artificial rinnenkarren system. A relation was looked for between the density of tributary channels and the average shape of the main channel, between the distance of tributary channels from each other and the shape of a given place of the main channel. The density and total length of the tributary channels on the lower and upper sections of the main channels being narrow at their lower end (11 pieces) and being wide at their lower end (10 pieces) of the rinnenkarren systems were calculated as well as their average proportional distance from the lower end of the main channel. The number of channel hollows was determined on the lower and upper sections of these main channels. It can be stated that the average shape of the main channel calculated to its total length depends on the density of the tributary channels and on the distance of tributary channels from each other. The main channel shape is smaller if less water flows on the floor for a long time because of the small density of the tributary channels and the great distance between the tributary channels. In this case, the channel deepens, but it does not widen. The width of the main channel depends on the number and location of the rivulets developing on channel-free relief. The main channel becomes narrow towards its lower end if the tributary rivulets are denser and longer on the upper part of the main rivulet developing on the channel-free, plain terrain and their distance is larger compared to the lower end. The channel hollows develop mainly at those places where the later developing tributary channels are hanging above the floor of the main channel. Thus, the former ones are younger than the latter ones. It can be stated that the morphology of the main channels (shape, channel hollows, and width changes of the main channel) is determined by the tributary channels (their number, location and age).