The present study aims to analyze the available Ocean Model Intercomparison Project Phase-2 (OMIP2) models from the Coupled Model Intercomparison Project Phase-6 (CMIP6) group in representing the seasonal mean variations of temperature and salinity. The skill of the models in representing seasonal mean biases of temperature and salinity is assessed with World Ocean Atlas and Argo observations over the Tropical Indian Ocean (TIO), Bay of Bengal (BoB), Arabian Sea (AS), and Southern Indian Ocean. It is identified that most of the individual models and multi-model mean of OMIP2 models exhibit a cold (surface) and warm (subsurface) temperature bias over the entire TIO. The salinity analysis reveals that most of the TIO (except the equatorial TIO) is dominated by saltier biases, especially predominant over the south AS and the western BoB. Vertical shear of horizontal currents (VSHC) and the Brunt-Väisälä frequency have been analyzed to understand the stability of the Ocean, suggesting that the VSHC contributes to the vertical mixing resulting in weak stratification, is mainly responsible for the persistence of surface cold and subsurface warm biases. In addition, Freshwater transport (FWT) is estimated at different straits, suggesting that FWT can modulate the salinity in the fresh region of TIO. The current study summarizes the improvement and necessity of ocean models to depict vertical hydrodynamic conditions for skillful seasonal forecasts accurately.

Pan et al. (2025) claim to uncover numerous timing and datum errors in the University of Hawai‘i Sea Level Center Research Quality dataset and declare that the dataset is “less reliable than previously assumed”. We argue that Pan et al. overstate the scientific implications of unresolved issues in the dataset and demonstrate that Pan et al. have failed to consider the extensive metadata accompanying the observations, which already document most of the issues Pan et al. claim to expose. We also clarify the nature of the Research Quality dataset, which is not guaranteed to be free of errors, and assert that when data and metadata are used together as intended, the dataset remains a reliable basis for scientific research.

A thermocline inversion model based on multi-head attention mechanism within a neural network framework is developed to estimate and analyze the ocean thermocline features, including depth (updepth and base), thickness, and intensity, in the western Pacific Ocean. This model employs Argo-derived thermocline product alongside various satellite remote sensing observations of ocean surface parameters, such as sea surface height, salinity, temperature, and wind. Specifically, three independent inversion models are executed using a dataset spanning the previous five years for training purposes, with the resulting model parameters used to estimate thermocline features in March, June, September, and December of 2016. The analysis reveals that thermocline updepth mainly located in the east of the Kuroshio extension area occurring in winter and spring; the seasonal distribution of the thermocline base is characterized by deeper depths at higher latitudes in the northern hemisphere during winter and spring, and from summer to winter in the southern hemisphere; the thermocline intensity in tropical regions is observed to be shallower yet stronger, exhibiting significant variations along the latitude with distinct seasonal changes. The seasonal distribution characteristics and variation trends of the updepth, base and intensity of the thermocline calculated by the multi-head attention neural network are generally consistent with the referenced Argo-derived thermocline dataset. Notably, the proposed intelligent inversion model for thermocline could also be utilized under condition of certain position with high flexibility, and exhibits faster convergence and greater accuracy compared to the classic Bi-LSTM model under comparable experimental conditions.

Abstract Long-term observations, which were collected in the Sea of Okhotsk coastal zone under open-water conditions and at times when the sea was ice-covered from mid-January to mid-March 2022, are interpreted. The project augments preceding work using three synchronously-recording, seafloor-mounted, pressure transducers sampling at 1 Hz to acquire time series of inshore wave oscillations. Units are deployed nearer to the shore and closer together than was done during the previous studies. Spectral analyses of open-water sea level oscillations perpendicular to the coast reveal a wide band of energy, suggestive of a propagating edge wave at about 5-min period with an offshore-to-inshore gain up to 1.5. Similar wave characteristics occur alongshore. Intriguingly, the peak edge wave period at 5 min migrates to 6–7 min when the sea becomes covered with ice, and narrower bands at periods from 0.5 to about 3 min emerge. Other period ranges also appear to be affected by the onset, presence and eventual disintegration of sea ice. Whilst a shift of the dominant edge wave period can be attributed to changes in spectral refraction arising from the incoming swells being low-pass filtered by sea ice, because of the observed tidal signal it is speculated in this case that the observed period adjustment and attendant expansion of bandwidth may also be associated with instability around the fundamental frequency. The potential for edge wave solitons to exist is explored.

Abstract Marine heatwaves (MHWs) pose a significant threat to marine ecosystems and economies. Predicting MHWs is essential for mitigating their impact, but remains a challenge. Despite considerable progress having been made in understanding the regional-scale drivers of MHWs, a significant knowledge gap remains when it comes to understanding the synoptic-scale processes associated with these events. In this study, we used self-organising maps to identify the synoptic-scale atmospheric and oceanic patterns associated with MHWs identified in four sub-regions of the Tasman Sea between 1985 and 2014. Our results reveal patterns associated with recurring, as well as distinct extreme warming events. We show that anomalous atmospheric influence is consistently present during MHWs and that the two most recurring patterns are linked to a La Niña climate phase. Distinct synoptic air-sea patterns are also identified in the 1997/98 El Niño event. Furthermore, we identify a ‘reservoir’ of warm subsurface temperatures from 2000-2014, during which MHW frequency increased two-fold. Importantly, we have identified patterns of persistent anomalous conditions before the onset of MHWs with timescales on the order of days for atmospheric conditions and weeks to months for oceanic conditions, providing valuable insight into MHW predictors. These findings highlight the importance of understanding synoptic-scale drivers of MHWs and timescales of recurring patterns for MHW prediction. The temporal variability observed in the lead-up to MHWs underscores the potential significance of factors such as surface-layer temperature and sea-level anomalies in capturing longer-term warming trends, likely influenced by sustained atmospheric stress and oceanic dynamics, whilst atmospheric conditions at onset precipitate the transition to the extreme warming thresholds.

Ocean wave-current interactions are important physical processes at the sea surface, which can potentially cause extreme sea states under certain conditions. Usually, such interactions are more notable in regions with strong waves and background currents. In this study, focusing on the Kuroshio Extension, we used buoy-measured and altimeter-derived wave data to determine variations in wave properties with the background currents. Statistically, the wave height can be underestimated (overestimated) by approximately 4% (3%) when the current and waves are in the opposite (same) direction. In regions with warm (cold) eddies, the wave height and wavelength inside the eddy are larger (smaller) than those outside by approximately 5% and 8% (4% and 4%), respectively, and the wave direction is deflected by 11° anticlockwise (clockwise). The wavenumber spectra of wave height and surface current speed are highly correlated with a power law of k− 2–k− 3 at scales of 20–200 km for swell-dominated cases. Additionally, the convergence and divergence of wave energy resulting from the current-induced refraction of swell are captured. From another perspective, the wave-induced Stokes drift calculated using the directional spectrum accounts for 54% of the reanalysis surface currents, and the accuracy of the estimated surface current can be improved by up to 14% by considering Stokes drift. This study provided quantitative analysis of observed surface wave variations in the Kuroshio Extension region from multiple perspectives.

Abstract Hydrodynamics and physical processes that occur at various length and time scales strongly influence coral reefs. Therefore, understanding the interactions between reefs, hydrodynamics and other physical processes is crucial for the maintenance and survival of reef systems. Coral reefs around the world are under increasing threat to global climate change, and additionally coral bleaching is a major concern for the health and survival of these reefs. Some marginal coral reefs are situated in areas where the complex ocean flow patterns interact with topographical features, providing possible refuges to rising ocean temperatures and coral bleaching. A prominent example is the Sodwana Bay coral reef system which has shown resilience to coral bleaching. This resilience has been attributed to cold water temperature anomalies that cause short-term temperature fluctuations on the reefs. This study explores hydrodynamics at various scales around the Sodwana Bay coral reefs and associated short-term temperature anomalies using a flexible mesh hydrodynamic model of the southwest region of the Indian Ocean, nested within a global ocean model. The nested hydrodynamic model better replicates the observed temperature anomalies when compared to the the reanalysis NEMO global ocean model. The higher model resolution around Sodwana results in less numerical mixing and smoothing of the temperature fields in the nearshore region when compared to the reanalysed NEMO global ocean model leading to a better replication of the local hydrodynamics around the Sodwana region. The anomalies investigated were associated with remote upwelling of cold water near the Delagoa Peninsula, followed by advection from the Delagoa Bight towards the Sodwana region. The separation of the strong intermittent southward stream from the Delagoa Peninsula is strongly linked to the upwelling at the Delagoa Peninsula. An analysis of the hydrodynamic patterns during the anomaly periods reveal that when the strong southward stream reattaches to the coastline, it typically does so south of Sodwana. The reattachment of the stream has an inertial effect and pushes the flow of water against the coastline which deflects the flow northwards up past Sodwana resulting in a northward current reversal along the Sodwana coastline which agrees with observed current reversals during the anomaly periods by insitu measurements taken on the Sodwana reefs. The model also revealed that local upwelling occurs within the Sodwana canyons during this event, making the water in the canyons colder than the surrounding water. When the locally upwelled water spreads over the reef system, the anomaly amplitude is enhanced by approximately 20 %.

Wind energy is a vital renewable energy source, and accurate Wind Speed Prediction (WSP) plays a key role in optimizing wind energy production and managing power grids effectively. However, predicting Wind Speed (WS) remains a significant challenge due to the inherently complex and dynamic behavior of wind flow. This paper introduces WindForecastX, an innovative approach that improves prediction accuracy by leveraging a dynamic unified ensemble learning model combined with advanced data assimilation techniques. The ability to accurately predict WS is vital for wind energy planning and monitoring. The accuracy of WSP has been limited because previous studies predominantly relied on data from a single location to develop models and predictions. The proposed WindForecastX model combines the strengths of ensemble learning and data assimilation techniques to enhance long-term WSPaccuracy. WindForecastX utilizes a Stacked Convolutional Neural Network (CNN) and bidirectional long short-term memory (BiLSTM) with a Data assimilation (SCBLSTM + DA) model, Adaptive Wind Speed Assimilation and Quality (AWAQ) incorporating WS observations from nearby locations. By leveraging these advanced techniques, including the Kalman filter, WindForecastX assimilates data from multiple sources to enhance the accuracy of WSP. To evaluate WindForecastX, we utilize real-world wind speed data collected from nine meteorological stations in the Tirunelveli district of Tamil Nadu, India. These stations are used for training and testing, with two stations designated as target stations for WSP. The results demonstrate that WindForecastX outperforms existing WSPmodels. Furthermore, WindForecastX exhibits reduced sensitivity to changes in the prediction time scale compared to standalone models, enhancing its reliability.

The magnitude of change in Chlorophyll-a (Chl-a) concentration in the Northern Arabian Sea (NAS) and northwestern Arabian Sea (NWAS), associated winds and Sea Surface Temperature (SST) that could potentially contribute to the fisheries management policies were investigated. The aforementioned parameters were examined from 2003 to 2021 by considering the area into 7 regions. Off Oman, monthly average winds delineate the presence of weak anticyclonic circulation during October and another one from February strengthening and migrating southwards till April. This feature is absent in May. For the first time, the warming in different regions was investigated and observed a temperature range from 0.4 °C to 0.8 °C with highest off Iran and least off central Oman. The Ekman Mass Transport (EMT) decreased with the highest off Pakistan (18.91%) and least off central Oman (0.55%). The decrease in Chl-a concentration was highest off Iran and least off Oman. The correlation of Chl-a with SST off Yemen was highest (-0.55) indicating that any change in SST will more conclusively influence Chl-a off Yemen than the rest of the regions. From 2003 to 2021 the wind speed off Oman increased up to 0.46 m/s while the other regions show decreased wind speed. The maximum correlation between winds and Chl-a was observed off Oman, indicating that changes in winds are more likely to affect Chl-a concentration in this region compared to other regions. The study statistically establishes the differential influence of SST and Sea Surface Winds (SSW) in the study area.

Observations of Japanese eel Anguilla japonica recruitment on Hokkaido’s coast in 2020 revealed a poleward shift of the species’ northern limit by several hundred kilometers. Field observations conducted from April to July 2021 in a river in southern Hokkaido, as reported in this study, identified for the first time the potential recruitment period of juvenile glass eels in Hokkaido, suggesting that the recruitment season may have commenced in May and concluded in July. The long-term trend of Japanese eel recruitment to Hokkaido was examined using a three-dimensional particle-tracking model. Virtual larvae were programmed to swim both horizontally and vertically, in addition to being transported by ocean currents, after their release near eastern Taiwan (Scenario 1) and northeastern Japan (Scenario 2). Scenario 1 showed increased recruitment in northern Japan and decreased recruitment in southern Japan during 2014–2023 compared to 1994–2003, which was attributed to the shift in the Kuroshio path. In Scenario 2, focusing on local processes near Hokkaido, the spatial variation in estimated glass eel recruitment exhibited patterns consistent with the natural variation in eel abundance observed across 95 rivers in southern Hokkaido in 2022, with higher recruitment in southeastern Hokkaido and lower recruitment in the Tsugaru Strait. Simulated recruitment trends from 1994 to 2023 indicated an increase in southeastern Hokkaido and a decrease in the Tsugaru Strait. The increased recruitment to southeastern Hokkaido was linked to the northward shifts of the Kuroshio/Kuroshio Extension and Oyashio currents, which weakened the southward currents in the confluence zone of the Kuroshio/Kuroshio Extension and Oyashio. In contrast, reduced recruitment in the Tsugaru Strait was associated with the strengthening of the east-flowing Tsugaru Current. These findings suggest that long-term fluctuations in ocean currents significantly influence the northern limit of anguillid eel habitats, highlighting the impact of changing oceanic conditions on their natural distribution.

This study introduces an open software Lagrangian particle tracking model designed for simulating the transport of microplastics (MPs), which incorporates crucial processes such as horizontal diffusion, beaching, resuspension, and biofouling. A sensitivity analysis for the parametrization of these processes was conducted on a regional scale – in the Gulf of Finland (GoF), the eastern Baltic Sea – employing very high-resolution hydrodynamic model output to drive the particle model. The sensitivity analysis underscores the impact of each process on the number of particles in the water column, sediments, beach areas, and at the domain boundary. The results indicate a significant impact of including or excluding a process and relatively high sensitivity of the parametrization on the simulated MP pathways. Stronger diffusion dispersed particles widely throughout the gulf and enhanced the export of the MPs out from the gulf. Beaching and biofouling were the major contributing factors to particle removal from the water column, while resuspension promoted settling in offshore areas. The number of beached particles rapidly increased during the wind-induced downwelling process. Scenario simulations, including parametrizations favoring or hindering MP transport, showed that a coincidence of several factors could lead to very diverse MP pathways. The analysis offers valuable insights, providing a foundation for tuning the model parameters to improve simulations with realistic loads in the future.

Signatures of Indian Ocean Dipole (IOD) over the southeastern Arabian Sea (SEAS) during two representative IOD years: negative IOD 2016 (nIOD16) and positive IOD 2019 (pIOD19) are studied. Sea Surface Temperature (SST) over the SEAS reveals negative anomalies during the Fall transition months of nIOD16, whereas positive SST anomalies are observed during both southwest monsoon season and Fall transition months in pIOD19. In contrast, during the first half of the year, noticeable positive SST anomalies are observed in nIOD16, whereas SST anomalies are weak in pIOD19. The cooling rate during the SWM season of nIOD16 was twice that during pIOD19 and lasted one additional month. To understand the aforementioned anomalous features in the SST evolution, a heat budget analysis is performed based on a combination of in-situ data, satellite measurements and ocean reanalysis. The enhanced cooling during the SWM season of nIOD16 was driven largely by the net heat flux. In addition, the penetrative component of shortwave radiation also favoured enhanced cooling during both the SWM season and Fall transition months of nIOD16, driven by a relatively small mixed layer depth (MLD). The negative anomaly in MLD during the SWM and Fall transition months of nIOD16 is then shown to be in line with corresponding positive anomalies in the Ekman pumping. In contrast, relatively less intense cooling during the SWM season of pIOD19 is attributed to small resultant net heat flux and warming by the residual term.

Abstract Monthly surface currents at 2 km resolution near the mouths of three U.S. east coast bays were obtained from high-frequency radars (Coastal Ocean Dynamics Application Radar, CODAR) during 2012–2024. The currents near these bays, the Chesapeake Bay (CB), the Delaware Bay (DB) and the New York Bay (NB) were analyzed to infer similarity and differences, as well as potential common forcing from regional and basin-scale factors. The contribution to flow variability from local and remote forcing is evaluated by comparing surface currents with (a) river discharges into each bay, (b) with local and regional winds, and (c) with the North Atlantic Oscillation (NAO). The results show that surface flow variability near the mouth of the bays is linked with all three driving factors. The three bays often show similar flow patterns not only of the seasonal cycle, but also during extreme weather events. For example, increased surface flow into the bays from the Atlantic Ocean is seen when hurricanes are observed offshore in the fall, and increased surface flow from the bays is seen during winter storms. During positive NAO phases, eastward flow from all three bays increased due to intensified westerly winds, while during negative NAO phases flow decreased with weakening winds in the region. River discharges into the bays increased during 2012–2019 but decreased during 2019–2024. This change in river discharge trend was especially large in the CB, resulting in a change in trends of the surface currents. Monthly currents of each bay are only weakly correlated with the monthly river flow (R ~ 0.2–0.3; P < 0.05), while the seasonal cycles of rivers and currents have higher correlations (R ~ 0.6–0.7). Local winds show high correlations with the monthly currents (R ~ 0.75) with the current direction ~ 45° to the right of the wind, as expected from Ekman theory. However, contributions to current variability from regional and remote factors cannot be ignored. The results demonstrate the complex nature of the currents near the mouth of bays since multiple drivers, including estuarine, coastal and open ocean dynamics contribute to the observed variability.

High-frequency sea level observations from long-term tide gauges are indispensable for a wide array of scientific inquires, including storm surges, ocean tides, and tsunamis. The University of Hawaii Sea Level Center (UHSLC) stands as a pivotal sea level database, offering meticulously curated, hourly sea level records across the global ocean through rigorous and time-consuming quality control measures. These records are widely recognized as reliable and suitable for scientific analysis. However, this consensus is challenged by our research, which examines 595 tide gauges from the UHSLC distributed worldwide. Utilizing a newly developed detection algorithm, we identify significant timing errors in 35 tide gauges from the UHSLC, with these errors potentially altering tidal amplitude estimates of main constituents by more than 10 mm. Notably, at specific tide gauges, such as Tanjong Pagar and Manila, abnormal tidal amplitude changes caused by timing errors can be larger than 50 mm. Most tide gauges with severe timing errors lie on the coasts of developing countries, suggesting inadequate management. Conversely, sea level observations from developed countries generally exhibit high quality, albeit with some potential timing errors identified. Beyond timing issues, datum shifts and other exceptions are also uncovered in the UHSLC database. Consequently, when analyzing ‘science-ready’ sea level observations, caution is warranted. Anomalous variations in sea levels and tidal parameters may not reflect physical phenomena but rather artifacts stemming from multifarious observation errors.

This study examines the impact of different turbidity products on the Aegean Sea surface physical characteristics, by performing twin-experiment simulations using a high-resolution regional ocean model. The turbidity products used include an in-situ based diffuse attenuation coefficient dataset at 490 nm (kd490, in m− 1) and a satellite derived kd490 product. Satellite turbidity products are broadly used in ocean simulations due to their spatiotemporal coverage and algorithm universality. Their validation and empirical components are trained mainly in phytoplankton driven regions and this may cause systematic differences in oligotrophic areas of variable optical properties’ composition. In the Aegean Sea, the in-situ based turbidity product accounts for the contribution of suspended particles in the solar heating profile, having further implications in the surface characteristics. The Aegean Sea upper-ocean thermohaline characteristics and general circulation patterns, reveal distinct differences between the twin-experiment simulations, showcasing mesoscale to locally induced impact of the turbidity variations. The turbidity impact on the air-sea interaction fluxes affects both thermodynamic processes i.e., solar radiation penetration and absorption in the water column, as well as dynamic processes i.e., momentum fluxes due to changes of the sea surface temperature and subsequently to the momentum drag coefficient. The Aegean Sea surface characteristics in the in-situ based turbidity product simulation, show a stronger decoupling between the North and the South Aegean Sea, when compared with the satellite derived turbidity product simulation. These results highlight the importance of incorporating more realistic turbidity products in ocean models, especially for optically complex regions such as the Aegean Sea.

Tropical cyclones (TCs) trigger vigorous upper ocean mixing, pumping surface warm water into the ocean subsurface. TC-induced ocean subsurface warming (OSWTC) may be transported meridionally and influence ocean heat distribution, particularly in the Northwestern Pacific where TCs are most active. Based on ocean reanalysis products, a Lagrangian Particle Tracking analysis is conducted by releasing virtual particles representing OSWTC in the Northwestern Pacific to investigate the major pathways of OSWTC in the ocean. It is found that more than 50% of total particles are primarily moved westward and then join the Kuroshio Current, flowing northeastward towards the Kuroshio Extension after one year of release. Meanwhile, other particles moved eastward along the Subtropical Countercurrent (~ 4%), westward across the Luzon Strait into the South China Sea (~ 8%), and southward along Mindanao Current (~ 4%). In addition, about a quarter of the OSWTC particles remain locally beneath the mixed layer after one year of release. These results indicate that the majority of OSWTC can contribute to the meridional heat transport in the Northwestern Pacific. Both background currents and westward propagation constrained by differential rotation (Beta effect) contribute to the movements. We also find that the meridional movements of OSWTC are sensitive to the formation latitudes. A more northward formation location of OSWTC corresponds to a more northward heat transport. Since the TC tracks move poleward with global warming, our results imply more northward TC-induced ocean heat transport in the future.

To understand the variation in dissolved organic matter (DOM) source characteristics in relation to the dynamics of physical and biological processes in marine environments, fluorescence spectroscopy was utilized as a measurement tool. Observations were carried out in the mixed layer waters of the eastern Arabian Sea (EAS), from coastal to offshore waters, spanning from south to north during the winter and spring inter-monsoon (WM and SIM) periods. Parallel factor analysis of excitation-emission matrix fluorescence spectra revealed protein (Tryptophan, T, and Tyrosine, B) and humic-like (UV-visible humic, C, and marine humic, M) components. T, M, and C components showed a significant negative correlation with salinity in the coastal waters of north EAS during the WM, attributed to external DOM sources. In contrast, offshore waters had in situ DOM, showing a positive correlation of Chl a with C and B due to convective mixing. During the SIM, in situ DOM sources prevailed, with a significant positive correlation of humic-like components with Chl a and heightened B component intensity. T was low in the south EAS during WM but increased in T and B during SIM, linked to reduced Bay of Bengal water intrusion via the west Indian coastal current (WICC). However, the Central EAS exhibited mixed characteristics compared to the north and south EAS, accumulating bacteria-derived and humic-like DOM, reflected in T: B and M: C ratios decreasing in WM and SIM, respectively.

Density currents are caused by the difference of potential energy among fluids of distinct densities. Turbidity currents correspond to a specific type of density current where there is a considerable concentration of suspended particles that are transported by the flow. This phenomenon is often seen as one of the main factors responsible for the generation of sedimentary deposits that form hydrocarbon reservoirs. Numerical modelling is one of the most used approaches to model this process, based on solving a set of partial differential equations employing numerical methods. This paper presents the main hydrodynamic and depositional characteristics of turbidity currents formed in an unbounded configuration. The SuLi code, used in the present work, solves the incompressible continuity and Navier-Stokes equations using the projection method and the Large Eddy Simulation methodology. A finite difference scheme was implemented in order to solve the linear advection-diffusion equation. The code verification was carried out based on results present in the literature, demonstrating that the implemented numerical schemes are adequate to satisfactorily reproduce the dynamics of turbidity currents. Significant results were obtained considering different particle sizes, allowing the development of a comprehensive methodology for estimating the run-out distance of the deposit and the sedimentation rate over time as a function of the characteristic diameter of the grain.

AbstractThe Mississippi Sound and Bight is a complex coastal system with shallow estuarine waters that are highly vulnerable to the effects of climate change and anthropogenic influences. In order to further our understanding of the system and provide natural resource managers and decision-makers with science-based guidance, a pre-operational coastal ocean forecast system has been developed using the Coupled Ocean Atmosphere Wave Sediment Transport Modeling System (COAWST). The COAWST application for Mississippi Bight (msbCOAWST) can be run in hindcast mode, pre-operational near real-time mode, or forecast mode and relies on other operational models including the National Water Model (NWM) for river forcing, the High Resolution Rapid Refresh model (HRRR) for atmospheric forcing, and the Navy Coastal Ocean Model (NCOM) for open ocean boundary forcing. msbCOAWST is being validated using data from a variety of in situ measurements that quantify coastal processes, including tides and water quality parameters (i.e. temperature and salinity). The highest model skill is obtained for temperature followed by water levels and salinity. msbCOAWST has been used to provide guidance for quantifying how freshwater influences derived from river diversion operations impact water quality in estuarine waters. While the model is initially developed to study water quality and circulation in pre-operational near real-time and forecast modes, it is currently being extended to include waves, sediment transport, and biogeochemistry and also linked with habitat suitability models and ecological models in hindcast mode so as to comprehensively reveal consequential environmental concerns due to the onset and persistence of hypoxia, seasonal and storm-induced waves with their associated impacts on the region’s fisheries and shellfisheries.

Using fine-scale measurements in the northwestern Sea of Japan, we estimated the vertical mixing parameters in the sea water column extended from the lower part of the thermocline downward to the near-bottom layer above the continental slope. The vertical scales of the turbulent patches were determined together with the turbulent dissipation rate and diapycnal diffusivity based on the conductivity, temperature, and depth data obtained by an Aqualog moored profiler from April through October 2015. The Thorpe-scale method was used to estimate the vertical mixing parameters as well as the vertical heat and salt fluxes. The enhanced vertical mixing, as well as enhanced downward heat flux and upward salt flux, occurred below the mixed layer despite strong density stratification. By comparing the turbulent dissipation rate and diapycnal diffusivity estimates derived via the Thorpe-scale method and the estimates of the same parameters obtained earlier by applying the finescale parameterization method to the same dataset in addition to the collocates of the current velocity measurements, the comparative accuracy evaluation of both methods was carried out. Finally, by compiling the vertical mixing data obtained by the Thorpe-scale method and the finescale parameterization approach, the generalized depth profile for the diapycnal diffusivity is presented for the depth range from 70 to 350 m.