scholarly journals A reduced-order modeling approach to represent subgrid-scale hydrological dynamics for land-surface simulations: application in a polygonal tundra landscape

2014 ◽  
Vol 7 (5) ◽  
pp. 2091-2105 ◽  
Author(s):  
G. S. H. Pau ◽  
G. Bisht ◽  
W. J. Riley

Abstract. Existing land surface models (LSMs) describe physical and biological processes that occur over a wide range of spatial and temporal scales. For example, biogeochemical and hydrological processes responsible for carbon (CO2, CH4) exchanges with the atmosphere range from the molecular scale (pore-scale O2 consumption) to tens of kilometers (vegetation distribution, river networks). Additionally, many processes within LSMs are nonlinearly coupled (e.g., methane production and soil moisture dynamics), and therefore simple linear upscaling techniques can result in large prediction error. In this paper we applied a reduced-order modeling (ROM) technique known as "proper orthogonal decomposition mapping method" that reconstructs temporally resolved fine-resolution solutions based on coarse-resolution solutions. We developed four different methods and applied them to four study sites in a polygonal tundra landscape near Barrow, Alaska. Coupled surface–subsurface isothermal simulations were performed for summer months (June–September) at fine (0.25 m) and coarse (8 m) horizontal resolutions. We used simulation results from three summer seasons (1998–2000) to build ROMs of the 4-D soil moisture field for the study sites individually (single-site) and aggregated (multi-site). The results indicate that the ROM produced a significant computational speedup (> 103) with very small relative approximation error (< 0.1%) for 2 validation years not used in training the ROM. We also demonstrate that our approach: (1) efficiently corrects for coarse-resolution model bias and (2) can be used for polygonal tundra sites not included in the training data set with relatively good accuracy (< 1.7% relative error), thereby allowing for the possibility of applying these ROMs across a much larger landscape. By coupling the ROMs constructed at different scales together hierarchically, this method has the potential to efficiently increase the resolution of land models for coupled climate simulations to spatial scales consistent with mechanistic physical process representation.

2014 ◽  
Vol 7 (2) ◽  
pp. 2125-2172 ◽  
Author(s):  
G. S. H. Pau ◽  
G. Bisht ◽  
W. J. Riley

Abstract. Existing land surface models (LSMs) describe physical and biological processes that occur over a wide range of spatial and temporal scales. For example, biogeochemical and hydrological processes responsible for carbon (CO2, CH4) exchanges with the atmosphere range from molecular scale (pore-scale O2 consumption) to tens of kilometer scale (vegetation distribution, river networks). Additionally, many processes within LSMs are nonlinearly coupled (e.g., methane production and soil moisture dynamics), and therefore simple linear upscaling techniques can result in large prediction error. In this paper we applied a particular reduced-order modeling (ROM) technique known as "Proper Orthogonal Decomposition mapping method" that reconstructs temporally-resolved fine-resolution solutions based on coarse-resolution solutions. We applied this technique to four study sites in a polygonal tundra landscape near Barrow, Alaska. Coupled surface-subsurface isothermal simulations were performed for summer months (June–September) at fine (0.25 m) and coarse (8 m) horizontal resolutions. We used simulation results from three summer seasons (1998–2000) to build ROMs of the 4-D soil moisture field for the four study sites individually (single-site) and aggregated (multi-site). The results indicate that the ROM produced a significant computational speedup (> 103) with very small relative approximation error (< 0.1%) for two validation years not used in training the ROM. We also demonstrated that our approach: (1) efficiently corrects for coarse-resolution model bias and (2) can be used for polygonal tundra sites not included in the training dataset with relatively good accuracy (< 1.5% relative error), thereby allowing for the possibility of applying these ROMs across a much larger landscape. This method has the potential to efficiently increase the resolution of land models for coupled climate simulations, allowing LSMs to be used at spatial scales consistent with mechanistic physical process representation.


2008 ◽  
Vol 88 (5) ◽  
pp. 761-774 ◽  
Author(s):  
J. A. P. Pollacco

Hydrological models require the determination of fitting parameters that are tedious and time consuming to acquire. A rapid alternative method of estimating the fitting parameters is to use pedotransfer functions. This paper proposes a reliable method to estimate soil moisture at -33 and -1500 kPa from soil texture and bulk density. This method reduces the saturated moisture content by multiplying it with two non-linear functions depending on sand and clay contents. The novel pedotransfer function has no restrictions on the range of the texture predictors and gives reasonable predictions for soils with bulk density that varies from 0.25 to 2.16 g cm-3. These pedotransfer functions require only five parameters for each pressure head. It is generally accepted that the introduction of organic matter as a predictor improves the outcomes; however it was found by using a porosity based pedotransfer model, using organic matter as a predictor only modestly improves the accuracy. The model was developed employing 18 559 samples from the IGBP-DIS soil data set for pedotransfer function development (Data and Information System of the International Geosphere Biosphere Programme) database that embodies all major soils across the United States of America. The function is reliable and performs well for a wide range of soils occurring in very dry to very wet climates. Climatical grouping of the IGBP-DIS soils was proposed (aquic, tropical, cryic, aridic), but the results show that only tropical soils require specific grouping. Among many other different non-climatical soil groups tested, only humic and vitric soils were found to require specific grouping. The reliability of the pedotransfer function was further demonstrated with an independent database from Northern Italy having heterogeneous soils, and was found to be comparable or better than the accuracy of other pedotransfer functions found in the literature. Key words: Pedotransfer functions, soil moisture, soil texture, bulk density, organic matter, grouping


2013 ◽  
Vol 17 (7) ◽  
pp. 2781-2796 ◽  
Author(s):  
S. Shukla ◽  
J. Sheffield ◽  
E. F. Wood ◽  
D. P. Lettenmaier

Abstract. Global seasonal hydrologic prediction is crucial to mitigating the impacts of droughts and floods, especially in the developing world. Hydrologic predictability at seasonal lead times (i.e., 1–6 months) comes from knowledge of initial hydrologic conditions (IHCs) and seasonal climate forecast skill (FS). In this study we quantify the contributions of two primary components of IHCs – soil moisture and snow water content – and FS (of precipitation and temperature) to seasonal hydrologic predictability globally on a relative basis throughout the year. We do so by conducting two model-based experiments using the variable infiltration capacity (VIC) macroscale hydrology model, one based on ensemble streamflow prediction (ESP) and another based on Reverse-ESP (Rev-ESP), both for a 47 yr re-forecast period (1961–2007). We compare cumulative runoff (CR), soil moisture (SM) and snow water equivalent (SWE) forecasts from each experiment with a VIC model-based reference data set (generated using observed atmospheric forcings) and estimate the ratio of root mean square error (RMSE) of both experiments for each forecast initialization date and lead time, to determine the relative contribution of IHCs and FS to the seasonal hydrologic predictability. We find that in general, the contributions of IHCs to seasonal hydrologic predictability is highest in the arid and snow-dominated climate (high latitude) regions of the Northern Hemisphere during forecast periods starting on 1 January and 1 October. In mid-latitude regions, such as the Western US, the influence of IHCs is greatest during the forecast period starting on 1 April. In the arid and warm temperate dry winter regions of the Southern Hemisphere, the IHCs dominate during forecast periods starting on 1 April and 1 July. In equatorial humid and monsoonal climate regions, the contribution of FS is generally higher than IHCs through most of the year. Based on our findings, we argue that despite the limited FS (mainly for precipitation) better estimates of the IHCs could lead to improvement in the current level of seasonal hydrologic forecast skill over many regions of the globe at least during some parts of the year.


2018 ◽  
Vol 22 (9) ◽  
pp. 4649-4665 ◽  
Author(s):  
Anouk I. Gevaert ◽  
Ted I. E. Veldkamp ◽  
Philip J. Ward

Abstract. Drought is a natural hazard that occurs at many temporal and spatial scales and has severe environmental and socioeconomic impacts across the globe. The impacts of drought change as drought evolves from precipitation deficits to deficits in soil moisture or streamflow. Here, we quantified the time taken for drought to propagate from meteorological drought to soil moisture drought and from meteorological drought to hydrological drought. We did this by cross-correlating the Standardized Precipitation Index (SPI) against standardized indices (SIs) of soil moisture, runoff, and streamflow from an ensemble of global hydrological models (GHMs) forced by a consistent meteorological dataset. Drought propagation is strongly related to climate types, occurring at sub-seasonal timescales in tropical climates and at up to multi-annual timescales in continental and arid climates. Winter droughts are usually related to longer SPI accumulation periods than summer droughts, especially in continental and tropical savanna climates. The difference between the seasons is likely due to winter snow cover in the former and distinct wet and dry seasons in the latter. Model structure appears to play an important role in model variability, as drought propagation to soil moisture drought is slower in land surface models (LSMs) than in global hydrological models, but propagation to hydrological drought is faster in land surface models than in global hydrological models. The propagation time from SPI to hydrological drought in the models was evaluated against observed data at 127 in situ streamflow stations. On average, errors between observed and modeled drought propagation timescales are small and the model ensemble mean is preferred over the use of a single model. Nevertheless, there is ample opportunity for improvement as substantial differences in drought propagation are found at 10 % of the study sites. A better understanding and representation of drought propagation in models may help improve seasonal drought forecasting as well as constrain drought variability under future climate scenarios.


2018 ◽  
Vol 22 (10) ◽  
pp. 5341-5356 ◽  
Author(s):  
Seyed Hamed Alemohammad ◽  
Jana Kolassa ◽  
Catherine Prigent ◽  
Filipe Aires ◽  
Pierre Gentine

Abstract. Characterizing soil moisture at spatiotemporal scales relevant to land surface processes (i.e., of the order of 1 km) is necessary in order to quantify its role in regional feedbacks between the land surface and the atmospheric boundary layer. Moreover, several applications such as agricultural management can benefit from soil moisture information at fine spatial scales. Soil moisture estimates from current satellite missions have a reasonably good temporal revisit over the globe (2–3-day repeat time); however, their finest spatial resolution is 9 km. NASA's Soil Moisture Active Passive (SMAP) satellite has estimated soil moisture at two different spatial scales of 36 and 9 km since April 2015. In this study, we develop a neural-network-based downscaling algorithm using SMAP observations and disaggregate soil moisture to 2.25 km spatial resolution. Our approach uses the mean monthly Normalized Differenced Vegetation Index (NDVI) as ancillary data to quantify the subpixel heterogeneity of soil moisture. Evaluation of the downscaled soil moisture estimates against in situ observations shows that their accuracy is better than or equal to the SMAP 9 km soil moisture estimates.


2019 ◽  
Author(s):  
Bouchra Ait Hssaine ◽  
Olivier Merlin ◽  
Jamal Ezzahar ◽  
Nitu Ojha ◽  
Salah Er-raki ◽  
...  

Abstract. Thermal-based two-source energy balance modeling is very useful for estimating the land evapotranspiration (ET) at a wide range of spatial and temporal scales. However, the land surface temperature (LST) is not sufficient for constraining simultaneously both soil and vegetation flux components in such a way that assumptions (on either the soil or the vegetation fluxes) are commonly required. To avoid such assumptions, a new energy balance model (TSEB-SM) was recently developed in Ait Hssaine et al. (2018a) to integrate the microwave-derived near-surface soil moisture (SM), in addition to the thermal-derived LST and vegetation cover fraction (fc). Whereas, TSEB-SM has been recently tested using in-situ measurements, the objective of this paper is to evaluate the performance of TSEB-SM in real-life using 1 km resolution MODIS (Moderate resolution imaging spectroradiometer) LST and fc data and the 1 km resolution SM data disaggregated from SMOS (Soil Moisture and Ocean Salinity) observations by using DisPATCh. The approach is applied during a four-year period (2014–2018) over a rainfed wheat field in the Tensift basin, central Morocco, during a four-year period (2014–2018). The field was seeded for the 2014–2015 (S1), 2016–2017 (S2) and 2017–2018 (S3) agricultural season, while it was not ploughed (remained as bare soil) during the 2015–2016 (B1) agricultural season. The mean retrieved values of (arss, brss) calculated for the entire study period using satellite data are (7.32, 4.58). The daily calibrated αPT ranges between 0 and 1.38 for both S1 and S2. Its temporal variability is mainly attributed to the rainfall distribution along the agricultural season. For S3, the daily retrieved αPT remains at a mostly constant value (∼ 0.7) throughout the study period, because of the lack of clear sky disaggregated SM and LST observations during this season. Compared to eddy covariance measurements, TSEB driven only by LST and fc data significantly overestimates latent heat fluxes for the four seasons. The overall mean bias values are 119, 94, 128 and 181 W/m2 for S1, S2, S3 and B1 respectively. In contrast, these errors are much reduced when using TSEB-SM (SM and LST combined data) with the mean bias values estimated as 39, 4, 7 and 62 W/m2 for S1, S2, S3 and B1 respectively.


2019 ◽  
Vol 11 (1) ◽  
pp. 101-110 ◽  
Author(s):  
James W. Roche ◽  
Robert Rice ◽  
Xiande Meng ◽  
Daniel R. Cayan ◽  
Michael D. Dettinger ◽  
...  

Abstract. We present hourly climate data to force land surface process models and assessments over the Merced and Tuolumne watersheds in the Sierra Nevada, California, for the water year 2010–2014 period. Climate data (38 stations) include temperature and humidity (23), precipitation (13), solar radiation (8), and wind speed and direction (8), spanning an elevation range of 333 to 2987 m. Each data set contains raw data as obtained from the source (Level 0), data that are serially continuous with noise and nonphysical points removed (Level 1), and, where possible, data that are gap filled using linear interpolation or regression with a nearby station record (Level 2). All stations chosen for this data set were known or documented to be regularly maintained and components checked and calibrated during the period. Additional time-series data included are available snow water equivalent records from automated stations (8) and manual snow courses (22), as well as distributed snow depth and co-located soil moisture measurements (2–6) from four locations spanning the rain–snow transition zone in the center of the domain. Spatial data layers pertinent to snowpack modeling in this data set are basin polygons and 100 m resolution rasters of elevation, vegetation type, forest canopy cover, tree height, transmissivity, and extinction coefficient. All data are available from online data repositories (https://doi.org/10.6071/M3FH3D).


2010 ◽  
Vol 2 (2) ◽  
Author(s):  
Diandong Ren

AbstractBased on a 2-layer land surface model, a rather general variational data assimilation framework for estimating model state variables is developed. The method minimizes the error of surface soil temperature predictions subject to constraints imposed by the prediction model. Retrieval experiments for soil prognostic variables are performed and the results verified against model simulated data as well as real observations for the Oklahoma Atmospheric Surface layer Instrumentation System (OASIS). The optimization scheme is robust with respect to a wide range of initial guess errors in surface soil temperature (as large as 30 K) and deep soil moisture (within the range between wilting point and saturation). When assimilating OASIS data, the scheme can reduce the initial guess error by more than 90%, while for Observing Simulation System Experiments (OSSEs), the initial guess error is usually reduced by over four orders of magnitude.Using synthetic data, the robustness of the retrieval scheme as related to information content of the data and the physical meaning of the adjoint variables and their use in sensitivity studies are investigated. Through sensitivity analysis, it is confirmed that the vegetation coverage and growth condition determine whether or not the optimally estimated initial soil moisture condition leads to an optimal estimation of the surface fluxes. This reconciles two recent studies.With the real data experiments, it is shown that observations during the daytime period are the most effective for the retrieval. Longer assimilation windows result in more accurate initial condition retrieval, underlining the importance of information quantity, especially for schemes assimilating noisy observations.


2021 ◽  
Author(s):  
Maria Piles ◽  
Roberto Fernandez-Moran ◽  
Luis Gómez-Chova ◽  
Gustau Camps-Valls ◽  
Dara Entekhabi ◽  
...  

&lt;p&gt;The Copernicus Imaging Microwave Radiometer (CIMR) mission is currently being developed as a High Priority Copernicus Mission to support the Integrated European Policy for the Arctic. Due to its measurement characteristics, CIMR has exciting capabilities to enable a unique set of land surface products and science applications at a global scale. These characteristics go beyond what previous microwave radiometers (e.g. AMSR series, SMAP and SMOS) provide, and therefore allow for entirely new approaches to the estimation of bio-geophysical products from brightness temperature observations. Most notably, CIMR channels (L-,C-,X-,Ka-,Ku-bands) are very well fit for the simultaneous retrieval of soil moisture and vegetation properties, like biomass and moisture of different plant components such as leaves, stems or trunks. Also, the distinct spatial resolution of each frequency band allows for the development of approaches to cascade information and obtain these properties at multiple spatial scales. From a temporal perspective, CIMR has a higher revisit time than previous L-band missions dedicated to soil moisture monitoring (about 1 day global, sub-daily at the poles). This improved temporal resolution could allow resolving critical time scales of water processes, which is relevant to better model and understand land-atmosphere exchanges and feedbacks. In this presentation, new opportunities for soil moisture remote sensing made possible by the CIMR mission, as well as synergies and cross-sensor opportunities will be discussed.&amp;#160;&amp;#160;&lt;/p&gt;


Sign in / Sign up

Export Citation Format

Share Document