Software:BAITSSS

From HandWiki
BAITSSS
Developer(s)Ramesh Dhungel and group
Written inPython (programming language), shell script, GDAL, numpy
Operating systemMicrosoft Windows
TypeEvapotranspiration modeling, irrigation simulation, surface temperature simulation, soil moisture simulation, Geographic information system

BAITSSS (Backward-Averaged Iterative Two-Source Surface temperature and energy balance Solution) is biophysical Evapotranspiration (ET) computer model that determines water use, primarily in agriculture landscape, using remote sensing-based information.[1] It was developed and refined by Ramesh Dhungel and the water resources group at University of Idaho's Kimberly Research and Extension Center since 2010. It has been used in different areas in the United States including Southern Idaho, Northern California, northwest Kansas , Texas , and Arizona.

History of development

BAITSSS originated from the research of Ramesh Dhungel, a graduate student at the University of Idaho,[2] who joined a project called "Producing and integrating time series of gridded evapotranspiration for irrigation management, hydrology and remote sensing applications" under professor Richard G. Allen.[3]

In 2012, the initial version of landscape model was developed using the Python IDLE environment using NARR weather data (~ 32 kilometers).[1] Dhungel submitted his PhD dissertation in 2014 where the model was called BATANS (backward averaged two source accelerated numerical solution).[1][2] The model was first published in Meteorological Applications journal in 2016 under the name BAITSSS as a framework to interpolate ET between the satellite overpass when thermal based surface temperature is unavailable.[1] The overall concept of backward averaging was introduced to expedite the convergence process of iteratively solved surface energy balance components which can be time-consuming and can frequently suffer non-convergence, especially in low wind speed.[1]

In 2017, the landscape BAITSSS model was scripted in Python shell, together with GDAL and NumPy libraries using NLDAS weather data (~ 12.5 kilometers).[1] The detailed independent model was evaluated against weighing lysimeter measured ET, infrared temperature (IRT) and net radiometer of drought‐tolerant corn and sorghum at Conservation and Production Research Laboratory in Bushland, Texas by group of scientists from USDA-ARS and Kansas State University between 2017 and 2020.[1] Some later development of BAITSSS includes physically based crop productivity components, i.e. biomass and crop yield computation.[1][4][5]

Rationale

Major Components of BAITSSS ET Model

The majority of remote sensing based instantaneous ET models use evaporative fraction (EF) or reference ET fraction (ETrF), similar to crop coefficients, for computing seasonal values, these models generally lack the soil water balance and Irrigation components in surface energy balance.[1] Other limiting factors is the dependence on thermal-based radiometric surface temperature, which is not always available at required temporal resolution and frequently obscured by factors such as cloud cover.[1][6] BAITSSS was developed to fill these gaps in remote sensing based models liberating the use of thermal-based radiometric surface temperature and to serve as a digital crop water tracker simulating high temporal (hourly or sub-hourly) and spatial resolution (30 meter) ET maps.[1][7][8] BAITSSS utilizes remote sensing based canopy formation information, i.e. estimation of seasonal variation of vegetation indices and senescence.[1]

Approach and model structure

Surface energy balance is one of the commonly utilized approaches to quantify ET (latent heat flux in terms of flux), where weather variables and vegetation Indices are the drivers of this process. BAITSSS adopts numerous equations to compute surface energy balance and resistances where primarily are from Javis, 1976,[9] Choudhury and Monteith, 1988,[10] and aerodynamic methods or flux-gradient relationship equations[11][12] with stability functions associated with Monin–Obukhov similarity theory.

Underlying fundamental equations of surface energy balance

Latent heat flux (LE)

The aerodynamic or flux-gradient equations of latent heat flux in BAITSSS are shown below. [math]\displaystyle{ e^{o}_c }[/math] is saturation vapor pressure at the canopy and [math]\displaystyle{ e^{o}_s }[/math] is for soil, [math]\displaystyle{ e_a }[/math] is ambient vapor pressure, rac is bulk boundary layer resistance of vegetative elements in the canopy, rah is aerodynamic resistance between zero plane displacement (d) + roughness length of momentum (zom) and measurement height (z) of wind speed, ras is the aerodynamic resistance between the substrate and canopy height (d +zom), and rss is soil surface resistance.[1]

[math]\displaystyle{ LE_{c} = \frac{\rho_a c_{p} }{{\gamma}} \bigl(\frac{e^{o}_c - e_a }{{r_{ac} + r_{ah} + r_{sc}}}\bigr)\ \And \ LE_{s} = \frac{\rho_a c_{p} }{{\gamma}} \bigl(\frac{e^{o}_s - e_a }{{r_{as} + r_{ah} + r_{ss}}}\bigr) }[/math]

The schematic of BAITSSS latent heat flux (LE) and sensible heat flux (H) as electrical analogy showing various resistances (soil surface resistance : rss and canopy resistance: rsc) and surface temperatures (canopy temperature: Tc and soil surface temperature: Ts).

Sensible heat flux (H) and surface temperature calculation[1]

The flux-gradient equations of sensible heat flux and surface temperature in BAITSSS are shown below.

[math]\displaystyle{ H_{c} = {\rho_a c_{p} } \bigl(\frac{T_c - T_a }{{r_{ah} + r_{ac}}}\bigr)\Longleftrightarrow T_{c} = \frac{H_{c}(r_{ah} + r_{ac}) }{{\rho_a c_{p}}} + {T_a } }[/math]

[math]\displaystyle{ H_{s} = {\rho_a c_{p} } \bigl(\frac{T_s - T_a }{{r_{ah} + r_{as}}}\bigr)\Longleftrightarrow T_{s} = \frac{H_{s}(r_{ah} + r_{as}) }{{\rho_a c_{p}}} + {T_a } }[/math]

Canopy resistance (rsc)

Typical Jarvis type-equation of rsc adopted in BAITSSS is shown below, Rc-min is the minimum value of rsc, LAI is leaf area index, fc is fraction of canopy cover, weighting functions representing plant response to solar radiation (F1), air temperature (F2), vapor pressure deficit (F3), and soil moisture (F4) each varying between 0 and 1.[1]

An illustration of convergence of aerodynamic resistance of surface energy balance from backward averaging (modified - green) compared to non-averaging (orange)

[math]\displaystyle{ r_{sc} = \frac{R_{c-min}}{\frac{LAI} {f_c} F_1 F_2 F_3 F_4 } }[/math]

Equations of soil water balance and irrigation decision

Standard soil water balance equations for soil surface and the root zone are implemented in BAITSSS for each time step, where irrigation decisions are based on the soil moisture at the root zone.[1]

Data

Input

File:Temperature simulation updated.tifET models, in general, need information about vegetation (physical properties and vegetation indices) and environment condition (weather data) to compute water use. Primary weather data requirements in BAITSSS are solar irradiance (Rs↓), wind speed (uz), air temperature (Ta), relative humidity (RH) or specific humidity (qa), and precipitation (P). Vegetation indices requirements in BAITSSS are leaf area index (LAI) and fractional canopy cover (fc), generally estimated from normalized difference vegetation index (NDVI). Automated BAITSSS [1] can compute ET throughout United States using National Oceanic and Atmospheric Administration (NOAA) weather data (i.e. hourly NLDAS: North American Land Data Assimilation system at 1/8 degree; ~ 12.5 kilometers), Vegetation indices those acquired by Landsat, and soil information from SSURGO.

Output

BAITSSS generates large numbers of variables (fluxes, resistances, and moisture) in gridded form in each time-step. The most commonly used outputs are evapotranspiration, evaporation, transpiration, soil moisture, irrigation amount, and surface temperature maps and time series analysis.

Model features

Time series of BAITSSS simulated daily cumulative plots of corn a) transpiration (T), b) evaporation (Ess), c) mean soil moisture at root zone (θroot), d) mean soil moisture at surface (θsur), e) evapotranspiration (ET), f) gridded precipitation (P), and simulated irrigation (Irr; bar plot) of sampled pixel at Sheridan 6 (SD-6) LEMA (100° 38′ 22″ W, 39° 21′ 38″ N) between May 10 and September 15, Kansas , United States . Shade represents 5-year maximum and minimum and the black line represents mean value.
Feature Description
Two-source energy balance BAITSSS is a two-source energy balance model (separate soil and canopy section) which is integrated by fraction of vegetation cover (fc) based on vegetation indices.
Two-layers soil water balance BAITSSS simulates soil surface moisture (θsur) and root zone moisture (θroot) layers are related to the dynamics of evaporative (Ess) and transpirative (T) flux. Capillary rise (CR) from the layer below root zone into the root zone layer is neglected. The soil moisture at both layers is restricted to field capacity (θfc).
Surface temperature BAITSSS iteratively solves surface temperature inverting flux-gradient equations of H at the soil surface (subscript s) (Ts) and canopy level (subscript c) (Tc) for each time step using continuous weather variables and surface roughness defined by vegetation Indices.
Ground heat flux of soil BAITSSS estimates ground heat flux (G) of soil surface based on sensible heat flux (Hs) or net radiation (Rn_s) of soil surface and neglects G on vegetated surface.
Transpiration Variable canopy conductance in terms of canopy resistance (rsc), based on the Jarvis-type algorithm is used to compute transpiration.
Evaporation Evaporation (Ess) in BAITSSS is computed based on soil resistance (rss) and soil water content in soil surface layer (upper 100 millimeters of soil water balance).
Irrrigation BAITSSS simulates irrigation (Irr) in agricultural landscapes[13][14] by mimicking a tipping-bucket approach (applied to surface as sprinkler or sub-surface layer as drip), using Management Allowed Depletion (MAD), and soil water content regimes at rooting depth (lower 100-2000 millimeters of soil layer).
Biomass and Yield BAITSSS computes above biomass from transpiration efficiency normalized by vapor pressure deficit and grain fraction by empirical function of biomass.
Digital agriculture-Simulated cumulative seasonal evapotranspiration (mm at 30 m spatial resolution) from hourly weather data from NLDAS and Vegetation Indices from Landsat using automated BAITSSS assuming 0.5 MAD between 10 May and 15 September 2013 for regulated groundwater management district; SD-6 LEMA, Kansas , United States (black circles, water rights shapes).

Agriculture system applications and recognition

BAITSSS was implemented to compute ET in southern Idaho for 2008, and in northern California for 2010.[1] It was used to calculate corn and sorghum ET in Bushland, Texas for 2016, and multiple crops in northwest Kansas for 2013–2017.[1][15][16][4] BAITSSS has been widely discussed among the peers around the world, including Bhattarai et al. in 2017 and Jones et al. in 2019.[17] United States Senate Committee on Agriculture, Nutrition and Forestry listed BAITSSS in its climate change report.[18] BAITSSS was also covered by articles in Open Access Government,[6][19] Landsat science team,[20] Grass & Grain magazine,[21] National Information Management & Support System (NIMSS), [22] terrestrial ecological models, [23] key research contribution related to sensible heat flux estimation and irrigation decision in remote sensing based ET models.[24][25]

In September 2019, the Northwest Kansas Groundwater Management District 4 (GMD 4) along with BAITSSS received national recognition from American Association for the Advancement of Science (AAAS).[26][27][15][28][29][30] AAAS highlighted 18 communities across the U.S. that are responding to climate change[31][32][33] including Sheridan County, Kansas to prolong the life of Ogallala Aquifer by minimizing water use where this aquifer is depleting rapidly due to extensive agricultural practices . AAAS discussed the development and use of intricate ET model BAITSSS and Dhungel's and other scientists efforts supporting effective use of water in Sheridan County, Kansas.[15]

Furthermore, Upper Republican Regional Advisory Committee of Kansas (June 2019)[16] and GMD 4[34] discussed possible benefit and utilization of BAITSSS for managing water use, educational purpose, and cost-share. A short story about Ogallala Aquifer Conservation effort from Kansas State University and GMD4 using ET model was published in Mother Earth News (April/May 2020),[35] and Progressive Crop Consultant.[36]

Example application

Groundwater and Irrigation

Averaged evapotranspiration hysteresis between normalized evapotranspiration (ET) and normalized net radiation (Rn) and normalized water vapor pressure deficit (VPD).

Dhungel et al., 2020,[1] combined with field crop scientists, systems analysts, and district water managers, applied BAITSSS at the district water management level focusing on seasonal ET and annual groundwater withdrawal rates at Sheridan 6 (SD-6) Local Enhanced Management Plan (LEMA) for five years period (2013-2017) in northwest, Kansas , United States . BAITSSS simulated irrigation was compared to reported irrigation as well as to infer deficit irrigation within water right management units (WRMU). In Kansas , groundwater pumping records are legal documents and maintained by the Kansas Division of Water Resources. The in-season water supply was compared to BAITSSS simulated ET as well-watered crop water condition.

Evapotranspiration Hysterisis and Advection

Daily evapotranspiration uncertainties. The lysimeter measurements were taken as references. Linear regressions are red line lines and one-to-one lines are black.

A study related to ET uncertainty associated with ET hysteresis (Vapor pressure and net radiation) were conducted using lysimeter, Eddy covariance (EC), and BAITSSS model (point-scale) in an advective environment of Bushland, Texas.[1] Results indicated that the pattern of hysteresis from BAITSSS closely followed the lysimeter and showed weak hysteresis related to net radiation when compared to EC. However, both lysimeter and BAITSSS showed strong hysteresis related to VPD when compared to EC.[citation needed]

Lettuce Evapotranspiration

A study related to lettuce evapotranspiration was conducted at Yuma, Arizona using BAITSSS between 2016 and 2020, where model simulated ET closely followed twelve eddy covariance sites [14]

Challenges and limitations

Simulation of hourly ET at 30 m spatial resolution for seasonal time scale is computationally challenging and data-intensive.[1][37] The low wind speed complicates the convergence of surface energy balance components as well.[1] The peer group Pan et al. in 2017 [14] and Dhungel et al., 2019 [1] pointed out the possible difficulty of parameterization and validations of these kinds of resistance based models. The simulated irrigation may vary than that actually applied in field.[1]

See also

  • METRIC, another model developed by University of Idaho that uses Landsat satellite data to compute and map evapotranspiration
  • SEBAL, uses the surface energy balance to estimate aspects of the hydrological cycle. SEBAL maps evapotranspiration, biomass growth, water deficit and soil moisture

References

  1. 1.00 1.01 1.02 1.03 1.04 1.05 1.06 1.07 1.08 1.09 1.10 1.11 1.12 1.13 1.14 1.15 1.16 1.17 1.18 1.19 1.20 1.21 1.22 1.23 1.24 Peer reviewed papers from Dhungel and others:
  2. 2.0 2.1 Dhungel, Ramesh (May 2014). "Time Integration of Evapotranspiration Using a Two Source Surface Energy Balance Model Using Narr Reanalysis Weather Data and Satellite Based Metric Data". https://digital.lib.uidaho.edu/digital/collection/etd/id/829/. 
  3. "Univ of Idaho submitted to: Producing and Integrating Time Series of Gridded Evapotranspiration for Irrigation Management, Hydrology and Remote Sensing Applications". Research, Education & Economics Information System (REEIS). United States Department of Agriculture. https://reeis.usda.gov/web/crisprojectpages/0226143-producing-and-integrating-time-series-of-gridded-evapotranspiration-for-irrigation-management-hydrology-and-remote-sensing-applications.html. 
  4. 4.0 4.1 "Effective agricultural water management strategy with biophysical evapotranspiration algorithm (BAITSSS)". https://conferences.k-state.edu/govwater/sessions/2016-faculty-staff-professional-poster-abstracts/. 
  5. "Agricultural Economics, Staff, Programs, and Publications – Kansas State University". https://www.ageconomics.k-state.edu/research/research-publications/2018%20Publications%20.pdf. 
  6. 6.0 6.1 Eccleston, Sally (2018-11-30). "Water for plant growth: The foundation of the global food supply and ecosystem services" (in en-GB). https://www.openaccessgovernment.org/water-for-plant-growth-the-foundation-of-the-global-food-supply-and-ecosystem-services/55166/. 
  7. "'It's all about water' Global Food Systems meeting" (in en-US). https://www.k-state.edu/today/students/announcement/?id=46653. 
  8. Next Generation Landscape Evapotranspiration Tools: Is It Feasible?.. 6 November 2018. https://scisoc.confex.com/scisoc/2018am/meetingapp.cgi/Paper/113852. Retrieved 2019-05-11. 
  9. Jarvis, P.G. (1976). The interpretation of the variations in leaf water potential and stomatal conductance found in canopies in the field.. OCLC 709369248. 
  10. CHOUDHURY, BJ; MONTEITH, JL (1988-01-15). "A four-layer model for the heat budget of homogeneous land surfaces". Quarterly Journal of the Royal Meteorological Society 114 (480): 373–398. doi:10.1256/smsqj.48005. ISSN 1477-870X. 
  11. Hatfield, J.L.; Baker, J.M.; Prueger, John H.; Kustas, William P. (2005), "Aerodynamic Methods for Estimating Turbulent Fluxes", Micrometeorology in Agricultural Systems (American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America), doi:10.2134/agronmonogr47.c18, ISBN 978-0-89118-268-9, http://digitalcommons.unl.edu/cgi/viewcontent.cgi?article=2399&context=usdaarsfacpub 
  12. Pagán, Brianna; Maes, Wouter; Gentine, Pierre; Martens, Brecht; Miralles, Diego (2019-02-18). "Exploring the Potential of Satellite Solar-Induced Fluorescence to Constrain Global Transpiration Estimates" (in en). Remote Sensing 11 (4): 413. doi:10.3390/rs11040413. ISSN 2072-4292. Bibcode2019RemS...11..413P. 
  13. He, Ruyan; Jin, Yufang; Kandelous, Maziar; Zaccaria, Daniele; Sanden, Blake; Snyder, Richard; Jiang, Jinbao; Hopmans, Jan (2017-05-05). "Evapotranspiration Estimate over an Almond Orchard Using Landsat Satellite Observations" (in en). Remote Sensing 9 (5): 436. doi:10.3390/rs9050436. ISSN 2072-4292. Bibcode2017RemS....9..436H. 
  14. 14.0 14.1 14.2 Dhungel, Ramesh; Anderson, Ray G.; French, Andrew N.; Saber, Mazin; Sanchez, Charles A.; Scudiero, Elia (2022-08-26). "Assessing evapotranspiration in a lettuce crop with a two-source energy balance model" (in en). Irrigation Science 41 (2): 183–196. doi:10.1007/s00271-022-00814-x. ISSN 1432-1319. https://doi.org/10.1007/s00271-022-00814-x. 
  15. 15.0 15.1 15.2 "Kansas Farmers Minimize Water Use as the Southern Great Plains Become More Arid" (in en-US). https://howwerespond.aaas.org/community-spotlight/kansas-farmers-minimize-water-use-as-the-southern-great-plains-become-more-arid/. 
  16. 16.0 16.1 "Upper Republican Regional Advisory Committee". https://kwo.ks.gov/docs/default-source/regional-advisory-committees/upper-republican-rac/upper-republican-rac-meeting-notes/upper-republican---june-2019.pdf?sfvrsn=6a028214_0. 
  17. Discussion among peers:
  18. "Peer-Reviewed Research on Climate Change by USDA Authors January 2017-August 2019". https://www.agriculture.senate.gov/imo/media/doc/USDA%20Climate%20Change%20Research%20Jan17%20Aug19_FINAL.pdf. 
  19. "News | Ogallala Aquifer Program". https://ogallala.tamu.edu/news/. 
  20. "Landsat Science Team Publications". https://www.usgs.gov/land-resources/nli/landsat/landsat-science-team-publications#t. 
  21. "Sheridan 6 LEMA a success story in water conservation and management". https://www.grassandgrain.com/archived-newspaper-editions/file/2019/10October/gg_10-01-19_sect_1.pdf.html. 
  22. "Soil, Water, and Environmental Physics to Sustain Agriculture and Natural Resources". https://www.nimss.org/projects/view/publications/18606. 
  23. Li, Sinan; Zhang, Li; Ma, Rui; Yan, Min; Tian, Xiangjun (2020-11-01). "Improved ET assimilation through incorporating SMAP soil moisture observations using a coupled process model: A study of U.S. arid and semiarid regions" (in en). Journal of Hydrology 590: 125402. doi:10.1016/j.jhydrol.2020.125402. ISSN 0022-1694. Bibcode2020JHyd..59025402L. 
  24. Mohan, M. M. Prakash; Kanchirapuzha, Rajitha; Varma, Murari R. R. (2020-10-15). "Review of approaches for the estimation of sensible heat flux in remote sensing-based evapotranspiration models". Journal of Applied Remote Sensing 14 (4): 041501. doi:10.1117/1.JRS.14.041501. ISSN 1931-3195. Bibcode2020JARS...14d1501M. 
  25. Zhang, Jingwen; Guan, Kaiyu; Peng, Bin; Jiang, Chongya; Zhou, Wang; Yang, Yi; Pan, Ming; Franz, Trenton E. et al. (2021-02-08). "Challenges and opportunities in precision irrigation decision-support systems for center pivots". Environmental Research Letters 16 (5): 053003. doi:10.1088/1748-9326/abe436. ISSN 1748-9326. Bibcode2021ERL....16e3003Z. 
  26. (in en) How We Respond – Communities and Scientists Taking Action on Climate Change, https://www.youtube.com/watch?v=QIxcTbn2tUM, retrieved 2019-11-14 
  27. "How We Respond: Stories of Community Response to Climate Change" (in en-US). https://howwerespond.aaas.org/. 
  28. "GMD 4". https://gmd4.org/. 
  29. "Kansas Rural Water Association > ONLINE RESOURCES > News Article". https://krwa.net/ONLINE-RESOURCES/News-Article/urlid/2206. 
  30. "Report spotlights Sheridan County climate effort" (in en). 2019-09-20. https://www.farmprogress.com/weather/report-spotlights-sheridan-county-climate-effort. 
  31. "'How We Respond' spotlights how US communities are addressing climate change impacts" (in en). https://www.eurekalert.org/pub_releases/2019-09/aaft-wr082819.php. 
  32. Hoy, Anne Q. (25 October 2019). "Science Magazine". Science 366 (6464): 436–437. doi:10.1126/science.366.6464.436. https://www.science.org/doi/pdf/10.1126/science.366.6464.436. 
  33. "How American communities are preparing for the impacts of climate change • Earth.com" (in en). https://www.earth.com/news/communities-preparing-climate-change/. 
  34. Kansas GMD 4 meeting minutes:
  35. "Cities Move Toward Renewable Power for Transit". 5 March 2020. https://www.motherearthnews.com/green-transportation/cities-move-toward-renweable-power-transit-zm0z2002znad/. 
  36. Scudiero, Elia (2021-07-08). "Can Artificial Intelligence Enhance the Profit and Environmental Sustainability of Agriculture?" (in en). https://progressivecrop.com/2021/07/can-artificial-intelligence-enhance-the-profit-and-environmental-sustainability-of-agriculture/. 
  37. "Point to Landscape-Scale: Traits in Agroecology" (in en-US). http://eesa.lbl.gov/event/point-to-landscape-scale-traits-in-agroecology/. 

External links