Beach groundwater and nearshore hydrodynamic data were collected during a field experiment along two dissipative beach transects on Galveston Island, Texas, in the fall of 2023. The monitored beaches serve as nesting habitat for the critically endangered Kemp’s ridley sea turtle. Conditions ranged from calm to stormy, with two storms occurring during the experiment, inundating the entire beach up to the dune toe. Collected hydrodynamic data include readings from pressure loggers submerged in the foreshore and mounted in groundwater wells in the backshore, data from two wave buoys about 1.5 km offshore, and GoPro timestacks of the instantaneous waterline (wave runup). Other collected data include bathymetry and topography surveys, subsurface temperature and moisture content readings, and sediment characteristics. This comprehensive dataset can be used to (1) study relevant beach inundation and groundwater processes, including their effect on the local ecosystem (e.g., repeated flooding of sea turtle nests), (2) study the propagation of nearshore hydrodynamic processes into the beach matrix and groundwater table, and (3) validate existing beach groundwater models.@en

Tropical-cyclone impacts can have devastating effects on the population, infrastructure, and natural habitats. However, predicting these impacts is difficult due to the inherent uncertainties in the storm track and intensity. In addition, due to computational constraints, both the relevant ocean physics and the uncertainties in meteorological forcing are only partly accounted for. This paper presents a new method, called the Tropical Cyclone Forecasting Framework (TC-FF), to probabilistically forecast compound flooding induced by tropical cyclones, considering uncertainties in track, forward speed, and wind speed and/or intensity. The open-source method accounts for all major relevant physical drivers, including tide, surge, and rainfall, and considers TC uncertainties through Gaussian error distributions and autoregressive techniques. The tool creates temporally and spatially varying wind fields to force a computationally efficient compound-flood model, allowing for the computation of probabilistic wind and flood hazard maps for any oceanic basin in the world as it does not require detailed information on the distribution of historical errors. A comparison of TC-FF and JTWC operational ensembles, both based on DeMaria et al. (2009), revealed minor differences of <10 %, suggesting that TC-FF can be employed as an alternative, for example, in data-scarce environments. The method was applied to Cyclone Idai in Mozambique. The underlying physical model showed reliable skill in terms of tidal propagation, reproducing the storm surge generation during landfall and flooding near the city of Beira (success index of 0.59). The method was successfully applied to forecasting the impact of Idai with different lead times. The case study analyzed needed at least 200 ensemble members to get reliable water levels and flood results 3 d before landfall (<1 % flood probability error and <20 cm sampling errors). Results showed the sensitivity of forecasting, especially with increasing lead times, highlighting the importance of accounting for cyclone variability in decision-making and risk management.@en

Low-lying, tropical, coral-reef-lined coastlines are becoming increasingly vulnerable to wave-driven flooding due to population growth, coral reef degradation, and sea-level rise. Early-warning systems (EWSs) are needed to enable coastal authorities to issue timely alerts and coordinate preparedness and evacuation measures for their coastal communities. At longer timescales, risk management and adaptation planning require robust assessments of future flooding hazard considering uncertainties. However, due to diversity in reef morphologies and complex reef hydrodynamics compared to sandy shorelines, there have been no robust analytical solutions for wave runup to allow for the development of large-scale coastal wave-driven flooding EWSs and risk assessment frameworks for reef-lined coasts. To address the need for fast, robust predictions of runup that account for the natural variability in coral reef morphologies, we constructed the BEWARE-2 (Broad-range Estimator of Wave Attack in Reef Environments) meta-process modeling system. We developed this meta-process model using a training dataset of hydrodynamics and wave runup computed by the XBeach Non-Hydrostatic process-based hydrodynamic model for 440 combinations of water level, wave height, and wave period with 195 representative reef profiles that encompass the natural diversity in real-world fringing coral reef systems. Through this innovation, BEWARE-2 can be applied in a larger range of coastal settings than meta-models that rely on a parametric description of the coral reef geometry. In the validation stage, the BEWARE-2 modeling system produced runup results that had a relative root mean square error of 13 % and relative bias of 5 % relative to runup simulated by XBeach Non-Hydrostatic for a large range of oceanographic forcing conditions and for diverse reef morphologies (root mean square error and bias 0.63 and 0.26 m, respectively, relative to mean simulated wave runup of 4.85 m). Incorporating parametric modifications in the modeling system to account for variations in reef roughness and beach slope allows for systematic errors (relative bias) in BEWARE-2 predictions to be reduced by a factor of 1.5–6.5 for relatively coarse or smooth reefs and mild or steep beach slopes. This prediction provided by the BEWARE-2 modeling system is faster by 4–5 orders of magnitude than the full, process-based hydrodynamic model and could therefore be integrated into large-scale EWSs for tropical, reef-lined coasts and used for large-scale flood risk assessments.@en

Sandy beach-dune systems make up a large part of coastal areas world wide. Their function as an eco-system as well as a protective barrier for human and natural habitat is under increased threat due to climate change. A thorough understanding of change processes at the sediment surface is essential to facilitate prediction of future development and management strategies to maintain their function. Especially slow and small scale processes happening over several days up to weeks at cm level, such as aeolian sand transport are difficult to identify and analyse. Permanent laser scanning (PLS) is a useful tool in the study and analysis of coastal processes as it captures a data representation of the evolution of the sediment surface over extended periods of time (up to several years) with high detail (at cm-dm level). The PLS data set considered for this study, consists of hourly acquired 3D point clouds representing the surface evolution of a section of the Dutch coast during three years. However, it is challenging to extract concrete information on specific change processes from the large and complex PLS data set. We use multiple hypothesis testing in order to reduce the PLS data set to a so-called inventory of trends, consisting of 12.8 million partial time series with associated rate of change and elevation. The inventory of trends proofs to be a suitable tool to identify natural processes such as storms and aeolian sand transport in our test area in the aeolian zone of a sandy beach-dune system on the Dutch coast. We identify these processes and provide a tool to derive summarising data from the complex PLS data set. We find that all partial time series identified as most likely representing aeolian sand transport, result in 1354 m³ of sand deposition in our study area over the course of three years. We also show a comparison with transects from JarKus data and find a correlation between anthropogenic activities and erosion in our test area with a correlation coefficient of 0.3.

@en

Climate change and human activity threaten sea turtle nesting beaches through increased flooding and erosion. Understanding the environmental characteristics that enable nesting can aid to preserve and expand these habitats. While numerous local studies exist, a comprehensive global analysis of environmental influences on the distribution of sea turtle nesting habitats remains largely unexplored. Here, we relate the distribution of global sea turtle nesting to 22 coastal indicators, spanning hydrodynamic, atmospheric, geophysical, habitat, and human processes. Using state-of-the-art global datasets and a novel 50-km-resolution hexagonal coastline grid (Coastgons), we employ machine learning to identify spatially homogeneous patterns in the indicators and correlate these to the occurrence of nesting grounds. Our findings suggest sea surface temperature, tidal range, extreme surges, and proximity to coral and seagrass habitats significantly influence global nesting distribution. Low tidal ranges and low extreme surges appear to be particularly favorable for individual species, likely due to reduced nest flooding. Other indicators, previously reported as influential (e.g., precipitation and wind speed), were not as important in our global-scale analysis. Finally, we identify new, potentially suitable nesting regions for each species. On average, 23 % of global coastal regions between - 39 ^∘ and 48 ^∘ latitude could be suitable for nesting, while only 7 % is currently used by turtles, showing that the realized niche is significantly smaller than the fundamental niche, and that there is potential for sea turtles to expand their nesting habitat. Our results help identify suitable nesting conditions, quantify potential hazards to global nesting habitats, and lay a foundation for nature-based solutions to preserve and potentially expand these habitats.

@en

The bimodal wave climate of the semi-protected shallow Gulf St Vincent in South Australia has been analyzed through a forty-year (1980-2020) wave hindcast and an investigation into the climatic drivers of wave climate anomalies is presented. The sea and swell partitions of the wave climate were modelled independently as well as using an integrated model with both partitions represented. The wave hindcast was validated against two wave buoys located off the coast of Adelaide’s metropolitan beaches and key wave parameter anomalies were calculated across the gulf. Teleconnections were investigated, and the Southern Annular Mode is found to have the strongest correlations to wave parameter anomalies while the Southern Oscillation Index and the Dipole Mode Index fluctuations are found to correlate seasonally with wave parameter anomalies. Projected future trends of these climate drivers from literature have been related to the teleconnections found in this study to inform future trends of bimodal wave conditions in the gulf. The Southern Annular Mode is projected to trend positive which will reduce wave height and the westerly component of waves in the gulf, while the Southern Oscillation Index is projected to become more variable in the future which will lead to more extreme winter and spring wave conditions. An understanding of these trends allows coastal managers to pre-emptively manage the impacts of waves on the coastline at a seasonal to annual basis and provides insight into future wave conditions beyond these time periods.@en

All around the world, deltas are among the most densely populated and heavily utilised regions, where crucial functions, such as freshwater availability and safety against flooding, strongly relate to the natural dynamics of the system. Therefore, a thorough understanding of the estuarine system is crucial, especially when developing nature-based solutions for safeguarding these essential functions for today’s society as well as future generations. To better understand the effect of different estuarine parameters on salt intrusion, an extensive sensitivity analysis has been executed based on an idealised estuary layout. The idealised estuary is parametrically designed using thirteen parameters that represent both boundary conditions and geometric features, such as river discharge and water depth. Subsequently, the Delft3D Flexible Mesh (DFM) model has been used to determine the salt intrusion, allowing the exploration of a wide range of estuary layouts. @en

Estuaries are among the most densely populated and heavily utilised regions in the world, where crucial functions – e.g., freshwater availability and water safety – strongly relate to the natural dynamics of the system. When developing nature-based solutions to safeguard these essential functions, a thorough understanding of estuarine dynamics is required. This study describes an elaborate sensitivity analysis on the salt intrusion length using an idealised estuary, which is parametrically designed using key estuary-scale parameters – e.g., river discharge and tidal flats – to cover a wide range of estuary classes. We were able to systematically investigate such a wide range of estuary classes due to the combination of (1) state-of-the-art hydrodynamic modelling software, (2) high performance computing, and (3) reduction and analysis techniques using machine learning. The results show that the extent of the estuarine salt intrusion length is largely determined by four estuarine features: (1) river discharge; (2) cross-sectional area (especially water depth); (3) tidal damping/amplification; and (4) tidal asymmetry. In general, the salt intrusion length shows clear correlations with (a combination of) estuary-scale parameters, which all put an upper limit on the salt intrusion length. These relations provide crucial insights for successful development of nature-based solutions to mitigate salt intrusion in estuarine environments.

@en

Changing (wind) climate might influence the magnitude, direction, and frequency of wave systems (Lobeto et al., 2021). However, in coastal engineering applications, generalized wave parameters are commonly used in climate change assessments with the risk of, for example, misrepresenting the nearshore transformation of wind-driven wave climates (Hegermiller et al., 2017). In consequence, these uncertainties in the nearshore (wind) climate will affect, amongst others, ship navigation, the implementation of marine renewable energy farms, the feasibility of coastal infrastructure and defences, or the efficiency of sandy coastal maintenance, and thus the decision-making of long-term, multidecadal coastal strategies (Rijksoverheid, 2013), especially when they are designed accounting for the Building with Nature concept (de Vriend et al., 2015). This study analyses the importance and application of considering multiple coexisting wave trains on the Dutch shoreface. @en

In recent years, coastal management has been facing new challenges: socio-economic growth and consequent climate change impose new boundary conditions pushing coastal systems towards unseen states. For adaptation and mitigation strategies as well as risk management, the resilience of systems to these projected changes must be tested and quantified using predictive tools, given the scarcity of observations. Process-based models, which limit the number of assumptions, are the preferred tools. However, these models are computationally expensive and therefore unattractive for global sensitivity and uncertainty analyses. Input and model reduction techniques, as well as behavioural empirical models, have been widely used to overcome these computational difficulties. In this paper, we propose a process-based hybrid workflow—that combines statistical and machine learning with a process-based numerical model—to provide sensitivity analyses on complex systems. As an example we explore salt intrusion in estuaries. The novelty of the method presented is the implementation of an adaptive sampling technique of numerical experiments with a process-based hydrodynamic model, and the training of a neural network to augment the set of numerical runs executed. The first uses predictive uncertainty to automatically explore the response of the complex system to varying environmental boundary conditions and geomorphological configurations. The second is trained to provide system responses around the sampled points. This exploration is closed by simulating the extremes in the output space as found by a genetic algorithm. This scheme is shown to be highly efficient in non-linear, heteroscedastic, and highly non-stationary systems.

@en

In this study four experiments were conducted to investigate uncertainty in future longshore sediment transport (LST) projections due to: working with continuous time series of CSIRO CMIP6-driven waves (experiment #1) or sliced time series of waves from CSIRO-CMIP6-Ws and CSIRO-CMIP5-Ws (experiment #2); different wave-model-parametrization pairs to generate wave projections (experiment #3); and the inclusion/exclusion of sea level rise (SLR) for wave transformation (experiment #4). For each experiment, a weighted ensemble consisting of offshore wave forcing conditions, a surrogate model for nearshore wave transformation and eight LST models was used. The results of experiment # 1 indicated that the annual LST rates obtained from a continuous time series of waves were influenced by climate variability acting on timescales of 20-30 years. Uncertainty decomposition clearly reveals that for near-future coastal planning, a large part of the uncertainty arises from model selection and natural variability of the system (e.g., on average, 4% scenario, 57% model, and 39% internal variability). For the far future, the total uncertainty consists of 25% scenario, 54% model and 21% internal variability. Experiment #2 indicates that CMIP6 driven wave climatology yield similar outcomes to CMIP5 driven wave climatology in that LST rates decrease along the study area’s coast by less than 10%. The results of experiment #3 indicate that intra- and inter-annual variability of LST rates are influenced by the parameterization schemes of the wave simulations. This can increase the range of uncertainty in the LST projections and at the same time can limit the robustness of the projections. The inclusion of SLR (experiment #4) in wave transformation, under SSP1-2.6 and SSP5-8.5 scenarios, yields only meagre changes in the LST projections, compared to the case no SLR. However, it is noted that future research on SLR influence should include potential changes in nearshore profile shapes.@en

In the coming decades, coastal flooding will become more frequent due to sea-level rise and potential changes in storms. To produce global storm surge projections from 1950 to 2050, we force the Global Tide and Surge Model with a ∼25-km resolution climate model ensemble from the Coupled Model Intercomparison Project Phase 6 High Resolution Model Intercomparison Project (HighResMIP). This is the first time that such a high-resolution ensemble is used to assess changes in future storm surges across the globe. We validate the present epoch (1985–2014) against the ERA5 climate reanalysis, which shows a good overall agreement. However, there is a clear spatial bias with generally a positive bias in coastal areas along semi-enclosed seas and negative bias in equatorial regions. Comparing the future epoch (2021–2050) against the historical epoch (1951–1980), we project ensemble-median changes up to 0.1 (or 20%) in the 1 in 10-year storm surge levels. These changes are not uniform across the globe with decreases along the coast of Mediterranean and northern Africa and southern Australia and increases along the south coast of Australia and Alaska. There are also increases along (parts) of the coasts of northern Caribbean, eastern Africa, China and the Korean peninsula, but with less agreement among the HighResMIP ensemble. Information resulting from this study can be used to inform broad-scale assessment of coastal impacts under future climate change.

@en

This study unravels 58-years (1961–2018) of wind-waves in the Mediterranean Sea (MS). A wave dataset was developed using the wave model WAVEWATCH III forced with the high-resolution (5.5 km) UERRA-MESCAN-SURFEX downscaled wind fields which better contain the imprint of the local geomorphology compared to other, coarser wind datasets used in previous studies. Thus, improving the reliability and characterization of the wind-wave climate in the basin. Validation results revealed a higher performance than previous datasets, particularly on the wave direction (θ_m), with a bias<1°. Climate variability at seasonal and interannual scales, wind-seas and swells distribution, and long-term trends in storminess and in the mean and extreme regimes were analysed. Results show a slight swell influence over the wind-sea in the hourly spectra at a large portion of the basin, excluding the wave generation areas. We detected that the western MS is the most storminess region with an average of three storms/year. Moreover, the anomalies of the seasonal mean wave direction relative to θ_m are large (∼60°), with opposing behaviours between the winter and summer. Finally, the long-term trends in the mean and extreme conditions and in storminess are mild with values reaching 6 cm/decade and less than 2% in the absolute value, respectively.

@en

The last years, capacity of vegetation to reduce wave impact is receiving considerable attention. To predict wave attenuation processes within vegetation fields reliable estimates of vegetation parameters are needed. This proves to be difficult for woody vegetation as it consists of complex branch structures, characterized by varying branch densities, diameters and angles. State of the art physical and numerical models effectively use a single value for the diameter, b _v and density, N of vegetation, which is unrepresentative for complex vegetation, such as trees. Trees can be better described by the projected frontal-surface area, A _v. Hence, this work compares methods to quantify the A _v in space for a pollard willow forest, and determines suitability of these methods for predicting wave attenuation using a spectral wave model (SWAN). We use data from manual measurements and Terrestrial Laser Scans (TLS), to estimate the vertical distribution of A _v; and data from large-scale flume experiments performed on a willow forest to verify model sensitivity to A _v inferences. As a baseline for comparison, tree models that describe the structure of the trees in various degrees of complexity are compiled. The most realistic tree model is used to quantify potential errors in TLS and basic manual measurements of N and b _v. An initial comparison shows that the TLS data underestimates A _v, which indicates that conducting manual measurements is more suitable to quantify a homogeneous forest. We found that the TLS suffers from shadowing effects (i.e., blockage of laser beams) and we recommend to apply a correction factor to improve its measurements. Furthermore, we identified the impact that the different methods to determine A _v have on the estimation of wave attenuation using SWAN; in addition we verified the model results with data from large-scale flume experiments performed on the willow forest. The modeled sensitivity tests indicate large differences in wave attenuation and, consequently, a wide range (0.94–1.70) of bulk drag coefficients, (Formula presented.), for the various methods applied. This shows the variation of outcome between measuring methods and highlights the importance of stating the selected method for reliable frontal-surface area estimations, and consequently for reliable wave attenuation predictions.

@en

Sandy beaches and dune systems have high recreational and ecological value, and they offer protection against flooding during storms. At the same time, these systems are very vulnerable to storm impacts. Process-based numerical models are presently used to assess the morphological changes of dune and beach systems during storms. However, such models come with high computational costs, hindering their use in real-life applications which demand many simulations and/or involve a large spatial-temporal domain. Here we design a novel meta-model to predict dune erosion volume (DEV) at the Dutch coast, based on artificial neural networks (ANNs), trained with cases from process-based modeling. First, we reduce an initial database of 1/41400 observed sandy profiles along the Dutch coastline to 100 representative typological coastal profiles (TCPs). Next, we synthesize a set of plausible extreme storm events, which reproduces the probability distributions and statistical dependencies of offshore wave and water level records. We choose 100 of these events to simulate the dune response of the 100 TCPs using the process-based model XBeach, resulting in 10 000 cases. Using these cases as training data, we design a two-phase meta-model, comprised of a classifying ANN (which predicts the occurrence (or not) of erosion) and a regression ANN (which gives a DEV prediction). Validation against a benchmark dataset created with XBeach and a sparse set of available dune erosion observations shows high prediction skill with a skill score of 0.82. The meta-model can predict post-storm DEV 103-104 times faster (depending on the duration of the storm) than running XBeach. Hence, this model may be integrated in early warning systems or allow coastal engineers and managers to upscale storm forcing to dune response investigations to large coastal areas with relative ease.

@en

Worldwide, communities are facing increasing flood risk, due to more frequent and intense hazards and rising exposure through more people living along coastlines and in flood plains. Nature-based Solutions (NbS), such as mangroves, and riparian forests, offer huge potential for adaptation and risk reduction. The capacity of trees and forests to attenuate waves and mitigate storm damages receives massive attention, especially after extreme storm events. However, application of forests in flood mitigation strategies remains limited to date, due to lack of real-scale measurements on the performance under extreme conditions. Experiments executed in a large-scale flume with a willow forest to dissipate waves show that trees are hardly damaged and strongly reduce wave and run-up heights, even when maximum wave heights are up to 2.5 m. It was observed for the first time that the surface area of the tree canopy is most relevant for wave attenuation and that the very flexible leaves limitedly add to effectiveness. Overall, the study shows that forests can play a significant role in reducing wave heights and run-up under extreme conditions. Currently, this potential is hardly used but may offer future benefits in achieving more adaptive levee designs.

@en

Tropical Cyclones (TCs) are singular storms causing intense wind, large waves, extreme water levels, and heavy rainfall. TCs prove every year to be one of the most destructive natural phenomena worldwide. The quantitative assessment of the hazards resulting from TCs (i.e., flooding and extreme winds) is challenging since satellite data are only available for recent decades, whereas older historical observations are incomplete and less accurate. In addition, long-term prediction through numerical weather forecasting is still limited. This often results in large uncertainties in the definition of TC hazards associated with events with longer return periods or in areas infrequently impacted by TCs. Even when this information is available, for example through statistical sampling of synthetic TC tracks, the numerical modelling of the associated hazards for all the different TC conditions can lead to computational costs which are often infeasible. Several methodologies that overcome the issues of accuracy and computational efficiency currently exist, but these are not generically applicable, and they tend to focus on specific areas only, for example where TCs typically make landfall. The main contribution of this paper is a novel methodology for the estimation and analysis of TC hydro-meteorological conditions and induced hazards. The method is generically applicable and maximizes accuracy while accounting for computational efficiency. Our approach identifies a smaller but representative set of TC tracks (RTCs) that preserves the information about extremes and the frequency of events of the larger population. The method is successfully applied and validated in a case study in the Bay of Bengal, using a set of synthetic TC tracks representing 1000 years of TC climate. For the best-performing configuration, the required number of scenarios and associated computational costs were reduced by 90% while maintaining accuracy in the simulated offshore storm surges, significant wave height, and windspeeds typically within 10% of the prediction based on the original full set of scenarios. This method is globally applicable and greatly improves the efficiency of TC-related hazard estimation, making it particularly valuable for areas with limited historical data.

@en

This study quantifies the uncertainties in the projected changes in potential longshore sediment transport (LST) rates along a non-straight coastline. Four main sources of uncertainty, including the choice of emission scenarios, Global Circulation Model-driven offshore wave datasets (GCM-Ws), LST models, and their non-linear interactions were addressed through two ensemble modelling frameworks. The first ensemble consisted of the offshore wave forcing conditions without any bias correction (i.e., wave parameters extracted from eight datasets of GCM-Ws for baseline period 1979–2005, and future period 2081–2100 under two emission scenarios), a hybrid wave transformation method, and eight LST models (i.e., four bulk formulae, four process-based models). The differentiating factor of the second ensemble was the application of bias correction to the GCM-Ws, using a hindcast dataset as the reference. All ensemble members were weighted according to their performance to reproduce the reference LST patterns for the baseline period. Additionally, the total uncertainty of the LST projections was decomposed into the main sources and their interactions using the ANOVA method. Finally, the robustness of the LST projections was checked. Comparison of the projected changes in LST rates obtained from two ensembles indicated that the bias correction could relatively reduce the ranges of the uncertainty in the LST projections. On the annual scale, the contribution of emission scenarios, GCM-Ws, LST models and non-linear interactions to the total uncertainty was about 10–20, 35–50, 5–15, and 30–35%, respectively. Overall, the weighted means of the ensembles reported a decrease in net annual mean LST rates (less than 10% under RCP 4.5, a 10–20% under RCP 8.5). However, no robust projected changes in LST rates on annual and seasonal scales were found, questioning any ultimate decision being made using the means of the projected changes.@en

Wave climatologies from historical and projected simulations of the ACCESS1.0, MIROC5 and CNRM-CM5 Global Circulation Models (GCM) were sourced from the Coordinated Ocean Wave Climate Project (COWCLIP) and downscaled using the SWAN wave model. Biases between GCM's historical simulations and a regional hindcast were assessed, and the two best-performing models (ACCESS1.0, MIROC5) had their projections analysed. A general increase in wave height and period was observed along the south/west, together with a decrease in (Formula presented.) along the north/east coasts. The projected near-term (NEA21C) period shows mostly a (Formula presented.) increase, whilst for the long-term (END21C) period, increased and decreased (Formula presented.) are present. The areas of statistically significant changes are larger in the END21C than in the NEA21C period. The wave direction change is counter-clockwise along the west and clockwise along the east coasts. This study is a first assessment of historical and projected GCM-forced waves along New Zealand and the database we generated can be of great value for renewable energy research, risk assessment and the mitigation of future coastal hazards.

@en

Although there have been many efforts in the literature to hindcast the patterns of longshore sediment transport (LST), they mainly disregarded uncertainty issues. Forcing datasets, wave transformation methods, and LST models are among the main sources of uncertainty in LST estimations. The combination of the aforementioned sources of uncertainty makes the estimation of LST patterns challenging for non-straight coastlines as the uncertainty ranges might vary from site to site. In this paper, a simple ensemble modeling framework was employed to investigate LST rate uncertainty at seven sites along a non-straight coastline (Gold Coast, Australia). The ensemble was formed by two different forcing datasets (i.e., integral parameters of total wave energy from two hindcast datasets of ERA5 and CAWCR), two different wave transformation methods (i.e., inclusion and exclusion of local wind effects), and eight LST models (i.e., bulk formulae and process-based models). Moreover, the relative importance of each source of uncertainty was ranked using the ANOVA-variance-based model. Finally, the weighted ensemble mean was used to investigate intra- and inter-annual variability of LST rates. The results showed that the range of uncertainty of LST rates for open coasts of Gold Coast is much higher than that of semi-sheltered coasts. On annual scale, for open coasts, 40% to 50% of total uncertainty was due to the choice of wave transformation methods, while for semi-sheltered coasts, it was 20%–30%. Moreover, almost for all sites, 30% to 50% of total uncertainty was controlled by the choice of LST models and the interaction of wave transformation methods and LST models. Although the weighted ensemble mean could provide an estimate of LST patterns along the coast, addressing the residual uncertainties (arising from other sources, discussed at the end of this paper), in future works, would help increase the certainty of the estimations.

@en