Modelling Groundwater Flow with MIKE SHE Using Conventional Climate Data and Satellite Data As Model Forcing in Haihe Plain, China
Total Page:16
File Type:pdf, Size:1020Kb
water Article Modelling Groundwater Flow with MIKE SHE Using Conventional Climate Data and Satellite Data as Model Forcing in Haihe Plain, China Yunqiao Shu, Hongjun Li and Yuping Lei * Key Laboratory of Agricultural Water Resources, Hebei Key Laboratory of Agricultural Water-Saving, Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050021, China; [email protected] (Y.S.); [email protected] (H.L.) * Correspondence: [email protected]; Tel.: +86-311-8588-5030 Received: 3 July 2018; Accepted: 17 September 2018; Published: 20 September 2018 Abstract: In North China Plain, accurate spatial and temporal ET and precipitation pattern is very important in the groundwater resource assessment. This study demonstrated the potential for modelling ET and groundwater processes using remote sensing data for distributed hydrological modelling with MIKE SHE codes in the Haihe Plain, China. The model was successfully validated against independent groundwater level measurements following the calibration period and the model also provided a reasonable match of the lysimeter measurements of ET. The remote sensing data included ET derived from global radiation products of Fengyun-2C geostationary meteorological satellite (FY-2C) and FY-2C precipitation products. The comparisons show that precipitation is a critical factor for the hydrological response and for the spatial distribution of ET and groundwater flow. FY-2C precipitation products has a spatial resolution of about 11 km, which thus adds more spatial variability to the most important driving variable. The ET map based on FY-2C data has a higher spatial variability than that map based on conventional data, which are caused by higher resolution of ground information. The groundwater level changes in the aquifer system are shown in the quite different spatial patterns under two models, which is affected by the significant difference between two types of precipitation. In the Haihe Plain, accurate spatial and temporal ET pattern is very important in the groundwater resource assessment that determines the recharge to the saturated zone. Keywords: distributed hydrological model; remote sensing; evapotranspiration; groundwater; spatial pattern 1. Introduction In general, the surface and groundwater resources and their spatio-temporal variation and use are vaguely known which hinders proper water resources managements. These applications require an accurate prediction of hydrological variables spatially and temporally. Distributed hydrological models can provide the desired variables at required resolution [1,2]. MIKE SHE is a distributed and physically-based hydrological model developed by DHI Group, Inc., which delivers integrated modelling of groundwater, surface water, recharge [3] and evapotranspiration and has been widely used to study various water resource management problem [4–7]. We applied MIKE SHE to the part area of the alluvial plain of Mt. Taihang with an intensive crop growing, China to examine the components of the groundwater balance and to assess the impacts of crop water use and crop rotations on groundwater regime. However, the distributed nature of these models not only requires extensive data amounts for driving the models and for the parameterization of the land surface and subsurface Water 2018, 10, 1295; doi:10.3390/w10101295 www.mdpi.com/journal/water Water 2018, 10, x FOR PEER REVIEW 2 of 19 to assess the impacts of crop water use and crop rotations on groundwater regime. However, the distributed nature of these models not only requires extensive data amounts for driving the Water 2018, 10, 1295 2 of 19 models and for the parameterization of the land surface and subsurface but also requires that they should be calibrated and validated on a distributed basis [5,8,9]. Apart from land surface butdata also (such requires as thattopographic, they should soil, be land calibrated use and and vegetation validated data, on a etc.), distributed a hydrological basis [5,8 ,model9]. Apart need from landthe surface input of data meteorological (such as topographic, data including soil, land precipitation use and vegetation and reference data, evapotranspiration etc.), a hydrological (ET). model needAn the accurate input of representation meteorological of data meteorol includingogical precipitation data in space and and reference time is evapotranspiration critical for reliable (ET). Anpredictions accurate representation of the hydrological of meteorological responses [10] data. in space and time is critical for reliable predictions of the hydrologicalIn order to obtain responses these [10 kinds]. of meteorological data, people are struggling to set up global coverageIn order of to observing obtain these stations kinds making of meteorological surface measurements. data, people However, are struggling there is a to general set up lack global coverageof meteorological of observing data stations in arid makingand semiarid surface regions. measurements. Generally, a However, county (or there a city) is has a general a public lack of meteorologicalmeteorological station data in in arid China and and semiarid the area regions. of a county Generally, or city aranges county approximately (or a city) has 400 a km public2 meteorologicalto 1000 km2 in station Haihe in Plain, China China. and theFigure area 1 ofshows a county the meteorological or city ranges stations approximately that we received 400 km2 to 1000daily km 2datain Haihe from Plain,the China China. Meteorological Figure1 shows Data the meteorologicalService Center. stations Compared that weto this, received RS has daily dataattractive from the advantage China Meteorological of an unprecedented Data Service spatial Center. and temporal Compared coverage to this, of remote critical sensing land sur (RS)face has attractiveand atmospheric advantage ofdata. an unprecedented Therefore RS spatialdata have and temporalbeen extensively coverage ofused critical in landgroundwater surface and atmosphericmodelling, data. especially Therefore in the RS region data have with been a weak extensively observation used infrastructure in groundwater [11 modelling,–14]. RS has especially often in thebeen region used with in ahydrological weak observation modelling infrastructure as driving [11 –data14]. RSor hasfor oftenparameterizat been usedion in hydrological[5,15–19]. modellingFurthermore, as driving actual data ET or foror parameterizationsoil moisture estimated [5,15–19]. from Furthermore, remote-sensing actual ET data or soil provide moisture estimatedopportunities from remote-sensing for assimilation data in hydrological provide opportunities models and for for assimilation filling spatio in-temporal hydrological gaps models for andmodel for filling calibration spatio-temporal and validation gaps for[13,20 model–22] calibration. and validation [13,20–22]. FigureFigure 1. Location 1. Location of studyof study area. area. The The boundary boundary conditions conditions of of the the hydrological hydrological model model for for the the area arearea defined are asdefined follows: as b1:follows: specified b1: gradientspecified along gradient Mt. Taihangalong Mt. ridge, Taihang b2: specified ridge, b2: gradient specified at the south-westerngradient at boundary,the south-western b3: specified boundary, hydraulic b3: specified head along hydraulic the Yellow head river, along b4: the specified Yellow hydraulic river, headb4: (0 specified m) along hydraulic the coast line,head b5: (0 m) no-flow along along the coast a streamline line, b5: towardno-flownorth. along a streamline toward north. In current study, the objective is to examine if remote-sensing data from the Fengyun-2C (FY-2C), a ChineseIn operationalcurrent study, geostationary the objective satellite, is to examine can improve if remote the estimation-sensing data of actual from ET the and Fengyun groundwater-2C level.(FY We-2C), compared a Chinese two operational distributed geostationary models based satellite, on RS inputscan improve and conventional the estimation data of for actual the HaiheET Plain,and China, groundwater where there level. is oneWe of compared the most depleted two distributed aquifers inmodels the world based for cropon RS growing inputs depends and on intensiveconventional irrigation data for [23 the,24 ].Haihe FY-2C Plain, products China, including where there the daily is one precipitation of the most depleted and global aquifers radiation for potential evapotranspiration estimation applied in the distributed hydrological model. Meanwhile, the conventional distributed hydrological model also built up. It is expected that model results different both from distributed conventional hydrological model and RS hydrological model since there are Water 2018, 10, 1295 3 of 19 important difference in deriving variables. In addition, the triangle method that is widely used to estimate the daily ET is introduce. In this study, we developed, calibrated and validated an integrated and distributed hydrological model for the Haihe Plain based on available conventional data and information. Based on this calibrated hydrological model, the geostationary satellite FY-2C products were used as driving force for the second step. We compared the simulations results of the two models spatially and temporally. Further, the simulation results were compared to estimates of daily actual ET under an improved triangle method that is widely used based on MODIS and FY-2C data [25–27]. 2. Materials and