The Laurentian Great Lakes significantly influence the climate of the Midwest and Northeast United States due to their vast thermal inertia, moisture source potential, and complex heat and moisture flux dynamics. This study presents a newly developed coupled lake–ice–atmosphere (CLIAv1) modeling system for the Great Lakes by coupling the National Aeronautics and Space Administration (NASA) Unified Weather Research and Forecasting (NU-WRF) regional climate model (RCM) with the three-dimensional (3D) Finite Volume Community Ocean Model (FVCOM) and investigates the impact of coupled dynamics on simulations of the Great Lakes' winter climate. By integrating 3D lake hydrodynamics, CLIAv1 demonstrates superior performance in reproducing observed lake surface temperatures (LSTs), ice cover distribution, and the vertical thermal structure of the Great Lakes compared to the NU-WRF model coupled with the default 1D Lake Ice Snow and Sediment Simulator (LISSS). CLIAv1 also enhances the simulation of over-lake atmospheric conditions, including air temperature, wind speed, and sensible and latent heat fluxes, underscoring the importance of resolving complex lake dynamics for reliable regional Earth system projections. More importantly, the key contribution of this study is the identification of critical physical processes that influence lake thermal structure and ice cover – processes that are missed by 1D lake models but effectively resolved by 3D lake models. Through process-oriented numerical experiments, we identify key 3D hydrodynamic processes – ice transport, heat advection, and shear production in turbulence – that explain the superiority of 3D lake models to 1D lake models, particularly in cold season performance and lake–atmosphere interactions. Critically, all three of these processes are dynamically linked to water currents – spatially and temporally evolving flow fields that are structurally absent in 1D models. This study aims to advance our understanding of the physical mechanisms that underlie the fundamental differences between 3D and 1D lake models in simulating key hydrodynamic processes during the winter season, and it offers generalized insights that are not constrained by specific model configurations.

This paper presents the composite drought indicator (CDI) that Jordanian, Lebanese, Moroccan, and Tunisian government agencies now produce monthly to support operational drought management decision making, and it describes their iterative co-development processes. The CDI is primarily intended to monitor agricultural and ecological drought on a seasonal time scale. It uses remote sensing and modelled data inputs, and it reflects anomalies in precipitation, vegetation, soil moisture, and evapotranspiration. Following quantitative and qualitative validation assessments, engagements with policymakers, and consideration of agencies’ technical and institutional capabilities and constraints, we made changes to CDI input data, modelling procedures, and integration to tailor the system for each national context. We summarize validation results, drought modelling challenges and how we overcame them through CDI improvements, and we describe the monthly CDI production process and outputs. Finally, we synthesize procedural and technical aspects of CDI development and reflect on the constraints we faced as well as trade-offs made to optimize the CDI for operational monitoring to support policy decision-making—including aspects of salience, credibility, and legitimacy—within each national context.

This paper presents an approach to automated model integration into the NASA Land Information System (LIS). LIS defines interfaces for integrating models, which do not often have standardized coupling interfaces. Additional software development is needed when implementing models. The automated approach simplifies these barriers. LIS standardized modeling interface (LSMI) is proposed as a prerequisite Fortran subroutine, which is directly callable by LIS and calls model physics. Data communication between LIS and a model can be connected through an LSMI. The LSMI arguments can be functionally categorized. Taking their specifications, implementation code can be generated based on code templates with programming logic and language elements incorporated. An Excel/VBA-based tool (LIS-MIT) has been developed for automating LIS model implementation. Since the development of LSMI requires basic skills, LIS-MIT can significantly lower the development difficulties and efforts. The proposed approach can be generalized and applied for integrating models into other modeling environments.

Using information theory, our study quantifies the importance of selected indicators for the U.S. Drought Monitor (USDM) maps. We use the technique of mutual information (MI) to measure the importance of any indicator to the USDM, and because MI is derived solely from the data, our findings are independent of any model structure (concep-tual, physically based, or empirical). We also compare these MIs against the drought representation effectiveness ratings in the North America Drought Indices and Indicators Assessment (NADIIA) survey for Ko center dot ppen climate zones. This reveals 1) agreement between some ratings and our MI values [high for example indicators like standardized precipitation evapo-transpiration index (SPEI)]; 2) some divergences (e.g., soil moisture has high ratings but near-zero MIs for ESA Climate Change Initiative (CCI) soil moisture in the Western United States, indicating the need of another remotely sensed soil moisture source); and 3) new insights into the importance of variables such as snow water equivalent (SWE) that are not included in sources like NADIIA. Further analysis of the MI results yields findings related to 1) hydrological mechanisms (summertime SWE domination during individual drought events through snowmelt into the water-scarce soil); 2) hydrocli-matic types (the top pair of inputs in the Western and non-Western regions are SPEIs and soil moistures, respectively); and 3) predictability (high for the California 2012-17 event, with longer-time scale indicators dominating). Finally, the high MIs between multiple indicators jointly and the USDM indicate potentially high drought forecasting accuracies achievable using only model-based inputs, and the potential for global drought monitoring using only remotely sensed inputs, espe-cially for locations having insufficient in situ observations.SIGNIFICANCE STATEMENT: Drought maps from the U.S. Drought Monitor and the Objective Short-and Long-Term Drought Indicator Blends and Blend Equivalents are integrated information sources of the different types of drought. Multiple indicators go into creation of these maps, yet it is usually not clear to both public and private stake-holders like local agencies and insurance companies about the importance of any indicator in any region and season to the drought maps. Our study provides such objective information to enable understanding the mechanism and type of drought occurring at a location, season, and possibly event of interest, as well as to potentially aid in better drought monitoring and forecasting using smaller custom sets of indicators.

The Modern-Era Retrospective Analysis for Research and Application version 2 (MERRA-2) is a well-established reanalysis dataset and is widely used for driving global-scale hydrological models. However, owing to its relatively coarse spatial resolution (0.5°), the capability of MERRA-2 is repeatedly challenged in regional-scale studies, especially for smaller areas of interest. In addition, the availability of in situ observation data is a pressing issue for generating meteorological forcing. We developed a grid-based high spatial (0.125°) and temporal (hourly) resolution meteorological forcing dataset, which can evaluate hydrological processes in South Korea using state-of-the-art meteorological observations from 1980 to 2020. The forcing dataset was created by combining Automated Synoptic Observing System (ASOS) in situ measurement data from the Korean Meteorological Administration and MERRA-2 reanalysis datasets. Five meteorological variables were provided in the ASOS-MERRA2 (precipitation, air temperature, surface pressure, specific humidity, and wind speed). The study demonstrates that the region-based and high spatial resolution of ASOS-MERRA2 is superior to the existing MERRA-2 with improvements of all five weather variables, for example, from 5.6 to 2.8 mm root mean square error of precipitation. The ASOS-MERRA2 was more capable of reducing the biases and root mean squared error by improving the coefficient of determination compared with MERRA-2 for all five variables. The newly developed ASOS-MERRA2 provides an opportunity to drive land surface models to evaluate the hydroclimatic conditions in South Korea.

In situ hydrologic monitoring over regions most susceptible to food insecurity can be a challenge in current times due to various socio-economic and political issues in combination with environmental factors such as ongoing famine or drought. Hydrologic monitoring and initializing forecasts based on remotely sensed and analyzed data can contribute significantly to early warning in such regions. Routine hydrologic forecasts, as provided by NASA's Hydrologic Forecasting and Analysis System (NHyFAS), are a recent addition to early warning systems. A custom instance of NHyFAS, termed FLDAS-Forecast, is used by FEWS NET's Land Data Assimilation System (FLDAS). The FLDAS-Forecast's dynamic forecasting component was originally set up with Goddard Earth Observing System (GEOS) forecast inputs and has been recently expanded with precipitation forecast forcing from the North American Multi-Model Ensemble (NMME). This paper describes the improvements in seasonal hydrologic forecasts produced with this updated system. Evaluations in this study focus on soil moisture across southern Africa's growing season. Soil moisture forecasts are benchmarked and evaluated relative to climatology-based forecasts and historic runs, which are driven by observation-based meteorological forcing fields, and they are verified with remotely sensed observations of soil moisture and vegetation. Through multiple deterministic and probabilistic skill assessments, we show that using the larger ensemble of NMME precipitation inputs in the forecast system results in higher quality hydrologic forecasts than are allowed by climatology- or GEOS-only-based forecasts. Further, the near-real-time NMME-based rootzone soil moisture forecasts were able to correctly predict developing drought conditions over southern Africa through late 2019 and into early 2020.

This article describes a new precipitation analysis algorithm developed by NASA for time-sensitive operations at the United States Air Force. Implemented as part of the Land Information System-a land modeling and data assimilation software framework-this NASA-Air Force Precipitation Analysis (NAFPA) combines numerical weather prediction model outputs with rain gauge measurements and satellite estimates to produce global, gridded 3-h accumulated precipitation fields at approximately 10-km resolution. Input observations are subjected to quality control checks before being used by the Bratseth analysis algorithm that converges to optimal interpolation. NAFPA assimilates up to 3.5 million observations without artificial data thinning or selection. To evaluate this new approach, a multiyear reanalysis is generated and intercompared with eight alternative precipitation products across the contiguous United States, Africa, and the mon-soon region of eastern Asia. NAFPA yields superior accuracy and correlation over low-latency (up to 14 h) alternatives (numerical weather prediction and satellite retrievals), and often outperforms high-latency (up to 3.5 months) products, although the details for the latter vary by region and product. The development of NAFPA offers a high-quality, near-real-time product for use in meteorological, land surface, and hydrological research and applications. SIGNIFICANCE STATEMENT: Precipitation is a key input to land modeling systems due to effects on soil moisture and other parts of the hydrologic cycle. It is also of interest to government decision-makers due to impacts on human activities. Here we present a new precipitation analysis based on available near-real-time data. By running the program for prior years and comparing with alternative products, we demonstrate that our analysis provides better accuracy and usually less bias than near-real-time satellite data alone, and better accuracy and correlation than data provided by numerical weather models. Our analysis is also competitive with other products created months after the fact, justifying confidence in using our analysis in near-real-time operations.

Representation of irrigation in Earth System Models has advanced over the past decade, yet large uncertainties persist in the effective simulation of irrigation practices, particularly over locations where the on-ground practices and climate impacts are less reliably known. Here we investigate the utility of assimilating remotely sensed vegetation data for improving irrigation water use and associated fluxes within a land surface model. We show that assimilating optical sensor-based leaf area index estimates significantly improves the simulation of irrigation water use when compared to the USGS ground reports. For heavily irrigated areas, assimilation improves the evaporative fluxes and gross primary production (GPP) simulations, with the median correlation increasing by 0.1-1.1 and 0.3-0.6, respectively, as compared to the reference datasets. Further, bias improvements in the range of 14-35 mm mo(-1) and 10-82 g m(-2) mo(-1) are obtained in evaporative fluxes and GPP as a result of incorporating vegetation constraints, respectively. These results demonstrate that the use of remotely sensed vegetation data is an effective, observation-informed, globally applicable approach for simulating irrigation and characterizing its impacts on water and carbon states.

This study describes the development of a new, near real-time (NRT) global operational snow analysis, called the US Air Force Snow and Ice Analysis (USAF-SI). USAF-SI is a timely replacement for USAF 557th Weather Wing (557 WW)’s Snow Depth Analysis Model (SNODEP), which provides a six-hourly global analysis of snow at 25-km resolution and is used by numerous other 557 WW programs and external customers. USAF-SI incorporates several major enhancements over SNODEP. USAF-SI employs additional ancillary datasets and updates the algorithms related to passive microwave-based snow depth (SD) retrieval, surface observation blending, and quality control. The spatial resolution of USAF-SI estimates is also enhanced from 25-km to 10-km. The quality of USAF-SI is evaluated by comparing it against independent in-situ measurements at ~2300 sites. USAF-SI demonstrates improved performance over SNODEP, particularly during the peak snow accumulation and spring snow melt periods. In addition, USAF-SI is compared with routinely available SD analysis such as Modern-Era Retrospective Analysis for Research and Applications, version 2, European Centre for Medium-Range Weather Forecasts ERA5-Land, and Canadian Meteorological Center. Overall, the NRT USAF-SI provides comparable skills to these high-latency products, while outperforming them in some of the regional evaluations.

Eric Wood's contributions to distributed modeling were partially motivated by a desire to test the REA hypothesis as well as by a desire to demonstrate the impact of remotely sensed data on hydrologic prediction. In this brief talk, I will review the advances in distributed modeling, such as high-resolution terrain, distributed hydrometeorological forcings and soil-vegetation parameters, high performance computing and communications, data assimilation, coupled land-atmosphere modeling, that laid the foundation for macroscale and ultimately "hyperresolution" modeling. These foundational advances exemplify the 3rd paradigm in hydrology and are moving us towards embracing a 4th paradigm in hydrology, where we enable a rigorous confrontation of our hypotheses embodied within our models with a range of data types across many locations and spatial-temporal scales.

A globally consistent ground validation method for remotely sensed precipitation products is crucial for building confidence in these products. This study develops a new methodology to validate the IMERG precipitation products through the use of SMAP soil moisture changes as a proxy for precipitation occurrence. Using a standard 2 3 2 contingency table method, preliminary results provide confidence in SMAP's ability to be utilized as a validation tool for IMERG as results are comparable to previous validation studies. However, the method allows for an overestimate of false alarm frequency due to light precipitation events that can evaporate before the subsequent SMAP overpass and changes in overpass-to-overpass SMAP soil moisture that are within the range of SMAP uncertainty. To counter these issues, a 3 x 3 contingency table is used to reduce noise and extract more signal from the detection method. Through the use of this novel approach, the validation method produces a global mean POD of 0.64 and global mean FAR of 0.40, the first global-scale ground validation skill scores for the IMERG products. Advancing the method to validate precipitation quantity and the development of a real-time validation for the IMERG Early product are the crucial next developments.

Abstract. From the Hindu Kush Mountains to the Registan desert, Afghanistan is a diverse landscape where droughts, floods, conflict, and economic market accessibility pose challenges for agricultural livelihoods and food security. The ability to remotely monitor environmental conditions is critical to support decision making for humanitarian assistance. The FEWS NET Land Data Assimilation System (FLDAS) global and Central Asia data streams described here combine meteorological reanalysis datasets and land surface models to generate routine estimates of snow-covered fraction, snow water equivalent, soil moisture, runoff and other variables representing the water and energy balance. This approach allows us to fill the gap created by the lack of in situ hydrologic data in the region. First, we describe the configuration of the FLDAS and the two resultant data streams: one, global, at ~1 month latency, provides monthly average outputs on a 10 km2 grid from 1982–present. The second data stream, Central Asia, at ~1 day latency, provides daily average outputs on a 1 km2 grid from 2001–present. We describe our verification of these data that are compared to other remotely sensed estimates as well as qualitative field reports. These data and value-added products (e.g., anomalies and interactive time series) are hosted by NASA and USGS data portals for public use. The global data stream with a longer record, is useful for exploring interannual variability, relationships with atmospheric-oceanic teleconnections (e.g., ENSO), trends over time, and monitoring drought. Meanwhile, the higher spatial resolution Central Asia data stream, with lower latency, is useful for simulating snow-hydrologic dynamics in complex topography for monitoring snowpack and flood risk.

The USAF Weather (AFW) supports a number of military and U.S. government agencies by providing authoritative weather analysis and forecast products for any location globally, including soil moisture analyses. The long history of supporting soil moisture products and partnering with other U.S. government agencies led to the partnering between the U.S. Air Force (USAF) and NASA Goddard Space Flight Center, resulting in a merger of those organizations' modeling systems, collaborative development of the Land Information System (LIS), and operational fielding of the system within the USAF 557th Weather Wing [557 WW; formerly, Headquarters Air Force Weather Agency (HQ AFWA)]. In 2009, the USAF implemented the NASA LIS and later made it the primary software system to generate global soil hydrology and energy budget products. The implementation of LIS delivered a significant upgrade over the existing Land Data Assimilation System (LDAS) the USAF operated, the Agriculture Meteorology (AGRMET) system. Implementation enabled the rapid integration of new LDAS technology into USAF operations, and led to a long-term NASA-USAF partnership resulting in continued development, integration, and implementation of new LIS capabilities. This paper documents both the history of the USAF Weather organization capabilities enabling the generation of soil moisture and other land surface analysis products, and describes the USAF-NASA partnership leading to the development of the merged LIS-AGRMET system. The article also presents a successful example of a mutually beneficial partnership that has enabled cutting-edge land analysis capabilities at the USAF, while transitioning NASA software and satellite data into USAF operations.

The Middle East and North Africa (MENA) region has experienced more frequent and severe drought events in recent decades, leading to increasingly pressing concerns over already strained food and water security. An effective drought monitoring and early warning system is thus critical to support risk mitigation and management by countries in the region. Here we investigate the potential for assimilation of leaf area index (LAI) and soil moisture observations to improve the representation of the overall hydrological and carbon cycles and drought by an advanced land surface model. The results reveal that assimilating soil moisture does not meaningfully improve model representation of the hydrological and biospheric processes for this region, but instead it degrades the simulation of the interannual variation in evapotranspiration (ET) and carbon fluxes, mainly due to model weaknesses in representing prognostic phenology. However, assimilating LAI leads to greater improvement, especially for transpiration and carbon fluxes, by constraining the timing of simulated vegetation growth response to evolving climate conditions. LAI assimilation also helps to correct for the erroneous interaction between the prognostic phenology and irrigation during summertime, effectively reducing a large positive bias in ET and carbon fluxes. Independently assimilating LAI or soil moisture alters the categorization of drought, with the differences being greater for more severe drought categories. We highlight the vegetation representation in response to changing land use and hydroclimate as one of the key processes to be captured for building a successful drought early warning system for the MENA region.

As human exposure to hydroclimatic extremes increase and the number of in situ precipitation observations declines, precipitation estimates, such as those provided by the Integrated Multisatellite Retrievals for Global Precipitation Measurement (GPM) (IMERG) mission, provide a critical source of information. Here, we present a new gauge-enhanced dataset [the Climate Hazards Center IMERG with Stations (CHIMES)] designed to support global crop and hydrologic modeling and monitoring. CHIMES enhances the IMERG Late Run product using an updated Climate Hazards Center (CHC) high-resolution climatology (CHPclim) and low-latency rain gauge observations. CHPclim differs from other products because it incorporates long-term averages of satellite precipitation, which increases CHP(clim')s fidelity in data-sparse areas with complex terrain. This fidelity translates into performance increases in unbiased IMERG(late) data, which we refer to as CHIME. This is augmented with gauge observations to produce CHIMES. The CHC's curated rain gauge archive contains valuable contributions from many countries. There are two versions of CHIMES: preliminary and final. The final product has more copious and better-curated station data. Every pentad and month, bias-adjusted IMERG(late) fields are combined with gauge observations to create pentadal and monthly CHIMESprelim and CHIMESfinal. Comparisons with pentadal, high-quality gridded station data show that IMERG(late) performs well (r = 0.75), but has some systematic biases which can be reduced. Monthly cross-validation results indicate that unbiasing increases the variance explained from 50% to 63% and decreases the mean absolute error from 48 to 39 mm month(-1). Gauge enhancement then increases the variance explained to 75%, reducing the mean absolute error to 27 mm month(-1).

This study presents an evaluation of the impact of vegetation conditions on a land surface model (LSM) simulation of agricultural drought. The Noah-MP LSM is used to simulate water and energy fluxes and states, which are transformed into drought categories using percentiles over the continental United States from 1979 to 2017. Leaf area index (LAI) observations are assimilated into the dynamic vegetation scheme of Noah-MP. A weekly operational drought monitor (the U.S. Drought Monitor) is used for the evaluation. The results show that LAI assimilation into Noah-MP's dynamic vegetation scheme improves the model's ability to represent drought, particularly over cropland areas. LAI assimilation improves the simulation of the drought category, detection of drought conditions, and reduces the instances of drought false alarms. The assimilation of LAI in these locations not only corrects model errors in the simulation of vegetation, but also can help to represent unmodeled physical processes such as irrigation toward improved simulation of agricultural drought.

Though coarse in spatial resolution, the nearly all weather measurements from passive microwave sensors can help in improving the spatio‐temporal coverage of optical and thermal infrared sensors for monitoring vegetation changes on the land surface. This study demonstrates the use of vegetation optical depth (VOD) retrievals from the Soil Moisture Active Passive mission for capturing the vegetation alterations from the recent 2019 to 2020 Australian bushfires and drought. The impact of vegetation disturbances on terrestrial water budget is examined by assimilating the VOD retrievals into a dynamic phenology model. The results demonstrate that assimilating VOD observations lead to improved simulation of evapotranspiration, runoff, and soil moisture states. The study also demonstrates that the vegetation changes from the 2019 to 2020 Australian drought and fires led to significant modifications in the partitioning of evaporative and runoff fluxes, resulting in increased bare soil evaporation, reduced transpiration, and higher runoff.

As Earth’s largest collection of fresh water, the Laurentian Great Lakes have enormous ecological and socio-economic value. Their basin has become a regional hotspot of climatic and limnological change, potentially threatening its vital natural resources. Consequentially, there is a need to assess the current state of climate models regarding their performance across the Great Lakes region and develop the next generation of high-resolution regional climate models to address complex limnological processes and lake-atmosphere interactions. In response to this need, the current paper focuses on the generation and analysis of a 20-member ensemble of 3-km National Aeronautics and Space Administration (NASA)-Unified Weather Research and Forecasting (NU-WRF) simulations for the 2014-2015 cold season. The study aims to identify the model’s strengths and weaknesses; optimal configuration for the region; and the impacts of different physics parameterizations, coupling to a 1D lake model, time-variant lake-surface temperatures, and spectral nudging. Several key biases are identified in the cold-season simulations for the Great Lakes region, including an atmospheric cold bias that is amplified by coupling to a 1D lake model but diminished by applying the Community Atmosphere Model radiation scheme and Morrison microphysics scheme; an excess precipitation bias; anomalously early initiation of fall lake turnover and subsequent cold lake bias; excessive and overly persistent lake ice cover; and insufficient evaporation over Lakes Superior and Huron. The research team is currently addressing these key limitations by coupling NU-WRF to a 3D lake model in support of the next generation of regional climate models for the critical Great Lakes Basin.

Rainfall retrieval algorithms for passive microwave radiometers often exploit the brightness temperature depression due to ice scattering at high-frequency channels (≥85 GHz) over land. This study presents an alternate method to estimate the daily rainfall amount using the emissivity temporal variation (i.e., Δ e ) under rain-free conditions at low-frequency channels (19, 24, and 37 GHz). Emissivity is derived from 10 passive microwave radiometers, including the Global Precipitation Measurement (GPM) Microwave Imager (GMI), the Advanced Microwave Scanning Radiometer 2 (AMSR2), three Special Sensor Microwave Imager/Sounders (SSMIS), the Advanced Technology Microwave Sounder (ATMS), and four Advanced Microwave Sounding Units-A (AMSU-A). Four different satellite combination schemes are used to derive the Δ e for daily rainfall estimates. They are all 10 satellites, 5 imagers, 6 satellites with very different equator crossing times, and GMI only. Results show that Δ e from all 10 satellites has the best performance with a correlation of 0.60 and RMSE of 6.52 mm, compared with the Integrated Multisatellite Retrievals for GPM (IMERG) Final run product. The 6-satellites scheme has comparable performance with the all-10-satellites scheme. The 5-imagers scheme performs noticeably worse with a correlation of 0.49 and RMSE of 7.28 mm, while the GMI-only scheme performs the worst with a correlation of 0.25 and RMSE of 11.36 mm. The inferior performance from the 5-imagers and GMI-only schemes can be explained by the much longer revisit time, which cannot accurately capture the emissivity temporal variation.

Estimating diffuse recharge of precipitation is fundamental to assessing groundwater sustainability. Diffuse recharge is also the process through which climate and climate change directly affect groundwater. In this study, we evaluated diffuse recharge over the conterminous United States simulated by a suite of land surface models (LSMs) that were forced using a common set of meteorological input data. Simulated annual recharge exhibited spatial patterns that were similar among the LSMs, with the highest values in the eastern United States and Pacific Northwest. However, the magnitudes of annual recharge varied significantly among the models and were associated with differences in simulated ET, runoff, and snow. Evaluation against two independent datasets did not answer the question of whether the ensemble mean performs the best, due to inconsistency between those datasets. The amplitude and timing of seasonal maximum recharge differed among the models, influenced strongly by model physics governing deep soil moisture drainage rates and, in cold regions, snowmelt. Evaluation using in situ soil moisture observations suggested that true recharge peaks 1–3 months later than simulated recharge, indicating systematic biases in simulating deep soil moisture. However, recharge from lateral flows and through preferential flows cannot be inferred from soil moisture data, and the seasonal cycle of simulated groundwater storage actually compared well with in situ groundwater observations. Long-term trends in recharge were not consistently correlated with either precipitation trends or temperature trends. This study highlights the need to employ dynamic flow models in LSMs, among other improvements, to enable more accurate simulation of recharge.