Regional Ocean Modeling System Research Articles

Integrated 3D numerical modeling during La Niña and El Niño events using Regional Ocean Modeling System (ROMS) in Makassar Strait: preliminary study

This research showed the difference of thermocline layer depth during La Niña and El Niño events in the Makassar Strait, the main passage connecting the Pacific and Indian Oceans. La Niña event represented by 2011 and El Niño event represented by 2015. Three-dimensional numerical modeling using Regional Ocean Modeling System (ROMS) is used to determine ocean dynamic characteristics. The tides and atmospheric factor used as generating forces for coupled Ocean-Atmosphere-Wave-Sediment simulation model. Tidal elevation verification results from January to February 2017 showed a high confidence with RMSE = 0.148 m and MAPE = 13.51 %. Ocean circulation showed that water flows from north to south of the Makassar Strait in both conditions, El Niño and La Niña, indicating a transport from the Pacific to the Indian Oceans during those time. Two-dimensional seasonal sea surface temperature (SST) was used to determine the variability of thermocline layer depth. In general, based on simulation results, the depth of thermocline during La Niña was deeper than that during El Niño. The upper threshold of thermocline layer during El Niño phenomenon was shallower (average 27.29 m) compare to the threshold during La Niña (average 39.25 m). So do the lower threshold, during El Niño was shallower (average 160.47 m) compare to the threshold during La Niña (average 165.70 m). The thermocline thickness during El Niño (133.19 m) was found to be thicker compared to La Niña (126.43 m).

The changes of water mass characteristics using 3-dimensional Regional Ocean Modeling System (ROMS) in Balikpapan bay, Indonesia

Balikpapan Bay is located in the East Borneo Province of Borneo, Indonesia and directly connected to the Makassar Strait and also has very important role in determining the dynamics of ocean current, heat and freshwater content in the Makassar Strait. Numerical model has been used to simulate the dynamic process such as the changes of water mass in the Balikpapan Bay. We have used the terrain following three-dimensional Regional Ocean Modeling System (ROMS). Tides, river discharges, and atmospheric forcing (surface fields) have been used as generating forces for the model. The model has been simulated for three months from 1 October 2012 to 1 January 2013. The model has been validated by computing the RMSE, MAPE and Willmott’s index of agreement (d) using water level observation at the sampling station in Semayang port, Balikpapan. The verification showed a good agreement between model prediction and field observation data with RMSE = 7.8 cm, MAPE = 14.3% and d = 0.995 (with a value of perfect agreement equals to 1.0). The results analysed using vertical profile of salinity and temperature from 22 sampling stations with 11 stations located at the inner part of Balikpapan Bay. While the other stations located outside the Balikpapan Bay were used to analyse watermass. Water circulation is mostly dominated by the forcing from tides. The model results showed that the outer part of Balikpapan Bay is saltier than the inner part of Balikpapan bay. The highest average temperature occurred along coastal areas with value ~31 °C. The distribution of vertical salt transport showed that the water masses tend to be stratified during neap tides and mixed during spring tides condition.

Twentieth Century Seawater δ18O Dynamics and Implications for Coral‐Based Climate Reconstruction

AbstractThe oxygen isotopic composition of tropical coral skeletons (δ18O) is a crucial source of information on past El Niño/Southern Oscillation behavior. Both temperature and the δ18O of the surrounding seawater (δ18Osw) affect coral δ18O; a linear proportionality with δ18Osw is often used to infer past salinity variations, but the degree to which dynamical influences on δ18Osw may affect that relationship is still unclear. Here we use the isotope‐enabled Regional Ocean Modeling System to investigate the dynamics of δ18Osw and salinity variations in different twentieth century climate regimes. The dominant modes of δ18Osw variability are the background trend and eastern/central Pacific El Niño, similar to salinity; likewise, budget analysis reveals a strong impact of ocean dynamics (both advection and vertical mixing/diffusion) on δ18Osw variations during El Niño and La Niña events. These dynamics lead to alterations in the δ18Osw:salinity relationship across the Pacific: the linear approximation is most accurate near the eastern edge of the western Pacific warm pool but with nonlinearities due to large δ18Osw excursions during El Niño. At other sites, the δ18Osw:salinity relationship has substantial scatter and explains relatively little δ18Osw variance. Strikingly, interannual variability can lead to site‐specific δ18Osw:salinity regression coefficient changes of up to 0.4‰/psu on multiyear time scales. This calls the reliability of coral‐based estimates of El Niño/Southern Oscillation extremes into question and highlights the need for increased dynamical understanding as well as sustained in situ observations of δ18Osw and salinity variability.

Running the gauntlet: Connectivity between natal and nursery areas for Pacific ocean perch (Sebastes alutus) in the Gulf of Alaska, as inferred from a biophysical individual-based model

Little is known regarding the importance of early-life transport and dispersion mechanisms in determining recruitment variability for Pacific ocean perch (POP) in the Gulf of Alaska (GOA). These mechanisms influence the degree of, and variability in, connectivity between offshore natal areas (parturition sites) and inshore demersal nursery habitats for young-of-the year juveniles, and may thus play an important role in the “gauntlet” that individuals must survive from parturition to recruitment. As a first attempt to assess interannual variability in connectivity between natal and nursery areas for POP in the GOA in a synthetic manner, we developed a coupled biophysical individual-based model (IBM) for POP early life history and dispersal with simple representations of active vertical movement, passive horizontal movement, growth, and settlement in appropriate nursery habitat to integrate known early-life traits with variability in environmental forcing. We used an oceanographic model for the GOA based on the Regional Ocean Modeling system (ROMS) to provide the underlying daily physical environment to force the IBM for 1996–2011 and simulated hundreds of thousands of individual POP from parturition along the shelf break to settlement in inshore demersal nursery habitats as young-of-the-year. We used the IBM results to assess spatial patterns of annual “maximum potential” connectivity between presumed natal and nursery areas at alongshore scales of ~ 150 km, as well as the interannual variability in these patterns.Results showed that, even in the absence of mortality, most (> 70%) individuals were unsuccessful in dispersing from presumed natal areas along the continental shelf break to inshore nursery areas. For those that were successful, connectivity was directed in a counterclockwise fashion (southeast to northwest) around the GOA following prevailing current patterns. Typical dispersion distances were on the order of 100's of km alongshore, much larger than those inferred from genetic sampling. Natal areas from which the highest fractions of successful individuals originated were in the southeast GOA, while the nursery areas most frequently reached by those successful individuals were in the central GOA. POP from natal areas in the western GOA were consistently exported from the system and likely contribute little to the GOA population, although they may contribute to populations in the Aleutian Islands and eastern Bering Sea. We also found that annual indices derived from the connectivity matrices were not very strongly related to any of a suite of basin- and regional-scale environmental indices, reflecting the overall complexity and scale of the pathways POP in the GOA may undertake during their early life stages and suggesting that multiple drivers operating at different spatial and temporal scales influence connectivity patterns. Finally, while our results indicate that interannual variability in physical transport may have a substantial impact on connectivity, we found little support for the hypothesis that this alone drives variability in POP recruitment.

Variability in oceanographic barriers to coral larval dispersal: Do currents shape biodiversity?

The global center of marine biodiversity is located in the western tropical Pacific in a region known as the “Coral Triangle” (CT). This region is also considered the most threatened of all coral reef regions, because multiple impacts, including rising temperatures and coral bleaching, have already caused high mortality of reef corals over large portions of the CT. Larval dispersal and recruitment play a critical role in reef recovery after such disturbances, but our understanding of reproductive connectivity between reefs is limited by a paucity of observations. Oceanographic modeling can provide an economical and efficient way to augment our understanding of reef connectivity, particularly over an area as large as the CT, where marine ecosystem management has become a priority. This work combines daily averaged surface current velocity and direction from a Regional Ocean Modeling System developed for the CT region (CT-ROMS) with a Lagrangian particle tracking tool (TRACMASS) to investigate the probability of larval transport between reefs for a typical broadcasting coral. A 47-year historical simulation (1960–2006) was used to analyze the potential connectivity, the physical drivers of larval transport, and its variability following bi-annual spawning events in April and September. Potential connectivity between reefs was highly variable from year to year, emphasizing the need for long simulations. The results suggest that although reefs in this region are highly self-seeded, comparatively rare long-distance dispersal events may play a vital role in shaping regional patterns of reef biodiversity and recovery following disturbance. The spatial pattern of coral “subpopulations,” which are based on the potential connectivity between reefs, agrees with observed regional-scale patterns of biodiversity, suggesting that the physical barriers to larval dispersal are a first-order driver of coral biodiversity in the CT region. These physical barriers persist through the 21st Century when the model is forced with the Community Earth System Model (CESM) RCP8.5 climate scenario, despite some regional changes in connectivity between reefs.

Mean circulation of the Mid-Atlantic Bight from a climatological data assimilative model

The along-shelf momentum balance of the Mid-Atlantic Bight (MAB) coastal ocean includes a significant contribution from the along-shelf gradient in sea level. This sea level tilt, of order 10−7, and other features of the mean sea level are not captured well in global mean dynamic topography (MDT) derived from hydrographic observations or satellite altimetry and gravity data, and is poorly represented in global and basin scale dynamical models. This is problematic for applications that would use coastal satellite altimeter data to estimate total water level above datum. We have produced a MDT for the MAB using the Regional Ocean Modeling System (ROMS) with 4-dimensional variational (4D-Var) data assimilation configured to obtain climatological annual and monthly mean results. The observations assimilated were a regional hydrographic climatology of temperature and salinity, and long-term mean velocity from HF-radar, shipboard ADCP, and current-meters. Assimilation adjusts the 3-dimensional ocean state, boundary conditions, and air-sea fluxes to minimize the model-data misfit. The assimilation of mean velocity data is vital to obtaining a realistic circulation result. The MDT exhibits a strong across-shelf sea level slope in geostrophic balance with the southwestward mean flow. The subtle along-shelf tilt is recovered and is relatively uniform throughout the MAB inside the 50 m isobath, but on the southern outer shelf it reverses sign and drives significant across-isobath flow, partially draining the southward mean transport. In the north, across-shelf flow is offshore in the surface and bottom Ekman layers, but largely balanced locally by inflow in the interior depth range.

Cohesive and mixed sediment in the Regional Ocean Modeling System (ROMS v3.6) implemented in the Coupled Ocean–Atmosphere–Wave–Sediment Transport Modeling System (COAWST r1234)

Abstract. We describe and demonstrate algorithms for treating cohesive and mixed sediment that have been added to the Regional Ocean Modeling System (ROMS version 3.6), as implemented in the Coupled Ocean–Atmosphere–Wave–Sediment Transport Modeling System (COAWST Subversion repository revision 1234). These include the following: floc dynamics (aggregation and disaggregation in the water column); changes in floc characteristics in the seabed; erosion and deposition of cohesive and mixed (combination of cohesive and non-cohesive) sediment; and biodiffusive mixing of bed sediment. These routines supplement existing non-cohesive sediment modules, thereby increasing our ability to model fine-grained and mixed-sediment environments. Additionally, we describe changes to the sediment bed layering scheme that improve the fidelity of the modeled stratigraphic record. Finally, we provide examples of these modules implemented in idealized test cases and a realistic application.

Dynamical significance of tides over the Bay of Bengal

Tides play a significant role in the ocean surface circulations and vertical mixing thereby influencing the Sea Surface Temperatures (SST) as well. This, in turn, plays an important role in the global circulation when used as a lower boundary condition in a global atmospheric general circulation model. Therefore in the present study, the dynamics of tides over the Bay of Bengal (BoB) is investigated through numerical simulations using a high resolution (1/12°) Regional Ocean Modeling System (ROMS). Based on statistical analysis it is observed that incorporation of explicit tidal forcing improves the model performance in simulating the basin averaged monthly surface circulation features by 64% compared to the simulation without tides. The model simulates also Mixed Layer Depth (MLD) and SST realistically. The energy exchange between tidal oscillations and eddies leads to redistribution of surface kinetic energy density with a net decrease of 0.012 J m-3 in the western Bay and a net increase of 0.007 J m-3 in the eastern Bay. The tidal forcing also affects the potential energy anomaly and vertical mixing thereby leading to a fall in monthly MLD over the BoB. The mixing due to tides leads to a subsequent reduction in monthly SST and a corresponding reduction in surface heat exchange. These results from the numerical simulation using ROMS reveal that tides have a significant influence over the air-sea heat exchange which is the most important parameter for prediction of Tropical Cyclone frequency and its future variability over the BoB.

Impact of offshore eddies on shelf circulation and river plumes of the Sofala Bank, Mozambique Channel

Abstract A high-resolution, two-way nested Regional Ocean Modeling System, forced with monthly climatologies, has been set up for the Sofala Bank and adjacent deeper ocean of the Mozambique Channel to investigate the role of offshore mesoscale eddies on the shelf circulation, hydrographic structures and river plumes. The model is shown in comparison with available observations and published studies. Most known oceanographic features are reproduced by our model. We applied Self-Organizing Maps and showed that offshore passing eddies, depending on their strength and proximity to the shelf, modulate the shelf circulation and river plume direction and spread. The presence of a strong cyclonic eddy close to the shelf induces northward surface shelf currents. In contrast, the presence of a strong anticyclonic eddy close to the shelf induces a strong southward current over most of the shelf, except off Beira. Our analyses confirm that the plume of the Zambezi River is bi-directional. The southward-directed plume patterns, opposite to the dominant northwards, occur in response to nearby offshore anticyclonic eddies (26% of occurrence). This behavior could have an influence on water dispersal, shelf ecosystems and important fisheries. Therefore, offshore mesoscale eddies should be taken into account when studying the ocean dynamics of the Sofala Bank.

Effect of mesoscale wind stress-SST coupling on the Kuroshio extension jet

Effect of mesoscale wind stress-SST coupling on the Kuroshio extension jet is studied using the Regional Ocean Modeling System. The mesoscale wind stress perturbation (τMS) is diagnostically determined from modelled mesoscale SST perturbation (SSTMS) by using their empirical relationship derived from corresponding observation. From comparing two experiments with and without the τMS feedback, it is found that the interactively represented τMS-SSTMS coupling can modulate the kinetic energy along the Kuroshio extension jet, with little effect on the Kuroshio pathway. Similar results are also obtained in three additional sensitivity experiments, which consider half strength of the τMS, and the momentum flux and heat flux effect induced by τMS, respectively. That means simply taking into account the τMS-SSTMS coupling has little effect on improving the simulation of the Kuroshio Current system.

Assimilation of Different SST Datasets to a Coastal Ocean Modeling System in the Yellow and East China Sea

Kwon, K.; Choi, B.-J., and Lee, S.-H., 2018. Assimilation of Different SST Datasets to a Coastal Ocean Modeling System in the Yellow and East China Sea. In: Shim, J.-S.; Chun, I., and Lim, H.S. (eds.), Proceedings from the International Coastal Symposium (ICS) 2018 (Busan, Republic of Korea). Journal of Coastal Research, Special Issue No. 85, pp. 1041–1045. Coconut Creek (Florida), ISSN 0749-0208.Global gridded gap-free sea surface temperature (SST) datasets were quantitatively evaluated using independent insitu datasets in the Yellow and East China Sea from September 2011 to February 2012. The assimilation effect of the four different SST datasets on a regional ocean modeling system was examined. The SST from the gridded SST datasets was compared with two ocean buoy SST data and bi-monthly routine hydrographic observation data in the Yellow Sea (YS). The root-mean-square differences (RMSDs) of the OISST, MGDSST, OSTIA, and MWIR SST datasets relative to the in-situ SST were 0.64, 0.53, 0.62, and 0.91°C, respectively. The RMSDs of SST were higher than 1°C for all datasets in the southeastern corner of the YS. Ocean circulation was simulated with a regional ocean model, and the SST data from the four SST datasets were assimilated. Assimilation of the SST to the ocean model improved the ocean current as well as temperature distribution. The daily mean simulated temperature was compared to the observed SST and vertical profiles of the temperature at the ocean buoys and hydrographic observation stations. After the assimilation of OISST, MGDSST, OSTIA and MWIR datasets, RMSDs of the simulated SST were 0.78, 0.92, 0.66 and 0.86°C, respectively, at the buoy stations. The RMSDs of the simulated subsurface temperature were 1.60, 1.67, 1.56 and 1.70°C, respectively, at the hydrographic observation stations. The simulated SSTs from the model with assimilation of the OSTIA and OISST datasets had relatively smaller RMSDs.

The Impact of Different Forcing Agents on the Residual Circulation in a Tropical Estuary (Baía de Todos os Santos, Brazil)

ABSTRACT Santana, R.; Teixeira, C., and Lessa, G., 2018. The impact of different forcing agents on the residual circulation in a tropical estuary (Baia de Todos os Santos, Brazil). This study investigated the forcing of the residual circulation in a tropical estuary, Baia de Todos os Santos, Brazil (BTS). A numerical model (Regional Ocean Modeling System), forced with climatological means, was used to investigate the seasonal roles of the tide, wind, net heat, surface water fluxes, and river discharge on the residual circulation. The tide is the main driver of the circulation and its residual flow is structured at the bay mouth with a net-landward flow in the channel center and a net-seaward flow on the shoulders. The addition of wind drag affected the surface circulation, forcing W-bound (N-bound) currents in the spring–summer (fall–winter) months and generating a positive sea-level slope toward the continent. The addition of heat and surface water fluxes to the simulation established an incipient gravit...

A modelling approach to assess the effects of atmospheric nitrogen deposition on the marine ecosystem in the Bohai Sea, China

Atmospheric deposited nitrogen (AD-N) approximates or exceeds riverine input in many coastal ecosystems, such as the Bohai Sea (BHS) which is one of the most eutrophic coastal waters in China. We construct a three-dimensional (3D) physical-biogeochemical model to understand the influence of atmospheric dissolved inorganic nitrogen (DIN) on the intra- and inter-annual variations of phytoplankton blooms and nutrient dynamics in the BHS. The biological component, which is coupled to the Regional Ocean Modelling System (ROMS), is a simple but widely used Nutrient-Phytoplankton-Zooplankton-Detritus (NPZD) model with eight state variables. The model simulation successfully reproduces the spatial and temporal variations of observed DIN and dissolved inorganic phosphorus (DIP) and the climatological features of phytoplankton biomasses (chlorophyll a), which confirms the major role of air-transported nutrients in controlling standing stocks and nutrient limitations. The modelling results show that the relative contributions of deposited nitrogen to the total DIN in the seawater for the Bohai, Laizhou and Liaodong Bays and Central Bohai Basin reach 84.8%, 49.3%, 37.4% and 44.4% on average, respectively. The relative contribution ratio is approximately 54% for the entire Bohai region, causing a 56.5% increase in the phytoplankton biomass on average. These results also indicate that the effects of AD-N on nutrient and phytoplankton dynamics vary widely in regional areas because of the uneven spatial distribution of nitrogen deposition fluxes and also partially because of the hydrodynamic conditions, shortwave radiation and water temperatures.

Model-based assessment of a Northwestern Tropical Pacific moored array to monitor intraseasonal variability

The Northwestern Tropical Pacific Ocean (NWTPO) moorings observing system, including 15 moorings, was established in 2013 to provide velocity profile data. Observing system simulation experiments (OSSEs) were carried out to assess the ability of the observation system to monitor intraseasonal variability in a pilot study, where ideal “mooring-observed” velocity was assimilated using Ensemble Optimal Interpolation (EnOI) based on the Regional Oceanic Modeling System (ROMS). Because errors between the control and “nature” runs have a mesoscale structure, a random ensemble derived from 20–90-day bandpass-filtered nine-year model outputs is proved to be more appropriate for the NWTPO mooring array assimilation than a random ensemble derived from a 30-day running mean. The simulation of the intraseasonal currents in the North Equatorial Current (NEC), North Equatorial Countercurrent (NECC), and Equatorial Undercurrent (EUC) areas can be improved by assimilating velocity profiles using a 20–90-day bandpass-filtered ensemble. The root mean square errors (RMSEs) of the intraseasonal zonal (U) and meridional velocity (V) above 500 m depth within the study area (between 0°N–18°N and 122°E–147°E) were reduced by 15.4% and 16.9%, respectively. Improvements in the downstream area of the NEC moorings transect were optimum where the RMSEs of the intraseasonal velocities above 500 m were reduced by more than 30%. Assimilating velocity profiles can have a positive impact on the simulation and forecast of thermohaline structure and sea level anomalies in the ocean.

4DVAR Data Assimilation with the Regional Ocean Modeling System (ROMS): Impact on the Water Mass Distributions in the Yellow Sea

In this study, we evaluate the performance of the recently developed incremental strong constraint 4-dimensional variational (4DVAR) data assimilation applied to the Yellow Sea (YS) using the Regional Ocean Modeling System (ROMS). Two assimilation experiments are compared: assimilating remote-sensed sea surface temperature (SST) and both the SST and in-situ profiles measured by shipboard CTD casts into a regional ocean modeling from January to December of 2011. By comparing the two assimilation experiments against a free-run without data assimilation, we investigate how the assimilation affects the hydrographic structures in the YS. Results indicate that the SST assimilation notably improves the model behavior at the surface when compared to the non-assimilative free-run. The SST assimilation also has an impact on the subsurface water structure in the eastern YS; however, the improvement is seasonally dependent, that is, the correction becomes more effective in winter than in summer. This is due to a strong stratification in summer that prevents the assimilation of SST from affecting the subsurface temperature. A significant improvement to the subsurface temperature is made when the in-situ profiles of temperature and salinity are assimilated, forming a tongue-shaped YS bottom cold water from the YS toward the southwestern seas of Jeju Island.

Intercomparison of Antarctic ice-shelf, ocean, and sea-ice interactions simulated by MetROMS-iceshelf and FESOM 1.4

Abstract. An increasing number of Southern Ocean models now include Antarctic ice-shelf cavities, and simulate thermodynamics at the ice-shelf/ocean interface. This adds another level of complexity to Southern Ocean simulations, as ice shelves interact directly with the ocean and indirectly with sea ice. Here, we present the first model intercomparison and evaluation of present-day ocean/sea-ice/ice-shelf interactions, as simulated by two models: a circumpolar Antarctic configuration of MetROMS (ROMS: Regional Ocean Modelling System coupled to CICE: Community Ice CodE) and the global model FESOM (Finite Element Sea-ice Ocean Model), where the latter is run at two different levels of horizontal resolution. From a circumpolar Antarctic perspective, we compare and evaluate simulated ice-shelf basal melting and sub-ice-shelf circulation, as well as sea-ice properties and Southern Ocean water mass characteristics as they influence the sub-ice-shelf processes. Despite their differing numerical methods, the two models produce broadly similar results and share similar biases in many cases. Both models reproduce many key features of observations but struggle to reproduce others, such as the high melt rates observed in the small warm-cavity ice shelves of the Amundsen and Bellingshausen seas. Several differences in model design show a particular influence on the simulations. For example, FESOM's greater topographic smoothing can alter the geometry of some ice-shelf cavities enough to affect their melt rates; this improves at higher resolution, since less smoothing is required. In the interior Southern Ocean, the vertical coordinate system affects the degree of water mass erosion due to spurious diapycnal mixing, with MetROMS' terrain-following coordinate leading to more erosion than FESOM's z coordinate. Finally, increased horizontal resolution in FESOM leads to higher basal melt rates for small ice shelves, through a combination of stronger circulation and small-scale intrusions of warm water from offshore.

Exploring the Potential Impact of Greenland Meltwater on Stratification, Photosynthetically Active Radiation, and Primary Production in the Labrador Sea

AbstractIn July 2012, the surface of the Greenland Ice Sheet (GrIS) melted to an extent unprecedented over the last 100 years; we questioned the potential for such an extreme melt event to impact marine phytoplankton offshore. We hypothesized that stratification from meltwater could reduce light limitation for phytoplankton, and used a suite of numerical models to quantify the impact for 2003–2012. Because much of the 2012 meltwater discharged from southern Greenland, our study focused on the southwestern and southeastern coasts of Greenland, and the Labrador Sea. A 1‐D phytoplankton model used output from a Regional Ocean Modeling System (ROMS) coupled with a Regional Climate Model and a hydrological model of meltwater from runoff sources on the ice sheet, peripheral glaciers, and tundra. ROMS was run with and without meltwater to test the sensitivity of phytoplankton photosynthetic rates to the meltwater input. With meltwater, the pycnocline was shallower during late summer and early fall and thus light limitation on photosynthesis was reduced. Averaged over all years, added meltwater had the potential to increase gross primary production by 3–12% in the summer (July–August), and 13–60% in the fall (September–October). This meltwater effect was amplified when light was more limiting, and thus was greatest in the fall, under cloudier conditions, with higher self‐shading, and with more light‐sensitive phytoplankton groups. As the GrIS melt is projected to increase, late summer primary production in this region has the potential to increase as well, which could constitute an important biosphere response to high‐latitude climate change.

Island Wakes in Shallow Water

ABSTRACTIsland wakes in shallow water are investigated using the Regional Ocean Modeling System (ROMS). In contrast to deep water where bottom stress can be neglected in island wakes, shallow water implies that inhomogeneity in the bottom stress plays an important role in the wake vorticity generation. A series of numerical experiments are conducted to investigate wake formation and evolution in shallow water. It is found that the vertical structures of shallow-water and deep-water wakes are significantly different because of the presence of a density frontal jet, which results from the interaction between stratification and bottom topography. The frontal jet reaches its maximum within the bottom boundary layer over the shelf, giving rise to vorticity. The potential vorticity (PV) balance analysis reveals that frictional and diapycnal processes play different roles in the PV anomalies. With the absence of lateral stress (i.e., a sea mountain case), the surface vorticity becomes much weaker than that in the presence of an island.

How biophysical interactions associated with sub- and mesoscale structures and migration behavior affect planktonic larvae of the spiny lobster in the Juan Fernández Ridge: A modeling approach

The Juan Fernández Ridge (JFR) is a chain of topographical elevations in the eastern South Pacific (∼33–35°S, 76–81.5°W). Rich in endemic marine species, this ridge is frequently affected by the arrival of mesoscale eddies originating in the coastal upwelling zone off central-southern Chile. The impacts of these interactions on the structure and dynamics of the JFR pelagic system have, however, not been addressed yet. The present model-based study is focused on the coupled influence of mesoscale-submesoscale processes and biological behavior (i.e., diel vertical migration) on the horizontal distribution of planktonic larvae of the spiny lobster (Jasus frontalis) around the JFR waters. Two case studies were selected from a hydrodynamic Regional Ocean Modeling System to characterize mesoscale and submesoscale structures and an Individual-based model (IBM) to simulate diel vertical migration (DVM) and its impact on the horizontal distribution and the patchiness level. DVM behavior of these larvae has not been clearly characterized, therefore, three types of vertical mechanisms were assessed on the IBM: (1) no migration (LG), (2) a short migration (0–50 m depth, DVM1), and (3) a long migration (10–200 m depth, DVM2). The influence of physical properties (eddy kinetic energy, stretching deformation and divergence) on larval aggregation within meso and submesoscale features was quantified. The patchiness index assessed for mesoscale and submesoscale structures showed higher values in the mesoscale than in the submesoscale. However, submesoscale structures revealed a higher accumulation of particles by unit of area. Both vertical migration mechanisms produced larger patchiness indices compared to the no migration experiment. DVM2 was the one that showed by far the largest aggregation of almost all the aggregation zones. Larval concentrations were highest in the submesoscale structures; these zones were characterized by low eddy kinetic energy, negative stretching deformation, and slight convergence. Stretching deformation flow appeared to be triggered by the eddy-eddy interactions and the Robinson Island barrier effect, and it likely promotes the aggregation of the spiny lobster larvae in the Juan Fernández system. These results highlighted the importance of the coupled effect of physical (mesoscale and submesoscale oceanographic features) and biological processes (DVM) in the generation of larval patchiness and concentration of spiny lobster larvae around the JFR, which could be key for their survival and retention in those waters.

Sardine (Sardina pilchardus) larval dispersal in the Iberian upwelling system, using coupled biophysical techniques

The European sardine (Sardina pilchardus) is the most important small pelagic fishery of the Western Iberia Upwelling Ecosystem (WIUE). Recently, recruitment of this species has declined due to changing environmental conditions. Furthermore, controversies exist regarding its population structure with barriers thought to exist between the Atlantic-Iberian Peninsula, Northern Africa, and the Mediterranean. Few studies have investigated the transport and dispersal of sardine eggs and larvae off Iberia and the subsequent impact on larval recruitment variability. Here, we examine these issues using a Regional Ocean Modeling System climatology (1989–2008) coupled to the Lagrangian transport model, Ichthyop. Using biological parameters from the literature, we conduct simulations that investigate the effects of spawning patchiness, diel vertical migration behaviors, and egg buoyancy on the transport and recruitment of virtual sardine ichthyoplankton on the continental shelf. We find that release area, release depth, and month of release all significantly affect recruitment. Patchiness has no effect and diel vertical migration causes slightly lower recruitment. Egg buoyancy effects are significant and act similarly to depth of release. As with other studies, we find that recruitment peaks vary by latitude, explained here by the seasonal variability of offshore transport. We find weak, continuous alongshore transport between release areas, though a large proportion of simulated ichthyoplankton transport north to the Cantabrian coast (up to 27%). We also show low level transport into Morocco (up to 1%) and the Mediterranean (up to 8%). The high proportion of local retention and low but consistent alongshore transport supports the idea of a series of metapopulations along this coast.