Simulating the Effects of Irrigation over the United States in a Land Surface Model Based on Satellite-Derived Agricultural Data

2010 ◽  
Vol 11 (1) ◽  
pp. 171-184 ◽  
Author(s):  
Mutlu Ozdogan ◽  
Matthew Rodell ◽  
Hiroko Kato Beaudoing ◽  
David L. Toll

Abstract A novel method is introduced for integrating satellite-derived irrigation data and high-resolution crop-type information into a land surface model (LSM). The objective is to improve the simulation of land surface states and fluxes through better representation of agricultural land use. Ultimately, this scheme could enable numerical weather prediction (NWP) models to capture land–atmosphere feedbacks in managed lands more accurately and thus improve forecast skill. Here, it is shown that the application of the new irrigation scheme over the continental United States significantly influences the surface water and energy balances by modulating the partitioning of water between the surface and the atmosphere. In this experiment, irrigation caused a 12% increase in evapotranspiration (QLE) and an equivalent reduction in the sensible heat flux (QH) averaged over all irrigated areas in the continental United States during the 2003 growing season. Local effects were more extreme: irrigation shifted more than 100 W m−2 from QH to QLE in many locations in California, eastern Idaho, southern Washington, and southern Colorado during peak crop growth. In these cases, the changes in ground heat flux (QG), net radiation (RNET), evapotranspiration (ET), runoff (R), and soil moisture (SM) were more than 3 W m−2, 20 W m−2, 5 mm day−1, 0.3 mm day−1, and 100 mm, respectively. These results are highly relevant to continental-to-global-scale water and energy cycle studies that, to date, have struggled to quantify the effects of agricultural management practices such as irrigation. On the basis of the results presented here, it is expected that better representation of managed lands will lead to improved weather and climate forecasting skill when the new irrigation scheme is incorporated into NWP models such as NOAA’s Global Forecast System (GFS).

2010 ◽  
Vol 11 (5) ◽  
pp. 1103-1122 ◽  
Author(s):  
Rolf H. Reichle ◽  
Sujay V. Kumar ◽  
Sarith P. P. Mahanama ◽  
Randal D. Koster ◽  
Q. Liu

Abstract Land surface (or “skin”) temperature (LST) lies at the heart of the surface energy balance and is a key variable in weather and climate models. In this research LST retrievals from the International Satellite Cloud Climatology Project (ISCCP) are assimilated into the Noah land surface model and Catchment land surface model (CLSM) using an ensemble-based, offline land data assimilation system. LST is described very differently in the two models. A priori scaling and dynamic bias estimation approaches are applied because satellite and model LSTs typically exhibit different mean values and variabilities. Performance is measured against 27 months of in situ measurements from the Coordinated Energy and Water Cycle Observations Project at 48 stations. LST estimates from Noah and CLSM without data assimilation (“open loop”) are comparable to each other and superior to ISCCP retrievals. For LST, the RMSE values are 4.9 K (CLSM), 5.5 K (Noah), and 7.6 K (ISCCP), and the anomaly correlation coefficients (R) are 0.61 (CLSM), 0.63 (Noah), and 0.52 (ISCCP). Assimilation of ISCCP retrievals provides modest yet statistically significant improvements (over an open loop, as indicated by nonoverlapping 95% confidence intervals) of up to 0.7 K in RMSE and 0.05 in the anomaly R. The skill of the latent and sensible heat flux estimates from the assimilation integrations is essentially identical to the corresponding open loop skill. Noah assimilation estimates of ground heat flux, however, can be significantly worse than open loop estimates. Provided the assimilation system is properly adapted to each land model, the benefits from the assimilation of LST retrievals are comparable for both models.


2016 ◽  
Author(s):  
Valentijn R. N. Pauwels ◽  
Edoardo Daly

Abstract. It is generally accepted that the ground heat flux accounts for a significant fraction of the surface energy balance. In land surface models, the ground heat flux is estimated through a numerical solution of the heat conduction equation. Recent research has shown that this approach introduces errors in the estimation of the energy balance. In this paper, we calibrate a land surface model using a numerical solution of the heat conduction equation with four different vertical spatial resolutions. It is found that the thermal conductivity is the most sensitive parameter to the spatial resolution. More importantly, the thermal conductivity values are directly related to the spatial resolution, thus rendering any physical interpretation of this value irrelevant. The numerical solution is then replaced by an analytical solution. The results of the numerical and analytical solutions are identical when fine spatial and temporal resolutions are used. However, when using resolutions that are typical for land surface models, significant differences are found. When using the analytical solution, the ground heat flux is directly calculated without calculating the soil temperature profile. The calculation of the temperature at each node in the soil profile is thus no longer required, unless the model contains parameters that depend on the soil temperature, which in this study is not the case. The calibration is repeated, and thermal conductivity values independent of the vertical spatial resolution are obtained. The main conclusion of this study is that care must be taken when interpreting land surface model results that have been obtained using numerical ground heat flux estimates. The use of exact analytical solutions is recommended.


2016 ◽  
Vol 20 (11) ◽  
pp. 4689-4706 ◽  
Author(s):  
Valentijn R. N. Pauwels ◽  
Edoardo Daly

Abstract. It is generally accepted that the ground heat flux accounts for a significant fraction of the surface energy balance. In land surface models, the ground heat flux is typically estimated through a numerical solution of the heat conduction equation. Recent research has shown that this approach introduces errors in the estimation of the energy balance. In this paper, we calibrate a land surface model using a numerical solution of the heat conduction equation with four different vertical spatial resolutions. It is found that the thermal conductivity is the most sensitive parameter to the spatial resolution. More importantly, the thermal conductivity values are directly related to the spatial resolution, thus rendering any physical interpretation of this value irrelevant. The numerical solution is then replaced by an analytical solution. The results of the numerical and analytical solutions are identical when fine spatial and temporal resolutions are used. However, when using resolutions that are typical of land surface models, significant differences are found. When using the analytical solution, the ground heat flux is directly calculated without calculating the soil temperature profile. The calculation of the temperature at each node in the soil profile is thus no longer required, unless the model contains parameters that depend on the soil temperature, which in this study is not the case. The calibration is repeated, and thermal conductivity values independent of the vertical spatial resolution are obtained. The main conclusion of this study is that care must be taken when interpreting land surface model results that have been obtained using numerical ground heat flux estimates. The use of exact analytical solutions, when available, is recommended.


2020 ◽  
Vol 21 (1) ◽  
pp. 59-71 ◽  
Author(s):  
Augusto Getirana ◽  
Matthew Rodell ◽  
Sujay Kumar ◽  
Hiroko Kato Beaudoing ◽  
Kristi Arsenault ◽  
...  

AbstractWe evaluate the impact of Gravity Recovery and Climate Experiment data assimilation (GRACE-DA) on seasonal hydrological forecast initialization over the United States, focusing on groundwater storage. GRACE-based terrestrial water storage (TWS) estimates are assimilated into a land surface model for the 2003–16 period. Three-month hindcast (i.e., forecast of past events) simulations are initialized using states from the reference (no data assimilation) and GRACE-DA runs. Differences between the two initial hydrological condition (IHC) sets are evaluated for two forecast techniques at 305 wells where depth to water table measurements are available. Results show that using GRACE-DA-based IHC improves seasonal groundwater forecast performance in terms of both RMSE and correlation. While most regions show improvement, degradation is common in the High Plains, where withdrawals for irrigation practices affect groundwater variability more strongly than the weather variability, which demonstrates the need for simulating such activities. These findings contribute to recent efforts toward an improved U.S. drought monitoring and forecast system.


2014 ◽  
Vol 7 (1) ◽  
pp. 361-386 ◽  
Author(s):  
D. N. Walters ◽  
K. D. Williams ◽  
I. A. Boutle ◽  
A. C. Bushell ◽  
J. M. Edwards ◽  
...  

Abstract. We describe Global Atmosphere 4.0 (GA4.0) and Global Land 4.0 (GL4.0): configurations of the Met Office Unified Model and JULES (Joint UK Land Environment Simulator) community land surface model developed for use in global and regional climate research and weather prediction activities. GA4.0 and GL4.0 are based on the previous GA3.0 and GL3.0 configurations, with the inclusion of developments made by the Met Office and its collaborators during its annual development cycle. This paper provides a comprehensive technical and scientific description of GA4.0 and GL4.0 as well as details of how these differ from their predecessors. We also present the results of some initial evaluations of their performance. Overall, performance is comparable with that of GA3.0/GL3.0; the updated configurations include improvements to the science of several parametrisation schemes, however, and will form a baseline for further ongoing development.


2020 ◽  
Author(s):  
Benjamin Fersch ◽  
Alfonso Senatore ◽  
Bianca Adler ◽  
Joël Arnault ◽  
Matthias Mauder ◽  
...  

<p>The land surface and the atmospheric boundary layer are closely intertwined with respect to the exchange of water, trace gases and energy. Nonlinear feedback and scale dependent mechanisms are obvious by observations and theories. Modeling instead is often narrowed to single compartments of the terrestrial system or bound to traditional viewpoints of definite scientific disciplines. Coupled terrestrial hydrometeorological modeling systems attempt to overcome these limitations to achieve a better integration of the processes relevant for regional climate studies and local area weather prediction. We examine the ability of the hydrologically enhanced version of the Weather Research and Forecasting Model (WRF-Hydro) to reproduce the regional water cycle by means of a two-way coupled approach and assess the impact of hydrological coupling with respect to a traditional regional atmospheric model setting. It includes the observation-based calibration of the hydrological model component (offline WRF-Hydro) and a comparison of the classic WRF and the fully coupled WRF-Hydro models both with identical calibrated parameter settings for the land surface model (Noah-MP). The simulations are evaluated based on extensive observations at the pre-Alpine Terrestrial Environmental Observatory (TERENO Pre-Alpine) for the Ammer (600 km²) and Rott (55 km²) river catchments in southern Germany, covering a five month period (Jun–Oct 2016).</p><p>The sensitivity of 7 land surface parameters is tested using the <em>Latin-Hypercube One-factor-At-a-Time</em> (LH-OAT) method and 6 sensitive parameters are subsequently optimized for 6 different subcatchments, using the Model-Independent <em>Parameter Estimation and Uncertainty Analysis software</em> (PEST).</p><p>The calibration of the offline WRF-Hydro leads to Nash-Sutcliffe efficiencies between 0.56 and 0.64 and volumetric efficiencies between 0.46 and 0.81 for the six subcatchments. The comparison of classic WRF and fully coupled WRF-Hydro shows only tiny alterations for radiation and precipitation but considerable changes for moisture- and energy fluxes. By comparison with TERENO Pre-Alpine observations, the fully coupled model slightly outperforms the classic WRF with respect to evapotranspiration, sensible and ground heat flux, near surface mixing ratio, temperature, and boundary layer profiles of air temperature. The subcatchment-based water budgets show uniformly directed variations for evapotranspiration, infiltration excess and percolation whereas soil moisture and precipitation change randomly.</p>


2014 ◽  
Vol 15 (3) ◽  
pp. 921-937 ◽  
Author(s):  
Donghai Zheng ◽  
Rogier van der Velde ◽  
Zhongbo Su ◽  
Martijn J. Booij ◽  
Arjen Y. Hoekstra ◽  
...  

ABSTRACT Current land surface models still have difficulties with producing reliable surface heat fluxes and skin temperature (Tsfc) estimates for high-altitude regions, which may be addressed via adequate parameterization of the roughness lengths for momentum (z0m) and heat (z0h) transfer. In this study, the performance of various z0h and z0m schemes developed for the Noah land surface model is assessed for a high-altitude site (3430 m) on the northeastern part of the Tibetan Plateau. Based on the in situ surface heat fluxes and profile measurements of wind and temperature, monthly variations of z0m and diurnal variations of z0h are derived through application of the Monin–Obukhov similarity theory. These derived values together with the measured heat fluxes are utilized to assess the performance of those z0m and z0h schemes for different seasons. The analyses show that the z0m dynamics are related to vegetation dynamics and soil water freeze–thaw state, which are reproduced satisfactorily with current z0m schemes. Further, it is demonstrated that the heat flux simulations are very sensitive to the diurnal variations of z0h. The newly developed z0h schemes all capture, at least over the sparse vegetated surfaces during the winter season, the observed diurnal variability much better than the original one. It should, however, be noted that for the dense vegetated surfaces during the spring and monsoon seasons, not all newly developed schemes perform consistently better than the original one. With the most promising schemes, the Noah simulated sensible heat flux, latent heat flux, Tsfc, and soil temperature improved for the monsoon season by about 29%, 79%, 75%, and 81%, respectively. In addition, the impact of Tsfc calculation and energy balance closure associated with measurement uncertainties on the above findings are discussed, and the selection of the appropriate z0h scheme for applications is addressed.


2018 ◽  
Vol 19 (12) ◽  
pp. 1917-1933 ◽  
Author(s):  
Li Fang ◽  
Xiwu Zhan ◽  
Christopher R. Hain ◽  
Jifu Yin ◽  
Jicheng Liu

Abstract Green vegetation fraction (GVF) plays a crucial role in the atmosphere–land water and energy exchanges. It is one of the essential parameters in the Noah land surface model (LSM) that serves as the land component of a number of operational numerical weather prediction models at the National Centers for Environmental Prediction (NCEP) of NOAA. The satellite GVF products used in NCEP models are derived from a simple linear conversion of either the normalized difference vegetation index (NDVI) from the Advanced Very High Resolution Radiometer (AVHRR) currently or the enhanced vegetation index (EVI) from the Visible Infrared Imaging Radiometer Suite (VIIRS) planned for the near future. Since the NDVI or EVI is a simple spectral index of vegetation cover, GVFs derived from them may lack the biophysical meaning required in the Noah LSM. Moreover, the NDVI- or EVI-based GVF data products may be systematically biased over densely vegetated regions resulting from the saturation issue associated with spectral vegetation indices. On the other hand, the GVF is physically related to the leaf area index (LAI), and thus it could be beneficial to derive GVF from LAI data products. In this paper, the EVI-based and the LAI-based GVF derivation methods are mathematically analyzed and are found to be significantly different from each other. Impacts of GVF differences on the Noah LSM simulations and on weather forecasts of the Weather Research and Forecasting (WRF) Model are further assessed. Results indicate that LAI-based GVF outperforms the EVI-based one when used in both the offline Noah LSM and WRF Model.


Atmosphere ◽  
2020 ◽  
Vol 11 (8) ◽  
pp. 815
Author(s):  
Marcelo Somos-Valenzuela ◽  
Francisco Manquehual-Cheuque

The use of numerical weather prediction (NWP) model to dynamically downscale coarse climate reanalysis data allows for the capture of processes that are influenced by land cover and topographic features. Climate reanalysis downscaling is useful for hydrology modeling, where catchment processes happen on a spatial scale that is not represented in reanalysis models. Selecting proper parameterization in the NWP for downscaling is crucial to downscale the climate variables of interest. In this work, we are interested in identifying at least one combination of physics in the Weather Research Forecast (WRF) model that performs well in our area of study that covers the Baker River Basin and the Northern Patagonian Icecap (NPI) in the south of Chile. We used ERA-Interim reanalysis data to run WRF in twenty-four different combinations of physics for three years in a nested domain of 22.5 and 4.5 km with 34 vertical levels. From more to less confident, we found that, for the planetary boundary layer (PBL), the best option is to use YSU; for the land surface model (LSM), the best option is the five-Layer Thermal, RRTM for longwave, Dudhia for short wave radiation, and Thompson for the microphysics. In general, the model did well for temperature (average, minimum, maximum) for most of the observation points and configurations. Precipitation was good, but just a few configurations stood out (i.e., conf-9 and conf-10). Surface pressure and Relative Humidity results were not good or bad, and it depends on the statistics with which we evaluate the time series (i.e., KGE or NSE). The results for wind speed were inferior; there was a warm bias in all of the stations. Once we identify the best configuration in our experiment, we run WRF for one year using ERA5 and FNL0832 climate reanalysis. Our results indicate that Era-interim provided better results for precipitation. In the case of temperature, FNL0832 gave better results; however, all of the models’ performances were good. Therefore, working with ERA-Interim seems the best option in this region with the physics selected. We did not experiment with changes in resolution, which may have improved results with ERA5 that has a better spatial and temporal resolution.


Sign in / Sign up

Export Citation Format

Share Document