Water Surface Elevation Research Articles

Hydraulics of Time-Variable Water Surface Slope in Rivers Observed by Satellite Altimetry

The ICESat-2 and SWOT satellite earth observation missions have provided highly accurate water surface slope (WSS) observations in global rivers for the first time. While water surface slope is expected to remain constant in time for approximately uniform flow conditions, we observe time varying water surface slope in many river reaches around the globe in the ICESat-2 record. Here, we investigate the causes of time variability of WSSs using simplified river hydraulic models based on the theory of steady, gradually varied flow. We identify bed slope or cross section shape changes, river confluences, flood waves, and backwater effects from lakes, reservoirs, or the ocean as the main non-uniform hydraulic situations in natural rivers that cause time changes of WSSs. We illustrate these phenomena at selected river sites around the world, using ICESat-2 data and river discharge estimates. The analysis shows that WSS observations from space can provide new insights into river hydraulics and can enable the estimation of river discharge from combined observations of water surface elevation and WSSs at sites with complex hydraulic characteristics.

Methodological evaluation of river discharges derived from remote sensing and land surface models

This study assesses river discharges derived using remote sensing and hydrologic modeling approaches throughout the CONUS. The remote sensing methods rely on total water storage anomalies (TWSA) from the GRACE satellite mission and water surface elevations from altimetry satellites (JASON-2/3, Sentinel-3). Surface and subsurface runoff from two Land Surface Models (NOAH, CLSM) are routed using the Hillslope River Routing model to determine discharge. The LSMs are part of NASA’s Global Land Data Assimilation System (GLDAS). Differences in key physical processes represented in each model, model forcings, and use of data assimilation provide an intriguing basis for comparison. Evaluation is performed using the Kling Gupta Efficiency and USGS stream gauges. Results highlight the effectiveness of both satellite-derived discharge methods, with altimetry generally performing well over a range of discharges and TWSA capturing mean flows. LSM-derived discharge performance varies based on hydroclimatic conditions and drainage areas, with NOAH generally outperforming CLSM. CLSM-derived discharges may be impacted by the use of data assimilation (GLDAS v2.2). Low correlation and high variability contribute to lower KGE values. GLDAS models tend to perform poorly in snow dominated, semi-arid and water-regulated systems where both the timing and magnitude of the simulated results are early and overestimated.

Solar powered integrated multi sensors to monitor inland lake water quality using statistical data fusion technique with Kalman filter

This study proposes a data-driven statistical model using multi sensor fusion and Kalman filtering for real-time water quality assessment in lakes. A recursive estimation technique, the Kalman Filter, is employed to handle uncertainties and enhance computational efficiency. The fusion process integrates data from sensors monitoring parameters like chlorophyll concentration, surface water elevation, temperature, and precipitation, producing Markov features to capture temporal transitions and environmental dynamics. Data synchronization and fusion are achieved through recursive KF methods, enabling real-time adaptive management in response to environmental fluctuations such as seasonal changes, precipitation (6–18%), and evaporation rates (1.2–11.9 mm/day). Over a 30-day evaluation period, the model accurately predicted chlorophyll concentrations, reaching 128 mg/m3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {mg}/\\hbox {m}^{3}$$\\end{document} in mid-level inflow regions (3.6 m water elevation) compared to 86 mg/m3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {mg}/\\hbox {m}^{3}$$\\end{document} in extreme inflow areas (5.5 m). The integration of Markov feature extraction and eigenvalue estimation enhanced prediction stability and sensitivity, with the KF maintaining computational efficiency at 7.8 ms per computation cycle. The model’s accuracy was validated by achieving a residual error of less than 0.05 with minimal noise interference. Overall, the system provides a resilient and precise framework for real-time lake water quality assessment, capable of handling multi-parameter uncertainties and dynamic environmental changes, thereby supporting informed decision-making for aquatic ecosystem management.

ExtendinG SUb-DAily River Discharge data over INdia (GUARDIAN).

River discharge information is crucial for various applications, but the measurement process often remains impeded by factors that hinder near real-time (NRT) data availability in India. However, leveraging telemetry-based water surface elevation (WSE) data across the country provides an opportunity to convert it into river discharge. This conversion is made possible by utilizing rating curves (RCs) derived from historical collocated measurements of WSE and discharge. In this study, NRT WSE from the Central Water Commission (CWC) flood portal is obtained via the web-scraping tool. Through the application of RCs, discharge data is extended across 210 gauging stations in India from the year 2020 to the present, encompassing sub-daily discharge during the non-monsoon and hourly discharge series during monsoon (June-September) across the Indian rivers. Annually, the study generated over 800,000 discharge data points for Indian rivers, accounting for more than 4000 discharge measurements per station annually. These comprehensive datasets provide valuable insights for water resource and flood management research, offering NRT access to WSE, and discharge, along with the local RCs.

Acoustic reflector remotely tunable by the acoustic radiation force

Specifically structured surface can be used to control acoustic waves reflection enabling wavefront shaping tailored for desired applications such as imaging or cloaking. With the development of metasurfaces, which enable changing the reflection of an incident acoustic wave at an interface using a sub-wavelength structure, it is possible to explore new frontiers of Snell-Descartes laws. However, the control over wave propagation is pre-established when using fixed structures, directing research towards metasurfaces reconfigurable in both time and space. This work proposes a metamaterial architecture using deformations of a fluid interface between two media through the acoustic radiation pressure (ARP). An analytical and numerical study is conducted to assess the phase shift on airborne acoustic wave induced by a unit cell featuring a Helmholtz resonator tuned by the water surface elevation induced by ARP. The possibility to program in real-time the phase shift is successfully tested experimentally for an incident wave of frequency 3400 Hz. The extension to a metasurface is introduced through the juxtaposition of 15 unit cell. A retroreflection effect is then numerically demonstrated, highlighting the effectiveness of the proposed approach. This concept of real-time and non-contact reconfigurable metasurface opens new possibilities for beam deflection, acoustic holography or information encoding.

Experimental and numerical study on focused wave generation using different energy modes

This paper reports an experimental and numerical study on the generation of freak waves related to the spatiotemporal focusing mechanism, where multiple water waves converge at a specific point in space and time to produce a large transient wave. A new energy distribution mode of constant-difference wave number (CDK) within the wave group is proposed and applied to both a physical and a numerical flume to generate focused waves, with primary attention given to the water surface elevation, underlying kinematics, and internal spectrum structure during the focused wave generation process. The classical experiment by Baldock et al. [“A laboratory study of nonlinear surface waves on water,” Philos. Trans. R. Soc., Ser. A 354, 649–676 (1996)], which used a constant-difference wave period within the wave group, is chosen for comparison. Spectral analysis shows that the CDK energy distribution mode can easily stimulate multiwave resonance in the focused wave group, where significant energy transfer from the main frequencies to higher harmonics is observed. Compared with the data from Baldock et al. [“A laboratory study of nonlinear surface waves on water,” Philos. Trans. R. Soc., Ser. A 354, 649–676 (1996)], the CDK wave group produces larger focused wave crests with higher propagation speeds and larger velocity gradients near the water surface due to strong multiwave resonance. The effects of the input spectrum, incident amplitude and initial water depth under different energy modes on focused-wave generation are analyzed. All the results indicate that the CDK energy distribution mode can produce larger focused waves compared to other energy distribution modes, which offers a novel approach for generating freak waves in the laboratory.

Data-Driven Calibration of SWOT’s Systematic Errors: First In-Flight Assessment

The SWOT satellite, carrying the KaRIN first wide-swath onboard altimeter, was launched in December 2022, and has now delivered more than a year of surface water elevation data over the ocean and inland lakes/rivers. These data are affected by systematic errors which constitute the dominant part of the error budget at scales larger than a few thousands of kilometers. Some strategies for their estimation and calibration were explored during the pre-launch studies with performance estimations. Now, based on the real data, we propose in this study to assess the systematic error budget with statistical methods relying on spectral and co-spectral analysis. From this assessment, suggesting very low error levels (below requirements), we propose the implementation of the calibration algorithms at Level-2 and Level-3 with a few minor adjustments justified by the error spectra. The calibrated products are then validated with usual CalVal metrics.

Downstream impacts of dam breach using HEC-RAS: a case of Budhigandaki concrete arch dam in central Nepal

Studies on concrete dam breach are limited compared to earthen and other types of dams. With an increase in the construction of concrete dams, particularly in the developing world, it is imperative to have a better understanding of the dam breach phenomena and the identification of the most influential breach parameters. This study aims to contribute to this gap by taking the case of the concrete arch dam proposed for the 1200 MW Budhigandaki Hydropower Project located in central Nepal. This study carries special significance for Nepal, primarily because of the increasing number of under construction and proposed large dams for water resources development in the country. We carry out dam breach analysis of the Budhigandaki dam using HEC-RAS 2D model to calculate the flood discharge peaks, time to peak, water surface elevation and the extent of inundation for two scenarios (with and without probable maximum flood) to estimate the damage on four downstream settlements. We carry out sensitivity analysis of the breach parameters on the flood magnitudes and severity. Results show that all the study locations lie in the high flood hazard zone. Flood peaks can reach as high as 286,000 m3s− 1 to 511,000 m3s− 1 in the considered settlements. The time to peak ranges from 11.3 to 17 h after the breach at these locations. We estimate that if a breach should happen, it would most likely inundate around 150,000 buildings, impact nearly 672,000 lives and flood 3,500 km of road downstream. Furthermore, dam breach elevation is found to be the most sensitive parameter to downstream floods. Hence, rather than structural measures, it is recommended that non-structural measures are implemented for minimizing the impacts of flood disasters at the study locations. The findings could be a useful reference for future dam projects in Nepal and other areas with similar hydrological and topographical conditions.

Physics-informed geostatistical approach (PIGA) for real-time fluvial flood modeling

This study presents the Physics-Informed Geostatistical Approach (PIGA), a novel methodology aimed at improving flood prediction accuracy in riverine systems during hurricane events. PIGA integrates reduced spatial correlations derived from fluvial hydrodynamic models with a reduced geostatistical approach (RGA) based on monitoring data to provide real-time estimations of water surface elevations (WSE) across river channels and reaches. Focusing on the Savannah River basin, an area prone to riverine flooding during hurricane seasons, we developed a 1D-2D HEC-RAS model as the physical model and compared its performance with PIGA for various hurricane events, including Hurricane Matthew, Irma, Dorian, and Idalia. Our findings reveal that PIGA consistently outperforms the HEC-RAS model in predicting WSE, particularly in areas with limited bathymetry data. Statistical analyses demonstrate significant improvements in prediction accuracy achieved by PIGA, emphasizing its potential as a reliable tool for flood prediction in data-deficient regions. Additionally, we explore the transferability of reduced spatial correlations across different hurricane events. Our results indicate that spatial correlations derived from one hurricane event can effectively predict WSE for other events, with Hurricane Irma’s spatial correlations proving to be the most versatile for our study site. However, we also observe variability in spatial correlations among different hurricane events, highlighting the importance of selecting the optimal alternative spatial reference for prediction. Furthermore, we highlight PIGA’s computational efficiency, with simulations completed in significantly less time compared to traditional physical models, rendering its potential to provide real-time prediction of WSE during flooding events.

Accuracy Assessment of Estimated River Water Surface Elevations from Landsat 8 and 9 Imagery among Twenty Water Indices

Accuracy Assessment of Estimated River Water Surface Elevations from Landsat 8 and 9 Imagery among Twenty Water Indices

Experimental investigation of flow and sediment transport on an equilibrium beach formed by surging breaker

Flow patterns and sediment transport induced by breaking waves on a beach profile have been studied in a laboratory wave flume. Two parallel experiments were conducted on the same profiled rigid bed and sediment (mobile) bed, each subjected to identical wave forcing. The rigid bed experiments include measurements of bed shear stress and water-surface elevation, while the sediment bed experiments involve the measurement of pore-water pressures. The beach profile features a beach step leading to the generation of a hydraulic jump at the end of the rundown phase and accordingly forms a backwash vortex. The bed shear stress measurements revealed that the maximum onshore directed mean bed shear stress occurs after the wave breaking whereas the offshore directed one occurs during rundown at the same section. The increase of the bed shear stress with respect to that in the approaching wave boundary layer during the wave breaking and hydraulic jump processes can be by as much as a factor of 13 and 4, respectively. The effect of the wave breaking on turbulence is also more pronounced by a factor of 2.2 than the hydraulic jump in terms of r.m.s values of the fluctuating component of the bed shear stress. The pore pressure measurements showed that the vortex generated during the wave breaking creates an impulsive, upward-directed pressure gradient force that can be up to 1.2 times the submerged weight of the sediment. The sediment transport occurs in the sheet flow regime in a large part of the beach profile. The turbulence generated by the backwash vortex and the hydraulic jump at the beach step directly impacts the local sediment suspension which serves as the source for the advected sediment in the swash zone. The results are compared with a previous study that used the same methodology but involved plunging waves and the initial bed profile.

Multi‐Satellite Data Assimilation for Large‐Scale Hydrological‐Hydrodynamic Prediction: Proof of Concept in the Amazon Basin

AbstractSatellite remote sensing enhances model predictions by providing insights into terrestrial and hydrological processes. While data assimilation techniques have proven promising, there is a lack of standardized and effective approaches for integrating multiple observations simultaneously. This study presents a novel assimilation framework, the multi‐observation local ensemble‐Kalman‐filter (MoLEnKF), designed to effectively integrate multiple variables, even at scales different than the model. Evaluation of MoLEnKF in the Amazon River basin includes assimilation experiments with remote sensing data only, including water surface elevation (WSE), terrestrial water storage (TWS), flood extent (FE), and soil moisture (SM). MoLEnKF demonstrates improvements in a scenario where regions lack in‐situ hydroclimatic records and when assuming uncertainties of large‐scale hydrologic‐hydrodynamic models. Assimilating WSE outperforms daily discharge and water‐level estimations, achieving 38% and 36% error reduction, respectively. However, the monthly evapotranspiration estimate achieves the greatest error reduction by assimilating SM with 11%. MoLEnKF always remains in second position in a ranking of error and uncertainty reduction, providing an intermediate condition, being able to holistically outperform univariate experiments. MoLEnKF also outperform state‐of‐the‐art models in many cases. This study suggests potential improvements, urging exploration of correlations between assimilated variables and adaptive localization methods based on seasonality. The flexibility and the elegant way of expressing the LEnKF equations by MoLEnKF facilitates their application with different types of variables, compatible with large‐scale hydrologic‐hydrodynamic models and missions such as SWOT. Its robustness ensures easy replicability worldwide, facilitating hydrological reanalysis and improved forecasting, establishing MoLEnKF as a valuable tool for the scientific community in hydrological research.

Simulation of flow characteristics in open channels with inflatable dam

Abstract Inflatable dams are structure cylindrical inflatable and deflated bile structures made of rubberized material placed on top of a rigid base and inflatable by air, water, or a combination of both. The present paper aims to measure the water surface elevation, and depth of flow in Shatt Al Ibrahim and comparison with design depth and study the flow state according to design discharge and also according to scarce discharge. The study effect of the height of the dam filled with water on upstream and downstream water depth and analyzed numerically by using HEC-RAS program one–one-dimensional flow in two states, steady flow and unsteady flow under different variables such as changing of discharge. Good results were obtained in two flow states. The results show that the proper solution for maintaining the required water levels is to construct inflatable dams, due to lower construction and maintenance costs, easy and simple of operation, simple and easy to inflate and deflate quickly to prevent upstream flooding. The results of the current study indicated an increase in the water surface height (WS) and water depth (Wd) across the channel after the implementation of inflatable dams compared to natural conditions. Specifically, the data indicate a significant increase in (WS) and (Wd) at “station 52”, where (WS) increased by 0.45 m with a 0.5 m dam and 1.43 m with a 1.5 m dam compared to the basic scenario. This confirms the adaptability of inflatable dams in adapting water levels to real-time requirements, underscoring their potential as a cost-effective solution to mitigate upstream floods and improve water resource management practices.

Understanding Salinity Intrusion and Residence Times in a Small-Scale Bar-Built Estuary under Drought Scenarios: The Maipo River Estuary, Central Chile

The Maipo River estuary is a low-inflow bar-built estuary that includes a protected wetland, which harbors a rich ecosystem. The estuary and wetland have been threatened by a persistent drought for more than a decade, which has resulted in greater salinity intrusion and increased residence times. Previous studies have described salinity and pollutants in estuaries; however, almost all have focused on deeper and/or wider estuaries with dimensions much larger than those of the small-scale Maipo River estuary. In this study, we used the numerical model FVCOM to simulate the dynamics of the Maipo River estuary under drought scenarios and explored the interactions between river discharge and tides in terms of saline intrusion and particle dispersal. The model was validated against observations collected during a field campaign near the river mouth. The simulations successfully reproduced the water surface elevation but underestimated salinity values, such that the vertical salinity structure observed in the field was not captured by the model in this shallow and morphologically complex estuary. Consequently, our model results provide qualitative insight related to salinity and baroclinic dynamics. Results of maximum saline intrusion showed an exponential decay with increasing river discharge, and the analysis of salinity intrusion time series revealed that droughts may cause permanent non-zero salinity levels in the estuary, potentially affecting ecological cycles. The incorporation of passive tracers showed that decreasing river discharge increases the residence time of particles by allowing the tracers to re-enter the estuary. Model results showed the formation of accumulation zones (hotspots) in the shallower zones of the estuary.

Extending the historical in situ gauge network over Indian river reaches utilizing the Surface Water and Ocean Topography (SWOT) mission

ABSTRACT We present a preliminary study on using the Surface Water and Ocean Topography (SWOT) mission to extend water surface elevation (WSE) measurements from a sparsely distributed in situ gauge network along a one-dimensional (1D) study reach within the Mahanadi River basin in India. We use the Centre National d’Etudes Spatiales (CNES) simulator to simulate 3 years of synthetic measurements of SWOT WSE, which are assimilated into a 1D stochastic model based on six independent configurations. For each configuration, the model is derived using historical measurements of in situ WSE to simulate daily WSE estimates every 200 m along the study reach. Assimilated outputs showed high Nash-Sutcliffe efficiency (0.697–0.811), high correlation (0.838–0.901) and low root mean square error (0.242–0.513 m) with respect to in situ validation data. Results indicate consistent performance of the SWOT-assimilated model framework and are relevant considering the declining availability of in situ monitoring stations.

Evaluation of repeat bathymetric surveys on water surface elevations on the Saint Clair River

The Great Lakes contain 20 percent of the world’s fresh surface water, a drinking source for 30- million people. The system is vital to manufacturing, largely due to efficient transportation of commodities and an abundance of natural resources. Understanding the water balance and predicting future water levels is important for many industries including commercial navigation, the tourist industry, hydropower and shoreline property owners, just to name a few. These predictions require an understanding of both the hydroclimate drivers across the Great Lakes as well as physical changes in the connecting channels that pass water between each lake. Since the International Upper Great Lakes Study began, the St. Clair River has been the topic of conveyance change investigations. This study tests the sensitivity of water surface elevations to small differences in bathymetric changes. Greater than 80 percent of the surveyed points had less than 0.24 m of change between any combination of survey years. Repeat bathymetric surveys collected in 2007, 2012 and 2021, along with hydrodynamic modeling, find changes in water surface elevations associated with bathymetric change are on the sub-centimeter scale. These changes are put into perspective and compared to known anthropogenic and naturally occurring conveyance changes. While these surveys represent a short duration in the observed water level record of the Great Lakes, they represent an important documentation of the geomorphic state of the St. Clair River. These changes can and should be compared to similar future surveys of the river over longer engineering and geomorphic timescales.

Nearshore Flow Dynamics Over Shore‐Oblique Bathymetric Features During Storm Wave Conditions

AbstractShore‐oblique bathymetric features occur around the world and have been statistically correlated with enhanced shoreline retreat on sandy beaches. However, the physical mechanisms that explain a causal relationship are not well understood. In this study, radar remote sensing observations and results from a phase‐resolved numerical model explore how complex morphology alters nearshore hydrodynamics. Observations at selected times during high‐energy storm events as well as a suite of idealized simulations indicate that shore‐oblique features induce strong spatial variations in the water surface elevation, wave breaking patterns, and mean current pathways. Re‐emergent offshore flows and longshore current accelerations occur near the shoreward apex of the oblique nearshore features. The results suggest that complex bathymetric morphology exerts a powerful control on nearshore hydrodynamics and increases the potential for enhanced cross‐shore and alongshore sediment transport.

Hydrodynamic and Sediment Transport Modeling at Lower Citarum and Cibeet River Confluence West Java

Flooding and changes in riverbed morphology often become issues near river confluences, as observed at the confluence of Citarum Hilir and Cibeet Rivers. This study aims to analyze the cross-sectional capacity and riverbed changes around the confluence of Citarum Hilir and Cibeet Rivers. This study utilizes 1D numerical modeling MIKE 11 with Hydrodynamic and Sediment Transport Modules. The hydrodynamic modeling inputs are river network and cross-sectional profiles. Hydrographs with return periods of 2, 5, 10, and 25 years are used as upstream boundaries, and a discharge curve is used as the downstream boundary of the model. The model is calibrated using water surface elevation measurements taken from September 13-16, 2021. The hydrodynamic modeling results indicate that Citarum Hilir and Cibeet Rivers cannot convey flood discharges with a 2-year return period. The sediment transport modeling results show no significant changes in the riverbed of Citarum Hilir before the confluence, but degradation of approximately 0.38 – 0.46 meters occurs after the confluence. Degradation happens in zones of maximum velocity where flow turbulence and shear stress increase, causing erosion of the Citarum Hilir riverbed. Degradation also occurs near the Cibeet River confluence due to flow acceleration at the outer bend of the river. Flooding and sediment transport are natural factors that can alter the morphology of river confluences, such as the confluence of Citarum Hilir and Cibeet Rivers.Keywords: Citarum river, Cibeet river, river confluence, hydrodinamic modeling, sedimen transport modeling

The Study of Levee Flood Control at Pengkol River - Meteseh, Tembalang Sub District - Semarang City

The Pengkol River, situated in Meteseh, Tembalang sub-district of Semarang City, underwent three significant flood occurrences from January to February 2023. An approach involving a structural flood control system, specifically a levee, was suggested as a solution to this issue. The objective of this research is to establish the arrangement and highest elevation of the levee flood control and evaluate its effectiveness in mitigating floods. To determine flood hydrographs, a hydrological analysis was carried out using the Rainfall-Runoff simulation with the SCS Curve Number method. Additionally, a hydraulic analysis was performed using Unsteady Flow simulation with HECRAS 1D/2D hydrodynamic models to define the river body, water profile, inundation, and levee design. The hydrologic analysis indicates that the flood discharge for the 50, 100, and 200-year return periods are 209.5 m³/sec, 225.1 m³/sec, and 240.6 m³/sec, respectively. Based on the hydraulic analysis, the maximum water levels resulting from these return periods' flood discharges are +35.62 m (Q50 years), +35.77 m (Q100 years), and +36.43 m (Q200 years). The Q200 years of return period was chosen for 2D modeling because it resembles documented flood occurrences. Using the HECRAS 2D Unsteady Flow model, it was found that before the levee implementation, the flooded area within the residential zone spanned 9462 m², with a peak water level of +37.3 m. With the levee application, using an existing layout with a total length of 230 m and a top levee level of +37.5 m, flooding was effectively prevented, reducing the maximum water surface elevation to +37.17 m. This demonstrates the levee's effectiveness in preventing floods.Keywords: levee, flood, HECRAS 1D 2D model, inundation. Unsteady flow.

Satellite-based water surface slope over a small mountain river in northern China

Satellite altimetry has become a crucial alternative for measuring inland water surface elevations apart from ground observations. Water surface slope (WSS), a fundamental parameter in river discharge simulation, also allowing for assimilating satellite and ground-observed water level data with different observation positions. There are various satellite data products offering precise measurements of river water levels and estimates of WSS on a global scale. However, satellites observation of water surface over small rivers, particularly those in mountainous areas with narrow river channels, continues to be difficult. This study focuses on evaluating the precision of the ICESat-2 ATL03 photon height data in estimating WSS throughout the mountainous river section of Yongding River, which flows across the Hebei Province and Beijing City in northern China. The ICESat-2 ATL03 data with minimum along-track sampling interval of 0.7 m makes it possible to accurately estimate the WSS for small river sections that ranged from 50 to 100 m in width. A histogram-based statistical approach was used for identifying the satellite lidar photon height located right over the water surface. The eleven pairs of satellite virtual stations selected for extracting river water surface elevation and estimating river WSS demonstrated an overall validation correlation coefficient of 0.98. The range of relative error in WSS estimation is between 0.13 % and 9.02 %. The results of this work have important implications for using satellite-based measurements to estimate critical hydrological parameters in small and mountainous river basins in the absence of ground observations.