A high resolution coupled hydrologic–hydraulic model (HiResFlood-UCI) for flash flood modeling

HiResFlood-UCI was developed by coupling the NWS's hydrologic model (HL-RDHM) with the hydraulic model (BreZo) for flash flood modeling at decameter resolutions. The coupled model uses HL-RDHM as a rainfall-runoff generator and replaces the routing scheme of HL-RDHM with the 2D hydraulic model (BreZo) in order to predict localized flood depths and velocities. A semi-automated technique of unstructured mesh generation was developed to cluster an adequate density of computational cells along river channels such that numerical errors are negligible compared with other sources of error, while ensuring that computational costs of the hydraulic model are kept to a bare minimum. HiResFlood-UCI was implemented for a watershed (ELDO2) in the DMIP2 experiment domain in Oklahoma. Using synthetic precipitation input, the model was tested for various components including HL-RDHM parameters (a priori versus calibrated), channel and floodplain Manning n values, DEM resolution (10 m versus 30 m) and computation mesh resolution (10 m+ versus 30 m+). Simulations with calibrated versus a priori parameters of HL-RDHM show that HiResFlood-UCI produces reasonable results with the a priori parameters from NWS. Sensitivities to hydraulic model resistance parameters, mesh resolution and DEM resolution are also identified, pointing to the importance of model calibration and validation for accurate prediction of localized flood intensities. HiResFlood-UCI performance was examined using 6 measured precipitation events as model input for model calibration and validation of the streamflow at the outlet. The Nash–Sutcliffe Efficiency (NSE) obtained ranges from 0.588 to 0.905. The model was also validated for the flooded map using USGS observed water level at an interior point. The predicted flood stage error is 0.82 m or less, based on a comparison to measured stage. Validation of stage and discharge predictions builds confidence in model predictions of flood extent and localized velocities, which are fundamental to reliable flash flood warning

Modification Methods For Soil And Water Assessment Tool (SWAT) Performance In Simulating Runoff And Sediment Of Watersheds In Cold Regions

Streamflow predication is an important task in water management studies. It is needed in the operation and optimization of water resources and flood control projects. The accuracy of these predictions has a great influence on the water resources management and decision making processes. Various models and tool packages have been developed for simulation and prediction of streamflow. Among them, the Soil and Water Assessment Tool (SWAT) is one of the most widely used models, which was originally developed to predict the impacts of land management on water, sediment and agricultural chemical yield in large watershed simulations. Results of the SWAT streamflow simulations have indicated that this tool has deficiencies in simulating the peaks in streamflow generated by snow melting processes in the cold regions. Since global temperature is projected to be increased and the phenomena will change the snow melting characteristics in the snow dominant areas, such as the time of first melt and rate of melting. This trend along with more precipitation will cause more flooding problems in these regions. To improve daily streamflow prediction in these regions, two methods were developed. Firstly, a method was performed by separation of winter and summer seasons simulated streamflow with subsequent validation conducted in two different seasons using Calibration Uncertainty Procedure (SWAT_CUP). It should be noted that sensitivity analysis was performed on each of the seasons separately. The second method was conducted based on coupling Artificial Neural Networks (ANNs ) with calibrated and validated results of SWAT_CUP without any separation of the seasons. The calibrated streamflow, precipitation, maximum temperature, minimum temperature, snow depth, wind speed, and relative humidity were used as inputs to the ANNs model. The results of both methods have indicated significant improvements in the simulated series. In comparison between these two methods, the operation of the second method is considered better than the first method. Although, the first method has shown improvement in the simulated results but there is still a difference between the peak streamflow and the measured streamflow by USGS (United State Geological Survey) stations. However, this difference was found diminished in the simulations using the second method. ANNs method have increased peak streamflow predication in about 70%. With this improvement, the weakness of the SWAT model in simulating sediment accumulation due to improper peak run off simulation was eliminated

NeuralHydrology -- Interpreting LSTMs in Hydrology

Despite the huge success of Long Short-Term Memory networks, their applications in environmental sciences are scarce. We argue that one reason is the difficulty to interpret the internals of trained networks. In this study, we look at the application of LSTMs for rainfall-runoff forecasting, one of the central tasks in the field of hydrology, in which the river discharge has to be predicted from meteorological observations. LSTMs are particularly well-suited for this problem since memory cells can represent dynamic reservoirs and storages, which are essential components in state-space modelling approaches of the hydrological system. On basis of two different catchments, one with snow influence and one without, we demonstrate how the trained model can be analyzed and interpreted. In the process, we show that the network internally learns to represent patterns that are consistent with our qualitative understanding of the hydrological system.Comment: Pre-print of published book chapter. See journal reference and DOI for more inf

Soil erosion in the Alps : causes and risk assessment

The issue of soil erosion in the Alps has long been neglected due to the low economic value of the agricultural land. However, soil stability is a key parameter which affects ecosystem services like slope stability, water budgets (drinking water reservoirs as well as flood prevention), vegetation productivity, ecosystem biodiversity and nutrient production. In alpine regions, spatial estimates on soil erosion are difficult to derive because the highly heterogeneous biogeophysical structure impedes measurement of soil erosion and the applicability of soil erosion models. However, remote sensing and geographic information system (GIS) methods allow for spatial estimation of soil erosion by direct detection of erosion features and supply of input data for soil erosion models. Thus, the main objective of this work is to address the problem of soil erosion risk assessment in the Alps on catchment scale with remote sensing and GIS tools. Regarding soil erosion processes the focus is on soil erosion by water (here sheet erosion) and gravity (here landslides). For these two processes we address i) the monitoring and mapping of the erosion features and related causal factors ii) soil erosion risk assessment with special emphasis on iii) the validation of existing models for alpine areas. All investigations were accomplished in the Urseren Valley (Central Swiss Alps) where the valley slopes are dramatically affected by sheet erosion and landslides. For landslides, a natural susceptibility of the catchment has been indicated by bivariate and multivariate statistical analysis. Geology, slope and stream density are the most significant static landslide causal factors. Static factors are here defined as factors that do not change their attributes during the considered time span of the study (45 years), e.g. geology, stream network. The occurrence of landslides might be significantly increased by the combined effects of global climate and land use change. Thus, our hypothesis is that more recent changes in land use and climate affected the spatial and temporal occurrence of landslides. The increase of the landslide area of 92% within 45 years in the study site confirmed our hypothesis. In order to identify the cause for the trend in landslide occurrence time-series of landslide causal factors were analysed. The analysis revealed increasing trends in the frequency and intensity of extreme rainfall events and stocking of pasture animals. These developments presumably enhanced landslide hazard. Moreover, changes in land-cover and land use were shown to have affected landslide occurrence. For instance, abandoned areas and areas with recently emerging shrub vegetation show very low landslide densities. Detailed spatial analysis of the land use with GIS and interviews with farmers confirmed the strong influence of the land use management practises on slope stability. The definite identification and quantification of the impact of these non-stationary landslide causal factors (dynamic factors) on the landslide trend was not possible due to the simultaneous change of several factors. The consideration of dynamic factors in statistical landslide susceptibility assessments is still unsolved. The latter may lead to erroneous model predictions, especially in times of dramatic environmental change. Thus, we evaluated the effect of dynamic landslide causal factors on the validity of landslide susceptibility maps for spatial and temporal predictions. For this purpose, a logistic regression model based on data of the year 2000 was set up. The resulting landslide susceptibility map was valid for spatial predictions. However, the model failed to predict the landslides that occurred in a subsequent event. In order to handle this weakness of statistic landslide modelling a multitemporal approach was developed. It is based on establishing logistic regression models for two points in time (here 1959 and 2000). Both models could correctly classify >70% of the independent spatial validation dataset. By subtracting the 1959 susceptibility map from the 2000 susceptibility map a deviation susceptibility map was obtained. Our interpretation was that these susceptibility deviations indicate the effect of dynamic causal factors on the landslide probability. The deviation map explained 85% of new independent landslides occurring after 2000. Thus, we believe it to be a suitable tool to add a time element to a susceptibility map pointing to areas with changing susceptibility due to recently changing environmental conditions or human interactions. In contrast to landslides that are a direct threat to buildings and infrastructure, sheet erosion attracts less attention because it is often an unseen process. Nonetheless, sheet erosion may account for a major proportion of soil loss. Soil loss by sheet erosion is related to high spatial variability, however, in contrast to arable fields for alpine grasslands erosion damages are long lasting and visible over longer time periods. A crucial erosion triggering parameter that can be derived from satellite imagery is fractional vegetation cover (FVC). Measurements of the radiogenic isotope Cs-137, which is a common tracer for soil erosion, confirm the importance of FVC for soil erosion yield in alpine areas. Linear spectral unmixing (LSU), mixture tuned matched filtering (MTMF) and the spectral index NDVI are applied for estimating fractional abundance of vegetation and bare soil. To account for the small scale heterogeneity of the alpine landscape very high resolved multispectral QuickBird imagery is used. The performance of LSU and MTMF for estimating percent vegetation cover is good (r²=0.85, r²=0.71 respectively). A poorer performance is achieved for bare soil (r²=0.28, r²=0.39 respectively) because compared to vegetation, bare soil has a less characteristic spectral signature in the wavelength domain detected by the QuickBird sensor. Apart from monitoring erosion controlling factors, quantification of soil erosion by applying soil erosion risk models is done. The performance of the two established models Universal Soil Loss Equation (USLE) and Pan-European Soil Erosion Risk Assessment (PESERA) for their suitability to model erosion for mountain environments is tested. Cs-137 is used to verify the resulting erosion rates from USLE and PESERA. PESERA yields no correlation to measured Cs-137 long term erosion rates and shows lower sensitivity to FVC. Thus, USLE is used to model the entire study site. The LSU-derived FVC map is used to adapt the C factor of the USLE. Compared to the low erosion rates computed with the former available low resolution dataset (1:25000) the satellite supported USLE map shows “hotspots” of soil erosion of up to 16 t ha-1 a-1. In general, Cs-137 in combination with the USLE is a very suitable method to assess soil erosion for larger areas, as both give estimates on long-term soil erosion. Especially for inaccessible alpine areas, GIS and remote sensing proved to be powerful tools that can be used for repetitive measurements of erosion features and causal factors. In times of global change it is of crucial importance to account for temporal developments. However, the evaluation of the applied soil erosion risk models revealed that the implementation of temporal aspects, such as varying climate, land use and vegetation cover is still insufficient. Thus, the proposed validation strategies (spatial, temporal and via Cs-137) are essential. Further case studies in alpine regions are needed to test the methods elaborated for the Urseren Valley. However, the presented approaches are promising with respect to improve the monitoring and identification of soil erosion risk areas in alpine regions

Assessment of climate change and development of data based prediction models of sediment yields in Upper Indus Basin

Hohe Raten von Sedimentflüssen und ihre Schätzungen in Flusseinzugsgebieten erfordern die Auswahl effizienter Quantifizierungsansätze mit einem besseren Verständnis der dominierten Faktoren, die den Erosionsprozess auf zeitlicher und räumlicher Ebene steuern. Die vorherige Bewertung von Einflussfaktoren wie Abflussvariation, Klima, Landschaft und Fließprozess ist hilfreich, um den geeigneten Modellierungsansatz zur Quantifizierung der Sedimenterträge zu entwickeln. Einer der schwächsten Aspekte bei der Quantifizierung der Sedimentfracht ist die Verwendung traditioneller Beziehung zwischen Strömungsgeschwindigkeit und Bodensatzlöschung (SRC), bei denen die hydrometeorologischen Schwankungen, Abflusserzeugungsprozesse wie Schneedecke, Schneeschmelzen, Eisschmelzen usw. nicht berücksichtigt werden können. In vielen Fällen führt die empirische Q-SSC Beziehung daher zu ungenauen Prognosen. Heute können datenbasierte Modelle mit künstlicher Intelligenz die Sedimentfracht präziser abschätzen. Die datenbasierten Modelle lernen aus den eingespeisten Datensätzen, indem sie bei komplexen Phänomenen wie dem Sedimenttransport die geeignete funktionale Beziehung zwischen dem Output und seinen Input-Variablen herstellen. In diesem Zusammenhang wurden die datenbasierten Modellierungsalgorithmen in der vorliegenden Forschungsarbeit am Lehrstuhl für Wasser- und Flussgebietsmanagement des Karlsruher Instituts für Technologie in Karlsruhe entwickelt, die zur Vorhersage von Sedimenten in oberen unteren Einzugsgebieten des oberen Indusbeckens von Pakistan (UIB) verwendet wurden. Die dieser Arbeit zugrunde liegende Methodik gliedert sich in vier Bearbeitungsschritte: (1) Vergleichende Bewertung der räumlichen Variabilität und der Trends von Abflüssen und Sedimentfrachten unter dem Einfluss des Klimawandels im oberen Indus-Becken (2) Anwendung von Soft-Computing-Modellen mit Eingabevektoren der schneedeckten Fläche zusätzlich zu hydro-klimatischen Daten zur Vorhersage der Sedimentfracht (3) Vorhersage der Sedimentfracht unter Verwendung der NDVI-Datensätze (Hydroclimate and Normalized Difference Vegetation Index) mit Soft-Computing-Modellen (4) Klimasignalisierung bei suspendierten Sedimentausträge aus Gletscher und Schnee dominierten Teileinzugsgebeiten im oberen Indus-Becken (UIB). Diese im UIB durchgeführte Analyse hat es ermöglicht, die dominiertenden Parameter wie Schneedecke und hydrologischen Prozesses besser zu und in eine verbesserte Prognose der Sedimentfrachten einfließen zu lassen. Die Analyse der Bewertung des Klimawandels von Flüssen und Sedimenten in schnee- und gletscherdominierten UIB von 13 Messstationen zeigt, dass sich die jährlichen Flüsse und suspendierten Sedimente am Hauptindus in Besham Qila stromaufwärts des Tarbela-Reservoirs im ausgeglichenen Zustand befinden. Jedoch, die jährlichen Konzentrationen suspendierter Sedimente (SSC) wurden signifikant gesenkt und lagen zwischen 18,56% und 28,20% pro Jahrzehnt in Gilgit an der Alam Bridge (von Schnee und Gletschern dominiertes Becken), Indus in Kachura und Brandu in Daggar (von weniger Niederschlag dominiertes Becken). Während der Sommerperiode war der SSC signifikant reduziert und lag zwischen 18,63% und 27,79% pro Jahrzehnt, zusammen mit den Flüssen in den Regionen Hindukush und West-Karakorum aufgrund von Anomalien des Klimawandels und im unteren Unterbecken mit Regen aufgrund der Niederschlagsreduzierung. Die SSC während der Wintersaison waren jedoch aufgrund der signifikanten Erwärmung der durchschnittlichen Lufttemperatur signifikant erhöht und lagen zwischen 20,08% und 40,72% pro Jahrzehnt. Die datenbasierte Modellierung im schnee und gletscherdominierten Gilgit Teilbecken unter Verwendung eines künstlichen neuronalen Netzwerks (ANN), eines adaptiven Neuro-Fuzzy-Logik-Inferenzsystems mit Gitterpartition (ANFIS-GP) und eines adaptiven Neuro-Fuzzy-Logik-Inferenzsystems mit subtraktivem Clustering (ANFIS) -SC), ein adaptives Neuro-Fuzzy-Logik- Inferenzsystem mit Fuzzy-C-Mittel-Clustering, multiplen adaptiven Regressionssplines (MARS) und Sedimentbewertungskurven (SRC) durchgeführt. Die Ergebnisse von Algorithmen für maschinelles Lernen zeigen, dass die Eingabekombination aus täglichen Abflüssen (Qt), Schneedeckenfläche (SCAt), Temperatur (Tt-1) und Evapotranspiration (Evapt-1) die Leistung der Sedimentvorhersagemodelle verbesserne. Nach dem Vergleich der Gesamtleistung der Modelle schnitt das ANN-Modell besser ab als die übrigen Modelle. Bei der Vorhersage der Sedimentfrachten in Spitzenzeiten lag die Vorhersage der ANN-, ANIS-FCM- und MARS-Modelle näher an den gemessenen Sedimentbelastungen. Das ANIS-FCM-Modell mit einem absoluten Gesamtfehler von 81,31% schnitt bei der Vorhersage der Spitzensedimente besser ab als ANN und MARS mit einem absoluten Gesamtfehler von 80,17% bzw. 80,16%. Die datenbasierte Modellierung der Sedimentfrachten im von Regen dominierten Brandu-Teilbecken wurde unter Verwendung von Datensätzen für Hydroklima und biophysikalische Eingaben durchgeführt, die aus Strömungen, Niederschlag, mittlerer Lufttemperatur und normalisiertem Differenzvegetationsindex (NDVI) bestehen. Die Ergebnisse von vier ANNs (Artificial Neural Networks) und drei ANFIS-Algorithmen (Adaptive Neuro-Fuzzy Logic Inference System) für das Brandu Teilnbecken haben gezeigt, dass der mittels Fernerkundung bestimmte NDVI als biophysikalische Parameter zusätzlich zu den Hydroklima-Parametern die Leistung das Modell nicht verbessert. Der ANFIS-GP schnitt in der Testphase besser ab als andere Modelle mit einer Eingangskombination aus Durchfluss und Niederschlag. ANN, eingebettet in Levenberg-Marquardt (ANN-LM) für den Zeitraum 1981-2010, schnitt jedoch am besten mit Eingabekombinationen aus Strömungen, Niederschlag und mittleren Lufttemperaturen ab. Die Ergebnisgenauigkeit R2 unter Verwendung des ANN-LM-Algorithmus verbesserte sich im Vergleich zur Sedimentbewertungskurve (SRC) um bis zu 28%. Es wurde gezeigt, dass für den unteren Teil der UIB-Flüsse Niederschlag und mittlere Lufttemperatur dominierende Faktoren für die Vorhersage von Sedimenterträgen sind und biophysikalische Parameter (NDVI) eine untergeordnete Rolle spielen. Die Modellierung zur Bewertung der Änderungen des SSC in schnee- und gletschergespeiste Gilgit- und Astore-Teilbecken wurde unter Verwendung des Temp-Index degree day modell durchgeführt. Die Ergebnisse des Mann-Kendall-Trendtests in den Flüssen Gilgit und Astore zeigten, dass der Anstieg des SSC während der Wintersaison auf die Erwärmung der mittleren Lufttemperatur, die Zunahme der Winterniederschläge und die Zunahme der Schneeschmelzen im Winter zurückzuführen ist. Während der Frühjahrssaison haben die Niederschlags- und Schneedeckenanteile im Gilgit-Unterbecken zugenommen, im Gegensatz zu seiner Verringerung im Astore-Unterbecken. Im Gilgit-Unterbecken war der SSC im Sommer aufgrund des kombinierten Effekts der Karakorum-Klimaanomalie und der vergrößerten Schneedecke signifikant reduziert. Die Reduzierung des Sommer-SSC im Gilgit Fluss ist auf die Abkühlung der Sommertemperatur und die Bedeckung der exponierten proglazialen Landschaft zurückzuführen, die auf erhöhten Schnee, verringerte Trümmerflüsse Trümmerflüsse und verringerte Schneeschmelzen von Trümmergletschern zurückzuführen sind. Im Gegensatz zum Gilgit River sind die SSC im Astore River im Sommer erhöht. Der Anstieg des SSC im Astore-Unterbecken ist auf die Verringerung des Frühlingsniederschlags und der Schneedecke, die Erwärmung der mittleren Sommerlufttemperatur und den Anstieg des effektiven Niederschlags zurückzuführen. Die Ergebnisse zeigen ferner eine Verschiebung der Dominanz von Gletscherschmelzen zu Schneeschmelzen im Gilgit-Unterbecken und von Schnee zu Niederschlägen im Astore-Unterbecken bei Sedimenteden Sedimentfrachten in UIB. Die vorliegende Forschungsarbeit zur Bewertung der klimabedingten Veränderungen des SSC und seiner Vorhersage sowohl in den oberen als auch in den unteren Teilbecken des UIB wird nützlich sein, um den Sedimenttransportprozess besser zu verstehen und aufbauen auf dem verbessertenProzessverständnis ein angepasstes Sedimentmanagement und angepasste Planungen der zukünftigen Wasserinfrastrukturen im UIB ableiten zu können

Fuels treatment and wildfire effects on runoff from Sierra Nevada mixed-conifer forests

We applied an eco-hydrologic model (Regional Hydro-Ecologic Simulation System [RHESSys]), constrained with spatially distributed field measurements, to assess the impacts of forest-fuel treatments and wildfire on hydrologic fluxes in two Sierra Nevada firesheds. Strategically placed fuels treatments were implemented during 2011–2012 in the upper American River in the central Sierra Nevada (43 km2) and in the upper Fresno River in the southern Sierra Nevada (24 km2). This study used the measured vegetation changes from mechanical treatments and modelled vegetation change from wildfire to determine impacts on the water balance. The well-constrained headwater model was transferred to larger catchments based on geologic and hydrologic similarities. Fuels treatments covered 18% of the American and 29% of the Lewis catchment. Averaged over the entire catchment, treatments in the wetter central Sierra Nevada resulted in a relatively light vegetation decrease (8%), leading to a 12% runoff increase, averaged over wet and dry years. Wildfire with and without forest treatments reduced vegetation by 38% and 50% and increased runoff by 55% and 67%, respectively. Treatments in the drier southern Sierra Nevada also reduced the spatially averaged vegetation by 8%, but the runoff response was limited to an increase of less than 3% compared with no treatment. Wildfire following treatments reduced vegetation by 40%, increasing runoff by 13%. Changes to catchment-scale water-balance simulations were more sensitive to canopy cover than to leaf area index, indicating that the pattern as well as amount of vegetation treatment is important to hydrologic response

Regional Scale Generalizations of Firn Thickness and Snow Accumulation in Southeast Alaska: an Applied Deep Learning Framework

High elevation mountain glaciers are referred to as the “water towers of the world” due to their ability to store water internally and release it later through melt and runoff. Much of the work to understand these complex environments has focused on cold based polar glaciers in high latitude regions, neglecting temperate glaciers and associated snowpack which are more responsive than polar systems to minor climatic changes. Advancements in the availability of open source data, deep learning theory, and computational efficiency has opened up new methods of “data hungry” modeling that may be well suited to the type of complex terrain-climate interactions common in these environments. Through the use of regression based deep neural networks (DNNs), this study explores the applicability of this method on a 400 MHz ground penetrating radar dataset collected on the Juneau Icefield between 2012 and 2021. The primary task of these networks is to find the statistical relationship between englacial layer thickness and a collection of topographic and climatological data with the goal of upscaling measurements from radar transects to the entire region. These models are separated into two primary tasks, the prediction of firn thickness, and the prediction of annual accumulation thickness. These tasks are approached using two different data structures, one that uses only a single year of data, and one that uses all available years of data. Models are validated using the standard “leave on glacier out (LOGO) and “leave one year out” (LOYO) holdout methods depending on the data structure. Results indicate that the available data for the Juneau Icefield is better suited for the prediction of annual accumulation thickness rather than firn thickness, with much better regional scale models produced by the multiyear data structure. Through this method, realistic simulations of icefield scale firn and snowpack distributions are generated and can be paired with existing models to output estimations of key environmental characteristics such as snow water equivalent and porosity

Applied Hydrological Modeling with the Use of Geoinformatics: Theory and Practice

Water resource management and catchment analysis are crucial aspects of the twenty-first century in hydrological and environmental sciences. Linked directly with studies and research about climate change effects in global resources (e.g., diminution of rainfall dynamic), as well as continuously growing extreme natural phenomena with catastrophic results (e.g., floods and erosion), hydrological modeling has become a key priority in modern academic research goals. On a national or lower administrative level, the need for coping with natural disasters—affecting mainly human life, property, local economy, infrastructure, etc.—and the need to design management plans and projects for sustainable exploitation of natural resources set hydrological modeling in high demand by government organizations and local authorities. Thus, hazard assessment and risk evaluation modeling have become a strategic aim and an extremely useful tool for stakeholders, decision-makers, and scientific community

Regional and global hydrology and water resources issues: The role of international and national programs

This paper presents an overview of water resources issues in the context of world population growth, climate change, and variability, and provides examples of how these issues affect local and regional water policy concerns. Also discussed is the associated research of the international scientific community in regard to physically-based modeling of the hydrological cycle, with special focus on the Global Energy and Water cycle EXperiment (GEWEX) Programme. The critical role of precipitation measurements for climate model accuracy is emphasized, with a review of several satellite methods and strategies for improving precipitation measurements. Finally, the impact of semiarid regions on global hydrologic issues is underscored with a review of research conducted by SAHRA, the National Science Foundation Science and Technology Center dedicated to Sustainability of semi-Add Hydrology and Riparian Areas

Application of artificial neural networks in hydrological modeling: A case study of runoff simulation of a Himalayan glacier basin

The simulation of runoff from a Himalayan Glacier basin using an Artificial Neural Network (ANN) is presented. The performance of the ANN model is found to be superior to the Energy Balance Model and the Multiple Regression model. The RMS Error is used as the figure of merit for judging the performance of the three models, and the RMS Error for the ANN model is the latest of the three models. The ANN is faster in learning and exhibits excellent system generalization characteristics