Spatial interpolation techniques play an important role in hydrology, as many point observations need to be interpolated to create continuous surfaces. Despite the availability of several tools and methods for interpolating data, not all of them work consistently for hydrologic applications. One of the techniques, the Laplace Equation, which is used in hydrology for creating flownets, has rarely been used for data interpolation. The objective of this study is to examine the efficiency of Laplace formulation (LF) in interpolating data used in hydrologic applications (hydrologic data) and compare it with other widely used methods such as inverse distance weighting (IDW), natural neighbor, and ordinary kriging. The performance of LF interpolation with other methods is evaluated using quantitative measures, including root mean squared error (RMSE) and coefficient of determination (R2) for accuracy, visual assessment for surface quality, and computational cost for operational efficiency and speed. Data related to surface elevation, river bathymetry, precipitation, temperature, and soil moisture are used for different areas in the United States. RMSE and R2 results show that LF is comparable to other methods for accuracy. LF is easy to use as it requires fewer input parameters compared to inverse distance weighting (IDW) and Kriging. Computationally, LF is faster than other methods in terms of speed when the datasets are not large. Overall, LF offers a robust alternative to existing methods for interpolating various hydrologic data. Further work is required to improve its computational efficiency.
more »
« less
This content will become publicly available on June 1, 2027
Comparative Evaluation of Machine Learning Models for Satellite Chlorophyll-a Gap Reconstruction in the Chesapeake Bay
Harmful algal blooms (HABs) are increasing in frequency in the Chesapeake Bay, posing risks to marine ecosystems, water quality, and public health. Chlorophyll-a (Chl-a) is a widely used indicator of algal biomass, and satellite observations such as Sentinel-3 Ocean and Land Color Instrument (OLCI) enable large-scale monitoring of bloom dynamics. However, cloud cover and atmospheric interference frequently introduce missing pixels in daily satellite products, reducing temporal continuity and limiting monitoring reliability. Satellite-derived chlorophyll-a (Chl-a) data exhibit substantial missingness, with daily pixel gaps ranging from approximately 52.30% to 100% (mean ≈ 88.95%). This study evaluates spatial interpolation, EOF-based, supervised machine-learning, deep-learning, and convolutional autoencoder approaches for reconstructing missing Chl-a values. Sentinel-3 OLCI Chl-a data from 2023–2024 were used for model training, while data from 2025 served as a temporally independent test set to avoid spatiotemporal leakage. To simulate cloud-induced data gaps, artificial missingness scenarios ranging from 50% to 90% were applied for the Inverse Distance Weighting (IDW) and Data Interpolating Empirical Orthogonal Functions (DINEOF) baseline approaches, while machine-learning, deep-learning, and convolutional autoencoder models were evaluated using real satellite-derived missing observations. The evaluated models include IDW, DINEOF, K-Nearest Neighbors (KNN), Random Forest (RF), Extra Trees (ET), XGBoost, a Long Short-Term Memory (LSTM) network, and a Temporal Data Interpolating Convolutional Autoencoder (Temporal DINCAE). Model performance was assessed using Root Mean Square Error (RMSE), Mean Absolute Error (MAE), prediction bias, and the coefficient of determination (R2). Results indicate that tree-based ensemble models outperform spatial interpolation and EOF-based methods, with XGBoost achieving the best overall performance (R2 ≈ 0.86; RMSE ≈ 9.61 mg m−3). The LSTM model achieved lower prediction errors (RMSE ≈ 5.87 mg m−3; MAE ≈ 2.16 mg m−3), highlighting the benefit of incorporating temporal dependencies, although with slightly reduced variance capture. The convolutional autoencoder-based Temporal DINCAE model achieved strong reconstruction performance (R2 ≈ 0.84; RMSE ≈ 11.15 mg m−3). Uncertainty quantification shows that Extra Trees tends to underestimate uncertainty with narrower prediction intervals, whereas XGBoost provides better-calibrated but wider intervals.
more »
« less
- Award ID(s):
- 2228603
- PAR ID:
- 10708515
- Publisher / Repository:
- MDPI
- Date Published:
- Journal Name:
- Remote Sensing
- Volume:
- 18
- Issue:
- 11
- ISSN:
- 2072-4292
- Page Range / eLocation ID:
- 1736
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Algal blooms in freshwater ecosystems can negatively impact aquatic and human health. Satellite remote sensing of chlorophyll a (Chl-a) is often used to help determine the severity of algal blooms. However, satellite revisit flyover schedules may not match the erratic nature of algal blooms. Studies have paired satellite and ground-based data that were not collected on the same day, assuming Chl-a concentrations did not change significantly by the flyover date. We determined the effects of an increasing time window between satellite overpass dates and field-based collection of Chl-a on algorithms for Landsat 5, Landsat 8, and Sentinel-2, using 14 years (2006–2020) of Chl-a data from 10 Oklahoma reservoirs. Multiple regression models were built, and selected statistics were used to rank the time windows. The Sentinel-2 results showed strong relationships between Chl-a and satellite data collected up to a ±5-day window. The strength of these relationships decreased beyond a ±3-day time window for Landsat 8 and a ±1-day time window for Landsat 5. Our results suggest that the time window between field sampling and satellite overpass can impact the use of satellite data for Chl-a monitoring in reservoirs. Furthermore, longer time windows can be used with higher resolution (spatial, spectral) satellites.more » « less
-
Chlorophyll-a (chl-a) is a key indicator of water quality in coastal and estuarine ecosystems. In regions like the Chesapeake Bay, elevated chl-a levels often signal the presence of harmful algal blooms and biomass accumulation. However, satellite-based chl-a observations are frequently obscured by cloud cover, limiting their use in continuous coastal monitoring. Spatial interpolation models, which estimate values at unsampled locations based on existing data, offer one solution to resolve this scarcity. Existing research has examined the effectiveness of different interpolation models for the bay’s salinity and temperature; however, few studies have investigated the application of these models for chl-a. In this study, we evaluate the performance of three kriging-based models—universal kriging (UK), ordinary kriging (OK), and empirical bayesian kriging (EBK)—and three machine learning-based models—k-Nearest Neighbor (KNN), Extra Trees (ET), and XGBoost—for interpolating chl-a concentrations. Using over 5 million remotely sensed observations across nine days in early 2025, we find that EBK exhibits the best performance among kriging-based models, while ET outperforms all other machine learning models. Our results also demonstrate that while kriging-based models outperform machine learning models in data-rich conditions, machine learning models are more adaptable and accurate for data-sparse conditions. Since data availability varies significantly from day-to-day based on cloud cover, these findings suggest that no one model is universally optimal. Integrating both approaches may offer a hybrid framework for improving the continuity and reliability of chl-a monitoring in coastal regions.more » « less
-
Harmful algal blooms (HAB) have become an increasingly pressing environmental issue in the Chesapeake Bay Region, posing a threat to surrounding ecosystems and human health through releasing toxins or sheer biomass. As the nation’s largest estuary, Chesapeake Bay is not only a center of agricultural activity such as shellfish farming but also supports a wide range of diverse habitats. Consequently, the prediction of algal blooms is essential to protecting water quality. Satellite data collection of chlorophyll-a (Chl-a) values has allowed for accessible and consistent collection of data to train prediction models; however, data is often obscured by cloud cover and other environmental factors. A lack of complete and continuous data has compromised the accuracy of predictions and thus it is essential to use machine learning for the imputation of data gaps. We used a full-coverage dataset from the year 2025 to compare the performance of Random Forest (RF) and K-Nearest Neighbors (KNN) models for imputing Chl-a values in the Chesapeake Bay Region. The best model was found to be RF with a coefficient of determination (R2) of 0.916, a notable improvement over KNN which yielded an R2 value of 0.588 when 30.46% of data was missing. The accuracy of the RF models was enhanced through iterative imputation to improve the accuracy of RF from 0.85 initially, to 0.916. Not only was the RF model more accurate, but it also exhibited a significantly lower runtime that proved essential when processing large amounts of satellite data. This study provides valuable insight into methods of adjusting RF models to better improve the quality of Chl-a data to advise water quality management.more » « less
-
The size and distribution of Phytoplankton populations are indicators of the ecological status of a water body. The chlorophyll-a (Chl-a) concentration is estimated as a proxy for the distribution of phytoplankton biomass. Remote sensing is the only practical method for the synoptic assessment of Chl-a at large spatial and temporal scales. Long-term records of ocean color data from the MODIS Aqua Sensor have proven inadequate to assess Chl-a due to the lack of a robust ocean color algorithm. Chl-a estimation in shallow and coastal water bodies has been a challenge and existing operational algorithms are only suitable for deeper water bodies. In this study, the Ocean Color 3M (OC3M) derived Chl-a concentrations were compared with observed data to assess the performance of the OC3M algorithm. Subsequently, a regression analysis between in situ Chl-a and remote sensing reflectance was performed to obtain a green-red band algorithm for coastal (case 2) water. The OC3M algorithm yielded an accurate estimate of Chl-a for deep ocean (case 1) water (RMSE = 0.007, r2 = 0.518, p < 0.001), but failed to perform well in the coastal (case 2) water of Chesapeake Bay (RMSE = 23.217, r2 = 0.009, p = 0.356). The algorithm developed in this study predicted Chl-a more accurately in Chesapeake Bay (RMSE = 4.924, r2 = 0.444, p < 0.001) than the OC3M algorithm. The study indicates a maximum band ratio formulation using green and red bands could improve the satellite estimation of Chl-a in coastal waters.more » « less
An official website of the United States government
