DOWNSLOPE EROSION PROCESS UNDER UPSLOPE RUNOFF AND SEDIMENT USING A DUAL-BOX SYSTEM

1 INTRODUCTION Soil erosion at the hilly-gully region of the Loess Plateau has obvious vertical erosion zones from watershed boundary to gully edge, i.e., sheet erosion-dominated zone, rill erosion-dominated zone and shallow gully erosion-dominated zone, from top to bottom (Chen et al., 1988). Meanwhile, upslope runoff and sediment have a significant impact on the downslope erosion process. But with the limits of research methods,there is not much data to quantify upslope runoff and sedi…     

风垂直切变对中尺度地形对流降水影响的研究

针对长江中下游中尺度地形特点以及暴雨过程发生发展期间风垂直切变的主要观测特征,设计了一系列中尺度地形的三维理想数值试验,分析了干大气地形流和重力波特征,探讨了条件不稳定湿大气地形对流降水的模态分布,在此基础上研究了圆形、直线风垂直切变和切变厚度对中尺度地形对流降水强度和模态分布的影响.结果发现:在 F_r≈1的干大气条件下,气流遇到地形后分支、绕流和爬升现象同时存在,地形激发的重力波在水平和垂直方向上传播,其在迎风坡、背风坡、地形上游和下游的振幅不同,并组织出不同强度的垂直上升运动.在F_r > 1的条件不稳定湿大气下,地形对流降水主要存在三种模态,即迎风坡和背风坡准静止对流降水以及地形下游移动性对流降水,地形对流降水的形成与重力波在低层组织的上升运动密切相关.风垂直切变对地形对流降水的强度和模态分布有重要作用,其中圆形风垂直切变(风随高度旋转)不仅影响地形下游对流降水系统的移动方向,而且影响迎风坡和背风坡山脚处对流降水中心的分布和强度;直线风垂直切变(风随高度无旋转)主要影响地形对流降水的移动速度和强度.风随高度自下而上顺(逆)时针旋转,地形对流系统向下游传播时向右(左)偏移.风垂直切变主要通过影响地形重力波的结构和传播以及对流系统的形成、移动方向和速度,来影响地形对流降水的模态分布,其中对流层中低层的风垂直切变对地形对流降水强度和模态分布有重要影响.     

RILL EROSION PROCESS AND RILL FLOW HYDRAULIC PARAMETERS

In the rill erosion process, run-on water and sediment from upslope areas, and rill flow hydraulic parameters have significant effects on sediment detachment and transport. However, there is a lack of data to quantify the effects of run-on water and sediment and rill flow hydraulic parameters on rill erosion process at steep hillslopes, especially in the Loess Plateau of China. A dual-box system, consisting of a 2-m-long feeder box and a 5-m-long test box with 26.8% slope gradient was used to quantify the effects of upslope runoff and sediment, and of rill flow hydraulic parameters on the rill erosion process. The results showed that detachment-transport was dominated in rill erosion processes; upslope runoff always caused the net rill detachment at the downslope rill flow channel, and the net rill detachment caused by upslope runoff increased with a decrease of runoff sediment concentration from the feeder box or an increase of rainfall intensity. Upslope runoff discharging into the rill flow channel or an increase of rainfall intensity caused the rill flow to shift from a stratum flow into a turbulent flow. Upslope runoff had an important effect on rill flow hydraulic parameters, such as rill flow velocity, hydraulic radius, Reynolds number, Froude number and the Darcy-Weisbach resistance coefficient. The net rill detachment caused by upslope runoff increased as the relative increments of rill flow velocity, Reynolds number and Froude number caused by upslope runoff increased. In contrast, the net rill detachment decreased with an increase of the relative decrement of the Darcy-Weisbach resistance coefficient caused by upslope runoff. These findings will help to improve the understanding of the effects of run-on water and sediment on the erosion process and to find control strategies to minimize the impact of run-on water.     

中国东部层积云发展过程中云微物理特征的演变

基于2007—2010年的CloudSat卫星观测数据,以云层液态水路径为指标将层积云的发展过程划分为五个阶段,对比研究了中国东部降水与非降水层积云发展过程中云微物理特征和云微物理机制的演变,并分析了其海陆差异.研究表明:非降水层积云中,云滴增长主要通过凝结过程完成,但云滴的凝结增长有限,难以形成降水,在非降水层积云发展的旺盛阶段,云层中上部云滴发生较弱的碰并过程.降水层积云中云滴碰并增长活跃,当云层液态水路径小于500 g·m~(-2)时,云滴在从云顶下落至云底的过程中持续碰并,并在云底附近出现云水向雨水的转化;当降水层积云液态水路径超过500 g·m~(-2)时,云滴碰并增长主要发生在云层上部,在云层中部,云液态水含量、液态粒子数浓度和液态粒子有效半径达到最大,云水向雨水的转化最为活跃.层积云微物理特征的海陆差异主要是由海陆上空气溶胶浓度和云中上升气流强度不同导致的.在非降水层积云中下部,陆地丰富的气溶胶为云滴凝结增长提供了充足的云凝结核,因而云微物理量的量值在陆地上空更大,而在云层中上部,云滴凝结增长达到极限,海洋充足的水汽输送使云微物理量的量值在海洋上空更大.当降水层积云液态水路径大于500 g·m~(-2)时,陆地层积云中更强的上升气流使大量云滴在云层中上部累积滞留,云滴碰并增长活跃,云层中上部云微物理量的量值在陆地上空更大.     

Reducing the biases in shortwave cloud radiative forcing in tropical and subtropical regions from the perspective of boundary layer processes

Biases in shortwave cloud radiative forcing (SWCF), which cause overestimates in tropical regions and underestimates in subtropical marine stratocumulus regions, are common in many climate models. Here, two boundary layer processes are investigated in the atmospheric model GAMIL2, entrainment at the top of the boundary layer and longwave radiative cooling at the top of stratocumulus clouds, in order to reduce biases and reveal the mechanisms underlying these processes. Our results show that including the entrainment process in the model can reduce negative SWCF biases in most tropical regions but increases positive SWCF biases in subtropical marine stratocumulus regions. This occurs because entrainment reduces the low-level cloud fraction and its cloud liquid water content by suppressing the vertical turbulent diffusion in the boundary layer and decreasing the relative humidity when warm and dry free atmosphere is entrained in the boundary layer. Longwave radiative cooling at the top of stratocumulus clouds can enhance turbulent diffusion within the stratocumulus-topped boundary layer. When combined with the entrainment process, longwave radiative cooling reduces the positive SWCF biases in subtropical marine stratocumulus regions that are observed using the entrainment process alone. The incorporation of these two boundary layer processes improves the simulated SWCF in tropical and subtropical regions in GAMIL2.     

Cumulus dynamics: Local compensating subsidence and its implications for cumulus parameterization

Observations of air flow in and around convective clouds are summarized and discussed in light of the requirements for parameterization of midlatitude convection. Both theory and observation indicate that a substantial portion of the compensating subsidence occurs as dry downdrafts in the immediate vicinity of convective clouds, which suggests that an additional physical mechanism is necessary in linking convection to the large scale. A conceptual three-dimensional model of midlatitude organized convection along with its implications for parameterization, particularly in mesoscale modeling, is presented. The effects of moist downdrafts are also considered and some differences between midlatitude and tropical convection are discussed.     

Characteristics of stemflow for tall stewartia (Stewartia monadelpha) growing on a hillslope

The characteristics of stemflow were observed in a tall stewartia (Stewartia monadelpha) deciduous forest on a hillslope in central Japan, revealing new findings for a previously unreported type of deciduous forest. Using 2-year observations of 250 rainfall events, we analyzed seasonal and spatial variations in stemflow for several trees, and applied additional data sets of throughfall and plant area index (PAI) to produce a rough estimate of seasonal variations in rainfall redistribution processes and canopy architecture for a single tree. Compared to previous findings for other deciduous tree species, the ratios of throughfall, stemflow, and interception to open-area rainfall obviously varied with PAI changes for tall stewartia. Meteorological conditions of rainfall amount, rainfall intensity, wind speed, and wind direction had little effect on stemflow generation, which was mainly affected by variation in canopy architecture. Three novel characteristics of stemflow were identified for several tall stewartia trees. First, the yearly stemflow ratio at the forest-stand level for tall stewartia (12%) was high compared to previous findings on beech and oak stands, indicating tall stewartia has considerably high potential to generate a great amount of stemflow. Second, stemflow tended to be 1.3–2.0 times greater in the leafed period than in the leafless period. Third, the amount of stemflow was 12–132 times greater on the downslope side of the stem than on the upslope side. It likely caused by the uneven area between the upslope and downslope sides of the canopy and by asymmetrical stemflow pathways between the upslope and downslope sides of the trunk due to downslope tilting of the tree trunk.     

Modelling subsurface flow conditions in a salinized catchment in south‐western Australia,with a view to improving management practices

Finite element modelling of the saturated–unsaturated surface–subsurface flow mechanisms operative in a small salinized catchment in south‐western Australia was used to help define the flow system and explain the causes of waterlogging and salinization there. Data available at the site from a previous study were used to obtain a first approximation to the flow system. Altering the properties of some of the strata gave a closer calibration. It was found that the modelled saturated hydraulic conductivity of the B horizon in the duplex soil zone needed to be at least an order of magnitude lower than that measured in order to reproduce the perching conditions observed in the field. Also, the model indicated the influence of a doleritic dyke, whose presence was confirmed by field measurement. Our analysis showed that there were two main flow systems operating in the hillslope. The first, and most dominant, was the recharge occurring through the upslope gradational soil zone and percolating down to both the deeply weathered regolith and the basal aquifer. The second flow system is an unsaturated flow system operating in the high permeability A horizon in the downslope duplex soil zone. The first system is primarily responsible for the saline seepage zone in the valley bottom. The second contributes to the waterlogging and perching occurring upslope of the seepage zone. Vertical flow through the higher permeability B horizon in the gradational soil zone in the upper slopes is a major contributor of recharge. Recharge by flow through macropores occurs where, but only where, perched aquifers develop and allow the macropores to be activated. Areas with perched aquifers occurred in downslope locations and near a doleritic dyke located upslope. Thus, the area where macropore recharge occurred was not large. The recharge rate required to maintain the piezometric levels at present values is only about 30 mm/yr (about 5% of the annual rainfall). The piezometric levels under the upper part of the catchment varied greatly with only small changes in recharge rate. A 50% reduction in recharge rate had the effect of reducing the length of the seepage zone at the end of winter by 40%. Changes in recharge rate had little effect on the extent of the perched aquifer at the end of winter. Deep‐rooted perennial forages, shrubs or trees on the gradational soil zone in the upper part of the catchment and on the zones upslope of geological barriers to flow would be required to reduce the recharge and to allow for rehabilitation of the saline valley floor. Waterlogging associated with the perched water table in the bottom part of the catchment would be best addressed by tree plantations located just upslope of the salinized zone in the valley floor. Copyright © 1999 John Wiley & Sons, Ltd.     

Investigation of the impacts of vegetation distribution and evaporative cooling on synthetic urban daytime climate using a coupled LES—LSM model

Local flow properties and regional weather or climate are strongly affected by land‐atmosphere interactions of momentum and scalars within the daytime convective boundary layer (CBL). In this study, we investigate the impact of green space scale on the daytime atmospheric boundary layer (ABL) over a synthetic urban domain using a recently developed large‐eddy simulation‐land surface model (LES–LSM) framework. With the use of realistic soundings as initial conditions, a series of numerical experiments over synthetic urban surfaces with varied scale of vegetated area is performed. Simulated micrometeorological properties, surface fluxes, basic CBL characteristics, and cloud distribution are analysed. The results show reference‐level air potential temperature and specific humidity as well as surface fluxes over green space are significantly affected by the scale of green space in the urban domain. The surface organization due to vegetated area scale also has impacts on horizontally averaged scalar and momentum profiles; however, the magnitude in this study is smaller than the results of a previous study using a set of offline surface fluxes as the lower boundary condition for LES. In addition, even though this study only performs a daytime diurnal cycle, the impact of green space scale on cloud distribution in simulations is significant. The cases with more organized green space yield lower‐elevated cumulus cloud and larger‐cloud cover fraction, which impacts the energy budget at the top of boundary layer and, in turn, could lead to additional surface cooling with respect to longer‐term weather and climate. Copyright © 2010 John Wiley & Sons, Ltd.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

The pattern of surface wind flow in Antarctica

Summary The main characteristics of surface winds are tabulated for 34 Antarctic stations. Using these data, supplemented by traverse records, the average wind flow is interpolated for each region and presented as a map showing the pattern of surface wind flow for the whole continent. Attention is focused on the flow in relation to surface contours. Statistics are presented for surface slope, wind speed, temperature, seasonal variations of speed and temperature, diurnal variation (including power spectra) of the wind speed and times of maximum and minimum speed at coastal and inland stations, wind frequency versus direction, the occurrence of calms, the deviation of the plateau wind from the downslope direction, the wind direction near the front of ice shelves, the proportion of cloud cover, and wind chill factors. In all cases data are grouped according to the environs of the stations in an attempt to isolate systematic differences depending on location: coastal stations near the foot of the ice slope and fully exposed to katabatic flow, coastal stations on offshore islands, coastal stations on peninsulas, coastal stations on extensive rock areas, ice shelf stations and inland stations.     

Hydrologic connectivity at the hillslope scale through intra-snowpack flow paths during snowmelt

The seasonal snowmelt period is a critical component of the hydrologic cycle for many mountainous areas. Changes in the timing and rate of snowmelt as a result of physical hydrologic flow paths, such as longitudinal intra-snowpack flow paths, can have strong implications on the partitioning of meltwater amongst streamflow, groundwater recharge, and soil moisture storage. However, intra-snowpack flow paths are highly spatially and temporally variable and thus difficult to observe. This study utilizes new methods to non-destructively observe spatio-temporal changes in the liquid water content of snow in combination with plot experiments to address the research question: What is the scale of influence that intra-snowpack flow paths have on the downslope movement of liquid water during snowmelt across an elevational gradient? This research took place in northern Colorado with study plots spanning from the rain-snow transition zone up to the high alpine. Results indicate an increasing scale of influence from intra-snowpack flow paths with elevation, showing higher hillslope connectivity producing larger intra-snowpack contributing areas for meltwater accumulation, quantified as the upslope contributing area required to produce observed changes in liquid water content from melt rate estimates. The total effective intra-snowpack contributing area of accumulating liquid water was found to be 17, 6, and 0 m² for the above tree line, near tree line, and below tree line plots, respectively. Dye tracer experiments show capillary and permeability barriers result in increased number and thickness of intra-snowpack flow paths at higher elevations. We additionally utilized aerial photogrammetry in combination with ground penetrating radar surveys to investigate the role of this hydrologic process at the small watershed scale. Results here indicate that intra-snowpack flow paths have influence beyond the plot scale, impacting the storage and transmission of liquid water within the snowpack at the small watershed scale.     

Multitemporal terrestrial laser scanning point clouds for thaw subsidence observation at Arctic permafrost monitoring sites

This paper investigates different methods for quantifying thaw subsidence using terrestrial laser scanning (TLS) point clouds. Thaw subsidence is a slow (millimetre to centimetre per year) vertical displacement of the ground surface common in ice-rich permafrost-underlain landscapes. It is difficult to quantify thaw subsidence in tundra areas as they often lack stable reference frames. Also, there is no solid ground surface to serve as a basis for elevation measurements, due to a continuous moss–lichen cover. We investigate how an expert-driven method improves the accuracy of benchmark measurements at discrete locations within two sites using multitemporal TLS data of a 1-year period. Our method aggregates multiple experts’ determination of the ground surface in 3D point clouds, collected in a web-based tool. We then compare this to the performance of a fully automated ground surface determination method. Lastly, we quantify ground surface displacement by directly computing multitemporal point cloud distances, thereby extending thaw subsidence observation to an area-based assessment. Using the expert-driven quantification as reference, we validate the other methods, including in-situ benchmark measurements from a conventional field survey. This study demonstrates that quantifying the ground surface using 3D point clouds is more accurate than the field survey method. The expert-driven method achieves an accuracy of 0.1 ± 0.1 cm. Compared to this, in-situ benchmark measurements by single surveyors yield an accuracy of 0.4 ± 1.5 cm. This difference between the two methods is important, considering an observed displacement of 1.4 cm at the sites. Thaw subsidence quantification with the fully automatic benchmark-based method achieves an accuracy of 0.2 ± 0.5 cm and direct point cloud distance computation an accuracy of 0.2 ± 0.9 cm. The range in accuracy is largely influenced by properties of vegetation structure at locations within the sites. The developed methods enable a link of automated quantification and expert judgement for transparent long-term monitoring of permafrost subsidence. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Verification of an ENSO-Based Long-Range Prediction of Anomalous Weather Conditions During the Vancouver 2010 Olympics and Paralympics

A brief review of the anomalous weather conditions during the Vancouver 2010 Winter Olympic and Paralympic Games and the efforts to predict these anomalies based on some preceding El Niño–Southern Oscillation (ENSO) signals are presented. It is shown that the Olympic Games were held under extraordinarily warm conditions in February 2010, with monthly mean temperature anomalies of +2.2 °C in Vancouver and +2.8 °C in Whistler, ranking respectively as the highest and the second highest in the past 30 years (1981–2010). The warm conditions continued, but became less anomalous, in March 2010 for the Paralympic Games. While the precipitation amounts in the area remained near normal through this winter, the lack of snow due to warm conditions created numerous media headlines and practical problems for the alpine competitions. A statistical model was developed on the premise that February and March temperatures in the Vancouver area could be predicted using an ENSO signal with considerable lead time. This model successfully predicted the warmer-than-normal, lower-snowfall conditions for the Vancouver 2010 Winter Olympics and Paralympics.     

偏振激光雷达探测的北京地区云的统计特性分析

本文给出了CALIPSO卫星偏振激光雷达自2006年6月至2010年3月间探测的北京地区大气云层的统计特性,并根据理论研究结果对激光雷达探测的云进行了分类,当云的温度低于0℃时发现了过冷水的存在.统计结果表明:温度在0℃以上时,大气中的云大部分为水云;温度在0℃至-25℃之间时,大气中会有过冷水和水平导向的片状晶体的冰...     

An improved method for single flow direction calculation in grid digital elevation models

This paper presents improvements to the global D8 (GD8) method for calculating single flow directions in a grid digital elevation model. Flow directions computed from grid digital elevation models serve as the foundation for much of the analysis and modeling of hydrological processes that are driven by topographic gradients. The literature includes both single flow direction methods, where flow goes to only one downslope cell, and multiple flow direction methods that apportion flow among multiple downslope cells. Among single flow direction methods, the standard D8 method, in which the flow direction is set based on the steepest local slope, results in bias on surfaces that do not align with the grid directions. Efforts to address this problem have led to the development of extended methods that account for elevation values further upslope in determining flow directions. We have identified discrepancies in one such method, GD8, and have examined ways to resolve these discrepancies. An improvement to GD8, named iGD8, is presented that allows replacing a reference cell from which path deviations are accumulated and that considers horizontal path deviation rather than global slope as a flow direction criterion. The improved method is found to be effective in resolving the problems encountered with GD8 and to be more efficient than a previously proposed alternative method (least transversal deviation (LTD) based D8, namely D8‐LTD) that uses recursive searching for the largest upstream area when multiple flow paths converge. The proposed improved GD8 method offers the opportunity for improved analysis and modeling of topographically driven hydrological processes by providing better foundational flow directions for these analyses.     

Dead level contours and infiltration pits for risk mitigation in smallholder cropping systems of southern Zimbabwe

The persistent droughts, dry spells, and chronic food insecurity in semi-arid areas necessitate the introduction of more robust rainwater harvesting and soil water management technologies. The study reported here was conducted to assess the influence of dead level contours and infiltration pits on in-field soil water dynamics over two growing seasons. A transect consisting of six access tubes, spaced at 5 m interval, was established across each dead level contour with or without an infiltration pit before the onset of the rains. Two access tubes were installed upslope of the contour while four tubes were installed on the downslope side. Dead level contours with infiltration pits captured more rainwater than dead level contours only resulting in more lateral soil water movement. Significant lateral soil water movement was detected at 3 m downslope following rainfall events of 60–70 mm/day. The 0.2–0.6 m soil layer benefited more from the lateral soil water movement at all the farms. Our results suggest that dead level contours have to be constructed at 3–8 m spacing for crops to benefit from the captured rainwater. It is probably worth exploring strip cropping of food and fodder crops on the downslope of the dead level contours and infiltration pits using the current design of these between-field structures. With the advent of in situ rainwater harvesting techniques included in some conservation agriculture practices it will benefit smallholder cropping systems in semi-arid areas if these between-field structures are promoted concurrently with other sustainable land management systems such as conservation agriculture.     

Representing the Sensitivity of Convective Cloud Systems to Tropospheric Humidity in General Circulation Models

Convective cloud variability on many times scales can be viewed as having three major components: a suppressed phase of shallow and congestus clouds, a disturbed phase of deep convective clouds, and a mature phase of transition to stratiform upper-level clouds. Cumulus parameterization development has focused primarily on the second phase until recently. Consequently, many parameterizations are not sufficiently sensitive to variations in tropospheric humidity. This shortcoming may affect global climate model simulations of climate sensitivity to external forcings, the continental diurnal cycle of clouds and precipitation, and intraseasonal precipitation variability. The lack of sensitivity can be traced in part to underestimated entrainment of environmental air into rising convective clouds and insufficient evaporation of rain into the environment. As a result, the parameterizations produce deep convection too easily while stabilizing the environment too quickly to allow the effects of convective mesoscale organization to occur. Recent versions of some models have increased their sensitivity to tropospheric humidity and improved some aspects of their variability, but a parameterization of mesoscale organization is still absent from most models. Evidence about the effect of these uncertainties on climate change projections suggests that climate modelers should make improved simulation of high and convective clouds as high a priority as better representations of low clouds.     

A physically based watershed partitioning method

A new contour-based watershed partitioning method, relying on tracking flow lines moving from upslope to downslope, is presented. The model works by using only contour data and consists of three principal steps: 1) Ideal Drainage Network tracking; 2) hillslope and channel elements building; and 3) actual contour-based element network construction. It automatically detects and manages sinks, peaks and saddles and calculates element attributes taking into account the non-planar shape of the elements. Watershed, subwatershed and draining area at any point can be automatically detected by using network elements interconnectivity. It is also possible to determine the Actual Drainage Network.