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

Robust and reliable models are needed to understand how coastlines will evolve over the coming decades, driven by both natural variability and climate change. This study evaluated how accurately five popular ‘reduced-complexity’ models replicate multi-decadal shoreline change at Narrabeen-Collaroy Beach, a sandy embayment in Sydney, Australia. Measured shoreline positions derived from approximately monthly field surveys were used for 20-year calibration and 20-year validation periods. The models performed similarly on average but with large variability between transects. The set-up of several models was modified to compensate for their sensitivity to imperfect input wave data, and further site-specific improvements were identified. Capturing interannual to decadal-scale variability in cross-shore and longshore dynamics at this site was challenging for all five models. Models appeared to aggregate key processes at this timescale into parameter values rather than representing them directly. This suggests time-varying parameters or changes to model structure may be necessary for decadal-scale simulations.

@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

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

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

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

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

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

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

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

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

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

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

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

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

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

The Earth's landscapes are shaped by processes eroding, transporting and depositing material over various timespans and spatial scales. To understand these surface activities and mitigate potential hazards they inflict (e.g., the landward movement of a shoreline), knowledge is needed on the occurrences and impact of these activities. Near-continuous terrestrial laser scanning enables the acquisition of large datasets of surface morphology, represented as three-dimensional point cloud time series. Exploiting the full potential of this large amount of data, by extracting and characterizing different types of surface activities, is challenging. In this research we use a time series of 2,942 point clouds obtained over a sandy beach in The Netherlands. We investigate automated methods to extract individual surface activities present in this dataset and cluster them into groups to characterize different types of surface activities. We show that, first extracting 2,021 spatiotemporal segments of surface activity using an object detection algorithm, and second, clustering these segments with a Self-organizing Map (SOM) in combination with hierarchical clustering, allows for the unsupervised identification and characterization of different types of surface activities present on a sandy beach. The SOM enables us to find events displaying certain type of surface activity, while it also enables the identification of subtle differences between different events belonging to one specific surface activity. Hierarchical clustering then allows us to find and characterize broader groups of surface activity, even if the same type of activity occurs at different points in space or time.

@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

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