Projecting Future TOC in Data-Scarce Agricultural Reservoirs: A WGAN-Enhanced Framework Revealing the Importance of Dynamic Variability

Projecting Future TOC in Data-Scarce Agricultural Reservoirs: A WGAN-Enhanced Framework Revealing the Importance of Dynamic Variability | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Projecting Future TOC in Data-Scarce Agricultural Reservoirs: A WGAN-Enhanced Framework Revealing the Importance of Dynamic Variability Christian Joseph Siose, Jeongho Han, Kyung-Sook Choi, Byoung Han Choi, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8194601/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 21 Feb, 2026 Read the published version in Water Resources Management → Version 1 posted 5 You are reading this latest preprint version Abstract Climate change poses a significant threat to water quality in agricultural reservoirs. However, projecting future changes using machine learning (ML) is often hampered by the limited availability of long-term observational data. This study develops a robust framework to project future Total Organic Carbon (TOC) concentrations by combining a Wasserstein Generative Adversarial Network (WGAN) for data augmentation with high-performance ML models. We applied this framework to two contrasting South Korean reservoirs (one pristine, one polluted). Historical data (n = 60) for each was augmented using WGAN, and a suite of nine ML models was optimized and evaluated. The best-performing model, M5 Model Tree with Stochastic Gradient Boosting (M5-SGB), was then used to project TOC from 2025–2100 using a 12 Global Climate Model (GCM) ensemble from CMIP6 under SSP2-4.5 and SSP5-8.5 scenarios. Notably, the effectiveness of WGAN augmentation was model-dependent, significantly boosting the performance of certain algorithms while having a negligible effect on others. Future projections revealed a statistically significant increasing TOC trend for both reservoirs. However, the response was site-specific: the pristine but temporally unstable reservoir showed vulnerability even under the moderate-emissions scenario, while the polluted but stable reservoir only exhibited a significant trend under the high-emissions scenario. This study provides a practical framework for water quality prediction in data-limited contexts. The findings demonstrate that a reservoir's vulnerability to climate change is critically linked to its dynamic variability, not just its static water quality grade. This offers crucial insights for designing more effective, risk-based water management and monitoring strategies. Reservoir water quality Total organic carbon Climate change Machine learning model WGAN Data augmentation Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 1 Introduction Climate change is one of the most critical global challenges, exerting profound and far-reaching consequences on both natural and human systems. As global temperatures continue to rise, precipitation patterns shift, and extreme weather events become more frequent, the impact on freshwater resources is becoming increasingly evident (Delpla et al. 2009 ). These disruptions not only threaten the water availability and accessibility but also degrade water quality, posing significant risks to ecosystems, public health, and socio-economic development (Tranvik et al. 2009 ). Consequently, there is a growing urgency to understand the intricate relationship between climate change and water quality and to develop effective mitigation strategies. Reservoirs serve as crucial infrastructures for water storage and distribution, providing resources for domestic, agricultural, and industrial purposes, as well as supporting flood control and ecosystem services (Jeppesen et al. 2014 ). However, climate change exerts considerable stress on reservoir, affecting both inflow patterns and water quality dynamics (Whitehead et al. 2009 ). These pressures highlight the urgent need for innovative management strategies and reliable predictive frameworks to anticipate water quality variations under climate change scenarios. Traditionally, water quality has been assessed using parameters such as biochemical oxygen demand (BOD) and chemical oxygen demand (COD). Although these parameters are widely used, they can be influenced by specific organic or chemical conditions in aquatic environments. In contrast, total organic carbon (TOC), which represents the total amount of carbon bound in organic compounds dissolved or suspended in water, is increasingly recognized as a more comprehensive indicator of water quality (Aguilar-Torrejón et al. 2023 ). TOC quantifies the overall organic carbon content without depending on particular oxidation pathways. Thus, TOC provide a more stable and comprehensive measure of organic pollution that can better capture various forms of organic matter in a single parameter (Oh et al. 2024 ). Previous studies have also shown that higher summer temperatures and greater seasonal temperature variability are associated with higher TOC concentrations in lake sediments, whereas reduced seasonality is associated with lower TOC (Zhou et al. 2025 ). Similarly, in catchments, increases in rainfall intensity and duration drive higher TOC concentrations and export in runoff, with heavy rainfall events generating TOC fluxes several times greater than during lower-intensity storms (Delpla et al. 2011 ). Given these findings, many process-based or physically-based hydrological and water quality models have been developed and widely used to predict changes in water quality considering climate dynamics. However, their limited ability to capture complex, non-linear relationships among environmental variables hinder accurate predictions (Yan et al. 2024 ). Compared to these traditional approaches, machine learning (ML) models are capable of identifying intricate patterns in large datasets and offer a promising alternative for water quality prediction (Nallakaruppan et al. 2024 ). Despite ML models superiority in handling non-linear relationships, the application of ML to water quality prediction often faces challenges due to sparse or irregular datasets compared to those available for hydrological modeling. Given that ML are inherently data-driven, insufficient data for developing ML models can significantly undermine their predictive performance (Zheng et al. 2025 ). To address this constraint, data augmentation techniques, such as Wasserstein Generative Adversarial Networks (WGAN), has been applied to generate synthetic data and enhance model performance (Arjovsky et al. 2017 ). WGAN is a type of generative model that can learn the underlying distribution of the original data and generate realistic synthetic data. In this context, the present study aims to develop a comprehensive framework to project the impact of climate change on reservoir water quality, particularly TOC, by integrating ML models with data augmentation approach. The specific objectives are (1) to develop and compare diverse ML models for TOC prediction; (2) to apply WGAN-based data augmentation to enhance small water quality datasets; and (3) to evaluate future reservoir water quality trends based on Coupled Model Intercomparison Project Phase 6 (CMIP6) global climate models (GCMs), using the best-performing ML models. By achieving these objectives, this study will contribute to improving water quality management strategies for reservoirs, offering insight into climate-resilient infrastructure planning and providing a foundation for sustainable resource management. 2 Materials and Methods The overall research workflow is depicted in Fig. 1 . First, observational data and CMIP6 climate projections were collected and pre-processed to construct the input dataset. Next, a WGAN-based data augmentation technique was applied to address the limitations of a small dataset. Subsequently, Optuna was used to optimize and validate nine ML models. The best-performing model was then selected to project future TOC concentrations for the period 2025–2100 and to analyze the corresponding trends. The detailed procedures and data sources for each step are described in the following sections. 2.1 Study Area This study focused on two representative agricultural reservoirs in South Korea: Sacheon Reservoir in Gangneung and Docheon Reservoir in Yeongdeok (Fig. 2 ). While both are situated in the eastern coastal region of the Korea Peninsula, they were selected to represent contrasting environmental conditions. Sacheon Reservoir is located in Gangwon Province, at approximately 37°47' N, 128°53' E (Fig. 2 a). It has a total storage capacity of 2.1 × 10 6 m 3 and a watershed area of 2,280 ha. This reservoir, managed by the Korea Rural Community Corporation (KRC), is characterized by high water quality, consistently meeting the criteria for Grade A under South Korea's agricultural water quality standards, indicating its suitability for all agricultural purposes without prior treatment. The region is defined by a mountainous topography that contributes to higher precipitation and cooler temperature. Docheon Reservoir is located in Gyeongsangbuk Province at approximately 36°31′ N, 129°35′ E (Fig. 2 b), with a total storage capacity of 1.1 × 10 6 m 3 and a watershed area of 1,380 ha. Also managed by the KRC, this reservoir's water quality is generally poor, often classified as Grade E, indicating the need for significant treatment before agricultural use. The area has a relatively warmer climate with more distinct seasonal variations in rainfall compared to the Sacheon watershed. These two sites were selected based on four key criteria: (1) their roles as typical agricultural reservoirs in South Korea; (2) the availability of comprehensive water quality data; (3) their contrasting water quality statuses (high vs. low); and (4) their regional climate conditions. 2.2 Data Acquisition and Preprocessing The target variable for this study was TOC, with historical data obtained from the Rural Agricultural Water Resource Information System (RAWRIS 2024). The dataset for each reservoir spanned from 2011 to 2024 and consisted of 60 observations, collected at monthly or bimonthly intervals. To quantify the temporal variability of this historical data, the coefficient of variation (CV) was calculated for the TOC time series of each reservoir. Sacheon Reservoir exhibited greater instability (CV = 34.59%) compared to the more stable Docheon Reservoir (CV = 21.72%). Meteorological variables, including daily precipitation ( Prcp ), minimum temperature ( Tmin ), maximum temperature ( Tmax ), and relative humidity ( Rhum ), were collected from the Korea Meteorological Administration (KMA) for the North Gangneung and Yeongdok stations (Fig. 2 ). To address the temporal discrepancy between the daily meteorological data and the less frequent TOC measurements, monthly averages of each meteorological variable were calculated to server as input features for the TOC prediction models. For the climate change impact assessment, future climate projection data were retrieved from the Earth System Grid Federation node under the CMIP6 framework. We selected two contrasting scenarios: the intermediate-emissions pathways (SSP2-4.5) and the high-emissions, fossil-fueled development pathway (SSP5-8.5). SSP2-4.5 assumes radiative forcing stabilize at 4.5 W/m 2 by 2100, whereas SSP5-8.5 represents a worst-case scenario with forcing reaching 8.5 W/m 2 (Tebaldi et al. 2021 ). Analyzing these two contrasting scenarios was intended to capture a plausible range of future climate impacts on TOC. To account for inter-model uncertainty, outputs from 12 GCMs from CMIP6 were incorporated into our analysis (refer to Table S1 in the Supplementary Material). These GCMs were selected based on their established performance in simulating East Asian climate patterns and the availability of the required variables. This multi-model ensemble approach was designed to enhance the robustness of our future TOC projections by considering a wide spectrum of potential climate futures. 2.3 WGAN-based Data Augmentation and Validation Given that ML model performance relies heavily on data volume, the 60 observations available for each reservoir were deemed insufficient for robust training. To address this limitation, we augmented the dataset using a WGAN. A Generative Adversarial Network (GAN) (Goodfellow et al. 2014 ) is a deep learning architecture consisting of two competing neural networks: a generator that creates synthetic data and a discriminator that distinguishes real data from synthetic counterparts. While conventional GANs often suffer from training instability, WGAN utilizes the Wasserstein distance metric to mitigate these issues and improve convergence (Arjovsky et al. 2017 ; Gulrajani et al. 2017 ). Despite its advantages, universally established guidelines for applying WGAN to environmental data augmentation remain scarce. Drawing upon precedents in related fields, we established our methodological parameters. Previous studies demonstrated the suitability of a learning rate of 0.00005 and the feasibility of 3000 training epochs for environmental prediction tasks (Lee et al. 2021 ; Peng et al. 2022 ). Accordingly, the present study adopted these values. The original dataset was augmented fivefold, resulting in a 5:1 synthetic-to-real data ratio for subsequent model training. To validate the fidelity of the generated data, we assessed the distributional similarity between the synthetic and real datasets using three complementary nonparametric tests: the Kolmogorov–Smirnov (KS) test, which is adept at detecting differences in the central region of distributions; the Anderson–Darling (AD) test, which assigns greater weight to the tails and is thus more sensitive to discrepancies in extreme values; and the Energy Distance (ED) test, which provides a single, robust measure of divergence between multivariate samples (Kolmogorov A. 1933 ; Anderson and Darling 1952 ; Rizzo and Székely 2016 ). 2.4 Predictive Modeling Framework 2.4.1 Machine Learning Models To identify the most effective algorithm for TOC prediction, nine distinct ML models were evaluated. The comprehensive suite, selected to cover diverse modeling principles, included the decision tree (DT), random forest (RF), gradient boosting (GB), eXtreme gradient boosting (XGB), support vector regression (SVR), ridge regression (RR), LassoLars (LL), the M5 Model Tree (M5) and M5 Model Tree with Stochastic Gradient Boosting (M5-SGB). For more comprehensive details on the theoretical background of each algorithm, readers are referred to the original studies cited below. DT organizes data into a tree-like structure for regression, offering high interpretability but requiring hyperparameter tuning to prevent overfitting (Breiman et al. 1984 ; Mingers 1989 ). RF is an ensemble model that constructs multiple decision tree on a random subset and aggregates their predictions, thereby reducing overfitting and capturing complex nonlinear relationships (Breiman 2001 ; Probst et al. 2019 ). GB sequentially combines weak learners (typically trees) into a strong predictive model by iteratively training new learners on the residual errors of the previous ones (Friedman 2001 ; Natekin and Knoll 2013 ). XGB is a highly efficient and scalable implementation of GB that incorporates regularization to control overfitting and leverages parallel processing for computational speed (Chen and Guestrin 2016 ). SVR is a kernel-based method that maps features into a higher-dimensional space to model nonlinear relationships, fitting a function within a specified error margin (Lee et al. 2005 ; Awad and Khanna 2015 ). Regularized Linear Models such as LL and Ridge Regression (RR) were also included. LL uses L1 regularization for feature selection by shrinking some coefficients to zero, while RR uses L2 regularization to handle multicollinearity by reducing coefficient variance (Hoerl and Kennard 1970 ; Tibshirani 1996 ). M5 is a hybrid algorithm that places linear regression models at the leaf nodes of decision trees, enabling it to capture localized linear trends and perform limited extrapolation beyond the training data range, addressing a key limitation of standard decision trees, which cannot predict values outside the observed data (Quinlan Basser 1992 ). M5-SGB enhances the M5 model by integrating it with Stochastic Gradient Boosting to iteratively refine weak learners, thereby reducing bias and variance. This improved robustness enhances the model's potential for extrapolation, making it better suited for predicting novel, non-stationary climate-TOC relationships (Friedman 2002 ; Sattari et al. 2018 ). 2.4.2 Model Training, Optimization, and Application The historical observation data were partitioned into training (80%) and test (20%) sets. Synthetic data generated by a WGAN were added exclusively to the training sets to enhance model learning. Subsequently, the key hyperparameters for each of the nine machine-learning models were optimized using Optuna, an automated hyperparameter optimization framework (Akiba et al. 2019 ). The search space for each hyperparameter (see Table S2 in the Supplementary Material) was defined based on established practices in environmental modeling (Sipper 2022 ). The model demonstrating the best performance on the test set was selected as the final model. This optimized model was then used to project future TOC concentrations throughout 2100, driven by climate projections from 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios. 2.4.3 Performance Evaluation of ML Models ML model performance was evaluated using two widely accepted metrics in hydrological modeling: the Nash-Sutcliffe Model Efficiency (NSE) (Nash and Sutcliffe 1970 ) and the coefficient of determination (R 2 ) due to their effectiveness demonstrated by numerous studies (Moriasi et al. 2015 ; Bhatta et al. 2019 ; Alqahtani et al. 2022 ; Modi and Chintalacheruvu 2024 ). NSE is a normalized statistic that determines the relative magnitude of the residual variance compared to the measured data variance, ranging from negative infinity to 1. A NSE value of 1 indicates a perfect match, while a value of 0 suggests the model is only as accurate as the mean of the observed data. R 2 quantifies the proportion of the variance in the observed data that is predictable from the model, ranging from 0 to 1. 2.5 Trend Analysis of Future TOC Concentrations under Climate Change Scenarios This study used an ML-based approach to project reservoir TOC concentrations from 2025 to 2100 under SSP2-4.5 and SSP5-8.5 climate scenarios. To rigorously analyze these long-term projections and identify significant patterns, we adopted a comprehensive trend analysis framework. To detect the presence of monotonic trends, we applied the non-parametric Mann–Kendall (MK) test, which is robust to outliers and non-normal data distributions (Hirsch et al. 1982 ). The magnitude of any identified trend was then quantified using Sen's slope, a method that calculates the median of all possible pairwise slopes and is therefore resistant to the influence of extreme values (Sen 1968 ). Furthermore, to investigate seasonal dynamics, the projected time series were aggregated into seasonal means (spring, summer, autumn). This seasonal approach can uncover distinct trends that might be obscured in an annual analysis (Ha et al. 2022 ). Winter data were excluded from this seasonal analysis because monitoring during ice-covered months is typically absent. 3 Results 3.1 Distribution Similarity Analysis between Synthetic and Real data The distributional similarity between the synthetic and real data for both Sacheon and Docheon reservoirs were evaluated using three complementary nonparametric tests (Table 1 ). For all applied tests, the results were statistically non-significant (p ≥ 0.05), thus failing to reject the null hypothesis that both datasets originate from the same underlying distribution. This outcome provides strong evidence that WGAN-generated data faithfully reproduced the key statistical properties of the original observations. Specifically, the results from the KS and AD tests confirmed a high degree of similarity in both the central regions and the tails of the distributions, respectively. The ED test further corroborated these findings, indicating a strong overall proximity between two datasets. Collectively, these results statistically validate that the synthetic data are of sufficient fidelity to reliably augment the training set for the subsequent development of ML models. Table 1 Distribution similarity analysis between synthetic data and real data in each reservoir. Asterisks (*) denote that trends are statistically significant (p < 0.05) Reservoir Variables Kolmogorov-Smirnov Test Anderson-Darling Test Energy Distance Test Statistic p-value Statistic p-value Statistic p-value Sacheon Prcp 0.130 0.380 0.037 0.250 - - Tmax 0.111 0.575 -0.108 0.250 - - Tmin 0.099 0.716 -0.695 0.250 - - Rhum 0.079 0.907 -0.880 0.250 - - TOC 0.225 0.016 2.099 0.028 - - Overall variables - - - - 0.157 0.948 Docheon Prcp 0.104 0.653 -0.575 0.250 - - Tmax 0.104 0.653 0.491 0.250 - - Tmin 0.083 0.878 0.004 0.250 - - Rhum 0.125 0.422 0.002 0.250 - - TOC 0.100 0.702 0.541 0.250 - - Overall variables - - - - 0.257 0.597 3.2 Comparative Performance of Machine Learning Models In Sacheon Reservoir (Fig. 3 a), the M5-SGB model outperformed all other algorithms, achieving the highest performance metrics (NSE = 0.862, R² = 0.864). Following M5-SGB, the standalone M5 model also demonstrated strong performance. The evaluation revealed a nuanced impact of WGAN-based data augmentation, whereas boosting-based models such as XGB and GB showed improved performance with the augmented data, other models like RF, LL, RR, and SVR exhibited a marginal decrease. For Docheon Reservoir (Fig. 3 b), M5-SGB again proved to be the most accurate predictor, recording the highest NSE of 0.804 and an R² of 0.895. In this reservoir, most models, including the regularized linear models (LL and RR), showed slight performance enhancements with the augmented dataset. In summary, M5-SGB was consistently the superior model for TOC prediction in both reservoirs. However, the effectiveness of WGAN-based data augmentation was observed to be model-dependent. 3.3 Projection of TOC in 12 GCMs under Climate Change Scenarios The projected time series TOC concentrations for Sacheon Reservoir were analyzed using the MK test and Sen's slope (Fig. 4 a and Table 2 ). Under the high emissions SSP5-8.5 scenarios, results indicated a statistically significant increasing trend in TOC levels for all 12 GCMs. Under the intermediate SSP2-4.5 scenarios, 11 of the 12 GCMs also showed a significant increasing trend, with the FGOALS-g3 model being the sole exception. The seasonal analysis (Fig. 5 a) further revealed that projected summer TOC concentrations consistently exceeded those of spring and autumn, emphasizing the potential influence of intensified temperature fluctuations under climate change. In Docheon Reservoir (Fig. 4 b and Table 3 ), 8 of the 12 GCMs under the SSP5-8.5 scenarios exhibited a statistically significant increasing TOC trend. In contrast, under SSP2-4.5 scenario, 9 of the 12 models displayed no significant trend, with only MRI-ESM2-0, NESM3 and NorESM2-LM showing a significant upward trajectory. The seasonal analysis (Fig. 5 b) illustrated distinct periodic fluctuations, with particularly pronounced summer peaks. Table 2 Results of the MK test and Sen's slope analysis applied to the projected TOC time series (2025–2100) for Sacheon Reservoir. The analysis covers 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios. Asterisks (*) denote a statistically significant trend (p < 0.05) Model Scenario Trend Z-Score Sen's Slope ACCESS-CM2 SSP245 Increasing* 2.938 0.003 SSP585 Increasing* 5.996 0.005 ACCESS-ESM1-5 SSP245 Increasing* 5.243 0.005 SSP585 Increasing* 7.557 0.010 CMCC-ESM2 SSP245 Increasing* 6.472 0.006 SSP585 Increasing* 7.862 0.010 CanESM5 SSP245 Increasing* 4.274 0.005 SSP585 Increasing* 9.288 0.013 FGOALS-g3 SSP245 No trend 0.489 0.000 SSP585 Increasing* 3.512 0.003 KIOST-ESM SSP245 Increasing* 2.381 0.002 SSP585 Increasing* 3.763 0.003 MIROC6 SSP245 Increasing* 2.435 0.002 SSP585 Increasing* 6.292 0.007 MRI-ESM2-0 SSP245 Increasing* 2.651 0.003 SSP585 Increasing* 4.041 0.004 NESM3 SSP245 Increasing* 3.108 0.002 SSP585 Increasing* 3.009 0.002 NorESM2-LM SSP245 Increasing* 2.731 0.003 SSP585 Increasing* 6.552 0.007 NorESM2-MM SSP245 Increasing* 3.826 0.004 SSP585 Increasing* 4.050 0.004 TaiESM1 SSP245 Increasing* 5.216 0.004 SSP585 Increasing* 6.813 0.007 Table 3 Results of the MK test and Sen's slope analysis applied to the projected TOC time series (2025–2100) for Docheon Reservoir. The analysis covers 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios. Asterisks (*) denote a statistically significant trend (p < 0.05) Model Scenario Trend Z-Score Sen's Slope ACCESS-CM2 SSP245 No trend 0.525 0.000 SSP585 No trend -0.049 0.000 ACCESS-ESM1-5 SSP245 No trend -0.112 0.000 SSP585 No trend 1.135 0.001 CMCC-ESM2 SSP245 No trend 1.377 0.001 SSP585 Increasing* 5.099 0.005 CanESM5 SSP245 No trend 1.673 0.002 SSP585 Increasing* 6.283 0.007 FGOALS-g3 SSP245 No trend 0.229 0.000 SSP585 Increasing* 5.557 0.006 KIOST-ESM SSP245 No trend 1.709 0.003 SSP585 Increasing* 4.005 0.006 MIROC6 SSP245 No trend 1.180 0.001 SSP585 Increasing* 5.521 0.005 MRI-ESM2-0 SSP245 Increasing* 1.978 0.002 SSP585 Increasing* 2.238 0.002 NESM3 SSP245 Increasing* 3.252 0.004 SSP585 Increasing* 7.952 0.015 NorESM2-LM SSP245 Increasing* 2.256 0.002 SSP585 Increasing* 3.539 0.004 NorESM2-MM SSP245 No trend 1.404 0.001 SSP585 No trend 0.740 0.001 TaiESM1 SSP245 No trend -0.184 0.000 SSP585 No trend -1.601 -0.002 4 Discussion 4.1 Augmentation Effects on Model Performance The superior performance of the M5-SGB model across both reservoirs highlights the advantages of its hybrid architecture for predicting complex environmental variables like TOC. By integrating the piecewise linear functions of the M5 model tree with the iterative error-correction of Stochastic Gradient Boosting, M5-SGB effectively captures the non-linear relationships between climatic drivers and water quality. Unlike conventional tree-based models that are limited to their training range, the model's ability to perform localized linear regression at its leaf nodes likely provided a critical advantage in extrapolating future TOC concentrations under novel climate conditions. The successful application of WGAN for data augmentation aligns with findings from various domains. For instance, WGAN has been shown to generate synthetic data that closely matches real distributions in medical time series (Esteban et al. 2017 ) and electronic health records (Yoon et al. 2019 ). More specifically, within the environmental sciences, studies have reported that WGAN-based augmentation enhanced the predictive performance of models for aquatic ecosystem health (Lee et al. 2021 ) and improved streamflow prediction accuracy (López-Chacón et al. 2023 ). Our finding that WGAN provided a distinct performance boost to boosting algorithms like XGB and GB is consistent with this established potential. However, our results also introduce a critical nuance. The observation that this positive impact did not extend to inherently robust models like RF or regularized linear models suggests that data augmentation is not a universally beneficial strategy. Its utility is highly dependent on the chosen model architecture. This implies that a tailored approach is crucial, as synthetic data may offer significant benefits for data-sensitive algorithms while providing negligible or even slightly negative effects for models already well-suited to smaller datasets. This distinction is also pragmatically important. Implementing WGAN is computationally intensive and requires specialized expertise; therefore, understanding precisely which models benefit most ensures that this powerful technique is applied efficiently. Consequently, studies like the present one, which investigate the specific interactions between data augmentation methods and model architectures, are essential for developing effective and resource-conscious modeling strategies. 4.2 Contrasting Future TOC Projections between Reservoirs under Climate Change A central finding of this study is the divergent response of the two reservoirs to climate change; Sacheon Reservoir is projected to experience significant TOC increases under both moderate and high-emissions scenarios, whereas Docheon Reservoir's increasing trend is only significant under the high-emissions scenario. Sacheon Reservoir presents a pristine but unstable paradox. Despite its high-water quality (Grade A), it exhibits a high variability (CV = 34.59%), indicating considerable instability and sensitivity to irregular meteorological events or pollution inputs. This inherent vulnerability makes the reservoir highly susceptible to the systematic pressures of climate change, such as accelerated warming and intensified precipitation, which drive increased organic matter input from the watershed (Monteith et al. 2007 ; Fenner et al. 2021 ). Consequently, even the moderate climatic shift under SSP2-4.5 appears sufficient to trigger a significant rising TOC trend. Conversely, Docheon Reservoir can be characterized as polluted but stable, with poor baseline water quality but a much lower variability (CV = 21.72%). This suggests a more predictable system, possibly buffered by its already high organic load. A compelling hypothesis for its behavior is a threshold-based response. The watershed and reservoir system may effectively buffer TOC mobilization under moderate warming, but this capacity is overwhelmed once a climatic threshold is surpassed under the more extreme high-emissions conditions (Laudon et al. 2011 ). This would explain the lack of a significant trend in the SSP2-4.5 scenario and the abrupt emergence of one under SSP5-8.5. These long-term trends are further contextualized by the seasonal analysis, which showed pronounced summer peaks in TOC for both reservoirs. This is consistent with established mechanisms where rising temperatures enhance microbial activity and the decomposition of organic matter, a process expected to intensify under future warming (Sobek et al. 2007 ). These findings underscore that predicting climate change impacts on water quality requires a nuanced approach that considers not only the current state (e.g., water quality grade) but also the system's temporal variability and potential for non-linear responses. 4.3 Implications and Future Directions The projected increase in TOC concentrations carries significant practical implications for agricultural water management in South Korea. Rising TOC levels can increase the potential for summer algal blooms and elevate the costs and complexity of water treatment, particularly concerning the formation of disinfection byproducts. These findings present a direct challenge to water management bodies like the Korea Rural Community Corporation (KRC) and underscore the need for proactive, climate-adaptive watershed management strategies. Furthermore, this study's finding that the pristine Sacheon Reservoir is more unstable (high CV) than the polluted Docheon Reservoir suggests that current monitoring strategies may need re-evaluation. A water body's vulnerability to climate change may be better characterized by its dynamic variability than its current quality grade alone. Therefore, designing robust monitoring and management plans should consider not just the current state, but also the inherent instability of the system. While this study provides a robust framework for future TOC prediction, the analysis was confined to two reservoirs, and the generalization of these findings requires further research across a wider range of hydro-climatic conditions. Furthermore, all climate projections are subject to inherent uncertainty from the GCMs themselves, although this was mitigated by using a 12-model ensemble. Therefore, future research should aim to expand this modeling framework to a larger, more diverse set of reservoirs. Comparing WGAN with other generative models, such as Variational Autoencoders (VAEs), could also yield valuable insights into the most effective data augmentation techniques. Finally, incorporating additional predictor variables, such as satellite-derived land-use data or point-source pollution information, could further enhance the accuracy and explanatory power of the predictive models. 5 Conclusion This study successfully developed and validated a predictive framework combining WGAN-based data augmentation with a high performance M5-SGB model to project future TOC concentrations in data-scarce agricultural reservoirs. The ability to achieve satisfactory prediction accuracy using only readily available meteorological data is particularly encouraging. High-frequency water quality monitoring in reservoirs is often constrained by the logistical challenges of manual, boat-based sampling, making this modeling approach a valuable tool for supplementing sparse observational records. The framework proved effective, accurately capturing the complex relationships between climatic drivers and water quality. The projections reveal a significant increasing trend in TOC under future climate change, a trend that is particularly pronounced in the pristine but temporally unstable reservoir. This highlights the vulnerability of even healthy aquatic ecosystems to climatic shifts. Methodologically, this research demonstrated that the effectiveness of data augmentation is highly model-dependent, providing a critical insight for the practical application of machine learning in hydrology. In conclusion, this study provides a robust tool for water resource managers to anticipate future challenges. It underscores that effective climate adaptation strategies must consider not only the current state of water bodies but also their dynamic variability to build true resilience. Declarations Acknowledg e ments This work was supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through Intelligent Agricultural Infra Management for Climate Change Development Program, funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA) (RS-2025-02263904). Supplementary Information Supplementary data to this article is submitted at the same time. A uthor Contributions Conceptualization: Christian Joseph Siose, Jeongho Han, Kyung-Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Methodology: Jeongho Han, Kyung Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Formal analysis and investigation: Christian Joseph Siose, Jeongho Han; Writing—original draft preparation: Christian Joseph, Jeongho Han; Writing—review and editing: Kyung Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Funding acquisition: Kyoung Jae Lim; Supervision: Kyoung Jae Lim; Project administration: Kyoung Jae Lim. Data Availability The data that support the findings of this study are available from the corresponding author ( [email protected] ) upon reasonable request. Competing Interests The authors declare there is no conflict. References Aguilar-Torrejón JA, Balderas-Hernández P, Roa-Morales G, et al (2023) Relationship, importance, and development of analytical techniques: COD, BOD, and, TOC in water—An overview through time. SN Appl Sci 5 Akiba T, Sano S, Yanase T, et al (2019) Optuna: A Next-generation hyperparameter optimization framework. In: Proceedings of the ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. Association for Computing Machinery, pp 2623–2631 Alqahtani A, Shah MI, Aldrees A, Javed MF (2022) Comparative assessment of individual and ensemble machine learning models for efficient analysis of river water quality. Sustainability (Switzerland) 14: https://doi.org/10.3390/su14031183 Anderson TW, Darling DA (1952) asymptotic theory of certain “Goodness of Fit” criteria based on stochastic processes Arjovsky M, Chintala S, Bottou L (2017) Wasserstein Generative Adversarial Networks Awad M, Khanna R (2015) Support vector machines for classification. in: efficient learning machines. Apress, Berkeley, CA, pp 39–66 Bhatta B, Shrestha S, Shrestha PK, Talchabhadel R (2019) Evaluation and application of a SWAT model to assess the climate change impact on the hydrology of the Himalayan River Basin. Catena (Amst) 181: https://doi.org/10.1016/j.catena.2019.104082 Breiman, L., Friedman, J.H., Olshen, R.A., Stone, C.J., 1984. Classification and regression trees. Chapman & Hall/CRC. https://doi.org/10.1201/9781315139470. Breiman L (2001) Random forests. Mach Learn 45:5–32. https://doi.org/https://doi.org/10.1023/A:1010933404324 Chen T, Guestrin C (2016) XGBoost: A Scalable Tree Boosting System. https://doi.org/10.1145/2939672.2939785 Delpla I, Baurès E, Jung AV, Thomas O (2011) Impacts of rainfall events on runoff water quality in an agricultural environment in temperate areas. Science of the Total Environment 409:1683–1688. https://doi.org/10.1016/j.scitotenv.2011.01.033 Delpla I, Jung A-V, Baures E, et al (2009) Impacts of climate change on surface water quality in relation to drinking water production. Environ Int 35:1225–1233. https://doi.org/10.1016/j.envint.2009.07.001 Esteban C, Hyland SL, Rätsch G (2017) Real-valued (medical) time series generation with recurrent conditional gans. arXiv preprint arXiv:170602633 Fenner N, Meadham J, Jones T, et al (2021) Effects of Climate Change on Peatland Reservoirs: A DOC Perspective. Global Biogeochem Cycles 35: https://doi.org/10.1029/2021GB006992 Friedman JH (2001) 999 Reitz Lecture Greedy function approximation: A Gradient Boosting Machine 1 Friedman JH (2002) Stochastic gradient boosting. Comput Stat Data Anal 38:367–378. https://doi.org/10.1016/S0167-9473(01)00065-2 Goodfellow IJ, Pouget-Abadie J, Mirza M, et al (2014) Generative Adversarial Nets Gulrajani I, Ahmed F, Arjovsky M, et al (2017) Improved training of wasserstein gans. Adv Neural Inf Process Syst 30: Ha DW, Jung KY, Baek J, et al (2022) Trend Analysis Using Long-Term Monitoring Data of Water Quality at Churyeongcheon and Yocheon Basins. Sustainability (Switzerland) 14:. https://doi.org/10.3390/su14159770 Hirsch RM, Slack JR, Smith RA (1982) Techniques of trend analysis for monthly water quality data. Water Resour Res 18:107–121 Hoerl AE, Kennard RW (1970) Ridge regression: Biased estimation for nonorthogonal problems. Technometrics 12:55–67 Jeppesen E, Meerhoff M, Davidson TA, et al (2014) Climate change impacts on lakes: An integrated ecological perspective based on a multi-faceted approach, with special focus on shallow lakes. J Limnol 73:88–111. https://doi.org/10.4081/jlimnol.2014.844 Kolmogorov A. (1933) Sulla determinazione empirica di una legge didistribuzione. Giorn Dell’inst Ital Degli Att 4:89–91 Laudon H, Berggren M, Ågren A, et al (2011) Patterns and dynamics of dissolved organic carbon (DOC) in boreal streams: The role of processes, connectivity, and scaling. Ecosystems 14:880–893 Lee S, Kim J, Lee G, et al (2021) Prediction of aquatic ecosystem health indices through machine learning models using the wgan-based data augmentation method. Sustainability (Switzerland) 13: https://doi.org/10.3390/su131810435 Lee Y-J, Hsieh W-F, Huang C-M (2005) /spl epsi/-SSVR: a smooth support vector machine for /spl epsi/-insensitive regression. IEEE Trans Knowl Data Eng 17:678–685. https://doi.org/10.1109/TKDE.2005.77 López-Chacón SR, Salazar F, Bladé E (2023) Combining Synthetic and Observed Data to Enhance Machine Learning Model Performance for Streamflow Prediction. Water (Switzerland) 15: https://doi.org/10.3390/w15112020 Mingers, J., 1989. An empirical comparison of pruning methods for decision tree induction. Mach. Learn. 4, 227–243. https://doi.org/10.1023/A:1022604100933. Modi P, Chintalacheruvu MR (2024) Investigating River water quality assessment through non-parametric analysis: A case study of the Godavari River in India. Environmental Quality Management 33:239–264. https://doi.org/10.1002/tqem.22117 Monteith DT, Stoddard JL, Evans CD, et al (2007) Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature 450:537–540. https://doi.org/10.1038/nature06316 Moriasi DN, Gitau MW, Pai N, Daggupati P (2015) Hydrologic and water quality models: Performance measures and evaluation criteria. Trans ASABE 58:1763–1785. https://doi.org/10.13031/trans.58.10715 Nallakaruppan MK, Gangadevi E, Shri ML, et al (2024) Reliable water quality prediction and parametric analysis using explainable AI models. Sci Rep 14:7520. https://doi.org/10.1038/s41598-024-56775-y Nash JE, Sutcliffe J V (1970) River flow forecasting through conceptual models’ part i-a discussion of principles*. © North-Holland Publishing Co Natekin A, Knoll A (2013) Gradient boosting machines, a tutorial. Front Neurorobot 7: https://doi.org/10.3389/fnbot.2013.00021 Oh H, Park HY, Kim JI, et al (2024) Enhancing machine learning models for total organic carbon prediction by integrating geospatial parameters in river watersheds. Science of the Total Environment 943: https://doi.org/10.1016/j.scitotenv.2024.173743 Peng L, Li S, Sun H, Huang S (2022) A Pipe Ultrasonic Guided Wave Signal Generation Network Suitable for Data Enhancement in Deep Learning: US-WGAN. Energies (Basel) 15: https://doi.org/10.3390/en15186695 Probst P, Wright MN, Boulesteix AL (2019) Hyperparameters and tuning strategies for random forest. Wiley Interdiscip Rev Data Min Knowl Discov 9 Quinlan Basser JR (1992) Learning with continuous classes. World Scientiic Rizzo ML, Székely GJ (2016) Energy distance. Wiley Interdiscip Rev Comput Stat 8:27–38 Rural Agricultural Water Resource Information System (RAWRIS) (2024) Rural Agricultural Water Resource Information System. Available at: https://rawris.ekr.or.kr Accessed 11 Nov 2024. Sattari MT, Pal M, Mirabbasi R, Abraham J (2018) Ensemble of M5 Model Tree-Based Modelling of Sodium Adsorption Ratio Sen PK (1968) Estimates of the regression coefficient based on Kendall’s tau. J Am Stat Assoc 63:1379–1389 Sipper M (2022) High Per Parameter: A Large-Scale Study of Hyperparameter Tuning for Machine Learning Algorithms. Algorithms 15:315. https://doi.org/10.3390/a15090315 Sobek S, Tranvik LJ, Prairie YT, et al (2007) Patterns and regulation of dissolved organic carbon: An analysis of 7,500 widely distributed lakes. Limnol Oceanogr 52:1208–1219. https://doi.org/10.4319/lo.2007.52.3.1208 Tebaldi C, Debeire K, Eyring V, et al (2021) Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6. Earth System Dynamics 12:253–293. https://doi.org/10.5194/esd-12-253-2021 Tibshirani R (1996) Regression Shrinkage and Selection via the Lasso Tranvik LJ, Downing JA, Cotner JB, et al (2009) Lakes and reservoirs as regulators of carbon cycling and climate. Limnol Oceanogr 54:2298–2314. https://doi.org/10.4319/lo.2009.54.6_part_2.2298 Whitehead PG, Wilby RL, Battarbee RW, et al (2009) A review of the potential impacts of climate change on surface water quality. Hydrological Sciences Journal 54:101–123. https://doi.org/10.1623/hysj.54.1.101 Yan X, Zhang T, Du W, et al (2024) A comprehensive review of machine learning for water quality prediction over the past five years. J Mar Sci Eng 12 Yoon J, Jarrett D, Van der Schaar M (2019) Time-series generative adversarial networks. Adv Neural Inf Process Syst 32: Zheng Y, Zhang X, Zhou Y, et al (2025) Deep representation learning enables cross-basin water quality prediction under data-scarce conditions. NPJ Clean Water 8: https://doi.org/10.1038/s41545-025-00466-2 Zhou S, Long H, Chen W, et al (2025) Temperature seasonality regulates organic carbon burial in lake. Nature Communications 16: https://doi.org/10.1038/s41467-025-56399-4 Supplementary Files SupplementaryInformation.docx highlights.docx Cite Share Download PDF Status: Published Journal Publication published 21 Feb, 2026 Read the published version in Water Resources Management → Version 1 posted Editorial decision: Major revisions 03 Dec, 2025 Reviewers agreed at journal 28 Nov, 2025 Reviewers invited by journal 25 Nov, 2025 Editor assigned by journal 25 Nov, 2025 First submitted to journal 24 Nov, 2025 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8194601","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":550608673,"identity":"a0fe1819-54da-43b0-a02d-b2d5650f9ac2","order_by":0,"name":"Christian Joseph Siose","email":"","orcid":"","institution":"Kangwon National University","correspondingAuthor":false,"prefix":"","firstName":"Christian","middleName":"Joseph","lastName":"Siose","suffix":""},{"id":550608674,"identity":"5d276b98-777e-4e96-b722-371c75ce7b15","order_by":1,"name":"Jeongho Han","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABIElEQVRIiWNgGAWjYBACAygGEswHDD4cQJKSIKyFLaFwBrFaoAweg888xGgxZz97oOBDxR277ewNhpttztjlyUckH3vAUGPHIDn7AFYtlj15CYYzzjxL3tlzINk450ZyseGNtHQDhmPJDNJ8CdgddiDHwJi37XCywY2EY8Y5H5gTN87OMZNgYDvAIMeDwy/n3xgY/wVrSWz/bfGhHqrlHx4tN4C2MLYdtjO4kcxgzHDjcOJ8aaAWxrYDDNI4tbwxMOw5czjB4MwxBiDjeOIG+WdpEol9yTySPbgclmNm8KPisL3B8f4PBj+OVSfO7zl8TOLDNzs5iTPYtQABGyhuEhsQAQIkEhgYcDgLDJgfAAl7OFe+AZfCUTAKRsEoGKkAAFBPZhnE8geGAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0002-6735-7755","institution":"Kangwon National University","correspondingAuthor":true,"prefix":"","firstName":"Jeongho","middleName":"","lastName":"Han","suffix":""},{"id":550608675,"identity":"3445b433-5a83-49c5-ad83-3bef0e2c1c09","order_by":2,"name":"Kyung-Sook Choi","email":"","orcid":"","institution":"Kyungpook National University","correspondingAuthor":false,"prefix":"","firstName":"Kyung-Sook","middleName":"","lastName":"Choi","suffix":""},{"id":550608676,"identity":"ff07fb40-339e-4d19-86a3-fe35e717ce17","order_by":3,"name":"Byoung Han Choi","email":"","orcid":"","institution":"Korea Rural Community Corporation","correspondingAuthor":false,"prefix":"","firstName":"Byoung","middleName":"Han","lastName":"Choi","suffix":""},{"id":550608677,"identity":"ab2169ac-7f8f-4287-b689-54dc3691a533","order_by":4,"name":"Kyoung Jae Lim","email":"","orcid":"","institution":"Kangwon National University","correspondingAuthor":false,"prefix":"","firstName":"Kyoung","middleName":"Jae","lastName":"Lim","suffix":""}],"badges":[],"createdAt":"2025-11-24 15:20:21","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-8194601/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-8194601/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s11269-026-04503-8","type":"published","date":"2026-02-21T15:59:22+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":97140510,"identity":"fe20af71-5b61-4984-bf59-a884db4e4cdd","added_by":"auto","created_at":"2025-12-01 10:05:07","extension":"docx","order_by":0,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1726992,"visible":true,"origin":"","legend":"","description":"","filename":"Figures.docx","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/c936a92dad854278deb2ef97.docx"},{"id":97111253,"identity":"b24ed395-6774-4379-8fb4-5788acc7af9f","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":50604,"visible":true,"origin":"","legend":"","description":"","filename":"Table.docx","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/43c2306e855308d005d412a0.docx"},{"id":97111256,"identity":"d56dd034-25ee-461c-8a8e-ee06529f1d34","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"xml","order_by":3,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":9896,"visible":true,"origin":"","legend":"","description":"","filename":"warmWARMD2503568.xml","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/1e5fad253d37f545e3b703d1.xml"},{"id":97111261,"identity":"981e5004-5144-4175-a738-c720e57c919f","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"xml","order_by":4,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1118,"visible":true,"origin":"","legend":"","description":"","filename":"WARMD250356846690.go.xml","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/d3e2768acb985474fda17fdb.xml"},{"id":97111274,"identity":"ea41a289-6739-43d5-95c2-c13f8f7c358a","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"xml","order_by":5,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":832,"visible":true,"origin":"","legend":"","description":"","filename":"WARMD2503568Import.xml","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/71496cf3e50b0565a309f0a5.xml"},{"id":97111269,"identity":"7b57e5f0-3803-4f9a-95ab-dd2e98b71abd","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"xml","order_by":9,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":160657,"visible":true,"origin":"","legend":"","description":"","filename":"WARMD25035680enriched.xml","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/d436c439fe48b5232ca16417.xml"},{"id":97142877,"identity":"578f8b71-004f-476c-ab57-243056d38810","added_by":"auto","created_at":"2025-12-01 10:08:02","extension":"jpeg","order_by":10,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":212463,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/3e1b6412dab2a5063a767860.jpeg"},{"id":97111270,"identity":"510e012f-fdfc-495b-80d7-84dec1c676b8","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":11,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":405828,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage10.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/39293c14d443355d0412583d.jpeg"},{"id":97111266,"identity":"5b00e8e2-844c-4b99-a9a7-ead6535b15c6","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":12,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":442004,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/969fdba3363cc8a107f316e2.jpeg"},{"id":97140312,"identity":"7ef97e8c-e96c-40df-9289-c8508b681b72","added_by":"auto","created_at":"2025-12-01 10:04:35","extension":"jpeg","order_by":13,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":471956,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/66e2aa0f3d3daaa5e8a021c7.jpeg"},{"id":97111267,"identity":"e83ea3bf-c513-4c49-85e7-bbfe54bf20b5","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":14,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1078010,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage4.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/a0fa0aa0e9aa4bef0a877a91.jpeg"},{"id":97111275,"identity":"3bbea663-a109-4fda-8ea9-b1253c798a11","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":15,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":405828,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage10.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/8e000a20a88cd9b178172e51.jpeg"},{"id":97111271,"identity":"6b488086-bcec-4c7d-8d87-fc07bebe01f5","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":16,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":212463,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage1.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/26b8641f3539b1169c870189.jpeg"},{"id":97142224,"identity":"d23b6766-9c9c-40f7-b1ab-079bf749b94c","added_by":"auto","created_at":"2025-12-01 10:07:26","extension":"jpeg","order_by":17,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":442004,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage2.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/0bea3106fef20fadf851d7c0.jpeg"},{"id":97111265,"identity":"d36f9f4a-6874-4064-9079-7c0e44275460","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":18,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":471956,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage3.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/e288716d4df1a9603a1b047a.jpeg"},{"id":97111268,"identity":"ef2dc261-2773-4161-ad62-ea9cfb59656b","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"jpeg","order_by":19,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":1078010,"visible":true,"origin":"","legend":"","description":"","filename":"floatimage4.jpeg","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/06efeeca66b05d90e0dbd3c8.jpeg"},{"id":97111273,"identity":"2ee894ef-80e6-4a56-beeb-12fc0646a2b5","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"png","order_by":20,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":76268,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/036ea89472f661ef3a0ff746.png"},{"id":97111263,"identity":"fc309158-aab1-4de7-aea8-3715ed218a80","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"png","order_by":21,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":85225,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage10.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/dd61845e1831031a6384d45f.png"},{"id":97111278,"identity":"ccbe8be2-33d2-4e2c-b040-f8dfdf09e551","added_by":"auto","created_at":"2025-12-01 06:29:23","extension":"png","order_by":22,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":233056,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/bda9add7194c7782210fdc50.png"},{"id":97111286,"identity":"b405b2e5-72d9-4732-87a5-961c37de8bb7","added_by":"auto","created_at":"2025-12-01 06:29:23","extension":"png","order_by":23,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":101054,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/7bf4592b7d379659a5486fa4.png"},{"id":97142900,"identity":"4f14f4e0-b0f5-4f26-9615-91821d616888","added_by":"auto","created_at":"2025-12-01 10:08:04","extension":"png","order_by":24,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":215244,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/b24da7c106bc5da546c4119d.png"},{"id":97111279,"identity":"b76b092f-a9b9-49c5-83be-4cf614c2a305","added_by":"auto","created_at":"2025-12-01 06:29:23","extension":"png","order_by":25,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":85225,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage10.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/f8fd50a3aeba45eeacb04683.png"},{"id":97142864,"identity":"62711414-3fa4-4ccf-a8b0-49cb5f6f2550","added_by":"auto","created_at":"2025-12-01 10:08:00","extension":"png","order_by":26,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":76268,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/6feb962e7c8ea6f332e21285.png"},{"id":97111284,"identity":"56d97369-bf53-45fe-b381-e8ff17275026","added_by":"auto","created_at":"2025-12-01 06:29:23","extension":"png","order_by":27,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":233056,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/e81164380d659635293549af.png"},{"id":97142236,"identity":"20519b94-d372-4e84-b9e5-c824779110c2","added_by":"auto","created_at":"2025-12-01 10:07:26","extension":"png","order_by":28,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":101054,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/d858df8a6b68fa39f089d82e.png"},{"id":97111272,"identity":"92518d9a-8e5b-443e-ab2b-ed2f2fd169a5","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"png","order_by":29,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":215244,"visible":true,"origin":"","legend":"","description":"","filename":"Onlinefloatimage9.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/f0215e87f034183f32bd6028.png"},{"id":97111281,"identity":"9a83ed59-9b6f-4130-9ecc-acb9942a89b6","added_by":"auto","created_at":"2025-12-01 06:29:23","extension":"xml","order_by":30,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":160389,"visible":true,"origin":"","legend":"","description":"","filename":"WARMD25035680structuring.xml","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/a6b2c779b039e4e9595f05aa.xml"},{"id":97140508,"identity":"52140541-eb76-444a-9795-b42c2f467bcc","added_by":"auto","created_at":"2025-12-01 10:05:07","extension":"html","order_by":31,"title":"","display":"","copyAsset":false,"role":"acdc-reference","size":166685,"visible":true,"origin":"","legend":"","description":"","filename":"earlyproof.html","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/d0d8b0ba8472f342b231a1dc.html"},{"id":97111288,"identity":"a4de3cb4-62bb-4974-bab4-e4d2f19de61b","added_by":"auto","created_at":"2025-12-01 06:29:35","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":198836,"visible":true,"origin":"","legend":"\u003cp\u003eMethodological workflow for developing TOC prediction models and projecting future concentrations under climate change scenarios, detailing the key steps from data preprocessing to model application\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/bd2495baf4bc716dee1ec846.png"},{"id":97111258,"identity":"ce6af269-94d1-4e8d-988b-3f3b4edab8ee","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":493838,"visible":true,"origin":"","legend":"\u003cp\u003eStudy area in South Korea, showing the locations of (a) Sacheon and (b) Docheon Reservoirs\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/80b32e3a292dc06a397a97a9.png"},{"id":97140316,"identity":"bea71886-dbcf-4cc3-9319-7d2e08bcdb41","added_by":"auto","created_at":"2025-12-01 10:04:36","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":169236,"visible":true,"origin":"","legend":"\u003cp\u003ePerformance of ML models trained on the original and WGAN-augmented datasets for (a) Sacheon and (b) Docheon Reservoirs\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/1a4156e6e3fa31f67254a604.png"},{"id":97142823,"identity":"3a494514-e98f-4c09-bc91-dccc8f05b4d1","added_by":"auto","created_at":"2025-12-01 10:07:58","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":341647,"visible":true,"origin":"","legend":"\u003cp\u003eProjected annual TOC concentrations from 2025-2100 for (a) Sacheon and (b) Docheon Reservoirs, based on 12 GCMs. For each GCM, the faint line illustrates the projected annual values, while the bold line is the long-term trend estimated by Sen's slope. Trend significance is determined by the MK test (p \u0026lt; 0.05), as detailed in Tables 2 and 3\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/33845b442c023b26db6c54fd.png"},{"id":97111254,"identity":"e539153d-ccec-4ea4-9781-8aa833bcdc48","added_by":"auto","created_at":"2025-12-01 06:29:22","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":124200,"visible":true,"origin":"","legend":"\u003cp\u003eProjected seasonal TOC cycles for the near-term (2050-2055) and long-term (2080-2085) in (a) Sacheon and (b) Docheon Reservoirs. Lines represent the ensemble mean of 12 GCMs for each period.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/05ef82818198b1a6583eacdf.png"},{"id":103252647,"identity":"82a20e5e-732d-4227-9018-84b466435395","added_by":"auto","created_at":"2026-02-23 16:15:38","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2227871,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/9ee328ff-7c06-454a-aadc-2d5375a82446.pdf"},{"id":97141693,"identity":"abf3d336-8447-4659-b2c1-bf3697b50fd8","added_by":"auto","created_at":"2025-12-01 10:06:53","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":40358,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryInformation.docx","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/14563a42c3789fc2cb893679.docx"},{"id":97141259,"identity":"86e99981-db43-47f7-be9f-c0272fdc6b66","added_by":"auto","created_at":"2025-12-01 10:06:28","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":126067,"visible":true,"origin":"","legend":"","description":"","filename":"highlights.docx","url":"https://assets-eu.researchsquare.com/files/rs-8194601/v1/f4137afcf8508ed6fae3e5af.docx"}],"financialInterests":"","formattedTitle":"Projecting Future TOC in Data-Scarce Agricultural Reservoirs: A WGAN-Enhanced Framework Revealing the Importance of Dynamic Variability","fulltext":[{"header":"1 Introduction","content":"\u003cp\u003eClimate change is one of the most critical global challenges, exerting profound and far-reaching consequences on both natural and human systems. As global temperatures continue to rise, precipitation patterns shift, and extreme weather events become more frequent, the impact on freshwater resources is becoming increasingly evident (Delpla et al. \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). These disruptions not only threaten the water availability and accessibility but also degrade water quality, posing significant risks to ecosystems, public health, and socio-economic development (Tranvik et al. \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). Consequently, there is a growing urgency to understand the intricate relationship between climate change and water quality and to develop effective mitigation strategies.\u003c/p\u003e\u003cp\u003eReservoirs serve as crucial infrastructures for water storage and distribution, providing resources for domestic, agricultural, and industrial purposes, as well as supporting flood control and ecosystem services (Jeppesen et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). However, climate change exerts considerable stress on reservoir, affecting both inflow patterns and water quality dynamics (Whitehead et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). These pressures highlight the urgent need for innovative management strategies and reliable predictive frameworks to anticipate water quality variations under climate change scenarios.\u003c/p\u003e\u003cp\u003eTraditionally, water quality has been assessed using parameters such as biochemical oxygen demand (BOD) and chemical oxygen demand (COD). Although these parameters are widely used, they can be influenced by specific organic or chemical conditions in aquatic environments. In contrast, total organic carbon (TOC), which represents the total amount of carbon bound in organic compounds dissolved or suspended in water, is increasingly recognized as a more comprehensive indicator of water quality (Aguilar-Torrej\u0026oacute;n et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). TOC quantifies the overall organic carbon content without depending on particular oxidation pathways. Thus, TOC provide a more stable and comprehensive measure of organic pollution that can better capture various forms of organic matter in a single parameter (Oh et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e\u003cp\u003ePrevious studies have also shown that higher summer temperatures and greater seasonal temperature variability are associated with higher TOC concentrations in lake sediments, whereas reduced seasonality is associated with lower TOC (Zhou et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2025\u003c/span\u003e). Similarly, in catchments, increases in rainfall intensity and duration drive higher TOC concentrations and export in runoff, with heavy rainfall events generating TOC fluxes several times greater than during lower-intensity storms (Delpla et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Given these findings, many process-based or physically-based hydrological and water quality models have been developed and widely used to predict changes in water quality considering climate dynamics. However, their limited ability to capture complex, non-linear relationships among environmental variables hinder accurate predictions (Yan et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eCompared to these traditional approaches, machine learning (ML) models are capable of identifying intricate patterns in large datasets and offer a promising alternative for water quality prediction (Nallakaruppan et al. \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2024\u003c/span\u003e). Despite ML models superiority in handling non-linear relationships, the application of ML to water quality prediction often faces challenges due to sparse or irregular datasets compared to those available for hydrological modeling. Given that ML are inherently data-driven, insufficient data for developing ML models can significantly undermine their predictive performance (Zheng et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2025\u003c/span\u003e). To address this constraint, data augmentation techniques, such as Wasserstein Generative Adversarial Networks (WGAN), has been applied to generate synthetic data and enhance model performance (Arjovsky et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). WGAN is a type of generative model that can learn the underlying distribution of the original data and generate realistic synthetic data.\u003c/p\u003e\u003cp\u003eIn this context, the present study aims to develop a comprehensive framework to project the impact of climate change on reservoir water quality, particularly TOC, by integrating ML models with data augmentation approach. The specific objectives are (1) to develop and compare diverse ML models for TOC prediction; (2) to apply WGAN-based data augmentation to enhance small water quality datasets; and (3) to evaluate future reservoir water quality trends based on Coupled Model Intercomparison Project Phase 6 (CMIP6) global climate models (GCMs), using the best-performing ML models. By achieving these objectives, this study will contribute to improving water quality management strategies for reservoirs, offering insight into climate-resilient infrastructure planning and providing a foundation for sustainable resource management.\u003c/p\u003e"},{"header":"2 Materials and Methods","content":"\u003cp\u003eThe overall research workflow is depicted in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. First, observational data and CMIP6 climate projections were collected and pre-processed to construct the input dataset. Next, a WGAN-based data augmentation technique was applied to address the limitations of a small dataset. Subsequently, Optuna was used to optimize and validate nine ML models. The best-performing model was then selected to project future TOC concentrations for the period 2025\u0026ndash;2100 and to analyze the corresponding trends. The detailed procedures and data sources for each step are described in the following sections.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e\u003ch2\u003e2.1 Study Area\u003c/h2\u003e\u003cp\u003eThis study focused on two representative agricultural reservoirs in South Korea: Sacheon Reservoir in Gangneung and Docheon Reservoir in Yeongdeok (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). While both are situated in the eastern coastal region of the Korea Peninsula, they were selected to represent contrasting environmental conditions.\u003c/p\u003e\u003cp\u003eSacheon Reservoir is located in Gangwon Province, at approximately 37\u0026deg;47' N, 128\u0026deg;53' E (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). It has a total storage capacity of 2.1 \u0026times; 10\u003csup\u003e6\u003c/sup\u003e m\u003csup\u003e3\u003c/sup\u003e and a watershed area of 2,280 ha. This reservoir, managed by the Korea Rural Community Corporation (KRC), is characterized by high water quality, consistently meeting the criteria for Grade A under South Korea's agricultural water quality standards, indicating its suitability for all agricultural purposes without prior treatment. The region is defined by a mountainous topography that contributes to higher precipitation and cooler temperature.\u003c/p\u003e\u003cp\u003eDocheon Reservoir is located in Gyeongsangbuk Province at approximately 36\u0026deg;31\u0026prime; N, 129\u0026deg;35\u0026prime; E (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb), with a total storage capacity of 1.1 \u0026times; 10\u003csup\u003e6\u003c/sup\u003e m\u003csup\u003e3\u003c/sup\u003e and a watershed area of 1,380 ha. Also managed by the KRC, this reservoir's water quality is generally poor, often classified as Grade E, indicating the need for significant treatment before agricultural use. The area has a relatively warmer climate with more distinct seasonal variations in rainfall compared to the Sacheon watershed.\u003c/p\u003e\u003cp\u003eThese two sites were selected based on four key criteria: (1) their roles as typical agricultural reservoirs in South Korea; (2) the availability of comprehensive water quality data; (3) their contrasting water quality statuses (high vs. low); and (4) their regional climate conditions.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec4\" class=\"Section2\"\u003e\u003ch2\u003e2.2 Data Acquisition and Preprocessing\u003c/h2\u003e\u003cp\u003eThe target variable for this study was TOC, with historical data obtained from the Rural Agricultural Water Resource Information System (RAWRIS 2024). The dataset for each reservoir spanned from 2011 to 2024 and consisted of 60 observations, collected at monthly or bimonthly intervals.\u003c/p\u003e\u003cp\u003eTo quantify the temporal variability of this historical data, the coefficient of variation (CV) was calculated for the TOC time series of each reservoir. Sacheon Reservoir exhibited greater instability (CV\u0026thinsp;=\u0026thinsp;34.59%) compared to the more stable Docheon Reservoir (CV\u0026thinsp;=\u0026thinsp;21.72%).\u003c/p\u003e\u003cp\u003eMeteorological variables, including daily precipitation (\u003cem\u003ePrcp\u003c/em\u003e), minimum temperature (\u003cem\u003eTmin\u003c/em\u003e), maximum temperature (\u003cem\u003eTmax\u003c/em\u003e), and relative humidity (\u003cem\u003eRhum\u003c/em\u003e), were collected from the Korea Meteorological Administration (KMA) for the North Gangneung and Yeongdok stations (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). To address the temporal discrepancy between the daily meteorological data and the less frequent TOC measurements, monthly averages of each meteorological variable were calculated to server as input features for the TOC prediction models.\u003c/p\u003e\u003cp\u003eFor the climate change impact assessment, future climate projection data were retrieved from the Earth System Grid Federation node under the CMIP6 framework. We selected two contrasting scenarios: the intermediate-emissions pathways (SSP2-4.5) and the high-emissions, fossil-fueled development pathway (SSP5-8.5).\u003c/p\u003e\u003cp\u003eSSP2-4.5 assumes radiative forcing stabilize at 4.5 W/m\u003csup\u003e2\u003c/sup\u003e by 2100, whereas SSP5-8.5 represents a worst-case scenario with forcing reaching 8.5 W/m\u003csup\u003e2\u003c/sup\u003e (Tebaldi et al. \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Analyzing these two contrasting scenarios was intended to capture a plausible range of future climate impacts on TOC.\u003c/p\u003e\u003cp\u003eTo account for inter-model uncertainty, outputs from 12 GCMs from CMIP6 were incorporated into our analysis (refer to Table \u003cspan refid=\"MOESM1\" class=\"InternalRef\"\u003eS1\u003c/span\u003e in the Supplementary Material). These GCMs were selected based on their established performance in simulating East Asian climate patterns and the availability of the required variables. This multi-model ensemble approach was designed to enhance the robustness of our future TOC projections by considering a wide spectrum of potential climate futures.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec5\" class=\"Section2\"\u003e\u003ch2\u003e2.3 WGAN-based Data Augmentation and Validation\u003c/h2\u003e\u003cp\u003eGiven that ML model performance relies heavily on data volume, the 60 observations available for each reservoir were deemed insufficient for robust training. To address this limitation, we augmented the dataset using a WGAN. A Generative Adversarial Network (GAN) (Goodfellow et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) is a deep learning architecture consisting of two competing neural networks: a generator that creates synthetic data and a discriminator that distinguishes real data from synthetic counterparts. While conventional GANs often suffer from training instability, WGAN utilizes the Wasserstein distance metric to mitigate these issues and improve convergence (Arjovsky et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Gulrajani et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eDespite its advantages, universally established guidelines for applying WGAN to environmental data augmentation remain scarce. Drawing upon precedents in related fields, we established our methodological parameters. Previous studies demonstrated the suitability of a learning rate of 0.00005 and the feasibility of 3000 training epochs for environmental prediction tasks (Lee et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Peng et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Accordingly, the present study adopted these values. The original dataset was augmented fivefold, resulting in a 5:1 synthetic-to-real data ratio for subsequent model training.\u003c/p\u003e\u003cp\u003eTo validate the fidelity of the generated data, we assessed the distributional similarity between the synthetic and real datasets using three complementary nonparametric tests: the Kolmogorov\u0026ndash;Smirnov (KS) test, which is adept at detecting differences in the central region of distributions; the Anderson\u0026ndash;Darling (AD) test, which assigns greater weight to the tails and is thus more sensitive to discrepancies in extreme values; and the Energy Distance (ED) test, which provides a single, robust measure of divergence between multivariate samples (Kolmogorov A. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e1933\u003c/span\u003e; Anderson and Darling \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e1952\u003c/span\u003e; Rizzo and Sz\u0026eacute;kely \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2016\u003c/span\u003e).\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec6\" class=\"Section2\"\u003e\u003ch2\u003e2.4 Predictive Modeling Framework\u003c/h2\u003e\u003cdiv id=\"Sec7\" class=\"Section3\"\u003e\u003ch2\u003e2.4.1 Machine Learning Models\u003c/h2\u003e\u003cp\u003eTo identify the most effective algorithm for TOC prediction, nine distinct ML models were evaluated. The comprehensive suite, selected to cover diverse modeling principles, included the decision tree (DT), random forest (RF), gradient boosting (GB), eXtreme gradient boosting (XGB), support vector regression (SVR), ridge regression (RR), LassoLars (LL), the M5 Model Tree (M5) and M5 Model Tree with Stochastic Gradient Boosting (M5-SGB). For more comprehensive details on the theoretical background of each algorithm, readers are referred to the original studies cited below.\u003c/p\u003e\u003cp\u003eDT organizes data into a tree-like structure for regression, offering high interpretability but requiring hyperparameter tuning to prevent overfitting (Breiman et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e1984\u003c/span\u003e; Mingers \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e1989\u003c/span\u003e). RF is an ensemble model that constructs multiple decision tree on a random subset and aggregates their predictions, thereby reducing overfitting and capturing complex nonlinear relationships (Breiman \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2001\u003c/span\u003e; Probst et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). GB sequentially combines weak learners (typically trees) into a strong predictive model by iteratively training new learners on the residual errors of the previous ones (Friedman \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2001\u003c/span\u003e; Natekin and Knoll \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). XGB is a highly efficient and scalable implementation of GB that incorporates regularization to control overfitting and leverages parallel processing for computational speed (Chen and Guestrin \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2016\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eSVR is a kernel-based method that maps features into a higher-dimensional space to model nonlinear relationships, fitting a function within a specified error margin (Lee et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2005\u003c/span\u003e; Awad and Khanna \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Regularized Linear Models such as LL and Ridge Regression (RR) were also included. LL uses L1 regularization for feature selection by shrinking some coefficients to zero, while RR uses L2 regularization to handle multicollinearity by reducing coefficient variance (Hoerl and Kennard \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e1970\u003c/span\u003e; Tibshirani \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e1996\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eM5 is a hybrid algorithm that places linear regression models at the leaf nodes of decision trees, enabling it to capture localized linear trends and perform limited extrapolation beyond the training data range, addressing a key limitation of standard decision trees, which cannot predict values outside the observed data (Quinlan Basser \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e1992\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eM5-SGB enhances the M5 model by integrating it with Stochastic Gradient Boosting to iteratively refine weak learners, thereby reducing bias and variance. This improved robustness enhances the model's potential for extrapolation, making it better suited for predicting novel, non-stationary climate-TOC relationships (Friedman \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2002\u003c/span\u003e; Sattari et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec8\" class=\"Section3\"\u003e\u003ch2\u003e2.4.2 Model Training, Optimization, and Application\u003c/h2\u003e\u003cp\u003eThe historical observation data were partitioned into training (80%) and test (20%) sets. Synthetic data generated by a WGAN were added exclusively to the training sets to enhance model learning. Subsequently, the key hyperparameters for each of the nine machine-learning models were optimized using Optuna, an automated hyperparameter optimization framework (Akiba et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). The search space for each hyperparameter (see Table \u003cspan refid=\"MOESM2\" class=\"InternalRef\"\u003eS2\u003c/span\u003e in the Supplementary Material) was defined based on established practices in environmental modeling (Sipper \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). The model demonstrating the best performance on the test set was selected as the final model. This optimized model was then used to project future TOC concentrations throughout 2100, driven by climate projections from 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec9\" class=\"Section3\"\u003e\u003ch2\u003e2.4.3 Performance Evaluation of ML Models\u003c/h2\u003e\u003cp\u003eML model performance was evaluated using two widely accepted metrics in hydrological modeling: the Nash-Sutcliffe Model Efficiency (NSE) (Nash and Sutcliffe \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e1970\u003c/span\u003e) and the coefficient of determination (R\u003csup\u003e2\u003c/sup\u003e) due to their effectiveness demonstrated by numerous studies (Moriasi et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Bhatta et al. \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2019\u003c/span\u003e; Alqahtani et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Modi and Chintalacheruvu \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2024\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eNSE is a normalized statistic that determines the relative magnitude of the residual variance compared to the measured data variance, ranging from negative infinity to 1. A NSE value of 1 indicates a perfect match, while a value of 0 suggests the model is only as accurate as the mean of the observed data. R\u003csup\u003e2\u003c/sup\u003e quantifies the proportion of the variance in the observed data that is predictable from the model, ranging from 0 to 1.\u003c/p\u003e\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e\u003ch2\u003e2.5 Trend Analysis of Future TOC Concentrations under Climate Change Scenarios\u003c/h2\u003e\u003cp\u003eThis study used an ML-based approach to project reservoir TOC concentrations from 2025 to 2100 under SSP2-4.5 and SSP5-8.5 climate scenarios. To rigorously analyze these long-term projections and identify significant patterns, we adopted a comprehensive trend analysis framework.\u003c/p\u003e\u003cp\u003eTo detect the presence of monotonic trends, we applied the non-parametric Mann\u0026ndash;Kendall (MK) test, which is robust to outliers and non-normal data distributions (Hirsch et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e1982\u003c/span\u003e). The magnitude of any identified trend was then quantified using Sen's slope, a method that calculates the median of all possible pairwise slopes and is therefore resistant to the influence of extreme values (Sen \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e1968\u003c/span\u003e).\u003c/p\u003e\u003cp\u003eFurthermore, to investigate seasonal dynamics, the projected time series were aggregated into seasonal means (spring, summer, autumn). This seasonal approach can uncover distinct trends that might be obscured in an annual analysis (Ha et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Winter data were excluded from this seasonal analysis because monitoring during ice-covered months is typically absent.\u003c/p\u003e\u003c/div\u003e"},{"header":"3 Results","content":"\u003cdiv id=\"Sec12\" class=\"Section2\"\u003e\u003ch2\u003e3.1 Distribution Similarity Analysis between Synthetic and Real data\u003c/h2\u003e\u003cp\u003eThe distributional similarity between the synthetic and real data for both Sacheon and Docheon reservoirs were evaluated using three complementary nonparametric tests (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e). For all applied tests, the results were statistically non-significant (p\u0026thinsp;\u0026ge;\u0026thinsp;0.05), thus failing to reject the null hypothesis that both datasets originate from the same underlying distribution.\u003c/p\u003e\u003cp\u003eThis outcome provides strong evidence that WGAN-generated data faithfully reproduced the key statistical properties of the original observations. Specifically, the results from the KS and AD tests confirmed a high degree of similarity in both the central regions and the tails of the distributions, respectively. The ED test further corroborated these findings, indicating a strong overall proximity between two datasets. Collectively, these results statistically validate that the synthetic data are of sufficient fidelity to reliably augment the training set for the subsequent development of ML models.\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eDistribution similarity analysis between synthetic data and real data in each reservoir. Asterisks (*) denote that trends are statistically significant (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05)\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"8\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c8\" colnum=\"8\"\u003e\u003c/div\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eReservoir\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eVariables\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colspan=\"2\" nameend=\"c4\" namest=\"c3\"\u003e\u003cp\u003eKolmogorov-Smirnov\u003c/p\u003e\u003cp\u003eTest\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colspan=\"2\" nameend=\"c6\" namest=\"c5\"\u003e\u003cp\u003eAnderson-Darling\u003c/p\u003e\u003cp\u003eTest\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colspan=\"2\" nameend=\"c8\" namest=\"c7\"\u003e\u003cp\u003eEnergy Distance\u003c/p\u003e\u003cp\u003eTest\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eStatistic\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003ep-value\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003eStatistic\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003ep-value\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003eStatistic\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003ep-value\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"5\" rowspan=\"6\"\u003e\u003cp\u003eSacheon\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003ePrcp\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.130\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.380\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e0.037\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTmax\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.111\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.575\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-0.108\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTmin\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.099\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.716\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-0.695\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eRhum\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.079\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.907\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-0.880\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTOC\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.225\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.016\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e2.099\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.028\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eOverall variables\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e0.157\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e0.948\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"5\" rowspan=\"6\"\u003e\u003cp\u003eDocheon\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003ePrcp\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.104\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.653\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-0.575\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTmax\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.104\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.653\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e0.491\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTmin\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.083\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.878\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eRhum\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.125\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.422\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eTOC\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e0.100\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e0.702\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e0.541\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e0.250\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eOverall variables\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c4\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c5\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c6\"\u003e\u003cp\u003e-\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c7\"\u003e\u003cp\u003e0.257\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c8\"\u003e\u003cp\u003e0.597\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec13\" class=\"Section2\"\u003e\u003ch2\u003e3.2 Comparative Performance of Machine Learning Models\u003c/h2\u003e\u003cp\u003eIn Sacheon Reservoir (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea), the M5-SGB model outperformed all other algorithms, achieving the highest performance metrics (NSE\u0026thinsp;=\u0026thinsp;0.862, R\u0026sup2; = 0.864). Following M5-SGB, the standalone M5 model also demonstrated strong performance. The evaluation revealed a nuanced impact of WGAN-based data augmentation, whereas boosting-based models such as XGB and GB showed improved performance with the augmented data, other models like RF, LL, RR, and SVR exhibited a marginal decrease.\u003c/p\u003e\u003cp\u003eFor Docheon Reservoir (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb), M5-SGB again proved to be the most accurate predictor, recording the highest NSE of 0.804 and an R\u0026sup2; of 0.895. In this reservoir, most models, including the regularized linear models (LL and RR), showed slight performance enhancements with the augmented dataset.\u003c/p\u003e\u003cp\u003eIn summary, M5-SGB was consistently the superior model for TOC prediction in both reservoirs. However, the effectiveness of WGAN-based data augmentation was observed to be model-dependent.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec14\" class=\"Section2\"\u003e\u003ch2\u003e3.3 Projection of TOC in 12 GCMs under Climate Change Scenarios\u003c/h2\u003e\u003cp\u003eThe projected time series TOC concentrations for Sacheon Reservoir were analyzed using the MK test and Sen's slope (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea and Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Under the high emissions SSP5-8.5 scenarios, results indicated a statistically significant increasing trend in TOC levels for all 12 GCMs. Under the intermediate SSP2-4.5 scenarios, 11 of the 12 GCMs also showed a significant increasing trend, with the FGOALS-g3 model being the sole exception. The seasonal analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea) further revealed that projected summer TOC concentrations consistently exceeded those of spring and autumn, emphasizing the potential influence of intensified temperature fluctuations under climate change.\u003c/p\u003e\u003cp\u003eIn Docheon Reservoir (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb and Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e), 8 of the 12 GCMs under the SSP5-8.5 scenarios exhibited a statistically significant increasing TOC trend. In contrast, under SSP2-4.5 scenario, 9 of the 12 models displayed no significant trend, with only MRI-ESM2-0, NESM3 and NorESM2-LM showing a significant upward trajectory. The seasonal analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb) illustrated distinct periodic fluctuations, with particularly pronounced summer peaks.\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eResults of the MK test and Sen's slope analysis applied to the projected TOC time series (2025\u0026ndash;2100) for Sacheon Reservoir. The analysis covers 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios. Asterisks (*) denote a statistically significant trend (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05)\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"5\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eModel\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eScenario\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eTrend\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003eZ-Score\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003eSen's Slope\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eACCESS-CM2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.938\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.996\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.005\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eACCESS-ESM1-5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.243\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.005\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e7.557\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.010\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eCMCC-ESM2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e6.472\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.006\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e7.862\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.010\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eCanESM5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e4.274\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.005\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e9.288\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.013\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eFGOALS-g3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.489\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.512\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eKIOST-ESM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.381\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.763\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eMIROC6\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.435\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e6.292\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.007\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eMRI-ESM2-0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.651\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e4.041\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNESM3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.108\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.009\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNorESM2-LM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.731\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e6.552\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.007\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNorESM2-MM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.826\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e4.050\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eTaiESM1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.216\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e6.813\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.007\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003e\u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e\u003ccaption language=\"En\"\u003e\u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e\u003cdiv class=\"CaptionContent\"\u003e\u003cp\u003eResults of the MK test and Sen's slope analysis applied to the projected TOC time series (2025\u0026ndash;2100) for Docheon Reservoir. The analysis covers 12 GCMs under SSP2-4.5 and SSP5-8.5 scenarios. Asterisks (*) denote a statistically significant trend (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05)\u003c/p\u003e\u003c/div\u003e\u003c/caption\u003e\u003ccolgroup cols=\"5\"\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e\u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e\u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e\u003cthead\u003e\u003ctr\u003e\u003cth align=\"left\" colname=\"c1\"\u003e\u003cp\u003eModel\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c2\"\u003e\u003cp\u003eScenario\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c3\"\u003e\u003cp\u003eTrend\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c4\"\u003e\u003cp\u003eZ-Score\u003c/p\u003e\u003c/th\u003e\u003cth align=\"left\" colname=\"c5\"\u003e\u003cp\u003eSen's Slope\u003c/p\u003e\u003c/th\u003e\u003c/tr\u003e\u003c/thead\u003e\u003ctbody\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eACCESS-CM2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.525\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-0.049\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eACCESS-ESM1-5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-0.112\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.135\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.001\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eCMCC-ESM2\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.377\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.001\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.099\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.005\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eCanESM5\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.673\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e6.283\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.007\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eFGOALS-g3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.229\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.557\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.006\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eKIOST-ESM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.709\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.003\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e4.005\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.006\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eMIROC6\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.180\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.001\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e5.521\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.005\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eMRI-ESM2-0\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.978\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.238\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNESM3\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.252\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e7.952\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.015\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNorESM2-LM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e2.256\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eIncreasing*\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e3.539\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.004\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eNorESM2-MM\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e1.404\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.001\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e0.740\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.001\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c1\" morerows=\"1\" rowspan=\"2\"\u003e\u003cp\u003eTaiESM1\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP245\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-0.184\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e0.000\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003ctr\u003e\u003ctd align=\"left\" colname=\"c2\"\u003e\u003cp\u003eSSP585\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"left\" colname=\"c3\"\u003e\u003cp\u003eNo trend\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e\u003cp\u003e-1.601\u003c/p\u003e\u003c/td\u003e\u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e\u003cp\u003e-0.002\u003c/p\u003e\u003c/td\u003e\u003c/tr\u003e\u003c/tbody\u003e\u003c/colgroup\u003e\u003c/table\u003e\u003c/div\u003e\u003c/p\u003e\u003cp\u003e\u003c/p\u003e\u003c/div\u003e"},{"header":"4 Discussion","content":"\u003cdiv id=\"Sec16\" class=\"Section2\"\u003e\u003ch2\u003e4.1 Augmentation Effects on Model Performance\u003c/h2\u003e\u003cp\u003eThe superior performance of the M5-SGB model across both reservoirs highlights the advantages of its hybrid architecture for predicting complex environmental variables like TOC. By integrating the piecewise linear functions of the M5 model tree with the iterative error-correction of Stochastic Gradient Boosting, M5-SGB effectively captures the non-linear relationships between climatic drivers and water quality. Unlike conventional tree-based models that are limited to their training range, the model's ability to perform localized linear regression at its leaf nodes likely provided a critical advantage in extrapolating future TOC concentrations under novel climate conditions.\u003c/p\u003e\u003cp\u003eThe successful application of WGAN for data augmentation aligns with findings from various domains. For instance, WGAN has been shown to generate synthetic data that closely matches real distributions in medical time series (Esteban et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) and electronic health records (Yoon et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). More specifically, within the environmental sciences, studies have reported that WGAN-based augmentation enhanced the predictive performance of models for aquatic ecosystem health (Lee et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) and improved streamflow prediction accuracy (L\u0026oacute;pez-Chac\u0026oacute;n et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Our finding that WGAN provided a distinct performance boost to boosting algorithms like XGB and GB is consistent with this established potential.\u003c/p\u003e\u003cp\u003eHowever, our results also introduce a critical nuance. The observation that this positive impact did not extend to inherently robust models like RF or regularized linear models suggests that data augmentation is not a universally beneficial strategy. Its utility is highly dependent on the chosen model architecture. This implies that a tailored approach is crucial, as synthetic data may offer significant benefits for data-sensitive algorithms while providing negligible or even slightly negative effects for models already well-suited to smaller datasets. This distinction is also pragmatically important. Implementing WGAN is computationally intensive and requires specialized expertise; therefore, understanding precisely which models benefit most ensures that this powerful technique is applied efficiently. Consequently, studies like the present one, which investigate the specific interactions between data augmentation methods and model architectures, are essential for developing effective and resource-conscious modeling strategies.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec17\" class=\"Section2\"\u003e\u003ch2\u003e4.2 Contrasting Future TOC Projections between Reservoirs under Climate Change\u003c/h2\u003e\u003cp\u003eA central finding of this study is the divergent response of the two reservoirs to climate change; Sacheon Reservoir is projected to experience significant TOC increases under both moderate and high-emissions scenarios, whereas Docheon Reservoir's increasing trend is only significant under the high-emissions scenario.\u003c/p\u003e\u003cp\u003eSacheon Reservoir presents a pristine but unstable paradox. Despite its high-water quality (Grade A), it exhibits a high variability (CV\u0026thinsp;=\u0026thinsp;34.59%), indicating considerable instability and sensitivity to irregular meteorological events or pollution inputs. This inherent vulnerability makes the reservoir highly susceptible to the systematic pressures of climate change, such as accelerated warming and intensified precipitation, which drive increased organic matter input from the watershed (Monteith et al. \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2007\u003c/span\u003e; Fenner et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Consequently, even the moderate climatic shift under SSP2-4.5 appears sufficient to trigger a significant rising TOC trend.\u003c/p\u003e\u003cp\u003eConversely, Docheon Reservoir can be characterized as polluted but stable, with poor baseline water quality but a much lower variability (CV\u0026thinsp;=\u0026thinsp;21.72%). This suggests a more predictable system, possibly buffered by its already high organic load. A compelling hypothesis for its behavior is a threshold-based response. The watershed and reservoir system may effectively buffer TOC mobilization under moderate warming, but this capacity is overwhelmed once a climatic threshold is surpassed under the more extreme high-emissions conditions (Laudon et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). This would explain the lack of a significant trend in the SSP2-4.5 scenario and the abrupt emergence of one under SSP5-8.5.\u003c/p\u003e\u003cp\u003eThese long-term trends are further contextualized by the seasonal analysis, which showed pronounced summer peaks in TOC for both reservoirs. This is consistent with established mechanisms where rising temperatures enhance microbial activity and the decomposition of organic matter, a process expected to intensify under future warming (Sobek et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). These findings underscore that predicting climate change impacts on water quality requires a nuanced approach that considers not only the current state (e.g., water quality grade) but also the system's temporal variability and potential for non-linear responses.\u003c/p\u003e\u003c/div\u003e\u003cdiv id=\"Sec18\" class=\"Section2\"\u003e\u003ch2\u003e4.3 Implications and Future Directions\u003c/h2\u003e\u003cp\u003eThe projected increase in TOC concentrations carries significant practical implications for agricultural water management in South Korea. Rising TOC levels can increase the potential for summer algal blooms and elevate the costs and complexity of water treatment, particularly concerning the formation of disinfection byproducts. These findings present a direct challenge to water management bodies like the Korea Rural Community Corporation (KRC) and underscore the need for proactive, climate-adaptive watershed management strategies. Furthermore, this study's finding that the pristine Sacheon Reservoir is more unstable (high CV) than the polluted Docheon Reservoir suggests that current monitoring strategies may need re-evaluation. A water body's vulnerability to climate change may be better characterized by its dynamic variability than its current quality grade alone. Therefore, designing robust monitoring and management plans should consider not just the current state, but also the inherent instability of the system.\u003c/p\u003e\u003cp\u003eWhile this study provides a robust framework for future TOC prediction, the analysis was confined to two reservoirs, and the generalization of these findings requires further research across a wider range of hydro-climatic conditions. Furthermore, all climate projections are subject to inherent uncertainty from the GCMs themselves, although this was mitigated by using a 12-model ensemble. Therefore, future research should aim to expand this modeling framework to a larger, more diverse set of reservoirs. Comparing WGAN with other generative models, such as Variational Autoencoders (VAEs), could also yield valuable insights into the most effective data augmentation techniques. Finally, incorporating additional predictor variables, such as satellite-derived land-use data or point-source pollution information, could further enhance the accuracy and explanatory power of the predictive models.\u003c/p\u003e\u003c/div\u003e"},{"header":"5 Conclusion","content":"\u003cp\u003eThis study successfully developed and validated a predictive framework combining WGAN-based data augmentation with a high performance M5-SGB model to project future TOC concentrations in data-scarce agricultural reservoirs. The ability to achieve satisfactory prediction accuracy using only readily available meteorological data is particularly encouraging. High-frequency water quality monitoring in reservoirs is often constrained by the logistical challenges of manual, boat-based sampling, making this modeling approach a valuable tool for supplementing sparse observational records. The framework proved effective, accurately capturing the complex relationships between climatic drivers and water quality.\u003c/p\u003e\u003cp\u003eThe projections reveal a significant increasing trend in TOC under future climate change, a trend that is particularly pronounced in the pristine but temporally unstable reservoir. This highlights the vulnerability of even healthy aquatic ecosystems to climatic shifts. Methodologically, this research demonstrated that the effectiveness of data augmentation is highly model-dependent, providing a critical insight for the practical application of machine learning in hydrology.\u003c/p\u003e\u003cp\u003eIn conclusion, this study provides a robust tool for water resource managers to anticipate future challenges. It underscores that effective climate adaptation strategies must consider not only the current state of water bodies but also their dynamic variability to build true resilience.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledg\u003c/strong\u003e\u003cstrong\u003ee\u003c/strong\u003e\u003cstrong\u003ements\u003c/strong\u003eThis work was supported by Korea Institute of Planning and Evaluation for Technology in Food, Agriculture and Forestry (IPET) through Intelligent Agricultural Infra Management for Climate Change Development Program, funded by Ministry of Agriculture, Food and Rural Affairs (MAFRA) (RS-2025-02263904).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eSupplementary Information\u0026nbsp;\u003c/strong\u003eSupplementary data to this article is submitted at the same time.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eA\u003c/strong\u003e\u003cstrong\u003euthor Contributions\u0026nbsp;\u003c/strong\u003eConceptualization: Christian Joseph Siose, Jeongho Han, Kyung-Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Methodology: Jeongho Han, Kyung Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Formal analysis and investigation: Christian Joseph Siose, Jeongho Han; Writing\u0026mdash;original draft preparation: Christian Joseph, Jeongho Han; Writing\u0026mdash;review and editing: Kyung Sook Choi, Byoung Han Choi, Kyoung Jae Lim; Funding acquisition: Kyoung Jae Lim; Supervision: Kyoung Jae Lim; Project administration: Kyoung Jae Lim.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData Availability\u0026nbsp;\u003c/strong\u003eThe data that support the findings of this study are available from the corresponding author ([email protected]) upon reasonable request.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting Interests\u0026nbsp;\u003c/strong\u003eThe authors declare there is no conflict.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003eAguilar-Torrej\u0026oacute;n JA, Balderas-Hern\u0026aacute;ndez P, Roa-Morales G, et al (2023) Relationship, importance, and development of analytical techniques: COD, BOD, and, TOC in water\u0026mdash;An overview through time. SN Appl Sci 5\u003c/li\u003e\n \u003cli\u003eAkiba T, Sano S, Yanase T, et al (2019) Optuna: A Next-generation hyperparameter optimization framework. In: Proceedings of the ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. Association for Computing Machinery, pp 2623\u0026ndash;2631\u003c/li\u003e\n \u003cli\u003eAlqahtani A, Shah MI, Aldrees A, Javed MF (2022) Comparative assessment of individual and ensemble machine learning models for efficient analysis of river water quality. Sustainability (Switzerland) 14: https://doi.org/10.3390/su14031183\u003c/li\u003e\n \u003cli\u003eAnderson TW, Darling DA (1952) asymptotic theory of certain \u0026ldquo;Goodness of Fit\u0026rdquo; criteria based on stochastic processes\u003c/li\u003e\n \u003cli\u003eArjovsky M, Chintala S, Bottou L (2017) Wasserstein Generative Adversarial Networks\u003c/li\u003e\n \u003cli\u003eAwad M, Khanna R (2015) Support vector machines for classification. in: efficient learning machines. Apress, Berkeley, CA, pp 39\u0026ndash;66\u003c/li\u003e\n \u003cli\u003eBhatta B, Shrestha S, Shrestha PK, Talchabhadel R (2019) Evaluation and application of a SWAT model to assess the climate change impact on the hydrology of the Himalayan River Basin. Catena (Amst) 181: https://doi.org/10.1016/j.catena.2019.104082\u003c/li\u003e\n \u003cli\u003eBreiman, L., Friedman, J.H., Olshen, R.A., Stone, C.J., 1984. Classification and regression trees. Chapman \u0026amp; Hall/CRC. https://doi.org/10.1201/9781315139470.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eBreiman L (2001) Random forests. Mach Learn 45:5\u0026ndash;32. https://doi.org/https://doi.org/10.1023/A:1010933404324\u003c/li\u003e\n \u003cli\u003eChen T, Guestrin C (2016) XGBoost: A Scalable Tree Boosting System. https://doi.org/10.1145/2939672.2939785\u003c/li\u003e\n \u003cli\u003eDelpla I, Baur\u0026egrave;s E, Jung AV, Thomas O (2011) Impacts of rainfall events on runoff water quality in an agricultural environment in temperate areas. Science of the Total Environment 409:1683\u0026ndash;1688. https://doi.org/10.1016/j.scitotenv.2011.01.033\u003c/li\u003e\n \u003cli\u003eDelpla I, Jung A-V, Baures E, et al (2009) Impacts of climate change on surface water quality in relation to drinking water production. Environ Int 35:1225\u0026ndash;1233. https://doi.org/10.1016/j.envint.2009.07.001\u003c/li\u003e\n \u003cli\u003eEsteban C, Hyland SL, R\u0026auml;tsch G (2017) Real-valued (medical) time series generation with recurrent conditional gans. arXiv preprint arXiv:170602633\u003c/li\u003e\n \u003cli\u003eFenner N, Meadham J, Jones T, et al (2021) Effects of Climate Change on Peatland Reservoirs: A DOC Perspective. Global Biogeochem Cycles 35: https://doi.org/10.1029/2021GB006992\u003c/li\u003e\n \u003cli\u003eFriedman JH (2001) 999 Reitz Lecture Greedy function approximation: A Gradient Boosting Machine 1\u003c/li\u003e\n \u003cli\u003eFriedman JH (2002) Stochastic gradient boosting. Comput Stat Data Anal 38:367\u0026ndash;378. https://doi.org/10.1016/S0167-9473(01)00065-2\u003c/li\u003e\n \u003cli\u003eGoodfellow IJ, Pouget-Abadie J, Mirza M, et al (2014) Generative Adversarial Nets\u003c/li\u003e\n \u003cli\u003eGulrajani I, Ahmed F, Arjovsky M, et al (2017) Improved training of wasserstein gans. Adv Neural Inf Process Syst 30:\u003c/li\u003e\n \u003cli\u003eHa DW, Jung KY, Baek J, et al (2022) Trend Analysis Using Long-Term Monitoring Data of Water Quality at Churyeongcheon and Yocheon Basins. Sustainability (Switzerland) 14:. https://doi.org/10.3390/su14159770\u003c/li\u003e\n \u003cli\u003eHirsch RM, Slack JR, Smith RA (1982) Techniques of trend analysis for monthly water quality data. Water Resour Res 18:107\u0026ndash;121\u003c/li\u003e\n \u003cli\u003eHoerl AE, Kennard RW (1970) Ridge regression: Biased estimation for nonorthogonal problems. Technometrics 12:55\u0026ndash;67\u003c/li\u003e\n \u003cli\u003eJeppesen E, Meerhoff M, Davidson TA, et al (2014) Climate change impacts on lakes: An integrated ecological perspective based on a multi-faceted approach, with special focus on shallow lakes. J Limnol 73:88\u0026ndash;111. https://doi.org/10.4081/jlimnol.2014.844\u003c/li\u003e\n \u003cli\u003eKolmogorov A. (1933) Sulla determinazione empirica di una legge didistribuzione. Giorn Dell\u0026rsquo;inst Ital Degli Att 4:89\u0026ndash;91\u003c/li\u003e\n \u003cli\u003eLaudon H, Berggren M, \u0026Aring;gren A, et al (2011) Patterns and dynamics of dissolved organic carbon (DOC) in boreal streams: The role of processes, connectivity, and scaling. Ecosystems 14:880\u0026ndash;893\u003c/li\u003e\n \u003cli\u003eLee S, Kim J, Lee G, et al (2021) Prediction of aquatic ecosystem health indices through machine learning models using the wgan-based data augmentation method. Sustainability (Switzerland) 13: https://doi.org/10.3390/su131810435\u003c/li\u003e\n \u003cli\u003eLee Y-J, Hsieh W-F, Huang C-M (2005) /spl epsi/-SSVR: a smooth support vector machine for /spl epsi/-insensitive regression. IEEE Trans Knowl Data Eng 17:678\u0026ndash;685. https://doi.org/10.1109/TKDE.2005.77\u003c/li\u003e\n \u003cli\u003eL\u0026oacute;pez-Chac\u0026oacute;n SR, Salazar F, Blad\u0026eacute; E (2023) Combining Synthetic and Observed Data to Enhance Machine Learning Model Performance for Streamflow Prediction. Water (Switzerland) 15: https://doi.org/10.3390/w15112020\u003c/li\u003e\n \u003cli\u003eMingers, J., 1989. An empirical comparison of pruning methods for decision tree induction. Mach. Learn. 4, 227\u0026ndash;243. https://doi.org/10.1023/A:1022604100933.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eModi P, Chintalacheruvu MR (2024) Investigating River water quality assessment through non-parametric analysis: A case study of the Godavari River in India. Environmental Quality Management 33:239\u0026ndash;264. https://doi.org/10.1002/tqem.22117\u003c/li\u003e\n \u003cli\u003eMonteith DT, Stoddard JL, Evans CD, et al (2007) Dissolved organic carbon trends resulting from changes in atmospheric deposition chemistry. Nature 450:537\u0026ndash;540. https://doi.org/10.1038/nature06316\u003c/li\u003e\n \u003cli\u003eMoriasi DN, Gitau MW, Pai N, Daggupati P (2015) Hydrologic and water quality models: Performance measures and evaluation criteria. Trans ASABE 58:1763\u0026ndash;1785. https://doi.org/10.13031/trans.58.10715\u003c/li\u003e\n \u003cli\u003eNallakaruppan MK, Gangadevi E, Shri ML, et al (2024) Reliable water quality prediction and parametric analysis using explainable AI models. Sci Rep 14:7520. https://doi.org/10.1038/s41598-024-56775-y\u003c/li\u003e\n \u003cli\u003eNash JE, Sutcliffe J V (1970) River flow forecasting through conceptual models\u0026rsquo; part i-a discussion of principles*. \u0026copy; North-Holland Publishing Co\u003c/li\u003e\n \u003cli\u003eNatekin A, Knoll A (2013) Gradient boosting machines, a tutorial. Front Neurorobot 7: https://doi.org/10.3389/fnbot.2013.00021\u003c/li\u003e\n \u003cli\u003eOh H, Park HY, Kim JI, et al (2024) Enhancing machine learning models for total organic carbon prediction by integrating geospatial parameters in river watersheds. Science of the Total Environment 943: https://doi.org/10.1016/j.scitotenv.2024.173743\u003c/li\u003e\n \u003cli\u003ePeng L, Li S, Sun H, Huang S (2022) A Pipe Ultrasonic Guided Wave Signal Generation Network Suitable for Data Enhancement in Deep Learning: US-WGAN. Energies (Basel) 15: https://doi.org/10.3390/en15186695\u003c/li\u003e\n \u003cli\u003eProbst P, Wright MN, Boulesteix AL (2019) Hyperparameters and tuning strategies for random forest. Wiley Interdiscip Rev Data Min Knowl Discov 9\u003c/li\u003e\n \u003cli\u003eQuinlan Basser JR (1992) Learning with continuous classes. World Scientiic\u003c/li\u003e\n \u003cli\u003eRizzo ML, Sz\u0026eacute;kely GJ (2016) Energy distance. Wiley Interdiscip Rev Comput Stat 8:27\u0026ndash;38\u003c/li\u003e\n \u003cli\u003eRural Agricultural Water Resource Information System (RAWRIS) (2024) Rural Agricultural Water Resource Information System. Available at: https://rawris.ekr.or.kr Accessed 11 Nov 2024.\u003c/li\u003e\n \u003cli\u003eSattari MT, Pal M, Mirabbasi R, Abraham J (2018) Ensemble of M5 Model Tree-Based Modelling of Sodium Adsorption Ratio\u003c/li\u003e\n \u003cli\u003eSen PK (1968) Estimates of the regression coefficient based on Kendall\u0026rsquo;s tau. J Am Stat Assoc 63:1379\u0026ndash;1389\u003c/li\u003e\n \u003cli\u003eSipper M (2022) High Per Parameter: A Large-Scale Study of Hyperparameter Tuning for Machine Learning Algorithms. Algorithms 15:315. https://doi.org/10.3390/a15090315\u003c/li\u003e\n \u003cli\u003eSobek S, Tranvik LJ, Prairie YT, et al (2007) Patterns and regulation of dissolved organic carbon: An analysis of 7,500 widely distributed lakes. Limnol Oceanogr 52:1208\u0026ndash;1219. https://doi.org/10.4319/lo.2007.52.3.1208\u003c/li\u003e\n \u003cli\u003eTebaldi C, Debeire K, Eyring V, et al (2021) Climate model projections from the Scenario Model Intercomparison Project (ScenarioMIP) of CMIP6. Earth System Dynamics 12:253\u0026ndash;293. https://doi.org/10.5194/esd-12-253-2021\u003c/li\u003e\n \u003cli\u003eTibshirani R (1996) Regression Shrinkage and Selection via the Lasso\u003c/li\u003e\n \u003cli\u003eTranvik LJ, Downing JA, Cotner JB, et al (2009) Lakes and reservoirs as regulators of carbon cycling and climate. Limnol Oceanogr 54:2298\u0026ndash;2314. https://doi.org/10.4319/lo.2009.54.6_part_2.2298\u003c/li\u003e\n \u003cli\u003eWhitehead PG, Wilby RL, Battarbee RW, et al (2009) A review of the potential impacts of climate change on surface water quality. Hydrological Sciences Journal 54:101\u0026ndash;123. https://doi.org/10.1623/hysj.54.1.101\u003c/li\u003e\n \u003cli\u003eYan X, Zhang T, Du W, et al (2024) A comprehensive review of machine learning for water quality prediction over the past five years. J Mar Sci Eng 12\u003c/li\u003e\n \u003cli\u003eYoon J, Jarrett D, Van der Schaar M (2019) Time-series generative adversarial networks. Adv Neural Inf Process Syst 32:\u003c/li\u003e\n \u003cli\u003eZheng Y, Zhang X, Zhou Y, et al (2025) Deep representation learning enables cross-basin water quality prediction under data-scarce conditions. NPJ Clean Water 8: https://doi.org/10.1038/s41545-025-00466-2\u003c/li\u003e\n \u003cli\u003eZhou S, Long H, Chen W, et al (2025) Temperature seasonality regulates organic carbon burial in lake. Nature Communications 16: https://doi.org/10.1038/s41467-025-56399-4\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"water-resources-management","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"warm","sideBox":"Learn more about [Water Resources Management](https://www.springer.com/journal/11269)","snPcode":"11269","submissionUrl":"https://submission.nature.com/new-submission/11269/3","title":"Water Resources Management","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Reservoir water quality, Total organic carbon, Climate change, Machine learning model, WGAN, Data augmentation","lastPublishedDoi":"10.21203/rs.3.rs-8194601/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8194601/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eClimate change poses a significant threat to water quality in agricultural reservoirs. However, projecting future changes using machine learning (ML) is often hampered by the limited availability of long-term observational data. This study develops a robust framework to project future Total Organic Carbon (TOC) concentrations by combining a Wasserstein Generative Adversarial Network (WGAN) for data augmentation with high-performance ML models. We applied this framework to two contrasting South Korean reservoirs (one pristine, one polluted). Historical data (n\u0026thinsp;=\u0026thinsp;60) for each was augmented using WGAN, and a suite of nine ML models was optimized and evaluated. The best-performing model, M5 Model Tree with Stochastic Gradient Boosting (M5-SGB), was then used to project TOC from 2025\u0026ndash;2100 using a 12 Global Climate Model (GCM) ensemble from CMIP6 under SSP2-4.5 and SSP5-8.5 scenarios. Notably, the effectiveness of WGAN augmentation was model-dependent, significantly boosting the performance of certain algorithms while having a negligible effect on others. Future projections revealed a statistically significant increasing TOC trend for both reservoirs. However, the response was site-specific: the pristine but temporally unstable reservoir showed vulnerability even under the moderate-emissions scenario, while the polluted but stable reservoir only exhibited a significant trend under the high-emissions scenario. This study provides a practical framework for water quality prediction in data-limited contexts. The findings demonstrate that a reservoir's vulnerability to climate change is critically linked to its dynamic variability, not just its static water quality grade. This offers crucial insights for designing more effective, risk-based water management and monitoring strategies.\u003c/p\u003e","manuscriptTitle":"Projecting Future TOC in Data-Scarce Agricultural Reservoirs: A WGAN-Enhanced Framework Revealing the Importance of Dynamic Variability","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-12-01 06:29:17","doi":"10.21203/rs.3.rs-8194601/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revisions","date":"2025-12-03T05:39:00+00:00","index":"","fulltext":""},{"type":"reviewerAgreed","content":"","date":"2025-11-29T03:37:01+00:00","index":0,"fulltext":""},{"type":"reviewersInvited","content":"","date":"2025-11-25T13:31:35+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2025-11-25T05:14:49+00:00","index":"","fulltext":""},{"type":"submitted","content":"Water Resources Management","date":"2025-11-24T10:19:57+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"water-resources-management","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"warm","sideBox":"Learn more about [Water Resources Management](https://www.springer.com/journal/11269)","snPcode":"11269","submissionUrl":"https://submission.nature.com/new-submission/11269/3","title":"Water Resources Management","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"5da01454-cc63-4b4d-b3cb-1470f6c37744","owner":[],"postedDate":"December 1st, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2026-02-23T16:14:43+00:00","versionOfRecord":{"articleIdentity":"rs-8194601","link":"https://doi.org/10.1007/s11269-026-04503-8","journal":{"identity":"water-resources-management","isVorOnly":false,"title":"Water Resources Management"},"publishedOn":"2026-02-21 15:59:22","publishedOnDateReadable":"February 21st, 2026"},"versionCreatedAt":"2025-12-01 06:29:17","video":"","vorDoi":"10.1007/s11269-026-04503-8","vorDoiUrl":"https://doi.org/10.1007/s11269-026-04503-8","workflowStages":[]},"version":"v1","identity":"rs-8194601","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8194601","identity":"rs-8194601","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}