A comparison of hillslope drainage area estimation methods using high-resolution DEMs with implications for topographic studies of gullies

Topographic models provide a useful tool for understanding gully occurrence in the landscape but require reliable estimates of gully head drainage areas. Modern high-resolution topography data (collected using structure from motion photogrammetry or light detection and ranging) is increasingly used for topographic studies of gullies, but little work has been done to assess the variability of gully head drainage area estimates using different methods. This study evaluated alternative approaches to using high-resolution digital elevation models (DEMs) so that gully topographic models can be more readily applied to any area with suitably high-resolution data. Specifically, we investigated the impact of single- or multiple-direction flow routing algorithms, DEM hydrologic-enforcement procedures and spatial resolution on gully head drainage area estimation. We tested these methods on a 40 km² site centred on Weany Creek, a low-relief semi-arid landscape draining towards the Great Barrier Reef, Australia. Using a subroutine to separate gully heads into those with divergent or convergent flow patterns upslope, we found that divergent flow conditions occurred at half of 484 studied gullies. Drainage areas estimated by different flow routing algorithms were more variable in these divergent cases than for convergent cases. This variation caused a significant difference between topographic threshold parameters (slope b and intercept k) derived from single- or multiple-direction flow routing algorithms, respectively. Different methods of hydrologic enforcement (filling or breaching) also affected threshold analysis, resulting in estimates of the exponent b being ~188% higher if the DEM was filled than if breached. The testing of the methods to date indicates that a finer resolution (≤2 m) DEM and a multiple-direction flow routing algorithm achieve the most realistic drainage area estimates in low-relief landscapes. For Weany Creek we estimated threshold parameters k = 0.033 and b = 0.189, indicating that it is highly susceptible to gully erosion.     

The geomorphic legacy of water and erosion control structures in a semiarid rangeland watershed

Control over water supply and distribution is critical for agriculture in drylands where manipulating surface runoff often serves the dual purpose of erosion control. However, little is known of the geomorphic impacts and legacy effects of rangeland water manipulation infrastructure, especially if not maintained. This study investigated the geomorphic impacts of structures such as earthen berms, water control gates, and stock tanks, in a semiarid rangeland in the southwestern USA that is responding to both regional channel incision that was initiated over a century ago, and a more recent land use change that involved cattle removal and abandonment of structures. The functional condition of remnant structures was inventoried, mapped, and assessed using aerial imagery and lidar data. Headcut initiation, scour, and channel incision associated with compromised lateral channel berms, concrete water control structures, floodplain water spreader berms, and stock tanks were identified as threats to floodplains and associated habitat. Almost half of 27 identified lateral channel berms (48%) have been breached and 15% have experienced lateral scour; 18% of 218 shorter water spreader berms have been breached and 17% have experienced lateral scour. A relatively small number of 117 stock tanks (6%) are identified as structurally compromised based on analysis of aerial imagery, although many currently do not provide consistent water supplies. In some cases, the onset of localized disturbance is recent enough that opportunities for mitigation can be identified to alter the potentially damaging erosion trajectories that are ultimately driven by regional geomorphic instability. Understanding the effects of prior land use and remnant structures on channel and floodplain morphologic condition is critical because both current land management and future land use options are constrained by inherited land use legacy effects. Published 2017. This article is a U.S. Government work and is in the public domain in the USA     

Land use effects on climate in China as simulated by a regional climate model

A regional climate model (RegCM3) nested within ERA40 re-analyzed data is used to investigate the climate effects of land use change over China. Two 15-year simulations (1987―2001), one with current land use and the other with potential vegetation cover without human intervention, are conducted for a domain encompassing China. The climate impacts of land use change are assessed from the difference between the two simulations. Results show that the current land use (modified by anthropogenic ac- tivities) influences local climate as simulated by the model through the reinforcement of the monsoon circulation in both the winter and summer seasons and through changes of the surface energy budget. In winter, land use change leads to reduced precipitation and decreased surface air temperature south of the Yangtze River, and increased precipitation north of the Yangtze River. Land use change signifi- cantly affects summer climate in southern China, yielding increased precipitation over the region, de- creased temperature along the Yangtze River and increased temperature in the South China area (south-end of China). In summer, a reduction of precipitation over northern China and a temperature rise over Northwest China are also simulated. Both daily maximum and minimum temperatures are affected in the simulations. In general, the current land use in China leads to enhanced mean annual precipitation and decreased annual temperature over south China along with decreased precipitation over North China.     

Improvement of snowpack simulations in a regional climate model

To improve simulations of regional‐scale snow processes and related cold‐season hydroclimate, the Community Land Model version 3 (CLM3), developed by the National Center for Atmospheric Research (NCAR), was coupled with the Pennsylvania State University/NCAR fifth‐generation Mesoscale Model (MM5). CLM3 physically describes the mass and heat transfer within the snowpack using five snow layers that include liquid water and solid ice. The coupled MM5–CLM3 model performance was evaluated for the snowmelt season in the Columbia River Basin in the Pacific Northwestern United States using gridded temperature and precipitation observations, along with station observations. The results from MM5–CLM3 show a significant improvement in the SWE simulation, which has been underestimated in the original version of MM5 coupled with the Noah land‐surface model. One important cause for the underestimated SWE in Noah is its unrealistic land‐surface structure configuration where vegetation, snow and the topsoil layer are blended when snow is present. This study demonstrates the importance of the sheltering effects of the forest canopy on snow surface energy budgets, which is included in CLM3. Such effects are further seen in the simulations of surface air temperature and precipitation in regional weather and climate models such as MM5. In addition, the snow‐season surface albedo overestimated by MM5–Noah is now more accurately predicted by MM5–CLM3 using a more realistic albedo algorithm that intensifies the solar radiation absorption on the land surface, reducing the strong near‐surface cold bias in MM5–Noah. The cold bias is further alleviated due to a slower snowmelt rate in MM5–CLM3 during the early snowmelt stage, which is closer to observations than the comparable components of MM5–Noah. In addition, the over‐predicted precipitation in the Pacific Northwest as shown in MM5–Noah is significantly decreased in MM5–CLM3 due to the lower evaporation resulting from the longer snow duration. Copyright © 2011 John Wiley & Sons, Ltd.     

Concentrated flow erosion as a main source of sediments in Galicia,Spain

Knowledge of soil loss rates by water erosion under given climate, soil, topography, and management conditions is important for establishing soil conservation schemes. In Galicia, a region with Atlantic climatic conditions in Spain, field observations over the last decade indicate that interrill, rill and ephemeral gully erosion may be an important sediment source. The aim of this work was to assess concentrated erosion rates, describe types of rills and ephemeral gullies and determine their origin, evolution and importance as sediment sources. Soil surface state and concentrated flow erosion were surveyed on medium textured soils, developed over basic schists of the Ordenes Complex series (Coruña province, Spain) from 1997 to 2006. Soil surface state was characterized by crust development, tillage features and roughness degree. Soil erosion rate was directly measured in the field. Concentrated flow erosion took place mainly on seedbeds and recently tilled surfaces in late spring and by autumn or early winter. During the study period, erosion rates were highly variable and the following situations could be distinguished: (a) no incision or limited rill incision, i.e. below 2 Mg ha^?1 year^?1; (b) generalized rill and ephemeral gully incision in the class of mean values between 2·5 and 6·25 Mg ha^?1 year^?1, this was the most common erosion pattern; and (c) heavy erosion as observed during an extremely wet winter period, between October 2000 and February 2001, with erosion figures that may be about ten orders of magnitude higher, up to 55–60 Mg ha^?1 year^?1. Therefore, low values of soil losses are dominant, but also large values of rill and ephemeral gully erosion occurred during the study period. Copyright © 2009 John Wiley & Sons, Ltd.     

Diagnostic calibration and cross‐catchment transferability of a simple process‐consistent hydrologic model

The transferability of hydrologic models is of ever increasing importance for making improved hydrologic predictions and testing hypothesized hydrologic drivers. Here, we present an investigation into the variability and transferability of the recently introduced catchment connectivity model (Smith et al., 2013 ). The catchment connectivity model was developed following extensive experimental observations identifying the key drivers of streamflow in the Tenderfoot Creek Experimental Forest (Jencso et al., 2009 ; Jencso et al., 2010 ), with the goal of creating a simple model consistent with internal observations of catchment hydrologic connectivity patterns. The model was applied across seven catchments located within Tenderfoot Creek Experimental Forest to investigate spatial variability and transferability of model performance and parameterization. The results demonstrated that the model resulted in historically good fits (based on previous studies at the sites) to both the hydrograph and internal water table dynamics (corroborated with experimental observations). The impact of a priori parameter limits was also examined. It was observed that enforcing field‐based limits on model parameters resulted in slight reductions to streamflow hydrograph fits, but significant improvements to model process fidelity (as hydrologic connectivity), as well as moderate improvement in the transferability of model parameterizations from one catchment to the next. Copyright © 2016 John Wiley & Sons, Ltd.     

Impact of anticipated climate change on direct groundwater recharge in a humid tropical basin based on a simple conceptual model

A simple conceptual semi‐distributed modelling approach for assessing the impacts of climate change on direct groundwater recharge in a humid tropical river basin is investigated. The study area is the Chaliyar river basin in the state of Kerala, India. Many factors affecting future groundwater recharge include decrease or increase in precipitation and temperature regimes, coastal flooding, urbanization and changes in land use. The model is based on the water‐balance concept and links the atmospheric and hydrogeologic parameters to different hydrologic processes. It estimates daily water‐table fluctuation and is calibrated and validated using 10 years of data. Data for the first 6 years (2000 to 2005) is used for model calibration, and data for the remaining four years (2006 to 2009) is used for validation. For assessing the impact of predicted climate change on groundwater recharge during the period 2071–2100, temperature and precipitation data in two post climate change scenarios, A2 and B2, were predicted using the Regional Climate Model (RCM), PRECIS (Providing Regional Climates for Impact Studies). These data were then corrected for biases and used in a hydrologic model to predict groundwater recharge in the post climate change scenario. Due to lack of reliable data and proper knowledge as to the magnitude and extent of future climatic changes, it may not be possible to include all the possible effects quantitatively in groundwater recharge modelling. However, the study presents a scientific method to assess the impact of predicted climate change on groundwater recharge and would help engineers, hydrologists, administrators and planners to devise strategies for the efficient use as well as conservation of freshwater resources. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of a time-magnitude prediction model for earthquakes

In this paper we discuss the physical meaning of the magnitude-time model parameters for earthquake prediction. The gestation process for strong earthquake in all eleven seismic zones in China can be described by the magnitude-time prediction model using the computations of the parameters of the model. The average model parameter values for China are： b = 0.383, c=0.154, d = 0.035, B = 0.844, C = -0.209, and D = 0.188. The robustness of the model parameters is estimated from the variation in the minimum magnitude of the transformed data, the spatial extent, and the temporal period. Analysis of the spatial and temporal suitability of the model indicates that the computation unit size should be at least 4°× 4°for seismic zones in North China, at least 3° × 3° in Southwest and Northwest China, and the time period should be as long as possible.     

地震预报中震级-时间模型的应用

本文从震级-时间模型的表达式出发,讨论了模型参数所含的物理意义,通过计算中国大陆各地震带的震级-时间模型参数,认为只有其中11条地震带的强震孕震过程可以用该模型来表达,其平均参数为b=0.383,c=0.154,d=0.035,B=0.844,C=-0.209,D=0.188。通过变换资料的震级下限、时间域和空间域的计算结果,讨论了模型参数的鲁棒性。通过对模型应用空间范围的分析,在华北地震区,在利用震级—时间模型对未来的地震危险进行预测时,计算单元应不小于4°×4°,而在西北和西南地区,计算单元不要小于3°×3°。而通过对模型应用时间范围的分析,在应用模型进行预测时,资料时间域要尽量长一些。     

Application of a simple finite difference model for estimating evaporation from open water

In estimating the evaporation from open water, the challenge is to accurately quantify the change in heat stored in the water body. A simple finite difference model is described and a comparison made between measured values of water temperature and evaporation, from a reservoir in southeast England, and the values predicted by an equilibrium temperature model. The values predicted by the new model are in excellent agreement with the measurements and are closer to the measured values than those predicted by the equilibrium temperature model. The difference in performance is attributed to improved methods used for calculating the net radiation and the wind function. The simpler formulation of the finite difference model is considered to offset the disadvantage of the greater number of calculations required.     

Multiple objective evaluation of a simple phosphorus transfer model

A simple phosphorus (P) transfer model of the Welland catchment, UK, is evaluated against multiple objective functions using a Monte Carlo approach that combines calibration, identifiability, sensitivity and uncertainty analysis. The model is based on simple conceptual rainfall‐runoff and river routing components, combined with estimates of the daily non‐point source load derived from annual landuse‐based export coefficients, disaggregated as a function of the runoff. The model has limited data requirements, consistent with data availability, and is parsimoneous with respect to the number of parameters identified through inverse modelling. The best performing parameter sets capture the main aspects of the observed flow and total P (TP) concentrations and provide a suitable basis for a decision‐support tool. However, a trade‐off is evident between matching the observed flow peaks, flow recessions and TP concentrations simultaneously, highlighting some limitations of the model structure and/or calibration data. Model analysis indicates that daily non‐point source load cannot be described as a function of near‐surface runoff and land use alone, but that other influences, including seasonality, are important. However, further model development to improve performance is likely to introduce additional complexity (in terms of parameter numbers), and hence additional problems of parameter identifiability and output uncertainty, which in turn raises issues of the information content of the available data. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulation of regional land subsidence in the southern Yangtze Delta

Investigation of the deformation characteristics of individual hydrostratigraphic units is the key to construct a regional land subsidence model. All of 12 hydrostratigraphic units in the study area were discussed throughout. On the basis of the measured data of groundwater level, five kinds of changing patterns of groundwater level were deduced and the relationship between the deformation characteris-tics of aquifer units and the corresponding changing patterns of groundwater level was discussed. The study area is 1.7×104 km2, where the geological condition is complex. The changing patterns of groundwater level the hydrostratigraphic units have experienced vary from site to site and from time to time. Consequently, the deformation characteristics of units are sophisticated. An identical hydros-tratigraphic unit may exhibit different deformation characteristics, such as elasticity, elasto-plasticity, visco-elasticity, and visco-elasto-plasticity, at different sites or during different periods, not to mention the different units. The existing models are difficult to describe the complex visco-elasto-plastic con-stitutive law under the condition of land subsidence. So the Merchant’s model was modified to depict the visco-elasto-plastic behavior of units. Then a three-dimensional flow model with variable parame-ters and a vertical one-dimensional subsidence model were constructed and coupled. The coupled model was applied in simulating land subsidence in the southern Yangtze Delta and a satisfactory re-sult was obtained. The simulation results show that the new coupled model can depict the complex geological conditions and describe the developing process of land subsidence very well in the south-ern Yangtze Delta. The new model can surely be used to predict land subsidence in the future, which is very helpful to taking measurements to control land subsidence.     

Coupling the Xinanjiang model to a kinematic flow model based on digital drainage networks for flood forecasting

The Xinanjiang model, which is a conceptual rainfall‐runoff model and has been successfully and widely applied in humid and semi‐humid regions in China, is coupled by the physically based kinematic wave method based on a digital drainage network. The kinematic wave Xinanjiang model (KWXAJ) uses topography and land use data to simulate runoff and overland flow routing. For the modelling, the catchment is subdivided into numerous hillslopes and consists of a raster grid of flow vectors that define the water flow directions. The Xinanjiang model simulates the runoff yield in each grid cell, and the kinematic wave approach is then applied to a ranked raster network. The grid‐based rainfall‐runoff model was applied to simulate basin‐scale water discharge from an 805‐km² catchment of the Huaihe River, China. Rainfall and discharge records were available for the years 1984, 1985, 1987, 1998 and 1999. Eight flood events were used to calibrate the model's parameters and three other flood events were used to validate the grid‐based rainfall‐runoff model. A Manning's roughness via a linear flood depth relationship was suggested in this paper for improving flood forecasting. The calibration and validation results show that this model works well. A sensitivity analysis was further performed to evaluate the variation of topography (hillslopes) and land use parameters on catchment discharge. Copyright © 2009 John Wiley & Sons, Ltd.     

Multicriteria sensitivity analysis as a diagnostic tool for understanding model behaviour and characterizing model uncertainty

Complex hydrological models are being increasingly used nowadays for many purposes such as studying the impact of climate and land‐use change on water resources. However, building a high‐fidelity model, particularly at large scales, remains a challenging task, due to complexities in model functioning and behaviour and uncertainties in model structure, parameterization, and data. Global sensitivity analysis (GSA), which characterizes how the variation in the model response is attributed to variations in its input factors (e.g., parameters and forcing data), provides an opportunity to enhance the development and application of these complex models. In this paper, we advocate using GSA as an integral part of the modelling process by discussing its capabilities as a tool for diagnosing model structure and detecting potential defects, identifying influential factors, characterizing uncertainty, and selecting calibration parameters. Accordingly, we conduct a comprehensive GSA of a complex land surface–hydrology model, Modélisation Environmentale–Surface et Hydrologie (MESH), which combines the Canadian land surface scheme with a hydrological routing component, WATROUTE. Various GSA experiments are carried out using a new technique, called Variogram Analysis of Response Surfaces, for alternative hydroclimatic conditions in Canada using multiple criteria, various model configurations, and a full set of model parameters. Results from this study reveal that, in addition to different hydroclimatic conditions and SA criteria, model configurations can also have a major impact on the assessment of sensitivity. GSA can identify aspects of the model internal functioning that are counter‐intuitive and thus help the modeller to diagnose possible model deficiencies and make recommendations for improving development and application of the model. As a specific outcome of this work, a list of the most influential parameters for the MESH model is developed. This list, along with some specific recommendations, is expected to assist the wide community of MESH and Canadian land surface scheme users, to enhance their modelling applications.     

An integrated model to simulate sown area changes for major crops at a global scale

Dynamics of land use systems have attracted much attention from scientists around the world due to their ecological and socio-economic implications. An integrated model to dynamically simulate future changes in sown areas of four major crops (rice, maize, wheat and soybean) on a global scale is pre- sented. To do so, a crop choice model was developed on the basis of Multinomial Logit (Logit) model to model land users' decisions on crop choices among a set of available alternatives with using a crop utility function. A GIS-based Environmental Policy Integrated Climate (EPIC) model was adopted to simulate the crop yields under a given geophysical environment and farming management conditions, while the International Food Policy and Agricultural Simulation (IFPSIM) model was utilized to estimate crop price in the international market. The crop choice model was linked with the GIS-based EPIC model and the IFPSIM model through data exchange. This integrated model was then validated against the FAO statistical data in 2001-2003 and the Moderate Resolution Imaging Spectroradiometer (MODIS) global land cover product in 2001. Both validation approaches indicated reliability of the model for ad- dressing the dynamics in agricultural land use and its capability for long-term scenario analysis. Finally, the model application was designed to run over a time period of 30 a, taking the year 2000 as baseline. The model outcomes can help understand and explain the causes, locations and consequences of land use changes, and provide support for land use planning and policy making.     

Prediction of regional-scale groundwater recharge and nitrate storage in the vadose zone: A comparison between a global model and a regional model

Extensive nitrogen loads at the soil surface exceed plant uptake and soil biochemical capacity, and therefore lead to nitrogen accumulation in the deep vadose zone. Studies have shown that stored nitrogen in the vadose zone can eventually reach the water table and affect the quality of groundwater resources. Recently, global scale models have been implemented to quantify nitrate storage and nitrate travel time in the vadose zone. These global models are simplistic and relatively easy to implement and therefore facilitate analysis of the considered transport processes at a regional scale with no further requirements. However, the suitability of applying these models at a regional scale has not been tested. Here, we evaluate, for the first time, the performance and utility of global scale models at the regional scale. Applied to the Loess Plateau of China, we compare estimates of groundwater recharge and nitrate storage derived from global scale models with results from a regional scale approach utilizing the Richards and advection-dispersion equations. The estimated nitrate storage was compared to nitrate observations collected in the deep vadose zone (>50 m) at five sites across the Loess Plateau. Although both models predict similar spatial patterns of nitrate storage, the recharge fluxes were three times smaller and the nitrate storage was two times higher compared with the regional model. The results suggest that global scale models are a potentially useful screening tool, but require refinement for local scale applications.     

WRF-Hydro模型与新安江模型在陈河流域的应用对比

本文利用全球陆面数据同化系统与降雨观测数据,以陕西半湿润区陈河流域为研究对象,驱动WRF-Hydro模型,研究该模型的表现和适用性,并在结构、参数、输入输出和模拟结果方面与新安江模型对比.考虑到次表面层与实际包气带的区别,引入土层厚度乘子ZSOILFAC对前者进行等比缩放,发现其与新安江模型反推包气带的厚度有较好的一致性.研究表明:在陈河流域中WRF-Hydro计算步长须在建议值的基础上缩小; WRF-Hydro模型善于模拟洪水细节,新安江模型表现好且稳定;前者的径流深和洪峰合格率平于或略低于后者;在两个指标均合格的洪水中,前者平均均方根误差比后者小21.5%,但对于其他洪水,前者平均均方根误差比后者大56.2%; WRF-Hydro在洪水起涨时刻模拟较好,表现出其在中小流域应用的潜力.     

加速遗传算法在预应力连续箱梁桥有限元模型修正中的应用

针对预应力混凝土连续箱梁桥的通扬运河大桥,建立其有限元模型.基于评价模型修正效果的频率指标与振型相关系数指标、柔度矩阵指标和挠度指标,定义三个目标函数,利用环境激励的试验模态参数,通过实数编码加速遗传算法对有限元模型进行了修正.相对于基于频率指标与振型相关系数指标定义的目标函数,基于柔度矩阵指标或挠度指标定义的目标函数修正的有限元模型,能同时较准确地反映其动态和静态特性,也能实现其试验模态振型较精确的质量归一化.     

Enhanced modeling of latent heat flux from urban surfaces in the Noah/single-layer urban canopy coupled model

The numerical modeling of the impacts of urban buildings in mesoscale meteorological models has gradually improved in recent years. Correctly representing the latent heat flux from urban surfaces is a key issue in urban land-atmosphere coupling studies but is a common weakness in current urban canopy models. Using the surface energy balance data at a height of 140 m from a 325 m meteorological tower in Beijing, we conducted a 1-year continuous off-line simulation by using a coupled land surface model and a single-layer urban canopy model and found that this model has a relatively large systematic error for simulated latent heat flux. To improve the numerical method for modeling latent heat flux from urban surfaces, we combined observational analysis and urban land surface model to derive an oasis effect coefficient for urban green areas; to develop a temporal variation formula for water availability in urban impervious surfaces; and to specify a diurnal profile and the maximum values of anthropogenic latent heat release for four seasons. These results are directly incorporated into the urban land surface model to improve model performance. In addition, this method serves as a reference for studies in other urban areas.