Water Surface Elevation Research Articles

Computing River Discharge Using Water Surface Elevation Based on Deep Learning Networks

Accurately computing river discharge is crucial, but traditional computing methods are complex and need the assistance of many other hydraulic parameters. Therefore, it is of practical value to develop a convenient and effective auto-computation technique for river discharge. Water surface elevation is relatively easy to obtain and there is a strong relationship between river discharge and water surface elevation, which can be used to compute river discharge. Unlike previous usage of deep learning to predict short-term river discharge that need multiple parameters besides water level, this paper proved that deep learning has the potential to accurately compute long-term river discharge purely based on water level. It showed that the majority of relative errors on the test dataset were within ±5%, particularly it could operate continuously for almost one year with high precision without retraining. Then, we used BiGRU to compute river flow with different hyperparameters, and its best RMSE, NSE, MAE, and MAPE values were 256 m3/s, 0.9973, 207 m3/s, and 0.0336, respectively. With this data-driven based technology, it will be more convenient to obtain river discharge time series directly from local water surface elevation time series accurately in natural rivers, which is of practical value to water resources management and flood protection.

River Recreational Activity Vulnerability Assessment and the Hydraulic Index Proposal

In this study, a vulnerability index and hydraulic index for recreational activities were developed for assessments in riverside areas to provide water quality and hydraulic information to the public. These novel indices consist of several river recreation activities such as swimming, water skiing, canoeing, etc., and are calculated using hydraulic information. The hydraulic information is integrated with fuzzy synthetic evaluation, with parameters such as velocity, water surface elevation, and water surface width. Also, a water quality index was created using information integrated with parameters such as DO, pH, and chlorophyll a, and then these parameters were combined into the vulnerability index. The proposed vulnerability index and hydraulic index were applied to the Nakdong River, downstream of a large weir. The hydraulic index was also combined with the results from a two-dimensional flow model for the spatial representation of the index for the categorization of safe recreational acceptability levels in the river. The results showed that the calculated index was sufficient to reflect changes in the hydraulic parameters, shown with spatial data, with comparisons to the index calculation from the gauge site measured data, with differences ranging from 0.8% to 6.5%.

Seawater intrusion assessment along the Volturno River (Italy) via numerical modeling and spectral analysis

Surface and groundwater salinization are becoming a significant challenge to inland water quality, negatively affecting people and ecosystems in coastal areas. Even if rivers provide critical pathways for seawater intrusion, this salinization phenomenon has received relatively little attention compared to other salinization mechanisms. To assess the distribution of salinity along the final reach of the Volturno River (Italy), an entire hydrologic year was modeled using the HEC-RAS software. The model was fed with high resolution time-series measurements (time interval of 10 min) of water surface elevations at both river mouth and Cancello Arnone (a hydrometric station located 13 km inland). Field observations and remote sensed data were used to perform the hydrodynamic analysis. The model showed good performance indicators (R2 = 0.878, NSE = 0.870, and MAE = 0.037 m) and well caught hydrometric variation over the simulation period. The tidal component was affected by dissipation moving upstream and showed the capability to shape the salinity profile during dry periods. Whereas during wet periods, even if a strong tidal component is present, the profile is totally regulated by the river discharge. The analysis of the salinity distribution, modelled via the Water Quality module, revealed the massive contribution of the river discharge in limiting seawater intrusion. A correlation between intrusion events and hydrometric stages was established over twenty years (2002–2022), showing a consistent trend between intrusion occurrence and the surface water storage anomaly in the lower Volturno River calculated by Global Land Data Assimilation System (GLDAS) model. Although the 1D approach here used may lead to uncertainties in the reproduction of the involved hydrodynamic and salinization processes, the results are useful for the understanding of seawater intrusion in rivers, and may be utilized to study seawater intrusion in aquifers.

Study of Tiu Kulit Dam Reservoir Operational Pattern Simulation in Sumbawa Island

Tiu Kulit Reservoir plays a pivotal role in managing water resources in West Nusa Tenggara, Indonesia, serving as a vital source for raw water, irrigation, and livestock consumption in the neighboring regions. This study aims to assess the operational dynamics of Tiu Kulit Reservoir across varying hydrological conditions, namely dry, normal, and wet years. The analysis is based on inflow rates and equilibrium water surface elevations. To achieve the objectives, an investigation of the reservoir's operation pattern was conducted. Inflow rates and equilibrium water surface elevations were scrutinized to discern operational behaviors under different hydrological scenarios. The findings indicate distinct operational strategies for Tiu Kulit Reservoir. During dry years, the outflow rate is restricted to 1.20 m3/s at a water surface elevation of +53.07 m. Conversely, the reservoir's operation in normal and wet years allows for a maximum outflow rate aimed at satisfying water demand, set at +57.00 m. In instances of excessive rainfall, the study highlights the potential for overflow beyond the spillway threshold. The study underscores the significance of these results for effective water management and resource allocation in the region. Understanding the reservoir's behavior under various conditions enables informed decision-making and sustainable utilization of water resources. To build upon this study, it is recommended to implement continuous monitoring systems for both inflow and outflow rates. Furthermore, the installation of rainfall monitoring stations within the catchment area would enhance data accuracy. Regular evaluations of the reservoir's operational patterns are also advised, allowing for adaptive strategies in response to evolving climate conditions. Moreover, future research endeavors should focus on augmenting water supply to meet the demands of a growing population and livestock, while optimizing agricultural practices to boost productivity.

Interferometric Calibration Based on a Constrained Evolutionary Algorithm without Ground Control Points for a Tiangong-2 Interferometric Imaging Radar Altimeter

The interferometric imaging radar altimeter (InIRA), mounted on the Tiangong-2 space laboratory, utilizes a small incidence and a short interferometric baseline to achieve altimetry for wide swathes of ocean surface topography and inland water surface elevation. To obtain a high-precision digital elevation model (DEM), calibration of the interferometric system parameters is necessary. Because InIRA utilizes the small-incidence interference system design, serious coupling occurs between the interferometric parameters. Commonly used interferometric calibration methods tend to fall into the local optimal solution for InIRA. Because evolutionary algorithms have a stronger robustness and global search ability, they are better suited to handling the solution space structure under the coupling of complex interferometric parameters. This article establishes an interferometric calibration optimization model for InIRA by utilizing the relative flatness of the lake surface as an inequality constraint. Furthermore, an adaptive penalty coefficient constraint evolutionary algorithm is designed to solve the model. The proposed method was tested on actual InIRA data, and the results indicate that it efficiently adjusts interferometric parameters, enhancing the precision of measurements for Qinghai Lake elevation.

Development of flood inundation maps for the Chaliyar Basin, Kerala under climate change scenarios

Abstract Floods are one of the extreme events and widespread natural disasters that significantly affect the civil infrastructure and livelihoods of people. Recently, climate change has significantly altered the rainfall pattern and increased flood events worldwide, especially in India. Therefore, it has become essential to map potential flood inundation regions for various future extreme events to develop appropriate flood mitigation and management strategies. This study aims to develop flood inundation maps for different return periods under climate change scenarios for the Chaliyar basin, Kerala. The Hydrologic Engineering Center-Hydrologic Modelling System model was used to simulate streamflow under SSP2-4.5 and SSP5-8.5 scenarios. Later, flood inundation maps were developed for different return periods using the Hydrologic Engineering Center-River Analysis System model. It was observed that for the near future (2031–2040) and far future (2071–2080), simulated streamflow is higher for SSP5-8.5. However, the mid-future (2051–2060) resulted in a higher streamflow for SSP2-4.5 than the SSP5-8.5 scenario. A maximum of 19.52 m of water surface elevation occurred at Kizhupparamba during mid-future for SSP2-4.5, followed by 18.38 m of water surface elevation at Cheekode during the near future for SSP5-8.5, for 100-year return period events. This study showed that hydrologic and hydraulic models could be effectively combined for mapping the flood inundation areas.

The role of water management and river morphology on stranding pool formation

Fish can be stranded in pools during flow recessions that may occur naturally or as a result of dam operations. Past studies have recommended dam operation guidelines that only consider water surface elevation (WSE) recession rate to reduce fish stranding potential. Here, we holistically consider the role of flow, WSE recession rate, and river-floodplain topography under unregulated and regulated flow regimes to predict the potential of stranding pool formation. The proposed methodology leverages advances in two-dimensional (2D) surface hydraulic (HD) modeling and high-resolution topobathymetric LiDAR surveys. We applied this method to the South Fork Boise River (SFBR), Idaho, USA. We developed a 2D HD model supported by meter-resolution topobathymetric data to simulate channel and floodplain hydraulics and stranding pool formation for different flow scenarios.Comparison between modeled and fieldsurveyed stranding pool locations shows good match, supporting the use of this method to predict stranding pool potential. Results show that potential for stranding pool formation (number of pools and their surface area) increases with flow for the same flow recession rate. Contrary to expectations, our model predicts that unregulated flows have a considerably higher number of potential fish stranding pools than for regulated flows in the SFBR. The reach-average WSE recession rates were higher for unregulated (0.07 to 0.09 cm/h) than for regulated (0.02 to 0.04 cm/h) flow based on daily mean flows. Our study has shown the importance of river and floodplain morphology and hydraulics to predict fish stranding pool locations and analyze impacts of flow regulations on fish stranding. This approach may provide helpful information to identify potential fish stranding pools and mitigation and to manage for fish stranding issues in regulated and unregulated river systems.

Calibrating a hydrodynamic model using water surface elevation determined from ICESat-2 derived cross-section and Sentinel-2 retrieved sub-pixel river width

Accurate observations of river water surface elevation (WSE) are critical for hydrodynamic (HD) model calibration in poorly gauged regions. This study presents a new river WSE retrieval approach that combines ICESat-2 derived cross-sections and Sentinel-2 observed sub-pixel river widths to achieve an effective HD model calibration. The unknown channel Manning coefficient, river bed elevation and cross-sectional shape in the HD model are calibrated simultaneously through a Bayesian optimization algorithm. The complete workflow is implemented in a flood-prone reach in Yiluo River basin. Results show that the WSE retrieved by the new approach has a higher spatiotemporal resolution in narrower rivers than that of Sentinel-3 altimetry derived WSE, and the root-mean-square error (RMSE) varies from 0.25 to 0.59 m against in-situ WSE. The HD model calibrated by WSE retrieved by the new approach simulates the WSE with RMSE of 0.36 m and the flood extent with critical success index (CSI) of 0.60 in validation, outperforming the results of the HD model calibrated by Sentinel-3 altimetry derived WSE, with corresponding results of 0.49 m and 0.58, respectively. Additional case studies in the Maqu and Lingkou reaches further demonstrate the applicability of the WSE retrieval approach proposed in this study for HD model calibration in different types of rivers. This study highlights the potential of the proposed WSE retrieval approach in HD model calibration to improve flood forecasting in poorly gauged regions.

A hydraulic model of the Amur River informed by ICESat-2 elevation

ABSTRACT Accurate predictions of water surface elevation (WSE) in rivers at high spatial and temporal resolution are important for flood/drought risk assessment and flood/drought forecasting and management. WSE in a river is controlled by three main factors: discharge, riverbed geometry, and hydraulic roughness. In remote and poorly instrumented rivers, discharge and riverbed geometry are highly uncertain and WSE is therefore hard to predict. ICESat-2 laser altimetry provides accurate elevation transects across the river at very high spatial resolution (70 cm along track). This paper demonstrates how ICESat-2 elevation transects can be used to parameterize a basin-scale hydraulic model of a continental-scale river. The workflow is demonstrated for the transboundary Amur River in North-East Asia. Simulated WSE is subsequently validated against a large dataset of in situ and satellite altimetry observations, and we demonstrate that the model can reproduce available WSE observations throughout the basin with an accuracy of 1–2 m.

A momentum-based prediction model of backwater rise within and downstream of an emergent canopy in nonuniform flow

Vegetation is widely used to restore waterways and wetland ecosystems. However, emergent canopies create a nonuniform water surface, and predicting the backwater rise is necessary to evaluate the flood hazard. This paper investigates the water surface elevation from downstream to upstream of an emergent canopy in a steady nonuniform flow and predicts the upstream backwater rise. Laboratory flume experiments are conducted on rigid vegetation with varying coverage exposed to different flow rates. The backwater rise and the nonuniformity of vegetated flow are explored based on a momentum equation, and a general formula for backwater rise within and downstream of an emergent canopy is proposed for nonuniform flow. The results reveal that a unified phase diagram of backwater rise can be derived from the momentum model, which presents the variation of both the upstream and downstream water surface elevations and acquires a threshold to identify strong or weak nonuniformity in vegetated flow. The upstream water depth H1 exhibits a nonmonotonic variation, first decreasing and then increasing with water depth H2 at the trailing edge of the emergent canopy. Two partitions of subcritical flow and supercritical flow for the downstream water depth H3 are quantified by the experimental data. The proposed model can be used to accurately predict the backwater rise based on a dimensionless vegetated parameter Cveg and unit discharge q, as demonstrated by both the present and previous experimental results. This work improves the understanding and predictability of the nonuniformity and backwater rise generated by emergent canopies in natural or ecology-restored rivers.

Expressive fluxes over Amazon floodplain revealed by 2D hydrodynamic modelling

Water fluxes in the Amazon River floodplain affect hydrodynamic and ecological processes from local to global scales, but they remain poorly understood due to difficult accessibility and limited data. We characterized the hydrodynamics of eight floodplain units of the central Amazon River (40000 km2) using the 2D hydraulic model HEC-RAS. High resolution modelling (∼400 m) improved the representation of river and floodplain discharge, water surface elevation (77 cm accuracy) and flood extent (∼80% - high water period, ∼52% - low water period) compared to past modelling studies. Our results show that floodplain flows during floods are very intense with upstream inflow and downstream outflow of the floodplain units. These gross flows are much larger than the net flows, with values of up to 20% of the Amazon River discharge and residence time around 6 days during floods (several months during low water period). Water extent did not show strong interannual variability during floods as the volume stored in the floodplain did, possibly due to topographic constraints. Significant hysteresis in flood extent and volume, and active and storage zones on the floodplain highlight the complexity of floodplain hydrodynamics. Extreme floods strongly impacted the onset and duration of the flood by up to one month and, consequently on duration of high water renewal period with the river. Our characterization is important to assess the effects of extreme floods on riverine communities, understand nutrient and sediment variations in the floodplain, and characterize the export of water, sediment, and carbon flux to the ocean from the world's largest hydrological system.

Steady flow at a breached levee

The present study uses a laboratory scale experimental model and an analytical model to study the steady flow at a breached levee. Flow past a zero height side weir is considered for the purpose. Water surface elevation and the three-dimensional velocity field are measured by a point gauge and an Acoustic Doppler Velocimeter (Vectrino), respectively. Experimental results include the velocity field, free surface elevation, turbulent kinetic energy and bed shear stress. Empirical relations to predict the minimum water depth and discharge at the breach are obtained. An analytical model to predict the approximate flow field at the breach is also developed. The model parameters are calibrated using the experimental results. The approximate flow conditions predicted by the analytical model can be used for field conditions. The 17th Street Canal breach is employed as a case study to demonstrate the performance of the analytical model.

Modeling of water surface profile in non-prismatic compound channels

Abstract Estimating the water surface elevation of river systems is one of the most complicated tasks in formulating hydraulic models for flood control and floodplain management. Consequently, utilizing simulation models to calibrate and validate the experimental data is crucial. HEC-RAS is used to calibrate and verify the water surface profiles for various converging compound channels in this investigation. Based on experimental data for converging channels (θ = 5°, 9°, and 12.38°), two distinct flow regimes were evaluated for validation. The predicted water surface profiles for two relative depths (β = 0.25 and 0.30) follow the same variational pattern as the experimental findings and are slightly lower than the observed values. The MAPE for the simulated and experimental results is less than 3%, indicating the predicted HEC-RAS value performance and accuracy. Therefore, our findings imply that in the case of non-prismatic rivers, the proposed HEC-RAS models are reliable for predicting water surface profiles with a high generalization capacity and do not exhibit overtraining. However, the results demonstrated that numerous variables impacting the water surface profile should be carefully considered since this would increase the disparities between HEC-RAS and experimental data.

Bias in Flood Hazard Grid Aggregation

AbstractReducing flood risk through disaster planning and risk management requires accurate estimates of exposure, damage, casualties, and environmental impacts. Models can provide such information; however, computational or data constraints often lead to the construction of such models by aggregating high‐resolution flood hazard grids to a coarser resolution, the effect of which is poorly understood. Through the application of a novel spatial classification framework, we derive closed‐form solutions for the location (e.g., flood margins) and direction of bias from flood grid aggregation independent of any study region. These solutions show bias of some key metric will always be present in regions with marginal inundation; for example, inundation area will be positively biased when water depth grids are aggregated and volume will be negatively biased when water surface elevation grids are aggregated through averaging. In a separate computational analysis, we employ the same framework to a 2018 flood and successfully reproduce the findings of our study‐region‐independent derivation. Extending the investigation to the exposure of buildings, we find regions with marginal inundation are an order of magnitude more sensitive to aggregation errors, highlighting the importance of understanding such artifacts for flood risk modelers. Of the two aggregation routines considered, averaging water surface elevation grids better preserved flood depths at buildings than averaging of water depth grids. This work provides insight into, and recommendations for, aggregating grids used by flood risk models.

Scattering of Long Waves by Freely Oscillating Submerged Plates

Abstract We consider a horizontal, submerged plate in shallow water that is allowed to oscillate in the vertical direction due to the wave loads. The plate is attached to a linear spring and damper to control the oscillations. The focus of the study is on the transformation of the wave field by the submerged oscillating plate. To estimate energy scattering, wave reflection and transmission coefficients are determined from four wave gauges: two placed upwave and two placed downwave of the oscillating plate. The flow is governed by the nonlinear level I Green–Naghdi (GN) equations, coupled with the equations of the vertical oscillations of the plate. Time series of water surface elevation recorded at gauges upwave and downwave of the plate obtained by the GN model are compared with the available laboratory experiments and other data, and very good agreement is observed. Wave reflection and transmission coefficients are then determined for a range of involved parameters, including wave conditions (wavelength and wave height), initial submergence depth of the plate, plate length, and the spring-damper system attached to the plate. It is found that a submerged oscillating plate can have a remarkable effect on the wave field and that nonlinearity plays an important role in this wave–structure interaction problem. Discussion is provided on how the wave reflection and transmission coefficients vary with the wave conditions, plate characteristics, initial submergence depth, and spring-damper system properties.

A comparison of spatial and temporal drivers of stream metabolism

Abstract Stream metabolism provides insight into the functional processes that regulate trophic status, biomass, and nutrient cycling in streams. Comparative (spatial) analyses of streams generally do not account for shifting at‐site controls on stream metabolism over time and, thus, may not identify the primary controls on cross‐site differences in trophic status, biomass, and nutrient cycling over longer time scales. The spatial component included assessing environmental factors controlling stream metabolism (daily values) during a summer low‐flow period (August) at 17 wadable stream sites in four regions of the United States. Explanatory variables for the spatial analysis included discrete nutrients and physical parameters. The temporal component included assessing 12 of the 17 sites with more than 90 days of stream metabolism values, and continuous nitrate, photosynthetically active radiation (PAR), turbidity, maximum water temperature (°C) and stream stage (water surface elevation). Spatial analysis at the 17 sites during August (31 days) found that average gross primary production (GPP) ranged from 0.1 to 4 g O2 m−2 day−1 and ecosystem respiration (ER) from 0.37 to 6.4 g O2 m−2 day−1. Sites were primarily heterotrophic; however, some sites varied from autotrophic to heterotrophic daily. Multiple regression indicated that GPP was a function of orthophosphate, percent canopy cover, and the number of days since a high flow (R2 = 0.51), whereas ER was primarily a function of GPP (R2 = 0.40). Temporal patterns in daily stream metabolism were evaluated using at‐site autoregressive models at 12 streams with relatively complete, continuous records. PAR, nitrate and maximum water temperature were the dominant variables influencing GPP, and water depth and nitrate were the most common explanatory variables for ER models. The autoregressive components (i.e., GPP or ER for the previous day) had a strong influence on both GPP and ER except after spates, indicating that biomass may be the dominant control on metabolism rather than variability in environmental factors during stable‐flow periods. We found two dominant temporal metabolism regimes: a seasonal pattern with elevated spring metabolism followed by low stable summer metabolism and a temporally disturbed pattern which tended to occur in streams with frequent spates. These patterns were not indicated by differences in metabolism among streams during stable low‐flow conditions. Because streams are generally heterotrophic with only limited periods when production may be high, comparative analyses must estimate stream metabolism over a time scale representing at‐site spatial variability and account for controls on periods of high production to adequately characterise differences between streams.

Turbulent boundary layer profiles in airflow over young wind waves in co- and counter-wind water current

Understanding the air velocity profile over random and three-dimensional wind-waves is crucial for evaluating momentum and energy transfer between air and water. Unfortunately, determining accurate air velocity profiles in field conditions is nearly impossible. It is well accepted, though, that the airflow above water waves has a logarithmic profile usually expressed in terms of the effective roughness parameter. However, this parameter cannot be evaluated directly from the wave measurements at the water surface. For young wind waves, it was recently shown that the airflow over the waves maintains wall similarity. In this case, the airflow vertical velocity distribution can be described by the Nikuradse fully rough logarithmic profile, using the root mean square value of the water surface elevation as the equivalent sand grain roughness height. The existence of mean current in water in either co- or counter-wind direction may significantly modify the wave field, the vertical wind-velocity profile and thus the momentum and energy transfer from air to water. The effect of the water current on the spatially developing boundary layer over young wind waves and the wall similarity is examined. Combined laboratory measurements at several fetches of finely resolved mean air velocity profile above the water surface and of the characteristics of the wind-wave field are performed at multiple wind-forcing and mean water current conditions. The shear stress at the air–water interface estimated using two independent approaches is weakly dependent on current, while the resultant wave field differs significantly.

SaTSeaD: Satellite Triangulated Sea Depth Open-Source Bathymetry Module for NASA Ames Stereo Pipeline

We developed the first-ever bathymetric module for the NASA Ames Stereo Pipeline (ASP) open-source topographic software called Satellite Triangulated Sea Depth, or SaTSeaD, to derive nearshore bathymetry from stereo imagery. Correct bathymetry measurements depend on water surface elevation, and whereas previous methods considered the water surface horizontal, our bathymetric module accounts for the curvature of the Earth in the imagery. The process is semiautomatic, reliable, and repeatable, independent of any external bathymetry data eliminating user bias in selecting bathymetry calibration points, and it can generate a fully integrated and seamless topo-bathymetry digital elevation model (TBDEM) in the same coordinate system, comparable with the band-ratio method irrespective of the regression method used for the band-ratio algorithm. The ASP output can be improved by applying a camera bundle adjustment to minimize reprojection errors and by alignment to a more accurate topographic (above water) surface without any bathymetric input since the derived TBDEM is a rigid surface. These procedures can decrease bathymetry root mean square errors from 30 to 80 percent, depending on environmental conditions, the quality of satellite imagery, and the spectral band used (e.g., blue, green, or panchromatic).

Analysing the Impact of Land Subsidence on the Flooding Risk: Evaluation Through InSAR and Modelling

Floods greatly impact human settlements in flood risk areas, such as floodplains and coastal lowlands, following heavy rainfall. The Alto Guadalentin valley, an orogenic tectonic depression, experiences extreme flash floods and land subsidence due to groundwater withdrawal, rendering it one of Europe's fastest subsiding regions. In this study, we compared two 2D flood event models representing different land subsidence scenarios for 1992 and 2016. To determine the flooded area and water depth variations due to land subsidence, the Hydrologic Engineering Centre River Analysis System 2D (HEC-RAS 2D) model was used to simulate flood inundation by the Alto Guadalentin River and its tributaries. Synthetic aperture radar (SAR) satellite (ERS, ENVISAT, and Cosmo-SkyMED) images were employed, along with the interferometric SAR (InSAR) technique, to calculate the magnitude and spatial distribution of land subsidence. By analysing the accumulated subsidence distributions obtained from InSAR, the original topography of the valley in 1992 and 2016 was reconstructed. These digital surface models (DSMs) were then used to generate 2D hydraulic models, simulating flood scenarios in the unsteady mode. The results demonstrated significant changes in the water surface elevation over the 14-year period, with a 2.04 km2 increase in areas with water depths exceeding 0.7 m. These findings were utilized to create a flood risk map and assess the economic flood risk. The data highlight the crucial role of land subsidence in determining the inundation risk in the Alto Guadalentin valley, providing valuable insights for emergency management and civil protection against future potential flooding events.

Recalibration of a three-dimensional water quality model with a newly developed autocalibration toolkit (EFDC-ACT v1.0.0): how much improvement will be achieved with a wider hydrological variability?

Abstract. Autocalibration techniques have the potential to enhance the efficiency and accuracy of intricate process-based hydrodynamic and water quality models. In this study, we developed a new R-based autocalibration toolkit for the Environmental Fluid Dynamics Code (EFDC) and implemented it into the recalibration of the Yuqiao Reservoir Water Quality Model (YRWQM), with long-term observations from 2006 to 2015, including dry, normal, and wet years. The autocalibration toolkit facilitated recalibration and contributed to exploring how a model recalibrated with long-term observations performs more accurately and robustly. Previously, the original YRWQM was calibrated and validated with observations of dry years in 2006 and 2007, respectively. Compared to the original YRWQM, the recalibrated YRWQM performed just as well in water surface elevation, with a Kling–Gupta efficiency (KGE) of 0.99, and water temperature, with a KGE of 0.91, while performing better in modeling total phosphorus (TP), chlorophyll a (Chl a), and dissolved oxygen (DO), with KGEs of 0.10, 0.30, and 0.74, respectively. Furthermore, the KGEs improved by 43 %–202 % in modeling the TP–Chl a–DO process when compared to the models calibrated with only dry, normal, and wet years. The model calibrated in dry years overestimated DO concentrations, probably explained by the parameter of algal growth rate that increased by 84 %. The model calibrated in wet years performed poorly for Chl a, due to a 50 % reduction in the carbon-to-chlorophyll ratio, probably triggered by changes in the composition of the algal population. Our study suggests that calibrating process-based hydrodynamic and water quality models with long-term observations may be an important measure to improve the robustness of models under severe hydrological variability. The newly developed general automatic calibration toolkit and a possible hierarchical autocalibration strategy will also be a powerful tool for future complex model calibration.