PUBLICATIONS
Journal of Geophysical Research: Atmospheres
RESEARCH ARTICLE Influence of lateral subsurface flow and connectivity 10.1002/2015JD024067 on soil water storage in land surface modeling Key Points: Jonggun Kim1 and Binayak P. Mohanty1 • Investigate the influence of lateral fl subsurface ow connectivity on soil 1Department of Biological and Agricultural Engineering, Texas A&M University, College Station, Texas, USA water storage • Demonstrate subsurface flow variability effectively using lateral fl fi connectivity Abstract Lateral surface/subsurface ow and their connectivity play a signi cant role in redistributing soil • Develop a connectivity-based lateral water, which has a direct effect on biological, chemical, and geomorphological processes in the root zone subsurface flow algorithm in complex (~1 m). However, most of the land surface models neglect the horizontal exchanges of water at the grid or landscapes subgrid scales, focusing only on the vertical exchanges of water as one-dimensional process. To develop better hydrologic understanding and modeling capability in complex landscapes, in this study we added connectivity-based lateral subsurface flow algorithms in the Community Land Model. To demonstrate the Correspondence to: impact of lateral flow and connectivity on soil water storage we designed three cases including the following: B. P. Mohanty, [email protected] (1) with complex surface topography only, (2) with complex surface topography in upper soil layers and soil hydraulic properties with uniform anisotropy. and (3) with complex surface topography and soil hydraulic properties with spatially varying anisotropy. The connectivity was considered as an indicator for the variation Citation: Kim, J., and B. P. Mohanty (2016), of anisotropy in the case 3, which was created by wetness conditions or geophysical controls (e.g., soil type, Influence of lateral subsurface flow and normalized difference vegetation index, and topographic index). These cases were tested in two study sites (ER connectivity on soil water storage in 5 field and ER-sub watershed in Oklahoma) comparing to the field (gravimetric and remote sensing) soil land surface modeling, J. Geophys. Res. Atmos., 121, doi:10.1002/2015JD024067. moisture observations. Through the analysis of spatial patterns and temporal dynamics of soil moisture predictions from the study cases, surface topography was found to be a crucial control in demonstrating the Received 19 AUG 2015 variation of near surface soil moisture, but not significantly affected the subsurface flow in deeper soil layers. In Accepted 19 DEC 2015 addition, we observed the best performance in case 3 representing that the lateral connectivity can Accepted article online 23 DEC 2015 contribute effectively to quantify the anisotropy and redistributing soil water in the root zone. Hence, the approach with connectivity-based lateral subsurface flow was able to better characterize the spatially distributed patterns of subsurface flow and improve the simulation of the hydrologic cycle.
1. Introduction Lateral surface/subsurface flow is an important hydrologic process and a key component of the water budget. Through its direct impacts on soil moisture, it can affect water and energy fluxes at the land surface and influence the regional climate and water cycle [Gochis and Chen, 2003; Kumar, 2004]. Further, the lat- eral flow and its connectivity play significant role in redistributing soil water, which have a direct effect on biological, chemical, and geomorphologic processes in the root zone [Lu et al.,2011;Western et al., 2001]. In spite of the importance of lateral flow, most of the land surface models (LSMs: Community Land Model (CLM), Noah Land Surface Model (Noah LSM), Variable Infiltration Capacity, etc.) neglect the horizontal exchanges of water at the grid or subgrid scales, focusing only on the vertical exchanges of water as a one-dimensional process. Surface routing models (e.g., River Transport Model, RTM) are already included to reflect the lateral movement of surface water in land surface modeling, but the lateral subsurface flow is excluded because the models generally assume that lateral transfers of subsurface moisture are fairly marginal in soil water budgets of a regional scale. Recently, 3-D hydrological surface-subsurface models were developed by coupling LSMs with distributed hydrological models to account for interactions between atmospheric, hydrological, and ecological processes (CATHY/NoahMP [Niu et al., 2014] and PARFLOW/CLM [Maxwell and Miller, 2005]). Although these hydrological models include a process for the lateral subsurface flow, it still has limitations for deriving lateral hydraulic parameters (e.g., lateral hydraulic conductivity) that might be related to connected patterns of subsurface properties. Furthermore, spatial variability of soil moisture in the unsaturated zone cannot be described successfully without relevant understanding of how the subsurface flow is distributed or connected vertically or laterally in complex landscapes [Hatton, 1998; Zhang et al., 1999; Jana and Mohanty, 2012a, 2012b, 2012c; Shen et al., 2013].
©2015. American Geophysical Union. More realistic understanding of surface and subsurface water movement at large scales can be also All Rights Reserved. resolved through a hyperresolution land surface modeling that allows for better representation of spatially
KIM AND MOHANTY LATERAL SUBSURFACE FLOW AND CONNECTIVITY 1 Journal of Geophysical Research: Atmospheres 10.1002/2015JD024067
heterogeneous land surfaces [Wood et al., 2011]. Thus, the lateral subsurface flow should be accounted for in hydrological modeling, characterizing vertical and lateral flow components effectively in the unsatu- rated zone. Various studies have been conducted to account for the lateral flow in the unsaturated soil. Zaslavsky and Sinai [1981] explained a theory of unsaturated lateral flow with the major causes such as soil surface slope, anisotropy, and layering. Famiglietti and Wood [1994] developed a land surface modeling approach based on the TOPMODEL framework to address the impact of topographic configuration on soil moisture heteroge- neity at a watershed scale. They showed a significant role of the topographic control in development of soil moisture heterogeneity and improved the simulation of hydrologic cycle using the modeling approach. Chen and Kumar [2001] explored the role of the topographic control in the seasonal and interannual variations of energy and water balances using statistical moments of topographic wetness indices and observed an improvement of streamflow predictions. Gravity and gradients in matric potential are also critical mechan- isms in the unsaturated zone, causing soil water movements from high to low potential [McCord and Stephens, 1987; Jana and Mohanty, 2012a, 2012b, 2012c]. Water moving vertically through a heterogeneous soil profile can be influenced by the heterogeneity of soil hydraulic properties between soil layers, which can cause lateral flow at the interface [Zhu and Lin, 2009]. In process-based Soil-Vegetation-Atmosphere Transfer models, soil hydraulic properties (e.g., saturated soil water content, soil matric potential, and saturated hydraulic conductivity) are critical inputs to account for water movement in soil. The soil hydraulic proper- ties are normally derived using several empirical equations (e.g., van Genuchten, Cosby, and Clapp and Hornberger) according to soil texture. Among the soil properties, an estimation of lateral hydraulic conduc- tivity is more challenging because of the lack of available information. Thus, anisotropy has been used to derive the lateral hydraulic conductivities from the relationship between vertical and lateral permeability because soil behaves as an anisotropic medium which can cause lateral subsurface flow [Zaslavsky and Sinai, 1981; Wang et al., 2011]. In the previous studies related to soil anisotropy, statistical or empirical anisotropy ratios were used at various scales [Chen and Kumar, 2001; Kumar, 2004; Assouline and Or, 2006; Maxwell and Kollet, 2008]. However, available experimental data and information for the anisotropy ratio in unsaturated soils might be limited to be applied in heterogeneous landscapes of large land areas. In order to overcome the limitations, the anisotropy ratio can be derived by spatially distributed patterns of wetness condition or its dominant physical controls such as soil texture, vegetation (NDVI), and topographic index (TI) to characterize the spatial pattern of subsurface flow in the unsaturated zone [Chen and Kumar, 2001]. A hydrologic connectivity has been proposed to address not only hydrologic flow paths but also spatial pat- terns of soil moisture variability at a catchment scale [Western et al., 2001; Hwang et al., 2009; Gaur and Mohanty, 2013]. The lateral connectivity is critically important for representing connected pathways of runoff in the landscapes and understanding movements of surface/subsurface flow [Mueller et al., 2007; Smith et al., 2010]. Jencso et al. [2009] derived hydrologic connectivity between catchment landscapes and channel network to identify runoff source areas based on the topographic characteristics. Hwang et al. [2012] found significant relationships between annual hydrologic metrics (e.g., runoff and evapotranspiration (ET)) and hydrologic vegetation gradient used as an indicator for lateral hydrologic connectivity at a watershed scale. Lateral subsurface flow connectivity can be derived from spatially distributed patterns of wetness condition or dominant physical factors and used to quantify the spatially varied anisotropy ratios in heterogeneous landscapes. In this study, we explored the influences of lateral subsurface flow and its connectivity on soil water storage in the unsaturated zone using a land surface model (Community Land Model (CLM)). None of previous studies have considered spatially varying anisotropy ratios derived from lateral connectivity to consider the lateral subsurface flow in hydrological modeling. Thus, the objectives of this study are (1) to develop better hydrologic understanding and modeling cap- ability in complex landscapes using a connectivity-based lateral subsurface flow algorithm and (2) to demonstrate the subsurface flow variability effectively using spatially distributed patterns of root zone wetnessconditionsanditsphysicalcontrolsatfield and subwatershed scales. Although this study was focused on smaller-scale hydrological processes compared to large-scale climate models, it still can pro- vide insights for large-scale land surface modeling to enhance their capability. In this study, the concept of lateral flow was used for the unsaturated zone that can be governed by topography and gradients in matric potential.
KIM AND MOHANTY LATERAL SUBSURFACE FLOW AND CONNECTIVITY 2 Journal of Geophysical Research: Atmospheres 10.1002/2015JD024067
Figure 1. Study sites for (a) El Reno 5 (ER 5) matching the ESTAR remote sensing footprint with multidepth ground-based soil water measurements using truck- mounted Giddings probe (100 m spacing) and (b) El Reno subwatershed (ER-sub) in Oklahoma.
2. Methodology 2.1. Study Area El-Reno site 5 (ER 5: field scale) and El-Reno subwatershed (ER-sub: subwatershed scale) located in North Canadian River basin in Oklahoma were selected to evaluate the proposed approach in this study (Figure 1). The ER 5 site (area: 0.8 km × 0.8 km) is located within the ER-sub boundary (area: 27 km2). These sites have a sub- humid climate with an average annual rainfall of approximately 805 mm. Daily-mean maximum temperature is 34°C in July with annual-mean temperature of 15°C. The topography of the ER 5 is generally flat with average slopes less than 4.0%, while the ER-sub site has a variety of slopes from 11.0% to 0.001%. The ER 5 site has a native grass with 1 m root depth and mostly silty loam across the study domain. Vegetation in the ER-sub ranges from short and tall grasses (predominant) and forest in the north and central area to cropland in the south. Various soil types (e.g., silty loam (dominant), loam, and clay loam) are represented across the region. Our proposed approach was validated with daily in situ soil moisture (49 sampling points) measured in top 5 cm soil (18 June to 17 July) and in depths of 0–15, 15–30, 30–45, 45–60, and 60–90 cm (6–15 July) during the Southern Great Plains experiment 1997 (SGP97) [Mohanty et al., 2002] for the ER 5 site. Using a truck mounted Giddings probe, soil samples between the land surface and 90 cm depth were collected on a 7 × 7 square sampling grid (100 m spacing between sampling points) across the ER 5 field (Figure 1a). For the ER-sub site, we validated model predictions with Electronically Scanning Thin Array Radiometer (ESTAR) pixel-based (800 × 800 m) near surface soil moisture products [Jackson et al., 1999] obtained during Southern Great Plains Experiment 1997, SGP97 (18 June to 17 July) (Figure 1b).
KIM AND MOHANTY LATERAL SUBSURFACE FLOW AND CONNECTIVITY 3 Journal of Geophysical Research: Atmospheres 10.1002/2015JD024067
2.2. Description of Model Condition and Forcing Data Community Land Model (CLM) [Oleson et al., 2010] serves as the dynamic land surface model component of Community Earth System Model (CESM) [Oleson et al., 2010], which consists of various processes such as biogeophysics, hydrologic cycle, biogeochemistry, and dynamic vegetation. The model can be run in off-line mode with prescribed forcing data or in a mode fully coupled to CESM with output from Community Atmosphere Model [Collins et al., 2006], which is the atmospheric component of CESM. CLM simulates surface and subsurface runoff based on the simple TOPMODEL-based runoff model (SIMTOP) [Niu et al., 2005]. The model considers water table dynamics as the lower boundary using the SIMple Groundwater Model (SIMGM) [Niu et al., 2007]. Bare soil evaporation is simulated based on the Philip and De Vries [1957] diffusion model, and transpiration process uses an aerodynamic approach based on the Biosphere Atmosphere Transfer Scheme model [Dickinson et al., 1993] and a stomatal resistance from the LSM model [Bonan, 1996]. River Transport Model (RTM) is coupled to CLM for the runoff routing process over a domain [Oleson et al., 2010]. In this study, we used CLM4.0 and ran the model with RTM in off-line mode. The soil column in CLM consists of 10 soil layers with the thickness of 1.75, 2.76, 4.55, 7.5, 12.36, 20.38, 33.60, 55.39, 91.33, and 113.7 cm (total depth of 343 cm). Soil water flow in CLM is simulated by the modified one-dimensional (1-D) Richards’ equation [Zeng and Decker, 2009]. CLM has been enhanced to improve hydrological cycle (water balance), vegetation dynamics, and computational performance in the last decade. Nevertheless, the model still simplify the complex processes for the root zone soil hydrology considering only vertical flow using 1-D Richards’s equation. In this study, we modified soil water flow process including a lateral flow component in the unsaturated zone to improve the model perfor- mance (as described in section 2.3). We run the model in off-line mode with atmospheric forcing data (precipitation, temperature, specific humidity, wind speed, surface air pressure, and solar radiation) collected from North American Land Data Assimilation System (NLDAS), which were applied uniformly for the study sites. In this study, we gen- erated model input at spatial resolutions of 50 m and 100 m for the ER 5 site and the ER-sub, respectively. As required input data sets, land cover, soil types with depth, and topographic information were obtained from NLCD (National Land Cover Database), SSURGO (Soil Survey Geographic database), and NED (National Elevation Dataset), respectively. The bottom boundary condition of the model is decided with the water table dynamics calculated from aquifer water storage via the SIMGM [Niu and Yang, 2007], and then the model performed a spinning up to initialize the soil profile for the initial condition. In CLM, soil hydraulic properties are determined based on percentages of clay and sand using an empirical equation developed by Clapp and Hornberger [1978]. However, CLM tend to simulate the soil moisture lower than the observations in this study because the parameters estimated from the model input (per- centages of clay and sand) for the ER 5 site were deviated from the referenced parameter ranges (Clapp and Hornberger table) of silty loam soil (predominant in the ER 5 site). Thus, we adjusted the para- meters (trial and error) to satisfy the possible ranges of parameters and applied in CLM and modified CLM (section 2.3). 2.3. Lateral Subsurface Flow Process CLM (based on one-dimensional simulation) assumes that soil water drains only vertically to the water table, and there are no interactions between parallel soil columns. To improve the simplified subsurface flow process in the unsaturated zone by CLM, we modified the one-dimensional vertical soil water flow with three- dimensional flow based on Richards’s equation to consider the lateral subsurface flow in the model (Figure 1). The three-dimensional water flow can be expressed as follows,