Assessing different sources of uncertainty in hydrological projections of high and low flows: case study for Omerli Basin,Istanbul, Turkey

This study investigates the assessment of uncertainty contribution in projected changes of high and low flows from parameterization of a hydrological model and inputs of ensemble regional climate models (RCM). An ensemble of climate projections including 15 global circulation model (GCM)/RCM combinations and two bias corrections (change factor (CF) and bias correction in mean (BC)) was used to generate streamflow series for a reference and future period using the Hydrologiska Byråns Vattenbalansavdelning (HBV) model with the 25 best-fit parameter sets based on four objective functions. The occurrence time of high flows is also assessed through seasonality index calculation. Results indicated that the inputs of hydrological model from ensemble climate models accounts for greater contribution to the uncertainty related to projected changes in high flows comparing to the contribution from hydrological model parameterization. However, the uncertainty contribution is opposite for low flows, particularly for CF method. Both CF and BC increases the total mean variance of high and low flows. The variability in the occurrence time of high flows through RCMs is greater than the variability resulted from hydrological model parameters with and without statistical downscaling. The CF provides more accurate timing than BC and it shows the most pronounced changes in flood seasonality.     

基于OPAQ的城市空气质量预报系统研究

空气质量预测在国内的关注度日益提高,传统的空气质量预测系统通常运用数值化学传输模型,利用物理方程来计算污染物的扩散、沉降和化学反应。而化学传输模型的预测准确性很大程度上需要依赖详细的污染源排放信息和气象模型的输出结果。基于统计模型的OPAQ空气质量预报业务系统,采用人工神经网络算法,可预测各污染物的日均值或日最大值。并对北京空气质量预报的结果进行了评价,OPAQ空气质量预报业务系统对空气质量预测的准确性较高,能够利用较低的计算资源得到较为准确的预测结果。与数值预报相比,OPAQ空气质量预报业务系统不需要大量的基础数据作为输入,可弥补数值预报的不足,并成为数值预报的有力补充。     

Deterministic and Stochastic Simulations for Solid Waste Collection Systems – A SWIM Approach

The SWIM model is the first systems model in Australia that deals with integrated waste management systems. The main modelling approach adopted is simulation, which is based on both deterministic and stochastic models for collection systems.These models are described in this paper, after a number of modelling approaches are reviewed. An example of the application of the SWIM model is given, and planned extensions to the SWIM model are briefly outlined.     

ANN modelling of sediment concentration in the dynamic glacial environment of Gangotri in Himalaya

The present study explores for the first time the possibility of modelling sediment concentration with artificial neural networks (ANNs) at Gangotri, the source of Bhagirathi River in the Himalaya. Discharge, rainfall and temperature have been considered as the main controlling factors of variations in sediment concentration in the dynamic glacial environment of Gangotri. Fourteen feed forward neural networks with error back propagation algorithm have been created, trained and tested for prediction of sediment concentration. Seven models (T1-T7) have been trained and tested in the non-updating mode whereas remaining seven models (T1a-T7a) have been trained in the updating mode. The non-updating mode refers to the scenario where antecedent time (previous time step) values are not used as input to the model. In case of the updating mode, antecedent time values are used as network inputs. The inputs applied in the models are either the variables mentioned above as individual factors (single input networks) or a combination of them (multi-input networks). The suitability of employing antecedent time-step values as network inputs has hence been checked by comparative analysis of model performance in the two modes. The simple feed forward network has been improvised with a series parallel non-linear autoregressive with exogenous input (NARX) architecture wherein true values of sediment concentration have been fed as input during training. In the glacial scenario of Gangotri, maximum sediment movement takes place during the melt period (May–October). Hence, daily data of discharge, rainfall, temperature and sediment concentration for five consecutive melt periods (May–October, 2000–2004) have been used for modelling. High Coefficient of determination values [0.77–0.88] have been obtained between observed and ANN-predicted values of sediment concentration. The study has brought out relationships between variables that are not reflected in normal statistical analysis. A strong rainfall: sediment concentration and temperature: sediment concentration relationship is shown by the models which are not reflected in statistical correlation. It has also been observed that usage of antecedent time-step values as network inputs does not necessarily lead to improvement in model performance.     

Forecasting of cyanobacterial density in Torrão reservoir using artificial neural networks

The ability of general regression neural networks (GRNN) to forecast the density of cyanobacteria in the Torr?o reservoir (Tamega river, Portugal), in a period of 15 days, based on three years of collected physical and chemical data, was assessed. Several models were developed and 176 were selected based on their correlation values for the verification series. A time lag of 11 was used, equivalent to one sample (periods of 15 days in the summer and 30 days in the winter). Several combinations of the series were used. Input and output data collected from three depths of the reservoir were applied (surface, euphotic zone limit and bottom). The model that presented a higher average correlation value presented the correlations 0.991; 0.843; 0.978 for training, verification and test series. This model had the three series independent in time: first test series, then verification series and, finally, training series. Only six input variables were considered significant to the performance of this model: ammonia, phosphates, dissolved oxygen, water temperature, pH and water evaporation, physical and chemical parameters referring to the three depths of the reservoir. These variables are common to the next four best models produced and, although these included other input variables, their performance was not better than the selected best model.     

Optimal Sampling and Problematic Likelihood Functions in a Simple Population Model

Markov chains provide excellent statistical models for studying many natural phenomena that evolve with time. One particular class of continuous-time Markov chain, called birth–death processes, can be used for modelling population dynamics in fields such as ecology and microbiology. The challenge for the practitioner when fitting these models is to take measurements of a population size over time in order to estimate the model parameters, such as per capita birth and death rates. In many biological contexts, it is impractical to follow the fate of each individual in a population continuously in time, so the researcher is often limited to a fixed number of measurements of population size over the duration of the study. We show that, for a simple birth–death process, with positive Malthusian growth rate, subject to common practical constraints, there is an optimal schedule for measuring the population size that minimises the expected confidence region of the parameter estimates. Throughout our exposition of the optimal experimental design, we compare it to a simpler equidistant design, where the population is sampled at regular intervals. This is an experimental design worthy of comparison since it can represent a much simpler design to implement in practice. In order to find optimal experimental designs for our population model, we make use of a combination of useful statistical machinery. Firstly, we use a Gaussian diffusion approximation of the underlying discrete-state Markov process, which allows us to obtain analytical expressions for Fisher’s information matrix (FIM), which is crucial to optimising the experimental design. We also make use of the cross-entropy method of stochastic optimisation for the purpose of maximising the determinant of FIM to obtain the optimal experimental designs. Our results show that the optimal schedule devised by others for a simple model of population growth without death can be extended, for large populations, to the two-parameter model that incorporates both birth and death. For the simple birth–death process, we find that the likelihood surface is also problematic and poses serious problems for point estimation and easily defining confidence regions. A Bayesian approach to inference is proposed as a way in which these problems could be circumvented.     

Water quality trend analysis for the Karoon River in Iran

The Karoon River basin, with a basin area of 67,000 km², is located in the southern part of Iran. Monthly measurements of the discharge and the water quality variables have been monitored at the Gatvand and Khorramshahr stations of the Karoon River on a monthly basis for the period 1967–2005 and 1969–2005 for Gatvand and Khorramshahr stations, respectively. In this paper the time series of monthly values of water quality parameters and the discharge were analyzed using statistical methods and the existence of trends and the evaluation of the best fitted models were performed. The Kolmogorov–Smirnov test was used to select the theoretical distribution which best fitted the data. Simple regression was used to examine the concentration-time relationships. The concentration-time relationships showed better correlation in Khorramshahr station than that of Gatvand station. The exponential model expresses better concentration – time relationships in Khorramshahr station, but in Gatvand station the logarithmic model is more fitted. The correlation coefficients are positive for all of the variables in Khorramshahr station also in Gatvand station all of the variables are positive except magnesium (Mg²⁺), bicarbonates () and temporary hardness which shows a decreasing relationship. The logarithmic and the exponential models describe better the concentration-time relationships for two stations.     

Assessing land use impacts on flood processes in complex terrain by using GIS and modeling approach

A distributed hydrologic modeling and GIS approach is applied for the assessment of land use impact in the Steinsel sub-basin, Alzette, Grand-Duchy of Luxembourg. The assessment focuses on the runoff contributions from different land use classes and the potential impact of land use changes on runoff generation. The results show that the direct runoff from urban areas is dominant for a flood event compared with runoff from other land use areas in this catchment, and tends to increase for small floods and for the dry season floods, whereas the interflow from forested, pasture and agricultural field areas contributes to the recession flow. Significant variations in flood volume, peak discharge, time to the peak, etc., are found from the model simulation based on the three hypothetical land use change scenarios.     

Attributes of reliable long-term landscape-scale studies: Malpractice insurance for landscape ecologists

Monitoring long-term change in forested landscapes is an intimidating challenge with considerable practical, methodological, and theoretical limitations. Current field approaches used to assess vegetation change at the plot-to-stand scales and nationwide forest monitoring programs may not be appropriate at landscape scales. We emphasize that few vegetation monitoring programs (and, thus, study design models) are designed to detect spatial and temporal trends at landscape scales. Based primarily on advice from many sources, and trial and error, we identify 14 attributes of a reliable long-term landscape monitoring program: malpractice insurance for landscape ecologists. The attributes are to: secure long-term funding and commitment; develop flexible goals; refine objectives; pay adequate attention to information management; take an experimental approach to sampling design; obtain peer-review and statistical review of research proposals and publications; avoid bias in selection of long-term plot locations; insure adequate spatial replication; insure adequate temporal replication; synthesize retrospective, experimental, and related studies; blend theoretical and empirical models with the means to validate both; obtain periodic research program evaluation; integrate and synthesize with larger and smaller scale research, inventory, and monitoring programs; and develop an extensive outreach program. Using these 14 attributes as a guide, we describe one approach to assess the potential effect of global change on the vegetation of the Front Range of the Colorado Rockies. This self-evaluation helps identify strengthes and weaknesses in our program, and may serve the same role for other landscape ecologists in other programs.     

Estimation of Vegetation Coverage in Semi-arid Sandy Land Based on Multivariate Statistical Modeling Using Remote Sensing Data

The estimation of vegetation coverage is essential in the monitoring and management of arid and semi-arid sandy lands. But how to estimate vegetation coverage and monitor the environmental change at global and regional scales still remains to be further studied. Here, combined with field vegetation survey, multispectral remote sensing data were used to estimate coverage based on theoretical statistical modeling. First, the remote sensing data were processed and several groups of spectral variables were selected/proposed and calculated, and then statistically correlated to measured vegetation coverage. Both the single- and multiple-variable-based models were established and further analyzed. Among all single-variable-based models, that is based on Normalized Difference Vegetation Index showed the highest R (0.900) and R ² (0.810) as well as lowest standard estimate error (0.128024). Since the multiple-variable-based model using multiple stepwise regression analysis behaved much better, it was determined as the optimal model for local coverage estimation. Finally, the estimation was conducted based on the optimal model and the result was cross-validated. The coefficient of determination used for validation was 0.867 with a root-mean-squared error (RMSE) of 0.101. The large-scale estimation of vegetation coverage using statistical modeling based on remote sensing data can be helpful for the monitoring and controlling of desertification in arid and semi-arid regions. It could serve for regional ecological management which is of great significance.     

Development of geomorphologic instantaneous unit hydrograph for a large watershed

Hill torrents cause a lot of environmental and property damage in Pakistan every year. Proper assessment of direct runoff in the form of hill torrents is essential for protection of environment, property, and human life. In this paper, direct surface runoff hydrograph (DSRH) was derived for a large catchment using the geomorphologic instantaneous unit hydrograph concept. The catchment with hill torrent flows in semi-arid region of Pakistan was selected for this study. It was divided into series of linear cascades and hydrologic parameters required for Nash's conceptual model, and were estimated using geomorphology of the basin. Geomorphologic parameters were derived from satellite images of the basin and ERDAS and ArcGIS were used for data processing. Computer program was developed to systematically estimate the dynamic velocity, its related parameters by optimization and thereby to simulate the DSRH. The data regarding rainfall-runoff and satellite images were collected from Punjab Irrigation and Power Department, Pakistan. Model calibration and validation was made for 15 rainfall-runoff events. Ten events were used for calibration and five for validation. Model efficiency was found to be more than 90% and root mean square error to be about 5%. Impact of variation in model parameters (shape parameter and storage coefficient) on DSRH was investigated. For shape parameter, the number of linear cascades varied from 1 to 3 and it was found that the shaper parameter value of 3 produced the best DSRH. Various values of storage coefficient were used and it was observed that the value determined from geomorphology and the dynamic velocity produced the best results.     

Evaluation of the temporal scaling variability in forecasting ground-level ozone concentrations obtained from multiple linear regressions

Ozone is a highly unpredictable pollutant which severely affects living conditions in urban and surrounding areas in the Mediterranean basin. This secondary pollutant periodically reaches extremely high concentrations, damaging human health. Multiple linear regression has been widely used in previous works due to the fact that it is a simple and versatile method for forecasting ozone concentrations. However, these models usually prove their validity using fulfillment of statistical constraints, ignoring other intrinsic characteristics existing in the time series, such as the temporal scaling behavior and the data distribution over different time scales. In previous works, it has been demonstrated that observed ozone time series are of a multifractal nature, meaning that the data distribution can be described by using the multifractal spectrum. This work focuses on the capacity of a forecasting model to reproduce the scaling features existing in an observed time series when several chemical and meteorological explanatory variables are introduced following the stepwise procedure. A comparison between the observed spectrum and the simulated ones for each step is used to check which explanatory variables better reproduce the multifractal nature in real ozone time series. It has been confirmed that a model with few explanatory variables allows reproducing the multifractal nature in the simulated time series with an acceptable accuracy without compromising the values of the coefficient of determination and root-mean-squared error, which were used as performance indicators.     

Modeling Runoff Response to Land Cover and Rainfall Spatial Variability in Semi-Arid Watersheds

Hydrologic response is an integrated indicator of watershed condition, and significant changes in land cover may affect the overall health and function of a watershed. This paper describes a procedure for evaluating the effects of land cover change and rainfall spatial variability on watershed response. Two hydrologic models were applied on a small semi-arid watershed; one model is event-based with a one-minute time step (KINEROS), and the second is a continuous model with a daily time step (SWAT). The inputs to the models were derived from Geographic Information System (GIS) theme layers of USGS digital elevation models, the State Soil Geographic Database (STATSGO) and the Landsat-based North American Landscape Characterization classification (NALC) in conjunction with available literature and look up tables. Rainfall data from a network of 10 raingauges and historical stream flow data were used to calibrate runoff depth using the continuous hydrologic model from 1966 to 1974. No calibration was carried out for the event-based model, in which six storms from the same period were used in the calculation of runoff depth and peak runoff. The assumption on which much of this study is based is that land cover change and rainfall spatial variability affect the rainfall-runoff relationships on the watershed. To validate this assumption, simulations were carried out wherein the entire watershed was transformed from the 1972 NALC land cover, which consisted of a mixture of desertscrub and grassland, to a single uniform land cover type such as riparian, forest, oak woodland, mesquite woodland, desertscrub, grassland, urban, agriculture, and barren. This study demonstrates the feasibility of using widely available data sets for parameterizing hydrologic simulation models. The simulation results show that both models were able to characterize the runoff response of the watershed due to changes of land cover.     

Modeling the Hydroclimatic Disturbance of Soil Landscapes in the Southern Canadian Plains: The Problems of Scale and Place

The sensitivity of soil landscapes to climatic variability andhydroclimatic events can be expressed as a landscape change safety factor, the ratio of potential disturbance to resistance to change. The use of a geographic information system (GIS) enables the spatially-explicit modeling of landscape sensitivity, but also raises the risk of violating the characteristic scales of disturbance and resistance, because the GIS technically simplifies the extrapolation of models, and associated concepts, to landscapes and scales notrepresented by the digital data base. Embedding landscape sensitivity into hierarchy theory, the formal analysis of the hierarchical structure of complex systems, provides a conceptual framework for the transfer of models and variablesamong landscape scales. In the subhumid southern Canadian plains, major hydroclimatic events (strong winds, intense rain,rapid snow melt) cause much of the physical disturbance of soillandscapes and terrestrial ecosystems. Prolonged dry or wet weather influences the resistance of soil and vegetation to these events. The potential disturbance of soil landscapes therefore can be derived from the probabilities of extreme events and seasonal conditions, as recorded in instrumental and proxy climate records. This time series analysis can belinked to the modeling of landscape sensitivity by establishingthe probabilities of hydroclimatic events and climatic conditions which may exceed or lower the resistance of individual soil landscapes.     

Variable-source flood pulsing in a semi-arid transboundary watershed: the Chobe River, Botswana and Namibia

The Chobe River, characterized by an unusual flood pulsing regime and shared between Botswana and Namibia, lies at the heart of the world’s largest transfrontier conservation area (the Kavango–Zambezi Transfrontier Conservation Area). Significant ecological changes and vegetation conversions are occurring along its floodplains. Various scenarios for agricultural and urban water use are currently being proposed by the government of Botswana. However, the understanding of the river’s annual flow regime and timing of the relative contributions of water from three different sources is relatively poor. In light of past and future climate change and variability, this means that allocating water between ecological flows and economic and domestic uses will become increasingly challenging. We reconstruct the inundation history in this basin to help ease this challenge. This paper presents a spatiotemporal approach to estimate the contribution of water from various sources and the magnitude of changes in the flooding extent in the basin between 1985 and 2010. We used time series analysis of bimonthly NOAA AVHRR and NASA MODIS data and climatologic and hydrologic records to determine the flooding timing and extent. The results indicate that between 12 and 62 % of the basin is flooded on an annual basis and that the spatial extent of the flooding varies throughout the year as a function of the timing of peak discharge in two larger basins. A 30-year trend analysis indicates a consistent decline in the average monthly flooded area in the basin. The results may prove useful in future water utilization feasibility studies, in determining measures for protecting ecological flows and levels, and in ecosystem dynamics studies in the context of current and future climate change and variability.     

Statistical Simulation to Estimate Uncertain Behavioral Parameters of Hybrid Energy-Economy Models

In energy-economy modeling, new hybrid models attempt to combine the technological explicitness of bottom-up models with the macroeconomic feedbacks and statistically estimated behavioral parameters of top-down models. However, statistical estimation of behavioral parameters (portraying firm and household technology choices) with such models is challenged by the number of uncertain variables and the lack of historical data on technologies in terms of capital costs, operating costs, and market shares. Multiple combinations of parameter values might equally explain past technology choices. This paper reports on the application of a Bayesian statistical simulation approach for estimating the most likely values for these key behavioral parameters in order to best explain past technology choices and then simulate policies to influence future technology choices. The method included (1) data collection of key technology market shares, capital costs, and operating costs over the past; (2) backcasting a hybrid energy-economy model over a historical time period; and (3) the application of Markov chain Monte Carlo statistical simulation using the Metropolis–Hastings algorithm as a tool for estimating distributions for key parameters in the model. The results provide a means of indicating the uncertainty bounds around key behavioral parameters when generating forecasts of the effect of certain policies. However, the results also indicate that this approach may have limited applicability, given that future available technologies may differ substantially from past technologies and that it is difficult to separate the effects of parameter uncertainty from model structure uncertainty.     

A physically based method for prediction of runoff frequencies

A method is described for predicting changes in flood frequencies resulting from various causes affecting runoff in a watershed. The method utilizes the curve number and the traingular unit hydrograph procedures of the United States Department of Agriculture's Soil Conservation Service to transform rainfall frequencies into frequencies of peak discharge. The theory underlying the method is discussed. Application is illustrated by an example.     

常州地区植被排放VOC的估算研究

通过对2011年常州地区各类植物VOC排放因子，以及各类植被分布面积等数据统计分析，采用BEIS模型为参考的估算方法，建立起常州地区植被VOC的排放清单。结果表明，植被所排放VOC的变化规律既与植物本身有关又与气温和太阳辐射有关，区域内年植被VOC的总排放量为1．13×104 t。