Error assessment of grid-based diffuse solar radiation models

Two grid-based diffuse solar radiation models, ESRI’s Solar Analyst (SA) and Kumar’s model (KM), were assessed using a data-independent approach where each model’s numerical results of clear sky diffuse radiation on V/U-shaped surfaces were compared with analytical results derived using each model’s assumptions. SA and KM consistently underestimate and overestimate, respectively, diffuse radiation at daily, seasonal, and annual scale relative to the analytical results based on each model’s parameterizations. Overall, SA performs better than KM in modeling diffuse radiation at most timescales. While SA and KM have similar error in calculating diffuse radiation on a horizontal surface, SA models sky view factor much better than KM, with mean absolute relative differences of 0.76% and 17.02%, respectively. KM has a large error in sky view factor as it does not consider the shading effect from surrounding terrain. Sky view factor error in SA is small and use of more zenith divisions can further reduce the error. Based on our previous study, model performance on clear sky global solar radiation was also evaluated. Overall, KM performs better than SA in global radiation as KM performs better than SA in modeling direct radiation which is the major component of global radiation.     

Enabling point pattern analysis on spatial big data using cloud computing: optimizing and accelerating Ripley’s K function

Performing point pattern analysis using Ripley’s K function on point events of large size is computationally intensive as it involves massive point-wise comparisons, time-consuming edge effect correction weights calculation, and a large number of simulations. This article presented two strategies to optimize the algorithm for point pattern analysis using Ripley’s K function and utilized cloud computing to further accelerate the optimized algorithm. The first optimization sorted the points on their x and y coordinates and thus narrowed the scope of searching for neighboring points down to a rectangular area around each point in estimating K function. Using the actual study area in computing edge effect correction weights is essential to estimate an unbiased K function, but is very computationally intensive if the study area is of complex shape. The second optimization reused the previously computed weights to avoid repeating expensive weights calculation. The optimized algorithm was then parallelized using Open Multi-Processing (OpenMP) and hybrid Message Passing Interface (MPI)/OpenMP on the cloud computing platform. Performance testing showed that the optimizations effectively accelerated point pattern analysis using K function by a factor of 8 using both the sequential version and the OpenMP-parallel version of the optimized algorithm. While the OpenMP-based parallelization achieved good scalability with respect to the number of CPU cores utilized and the problem size, the hybrid MPI/OpenMP-based parallelization significantly shortened the time for estimating K function and performing simulations by utilizing computing resources on multiple computing nodes. Computational challenge imposed by point pattern analysis tasks on point events of large size involving a large number of simulations can be addressed by utilizing elastic, distributed cloud resources.     

西安市生态环境与经济发展的耦合关系

采用耦合协调模型对2000—2013年西安市生态环境系统和经济系统进行综合评价分析,并对两大系统的耦合关系及耦合度演变趋势进行研究。结果表明:生态环境系统综合水平呈现在波动中下降的演变趋势,经济系统综合发展水平呈现在波动中上升的演变趋势;西安市生态环境和经济发展二者协调性由低水平的经济滞后失调型发展为高水平的生态环境滞后协调型,并存在工业废气排量增加和不合理的产业结构等问题,因此,需加强污染物排放的控制与监管,实施绿色考核制度,坚持绿色生活方式,从转变经济增长方式和增长动力2个角度出发,从而实现区域经济社会可持续发展。     

A spatially distributed Clark’s unit hydrograph based hybrid hydrologic model (Distributed-Clark)

ABSTRACT

A hybrid hydrologic model (Distributed-Clark), which is a lumped conceptual and distributed feature model, was developed based on the combined concept of Clark’s unit hydrograph and its spatial decomposition methods, incorporating refined spatially variable flow dynamics to implement hydrological simulation for spatially distributed rainfall–runoff flow. In Distributed-Clark, the Soil Conservation Service (SCS) curve number method is utilized to estimate spatially distributed runoff depth and a set of separated unit hydrographs is used for runoff routing to obtain a direct runoff flow hydrograph. Case studies (four watersheds in the central part of the USA) using spatially distributed (Thiessen polygon-based) rainfall data of storm events were used to evaluate the model performance. Results demonstrate relatively good fit to observed streamflow, with a Nash-Sutcliffe efficiency (E_NS) of 0.84 and coefficient of determination (R²) of 0.86, as well as a better fit in comparison with outputs of spatially averaged rainfall data simulations for two models including HEC-HMS.     

Estimation of river discharge from non-trapezoidal open channel using QuickBird-2 satellite imagery/Utilisation des images satellites de Quickbird-2 pour le calcul des débits fluviaux en chenaux ouverts non-trapézoïdaux

Abstract

Abstract River discharge is traditionally acquired by measuring water stage and then converting the water stage to discharge by using a stage–discharge rating curve. The possibility of monitoring river discharge by satellite has not been adequately studied hitherto, because of the difficulty in making sufficiently precise measurements of the water surface. Since the successful launch of commercial satellites with very-high-resolution sensors, it has become possible to derive ground information from satellite data. To determine river discharge in a non-trapezoidal open channel, an efficient approach has been developed that uses mainly satellite data. The method, which focuses on the measurement of surface water width coupled with river width–stage and ?remote? stage–discharge rating curves, was applied to the Yangtze River (Changjiang) and an accurate estimate of river discharge was obtained. The method can be regarded as ancillary to traditional field measurement methods or other remote sensing methods.     

Temporal scaling in river flow: can it be chaotic? / L’invariance d’échelle de l’écoulement fluvial peut-elle être chaotique?

Abstract

Abstract Identification of the presence of scaling in the river flow process has been a challenging problem in hydrology. Studies conducted thus far have viewed this problem essentially from a stochastic perspective, because the river flow process has traditionally been assumed to be a result of a very large number of variables. However, recent studies employing nonlinear deterministic and chaotic dynamic concepts have reported that the river flow process could also be the outcome of a deterministic system with only a few dominant variables. In the wake of such reports, a preliminary attempt is made in this study to investigate the type of scaling behaviour in the river flow process (i.e. chaotic or stochastic). The investigation is limited only to temporal scaling. Flow data of three different scales (daily, 5-day and 7-day) observed in each of three rivers in the USA: the Kentucky River in Kentucky, the Merced River in California and the Stillaguamish River in Washington, are analysed. It is assumed that the dynamic behaviour of the river flow process at these individual scales provides clues about the scaling behaviour between these scales. The correlation dimension is used as an indicator to distinguish between chaotic and stochastic behaviours. The results are mixed with regard to the type of flow behaviour at individual scales and, hence, to the type of scaling behaviour, as some data sets show chaotic behaviour while others show stochastic behaviour. They suggest that characterization (chaotic or stochastic) of river flow should be a necessary first step in any scaling study, as it could provide important information on the appropriate approach for data transformation purposes.     

Stochastically conditioned flow paths and travel time distribution in highly heterogeneous synthesized aquifers / Chemins d’écoulement et distribution des temps de parcours au sein d’aquifères synthétiques hautement hétérogènes sous conditions stochastiques

Abstract

Abstract Characterization of heterogeneity at the field scale generally requires detailed aquifer properties such as transmissivity and hydraulic head. An accurate delineation of these properties is expensive and time consuming, and for many if not most groundwater systems, is not practical. As an alternative approach, stochastic representation of random fields is used and presented in this paper. Specifically, an iterative stochastic conditional simulation approach was applied to a hypothetical and highly heterogeneous pre-designed aquifer system. The approach is similar to the classical co-kriging technique; it uses a linear estimator that depends on the covariance functions of transmissivity (T), and hydraulic head (h), as well as their cross-covariances. A linearized flow equation along with a conditional random field generator constitutes the iterative process of the conditional simulation. One hundred equally likely realizations of transmissivity fields with pre-specified geostatistical parameters were generated, and conditioned to both limited transmissivity and head data. The successful implementation of the approach resulted in conditioned flow paths and travel-time distribution under different degrees of aquifer heterogeneity. This approach worked well for fields exhibiting small variances. However, for random fields exhibiting large variances (greater than 1.0), an iterative procedure was used. The results show that, as the variance of the ln[T] increases, the flow paths tend to diverge, resulting in a wide spectrum of flow conditions, with no direct discernable relationship between the degree of heterogeneity and travel time. The applied approach indicates that high errors may result when estimation of particle travel times in a heterogeneous medium is approximated by an equivalent homogeneous medium.     

Utility of high-resolution climate model simulations for water resources prediction over the Korean Peninsula: a sensitivity study / Utilité de simulations de modèle climatique à haute résolution pour la prévision des ressources en eau dans la Péninsule Coréenne: étude de sensibilité

Abstract

Abstract The utility of simulations of Global Climate Models (GCMs) for regional water resources prediction and management on the Korean Peninsula was assessed by a probabilistic measure. Global Climate Model simulations of an indicator variable (e.g. surface precipitation or temperature) were used for discriminating high vs low regional observations of a target variable (e.g. watershed precipitation or reservoir inflow). The formulation uses the significance probability of the Kolmogorov-Smirnov test for detecting differences between two distributions. High resolution Atmospheric Model Intercomparison Project-II (AMIP-II) type GCM simulations performed by the European Centre for Medium-Range Weather Forecasts (ECMWF) and AMIP-I type GCM simulations performed by the Korean Meteorological Research Institute (METRI) were used to obtain information for the indicator variables. Observed mean areal precipitation and temperature, and watershed-outlet discharge values for seven major river basins in Korea were used as the target variables. The results suggest that the use of the climate model nodal output from both climate models in the vicinity of the target basin with monthly resolution will be beneficial for water resources planning and management analysis that depends on watershed mean areal precipitation and temperature, and outlet discharge.     

CostMAP: an open-source software package for developing cost surfaces using a multi-scale search kernel

ABSTRACT

Cost surfaces are a crucial aspect of route optimization and least cost path (LCP) calculations and are used in awide range of disciplines including computer science, landscape ecology, and energy-infrastructure modeling. Linear features present akey weakness to traditional routing calculations along cost surfaces because they cannot identify whether moving from acell to its adjacent neighbors constitutes crossing alinear barrier (increased cost) or following acorridor (reduced cost). Following and avoiding linear features can drastically change predicted routes. We introduce an approach to address this adjacency issue using asearch kernel that identifies these critical barriers and corridors. We have built this approach into anew Java-based open-source software package– CostMAP (cost surface multi-layer aggregation program)– which calculates cost surfaces and cost networks using the search kernel. CostMAP allows users to input multiple GIS data layers and to set weights and rules for developing aweighted-cost network. We compare CostMAP performance with traditional cost surface approaches and show significant performance gains– both following corridors and avoiding barriers– by modeling the movement of alarge terrestrial animal– the Baird’s Tapir (Tapirus bairdii)– in amovement ecology framework and by modeling pipeline routing for carbon capture and storage (CCS).     

中国东北地壳结构的地震层析成像

中国东北地区自晚古生代晚期以来,受古亚洲洋、蒙古—鄂霍茨克洋、太平洋构造域的叠加作用,壳幔结构极为复杂。本文收集中国东北地区国家地震台网接收的100 980个P波和91 030个S波到时数据,采用地震走时层析成像方法获得了该地区地壳P波和S波速度结构,进而获得了泊松比结构,用以探讨复杂的地壳结构。成像结果显示：中国东北地区地壳地震波速度结构呈明显的横向不均匀性,不同构造单元和构造单元内部都存在不同程度的不均匀性。松辽盆地整体上浅层地壳以低速异常为主,尤其是S波速度,但部分区域分辨率较低,中下地壳存在较大范围高速异常,推测与太平洋俯冲、后撤导致的岩石圈拆沉和热物质上涌等动力学过程有关;北部的大兴安岭重力梯级带和长白山一线主要表现为低速异常,表明具有大范围的岩浆作用,广泛岩浆作用为固体矿产资源的形成提供热源或物源;长白山、五大连池等近代活动的火山下方部分地壳区域均表现了较强低地震波速度和高泊松比异常结构特征,表明仍存在活动的可能。