<p>Preserving and promoting the sustainable use of natural resources while stabilizing healthy ecosystems under rapid environmental changes is a tremendous challenge for the international community. Science-based strategies are imperative to maintain and improve Earth&#8217;s ecosystem. Our research is designed to improve predictive ability of managed ecosystems&#8217; responses to changing weather patterns and human management. Specifically, our research seeks to develop conservation plans to improve water quality in streams and lakes, while maintaining the economic sustainability of food production systems. Reducing pollution loading into aquatic systems help improve the water quality and enhance ecosystem sustainability. Non-point pollution sources are predominant factors in increasing pollution into the water bodies. Identifying the pollution sources is important to mitigate the impact. For this reason, the main objective of our study is to identify the &#8220;hot spots&#8221; and &#8220;hot moments&#8221; of excessive nitrogen and phosphorus leaching from managed landscapes in the midwestern United States.</p><p>We developed a simple lumped model with three parameters to simulate key water fluxes - surface and subsurface runoff, and evapotranspiration (ET) in the Maumee River Basin. We designed a machine learning algorithm to identify &#8220;hot moments&#8221; using nitrogen mass balance approach at watershed-scale. The simple model helps to link the relationship between applied fertilizer and retained nutrients in the soil that the heterogeneous landscape and land management influence. Nitrogen retained in the soil will be used as an output variable and connected with predictor variable ET. Relationships between crop yield and water use in crop growth (ET) could be interpreted in a simple empirical formulation where relative change in crop yield is related to the corresponding relative change in ET, which can be expressed as,</p><p>1&#8722;&#119884;<sub>&#119886;</sub>/ &#119884;<sub>&#119909;</sub>=&#119870;<sub>&#119910;</sub> (1&#8722; &#119864;&#119879;<sub>&#119909;</sub>/&#119864;&#119879;<sub>&#119886;</sub>)</p><p>where Yx and Ya are the maximum and actual yields, ETx and ETa are the maximum and actual evapotranspiration, and Ky is a yield response factor representing the effect of relative change in ET on crop yield. The developed algorithm will be trained, tested, and validated using the coupled water flux and crop yield models. We will then demonstrate how these relationships can be extended to complex watershed model simulations that account for key land management decisions, land use pattern, crop type, soil, and topographic variability. Ultimately, we hope our findings will enhance the knowledge related to the environmental policy and decision making.</p>