Regional Ocean Modeling System Research Articles

On the evidence of helico-spiralling recirculation within coherent cores of eddies using Lagrangian approach

Oceanic eddies exhibit remarkable coherence and longevity compared to other transient features in the surrounding flow. They possess the ability to transport properties over extensive distances while maintaining their material identity intact. The Lagrangian Coherent Structure (LCS) framework has proven effective in capturing these coherent eddies, where they display a solid-body-like rotation. Although various LCS approaches have been employed to investigate different facets of coherent eddies, a comprehensive understanding of their three-dimensional structures and internal dynamics remains elusive. This study aims to advance our comprehension of coherent eddies’ structural characteristics and delve into the precise nature of their internal dynamics by utilizing the Lagrangian Averaged Vorticity Deviation approach. Two eddies, one cyclonic and the other anti-cyclonic, were chosen from a high-resolution simulation carried out in the Bay of Bengal using the Regional Ocean Modeling System (ROMS). The findings unveil that these eddies have three-dimensional coherent cores resembling gently tapered cones that are broader at the surface and gradually narrow towards the bottom. Intriguingly, the dynamically coherent core of these eddies exhibits simultaneous upwelling and downwelling while maintaining their volumes during advection due to persistent material coherence. The three-dimensional trajectories followed by the fluid parcels inside the coherent core are helical. Their two-dimensional horizontal projections show alternating spiral bands of upwelling and downwelling which are the manifestations of Vortex Rossby Waves. These observations lead to a conceptual framework of a three-dimensional helico-spiralling recirculation pattern within the coherent cores of eddies.

Tropical Atmospheric Intraseasonal Oscillations Leading to Sea Level Extremes in Coastal Indonesia during Recent Decades

Abstract Coastal flooding induced by sea surface high extreme (HEX) events is an increasing risk to human society and infrastructure as both urban growth in coastal areas and anthropogenic sea level rise continue, especially for island nations like Indonesia. This paper investigates the role of atmospheric intraseasonal oscillations (ISOs), which are dominated by the Madden–Julian oscillation (MJO), in forcing HEXs on the coasts of Indonesia bordering the Indian Ocean. We use satellite altimetry data from 1993 to 2021 and tide gauge observations to detect HEXs, and modeling experiments using both the Regional Ocean Modeling System and a Bayesian dynamic linear model to understand the forcing and processes. We find that HEXs exhibit strong seasonality, with most events occurring during boreal winter (December–February) and spring (March–May) that are dominated by seasonal-to-decadal and intraseasonal variability respectively. In 32% of the 56 HEX events detected, the amplitude of ISO-induced sea level anomalies (SLAs) exceeds that of seasonal-to-decadal SLAs. Surface wind stress associated with atmospheric ISOs is the major forcing for intraseasonal SLAs, and both the remote westerly wind stress from the Indian Ocean equator and northwesterly longshore wind stress at the Indonesian coasts play important roles in driving the HEXs. The MJO is the dominant cause of ISO-dominated HEXs and its impact shows strong seasonal differences. Spring MJOs are associated with stronger convective anomalies over the eastern Indian Ocean equator that drive stronger zonal winds across the equatorial basin that lead to more HEX events compared to winter MJOs when the convection is shifted southward.

Interannual Variability of Marine Nitrogen Fixation in the Western Tropical Atlantic

AbstractN2 fixation is a central process of the marine nitrogen cycle, yet little is known about how this process varies from year‐to‐year. Here, we investigate this variability in the Western Tropical Atlantic (WTA), a region where N2 fixation is prevalent, fueled, in part, by the nutrient input from the Amazon River. We use hindcast simulations from 1983 through 2019 with the Regional Oceanic Modeling System (ROMS) coupled to the Biogeochemical Elemental Cycling (BEC) model expanded to include Diatom‐Diazotroph Assemblages (DDAs). Throughout the WTA, we find a substantial level of interannual variability of N2 fixation, altering it by up to 33%, and locally by up to nearly 60%. Part of this interannual variability is driven by variations in the Amazon River discharge, which during high discharge events leads to reduced rates in the upper parts of the plume and strongly enhanced rates in the downstream part. This dipole pattern is a consequence of the riverine inputs of phosphorus and the competition with non‐diazotrophs for this limiting resource. Another part of the N2 fixation variability is driven by the Atlantic Meridional Mode (AMM), and the El Niño‐Southern Oscillation (ENSO). These processes alter N2 fixation primarily through the supply of the limiting nutrients phosphorus and iron by vertical mixing, while the role of top‐down control through grazing is relatively limited in our model. The high, and so far not well recognized interannual N2 fixation variability can lead to erroneous extrapolation of field measurements and inaccuracies in the marine nitrogen budget, if unaccounted for.

Model Validation and Applications of Wave and Current Forecasts from the Hong Kong Observatory's Operational Marine Forecasting System

The Hong Kong Observatory has been running an Operational Marine Forecasting System (OMFS) adapted from the Regional Ocean Modelling System (ROMS) coupled with the WaveWatch III and SWAN wave models to provide wave, current and sea temperature forecasts up to 144 hours twice a day since December 2021. To facilitate users’ interpretation of model forecasts of significant wave height and current speed in coastal predictions and open seas which are of particular significance in high wind situations, model forecasts were validated against moored buoy observations and wave recorder measurements near the shores of Hong Kong and drifting buoy data over the South China Sea, as well as Mercator Ocean model reanalysis in 2022. The validation results showed that the wave forecasts generally agreed well with the buoy observations with coefficient of determination (R2) of around 0.7 and root-mean-square error (RMSE) of less than 0.2 m up to 72 hours ahead. The R2 for sea current forecasts ranged between 0.4 and 0.6, and the RMSE was around 8 to 11 cm/s in near shores up to T+144 forecast hours. Validation against drifting buoy demonstrated that the trend of current forecasts generally agreed well with the measurements. RMSE of surface current forecasts over open seas ranged from 19 cm/s for 24-hour forecast to around 30 cm/s for 144-hour forecast when compared against Mercator Ocean reanalysis. Results from the current downscaling approach could serve as a benchmark reference for HKO to enhance OMFS in the future. In this paper, applications of model forecasts in the provision of marine weather services in Hong Kong are also introduced.

Assessing the impact of subsurface temperature observations from fishing vessels on temperature and heat content estimates in shelf seas: a New Zealand case study using Observing System Simulation Experiments

We know that extremes in ocean temperature often extend below the surface, and when these extremes occur in shelf seas they can significantly impact ecosystems and fisheries. However, a key knowledge gap exists around the accuracy of model estimates of the ocean’s subsurface structure, particularly in continental shelf regions with complex circulation dynamics. It is well known that subsurface observations are crucial for the correct representation of the ocean’s subsurface structure in reanalyses and forecasts. While Argo floats sample the deep waters, subsurface observations of shelf seas are typically very sparse in time and space. A recent initiative to instrument fishing vessels and their equipment with temperature sensors has resulted in a step-change in the availability of in situ data in New Zealand’s shelf seas. In this study we use Observing System Simulation Experiments to quantify the impact of the recently implemented novel observing platform on the representation of temperature and ocean heat content around New Zealand. Using a Regional Ocean Modelling System configuration of the region with 4-Dimensional Variational Data Assimilation, we perform a series of data assimilating experiments to demonstrate the influence of subsurface temperature observations at two different densities and of different data assimilation configurations. The experiment period covers the 3 months during the onset of the 2017-2018 Tasman Sea Marine Heatwave. We show that assimilation of subsurface temperature observations in concert with surface observations results in improvements of 44% and 38% for bottom temperature and heat content in shelf regions (water depths< 400m), compared to improvements of 20% and 28% for surface-only observations. The improvement in ocean heat content estimates is sensitive to the choices of prior observation and background error covariances, highlighting the importance of the careful development of the assimilation system to optimize the way in which the observations inform the numerical model estimates.

Intercomparisons of Tracker v1.1 and four other ocean particle-tracking software packages in the Regional Ocean Modeling System

Abstract. Particle tracking is widely utilized to study transport features in a range of physical, chemical, and biological processes in oceanography. In this study, a new offline particle-tracking package, Tracker v1.1, is introduced, and its performance is evaluated in comparison to an online Eulerian dye, one online particle-tracking software package, and three offline particle-tracking software packages in a small, high-resolution model domain and a large coarser model domain. It was found that both particle and dye approaches give similar results across different model resolutions and domains when they were tracking the same water mass, as indicated by similar mean advection pathways and spatial distributions of dye and particles. The flexibility of offline particle tracking and its similarity against online dye and online particle tracking make it a useful tool to complement existing ocean circulation models. The new Tracker was shown to be a reliable particle-tracking package to complement the Regional Ocean Modeling System (ROMS) with the advantages of platform independence and speed improvements, especially in large model domains achieved by the nearest-neighbor search algorithm. Lastly, trade-offs of computational efficiency, modifiability, and ease of use that can influence the choice of which package to use are explored. The main value of the present study is that the different particle and dye tracking codes were all run on the same model output or within the model that generated the output. This allows some measure of intercomparison between the different tracking schemes, and we conclude that all choices that make each tracking package unique do not necessarily lead to very different results.

TYPHOON KARDING (NORU) STORM SURGE ANALYSIS USING THE COAWST MODELING SYSTEM

Abstract. The Philippines, frequently affected by typhoons, faces the hazard of storm surges. This study examined the Coupled-Ocean-Atmosphere-Wave-Sediment Transport Modeling System (COAWST) to simulate Typhoon Karding's September 2022 storm surge. COAWST integrates the Regional Ocean Modeling System (ROMS), Simulating WAves Nearshore (SWAN), and Weather Research and Forecasting (WRF) coupled using the Model Coupling Toolkit (MCT). Four setups were analyzed: i. ROMS only, ii. ROMS-SWAN, iii. ROMS-WRF, and iv. ROMS-SWAN-WRF, focusing on four variables: a. Surface air pressure, b. Wind speed, c. Free-surface elevation, and d. Significant wave height. Results show that the ROMS-WRF and ROMS-SWAN-WRF setups accurately simulated Typhoon Karding's track with minimal positional error and wind speed. However, the models overestimated the typhoon's minimum air pressure with p-biases of 7.74% (i and ii), 4.1% (iii), and 3.9% (iv), and RMSE values of 68.529 hPa (i and ii), 36.744 hPa (iii), and 36.789 hPa (iv). Additionally, water levels were underestimated, with RMSE ranging from 0.31 to 0.35 meters and p-biases from −72.56% to −154.89% at Baler, Aurora validation point. At the Real, Quezon validation point, RMSE and p-bias ranged from 0.30 to 0.34 meters and −84.80% to −166.44%, respectively. Nonetheless, the models were able to simulate the storm surge and significant wave height at Baler and Real points similar to recorded data, with setup iv performing best in storm surge simulation. In summary, COAWST may be employed for typhoon simulations, with coupling being able to increase accuracy.

Dynamical systems for remote validation of very high-resolution ocean models

This paper presents and investigates a novel methodology for validating high-resolution ocean models using satellite imagery. High-resolution ocean models provide detailed information in coastal areas where other available data products are too coarse. Models are usually fitted by comparing them with observations; However, accessing in situ data in all small coastal areas is not feasible, as in situ observations are scarce and obtained through dedicated ships or instruments in limited and selected regions. Our work aims to use alternative remote sensing information to overcome this challenge. The approach involves establishing connections between the satellite observations and the outcomes of various computational experiments carried out using the Regional Ocean Modeling System, which allows the selection of different parameters to run the ocean model. These choices are not fully determined a priori and each one produces distinct outputs, which are then linked to the images through dynamical systems objects. By defining a performance index, we are able to quantify which experiment provides a better representation of the local ocean state.

Assimilation of satellite swaths versus daily means of sea ice concentration in a regional coupled ocean–sea ice model

Abstract. Operational forecasting systems routinely assimilate daily means of sea ice concentration (SIC) from microwave radiometers in order to improve the accuracy of the forecasts. However, the temporal and spatial averaging of the individual satellite swaths into daily means of SIC entails two main drawbacks: (i) the spatial resolution of the original product is blurred (especially critical in periods with strong sub-daily sea ice movement), and (ii) the sub-daily frequency of passive microwave observations in the Arctic are not used, providing less temporal resolution in the data assimilation (DA) analysis and, therefore, in the forecast. Within the SIRANO (Sea Ice Retrievals and data Assimilation in NOrway) project, we investigate how challenges (i) and (ii) can be avoided by assimilating individual satellite swaths (level 3 uncollated) instead of daily means (level 3) of SIC. To do so, we use a regional configuration of the Barents Sea (2.5 km grid) based on the Regional Ocean Modeling System (ROMS) and the Los Alamos Sea Ice Model (CICE) together with the ensemble Kalman filter (EnKF) as the DA system. The assimilation of individual swaths significantly improves the EnKF analysis of SIC compared to the assimilation of daily means; the mean absolute difference (MAD) shows a 10 % improvement at the end of the assimilation period and a 7 % improvement at the end of the 7 d forecast period. This improvement is caused by better exploitation of the information provided by the SIC swath data, in terms of both spatial and temporal variance, compared to the case when the swaths are combined to form a daily mean before assimilation.

Dispersion and fate of methane emissions from cold seeps on Hikurangi Margin, New Zealand

The influence of cold seep methane on the surrounding benthos is well-documented but the fate of dissolved methane and its impact on water column biogeochemistry remains less understood. To address this, the distribution of dissolved methane was determined around three seeps on the south-east Hikurangi Margin, south-east of New Zealand, by combining data from discrete water column sampling and a towed methane sensor. Integrating this with bottom water current flow data in a dynamic Gerris model determined an annual methane flux of 3 x 105 kg at the main seep. This source was then applied in a Regional Ocean Modelling System (ROMS) simulation to visualize lateral transport of the dissolved methane plume, which dispersed over ∼100 km in bottom water within 1 year. Extrapolation of this approach to four other regional seeps identified a combined plume volume of 3,500 km3 and annual methane emission of 0.4–3.2 x 106 kg CH4 y-1. This suggests a regional methane flux of 1.1–10.9 x 107 kg CH4 y-1 for the entire Hikurangi Margin, which is lower than previous hydroacoustic estimates. Carbon stable isotope values in dissolved methane indicated that lateral mixing was the primary determinant of methane in bottom water, with potential methane oxidation rates orders of magnitude lower than the dilution rate. Calculations indicate that oxidation of the annual total methane emitted from the five seeps would not significantly alter bottom water dissolved carbon dioxide, oxygen or pH; however, superimposition of methane plumes from different seeps, which was evident in the ROMS simulation, may have localized impacts. These findings highlight the value of characterizing methane release from multiple seeps within a hydrodynamic model framework to determine the biogeochemical impact, climate feedbacks and connectivity of cold seeps on continental shelf margins.

Exploring temporal and spatial SST patterns and their impact in the Arabian Sea: Insights from the regional ocean modeling system

The Regional Ocean Modeling System (ROMS) is integrated across the Arabian Sea (AS) from 1992 to 2021 and covers 33°E to 80°E and 5°N to 32°N at high horizontal resolutions of 1/4°(∼25 km). The study demonstrates that the variability is controlled by a seasonal high-resolution setup, with superior performance for various seasons during the months of December-January-February (DJF), March-April-May (MAM), June-July-August-September (JJAS) and October–November (ON). The Sea Surface Temperature (SST) anomalies in the JJAS season werestudied. SST dominates the seasonal variability of the Indian summer monsoon (ISM). Our results show that the ROMS can simulate seasonal variability and its effects on the upper ocean properties over the AS. In addition, the impact of monthly variation is controlled. In recreating the spatio-temporal distribution of surface as well as subsurface hydrographic parameters such as surface and subsurface temperature, the simulation results are reasonable like observation and reanalysis.

Comparison of 4-dimensional variational and ensemble optimal interpolation data assimilation systems using a Regional Ocean Modeling System (v3.4) configuration of the eddy-dominated East Australian Current system

Abstract. Ocean models must be regularly updated through the assimilation of observations (data assimilation) in order to correctly represent the timing and locations of eddies. Since initial conditions play an important role in the quality of short-term ocean forecasts, an effective data assimilation scheme to produce accurate state estimates is key to improving prediction. Western boundary current regions, such as the East Australia Current system, are highly variable regions, making them particularly challenging to model and predict. This study assesses the performance of two ocean data assimilation systems in the East Australian Current system over a 2-year period. We compare the time-dependent 4-dimensional variational (4D-Var) data assimilation system with the more computationally efficient, time-independent ensemble optimal interpolation (EnOI) system, across a common modelling and observational framework. Both systems assimilate the same observations: satellite-derived sea surface height, sea surface temperature, vertical profiles of temperature and salinity (from Argo floats), and temperature profiles from expendable bathythermographs. We analyse both systems' performance against independent data that are withheld, allowing a thorough analysis of system performance. The 4D-Var system is 25 times more expensive but outperforms the EnOI system against both assimilated and independent observations at the surface and subsurface. For forecast horizons of 5 d, root-mean-squared forecast errors are 20 %–60 % higher for the EnOI system compared to the 4D-Var system. The 4D-Var system, which assimilates observations over 5 d windows, provides a smoother transition from the end of the forecast to the subsequent analysis field. The EnOI system displays elevated low-frequency (>1 d) surface-intensified variability in temperature and elevated kinetic energy at length scales less than 100 km at the beginning of the forecast windows. The 4D-Var system displays elevated energy in the near-inertial range throughout the water column, with the wavenumber kinetic energy spectra remaining unchanged upon assimilation. Overall, this comparison shows quantitatively that the 4D-Var system results in improved predictability as the analysis provides a smoother and more dynamically balanced fit between the observations and the model's time-evolving flow. This advocates the use of advanced, time-dependent data assimilation methods, particularly for highly variable oceanic regions, and motivates future work into further improving data assimilation schemes.

An interactive visualization and data portal tool (PALTIDE) for relative sea level and palaeotidal simulations of the northwest European shelf seas since the Last Glacial Maximum

ABSTRACTRelative sea level (RSL) predictions based on glacial isostatic adjustment (GIA) simulations and palaeotidal predictions generated by hydrodynamic models using GIA‐generated palaeotopographies are available in the published literature, and datasets are available via data repositories. However, these data are often difficult to extract for specific locations or timeslices, requiring users to request datasets from corresponding authors. To overcome the intractability of these data and to enable users to interrogate datasets themselves without requiring offline requests, we have developed PALTIDE, an online visualization tool with intuitive user interface accessible at https://shiny.bangor.ac.uk/paleotidal/. The model domain for this interactive visualization tool is the northwest European continental shelf, covering the period from the Last Glacial Maximum (LGM) to the present day, and is based on previous GIA simulations by Bradley and colleagues and hydrodynamic simulations using Regional Ocean Modelling System (ROMS) published by Ward and colleagues. The tool is developed in R and utilizes a number of packages including shiny and bslib for the frontend, and arrow, raster and the tidyverse for backend data processing. The tool enables visualizations and data downloads for RSL, tidal amplitude and tide‐dependent parameters for any location within the model domain over 1000‐year timesteps from the LGM to the present.

Subsurface temperature estimates from a Regional Ocean Modelling System (ROMS) reanalysis provide accurate coral heat stress indices across the Main Hawaiian Islands

As ocean temperatures continue to rise, coral bleaching events around the globe are becoming stronger and more frequent. High-resolution temperature data is therefore critical for monitoring reef conditions to identify indicators of heat stress. Satellite and in situ measurements have historically been relied upon to study the thermal tolerances of coral reefs, but these data are quite limited in their spatial and temporal coverage. Ocean circulation models could provide an alternative or complement to these limited data, but a thorough evaluation against in situ measurements has yet to be conducted in any Pacific Islands region. Here we compared subsurface temperature measurements around the nearshore Main Hawaiian Islands (MHI) from 2010 to 2017 with temperature predictions from an operational Regional Ocean Modeling System (ROMS) to evaluate the potential utility of this model as a tool for coral reef management. We found that overall, the ROMS reanalysis presents accurate subsurface temperature predictions across the nearshore MHI region and captures a significant amount of observed temperature variability. The model recreates several temperature metrics used to identify coral heat stress, including predicting the 2014 and 2015 bleaching events around Hawaiʻi during the summer and fall months of those years. The MHI ROMS simulation proves to be a useful tool for coral reef management in the absence of, or to supplement, subsurface and satellite measurements across Hawaiʻi and likely for other Pacific Island regions.

Multivariate Upstream Kuroshio Transport (UKT) Prediction and Targeted Observation Sensitive Area Identification of UKT Seasonal Reduction

Variation and seasonal reduction in the Upstream Kuroshio Transport (UKT) have important impacts on surrounding climate and oceanic circulation systems. Therefore, reliable UKT prediction is crucial. In this paper, we propose an intelligent UKT prediction model, KuroshioNet, which is firstly pre-trained with simulation data generated by the Regional Ocean Modeling System (ROMS) and then fine-tuned with reanalysis data of the Simple Ocean Data Assimilation (SODA). Operating at a five-day time resolution and a 0.5°spatial resolution, KuroshioNet has the capability to predict multivariate fields associated with upstream Kuroshio, including 3D variables like velocity, temperature, as well as salinity and 2D variables like sea surface height. Subsequently, the UKT is computed from the predicted fields. We evaluate and analyze the experimental results, which show that KuroshioNet has a lead time of 55 days for UKT prediction. In order to enhance the physical interpretability of KuroshioNet, we conduct an ablation experiment to evaluate the impact of each predictor on prediction skill. It demonstrates that selecting zonal velocity, meridional velocity, temperature, salinity, and SSH contributes to UKT prediction by KuroshioNet. Besides, by analyzing model performance and visualizing what the convolutional kernels learn, we find that KuroshioNet, which has learned from ROMS data, is capable of obtaining better initial performance and acquiring more active kernels to better learn the features in SODA data. Furthermore, we identify the targeted observation sensitive area of UKT seasonal reduction by KuroshioNet with the saliency map method, which is situated to the east of upstream kuroshio. The sensitive area is consistent with the result identified by numerical models and yields 38.1% improvement on prediction demonstrated by observing system simulation experiments.

Modeled Coastal‐Ocean Pathways of Land‐Sourced Contaminants in the Aftermath of Hurricane Florence

AbstractExtreme precipitation during Hurricane Florence, which made landfall in North Carolina in September 2018, led to breaches of hog waste lagoons, coal ash pits, and wastewater facilities. In the weeks following the storm, freshwater discharge carried pollutants, sediment, organic matter, and debris to the coastal ocean, contributing to beach closures, algae blooms, hypoxia, and other ecosystem impacts. Here, the ocean pathways of land‐sourced contaminants following Hurricane Florence are investigated using the Regional Ocean Modeling System (ROMS) with a river point source with fixed water properties from a hydrologic model (WRF‐Hydro) of the Cape Fear River Basin, North Carolina's largest watershed. Patterns of contaminant transport in the coastal ocean are quantified with a finite duration tracer release based on observed flooding of agricultural and industrial facilities. A suite of synthetic events also was simulated to investigate the sensitivity of the river plume transport pathways to river discharge and wind direction. The simulated Hurricane Florence discharge event led to westward (downcoast) transport of contaminants in a coastal current, along with intermittent storage and release of material in an offshore (bulge) or eastward (upcoast) region near the river mouth, modulated by alternating upwelling and downwelling winds. The river plume patterns led to a delayed onset and long duration of contaminants affecting beaches 100 km to the west, days to weeks after the storm. Maps of the onset and duration of hypothetical water quality hazards for a range of weather conditions may provide guidance to managers on the timing of swimming/shellfishing advisories and water quality sampling.

Enhancing the Assimilation of SWOT Simulated Observations Using a Multi-Scale 4DVAR Method in Regional Ocean Modeling System

This paper presents an innovative approach to enhance the assimilation of high-resolution simulated observations, specifically targeting Surface Water Ocean Topography (SWOT) Ka-band Radar Interferometer Sea Surface Height (SSH) products, within the Regional Ocean Modeling System (ROMS). Responding to the demand for improved assimilation techniques, we developed a multi-scale Four-Dimensional Variational Data Assimilation (4DVAR) system, building upon validated fine-scale correction capabilities from prior studies. The multi-scale strategy was extended to the ROMS-4DVAR system, providing a comprehensive solution for assimilating high-resolution observations. Leveraging the Observing System Simulation Experiment (OSSE) framework, we conducted a twin experiment comprising a nature run and a free run case. Subsequently, synthetic SWOT SSH measurements were decomposed, considering the model configuration resolution. These components, derived from dense SSH observations, were integrated into a two-step 4DVAR assimilation scheme. The first cycle targets large-scale features for model field correction, and the updated analysis serves as the background for the second assimilation step, addressing fine-scale observation components. Comparisons with the primitive ROMS-4DVAR using a single-scale scheme highlight the superiority of the multi-scale strategy in reducing gaps between the model and the SSH observations. The Root Mean Squared Error (RMSE) is halved, and the Mean Absolute Percentage Error (MAPE) decreases from 2.237% to 0.93%. The two-step assimilation procedure ensures comprehensive multi-scale updates in the SSH field simulation, enhancing fine-scale features in the analysis fields. The quantification of three-dimensional-model dynamic fields further validates the efficiency and superiority of the multi-scale 4DVAR approach, offering a robust methodology for assimilating high-resolution observations within the ROMS.

Iron “ore” nothing: benthic iron fluxes from the oxygen-deficient Santa Barbara Basin enhance phytoplankton productivity in surface waters

Abstract. The trace metal iron (Fe) is an essential micronutrient that controls phytoplankton productivity, which subsequently affects organic matter cycling with feedback on the cycling of macronutrients. Along the continental margin of the US West Coast, high benthic Fe release has been documented, in particular from deep anoxic basins in the Southern California Borderland. However, the influence of this Fe release on surface primary production remains poorly understood. In the present study from the Santa Barbara Basin, in situ benthic Fe fluxes were determined along a transect from shallow to deep sites in the basin. Fluxes ranged between 0.23 and 4.9 mmol m−2 d−1, representing some of the highest benthic Fe fluxes reported to date. To investigate the influence of benthic Fe release from the oxygen-deficient deep basin on surface phytoplankton production, we combined benthic flux measurements with numerical simulations using the Regional Ocean Modeling System coupled to the Biogeochemical Elemental Cycling (ROMS-BEC) model. For this purpose, we updated the model Fe flux parameterization to include the new benthic flux measurements from the Santa Barbara Basin. Our simulations suggest that benthic Fe fluxes enhance surface primary production, supporting a positive feedback on benthic Fe release by decreasing oxygen in bottom waters. However, a reduction in phytoplankton Fe limitation by enhanced benthic fluxes near the coast may be partially compensated for by increased nitrogen limitation further offshore, limiting the efficacy of this positive feedback.

Regional Comparison of Performance between EnKF and EnOI in the North Pacific

Abstract The North Pacific is divided into different regions based on ocean currents and sea surface temperature (SST) distribution. Data assimilation is a useful tool for generating accurate ocean estimates because of the limited availability of observational data. This study compared the performances of two data assimilation methods, ensemble optimal interpolation (EnOI) and ensemble Kalman filter (EnKF), in various North Pacific subregions using an ocean model configured with the Regional Ocean Modeling System (ROMS). Both methods assimilated spaceborne SST observations, and the simulation results varied by subregion. The study found that EnKF and EnOI methods performed better than the control model in all regions when compared against satellite SST. EnOI reproduced SST as well as EnKF and required fewer computational resources. However, EnOI performed worse than the control model at sea surface height (SSH) in the equatorial region, while EnKF’s performance improved. This was due to the crushed mean state in the EnOI, which used long-term historical data as an ensemble member. El Niño–Southern Oscillation at the equator drove substantial interannual variability that crushed the ensemble mean of SSH in the EnOI. It is crucial to use a suitable assimilation method for the target area, considering the regional properties of ocean variables. Otherwise, the performance of the assimilated model may be even worse than that of the control model. While EnKF is better suited for regions with high variability in ocean variables, EnOI requires fewer computational resources. Thus, it is crucial to use a suitable assimilation method for accurately predicting and understanding the dynamics of the North Pacific.

Kinetic Energy Cascade Induced by the Interaction of Mean Flow, Topography, and Mesoscale Eddies East of Taiwan: A Scale-to-Scale Analysis

Abstract The scale-to-scale kinetic energy (KE) cascade induced by the nonlinear interaction among topography, the Kuroshio, and mesoscale eddies is systematically investigated in the coarse-graining framework based on simulated data from the well-validated Regional Ocean Model System. The KE transfer exhibits inhomogeneous spatial and temporal distributions and varies with length scale. During current–topography interaction, the KE transfers downscale across larger scales and reversely across smaller scales with an inherent separation scale of 150 km northeast of Taiwan, resulting in a significant positive net KE flux for mesoscale motions. The transfer around Suao Ridge is consistently downscaled with significant seasonal variation that is stronger in summer and weaker in winter. South of Suao Ridge, the transfer is one order of magnitude weaker and changes greatly with time. The cyclonic (anticyclonic) eddy weakens (enhances) KE transfer in most study areas. In particular, the cyclonic eddy reverses the transfer direction around Suao Ridge. The anticyclonic eddy triggers a significant bidirectional transfer south of Suao Ridge. Analyses show that the special arc-shaped topographic feature and northwestward Kuroshio intrusion current are responsible for the nature of bidirectional KE transfer northeast of Taiwan. The direction of mean current relative to the topography gradient determines the Rossby number magnitude and the KE transfer direction. The large-scale circulation determines the transfer intensity by changing the horizontal shear and barotropic instabilities. The KE transfer caused by nonlinear dynamics contributes significantly to the total anticyclonic eddy-induced net KE flux changes. In particular, the inverse KE cascade plays a key role in net KE flux changes in mesoscale motions east of Taiwan. Significance Statement In this work, we clarify the scale-to-scale kinetic energy transfer east of Taiwan. Energy transfer varies in different subareas. It is bidirectional northeast of Taiwan, consistently downscale around Suao Ridge, and weak and variable south of Suao Ridge. These transfer characteristics are closely related to the direction of mean flow relative to the topography gradient. Mesoscale eddies can trigger bidirectional transfer and even reverse the transfer direction. The reverse cascade is as important as the forward cascade for the total mesoscale energy change. This study not only shows an example of KE cascade when strong currents pass through continental slope, sea ridge, and open ocean, but also is of great significance for studying the dynamical and ecological environment variation of marginal sea.