Ocean Modeling System Research Articles

Parameterization of the Vertical Mixing for the Luzon Undercurrent in the Northern Western Pacific Ocean

AbstractThe Luzon Undercurrent (LUC) is one of the most significant western boundary undercurrents in the northern western Pacific Ocean (WPO), essential for subsurface water transport and connecting subtropical–equatorial circulations. Over the last three decades, abundant observational progress has been made in revealing the basic features of the LUC. However, the dynamics and thus successful modeling of the LUC remain unresolved. In this work, we conducted a high‐resolution (3 km, 60 levels) numerical investigation of the WPO circulation using the China sea multi‐scale ocean modeling system. We paid particular attention to the vertical mixing, aiming to reasonably resolve the LUC by modifying the vertical mixing parameterization. Based on physics reasoning and experiments of physically based modeling, we designed an adaptive mixing scheme (AMS), which used a Munk‐like function in the thermocline for enhanced vertical mixing in the areas of small Richardson numbers. Using the AMS, we reproduced the two‐layer WPO circulations well and captured the inshore component of the LUC consistent with observations. Furthermore, we studied the dynamics of the LUC by analyzing the momentum balance and found that the LUC is primarily maintained by baroclinic pressure gradient force due to strong thermocline tilting near the western boundary. Enhanced vertical mixing in this highly sheared region crucially provides sufficient geostrophic support to sustain the LUC. Also, nonlinearity and submesoscale motions contribute positively to the LUC. This work advances physical understanding of the current‐undercurrent system and improves numerical capability in capturing the LUC and WPO circulations.

Implementation and evaluation of open boundary conditions for sea ice in a regional coupled ocean (ROMS) and sea ice (CICE) modeling system

Abstract. The Los Alamos Sea Ice Model (CICE) is used by several Earth system models where sea ice boundary conditions are not necessary, given their global scope. However, regional and local implementations of sea ice models require boundary conditions describing the time changes of the sea ice and snow being exchanged across the boundaries of the model domain. The physical detail of these boundary conditions regarding, for example, the usage of different sea ice thickness categories or the vertical resolution of thermodynamic properties, must be considered when matching them with the requirements of the sea ice model. Available satellite products do not include all required data. Therefore, the most straightforward way of getting sea ice boundary conditions is from a larger-scale model. The main goal of our study is to describe and evaluate the implementation of time-varying sea ice boundaries in the CICE model using two regional coupled ocean–sea ice models, both covering a large part of the Barents Sea and areas around Svalbard: the Barents-2.5 km​​​​​​​, implemented at the Norwegian Meteorological Institute (MET), and the Svalbard 4 km (S4K) model, implemented at the Norwegian Polar Institute (NPI). We use the TOPAZ4 model and a Pan-Arctic 4 km resolution model (A4) to generate the boundary conditions for the sea ice and the ocean. The Barents-2.5 km model is MET's main forecasting model for ocean state and sea ice in the Barents Sea. The S4K model covers a similar domain but it is used mainly for research purposes. Obtained results show significant improvements in the performance of the Barents-2.5 km model after the implementation of the time-varying boundary conditions. The performance of the S4K model in terms of sea ice and snow thickness is comparable to that of the TOPAZ4 system but with more accurate results regarding the oceanic component because of using ocean boundary conditions from the A4 model. The implementation of time-varying boundary conditions described in this study is similar regardless of the CICE versions used in different models. The main challenge remains the handling of data from larger models before its usage as boundary conditions for regional/local sea ice models, since mismatches between available model products from the former and specific requirements of the latter are expected, implying case-specific approaches and different assumptions. Ideally, model setups should be as similar as possible to allow a smoother transition from larger to smaller domains.

A New Norm for Seasonal Sea Ice Advance Predictability in the Chukchi Sea: Rising Influence of Ocean Heat Advection

Abstract Predictability of seasonal sea ice advance in the Chukchi Sea has been investigated in the context of ocean heat transport from the Bering Strait; however, the underlying physical processes have yet to be fully clarified. Using the Pan-Arctic Ice–Ocean Modeling and Assimilation System (PIOMAS) reanalysis product (1979–2016), we examined seasonal predictability of sea ice advance in early winter (November–December) and its source using canonical correlation analysis. It was found that 2-month leading (September–October) surface heat flux and ocean heat advection is the major predictor for interannual variability of sea ice advance. Surface heat flux is related to the atmospheric cooling process, which has influenced sea ice area in the southeastern Chukchi Sea particularly in the 1980s and 1990s. Anomalous surface heat flux is induced by strong northeasterly winds related to the east Pacific/North Pacific teleconnection pattern. Ocean heat advection, which is related to fluctuation of volume transport in the Bering Strait, leads to decrease in the sea ice area in the northwestern Chukchi Sea. Diagnostic analysis revealed that interannual variability of the Bering Strait volume transport is governed by arrested topographic waves (ATWs) forced by southeasterly wind stress along the shelf of the East Siberian Sea. The contribution of ocean heat flux to sea ice advance has increased since the 2000s; therefore, it is suggested that the major factor influencing interannual variability of sea ice advance in early winter has shifted from atmospheric cooling to ocean heat advection processes. Significance Statement Predictability of sea ice advance in the marginal Arctic seas in early winter is a crucial issue regarding future projections of the midlatitude winter climate and marine ecosystem. This study examined seasonal predictability of sea ice advance in the Chukchi Sea in early winter using a statistical technique and historical model simulation data. We identified that atmospheric cooling and ocean heat transport are the two main predictors of sea ice advance, and that the impact of the latter has become amplified since the 2000s. Our new finding suggests that the precise information on wind-driven ocean currents and temperatures is crucial for the skillful prediction of interannual variability of sea ice advance under present and future climatic regimes.

Evaluation of the Operational CMEMS and Coastal Downstream Ocean Forecasting Services During the Storm Gloria (January 2020)

Storm Gloria was the 10th named storm in Europe for the 2019–2020 winter season, and it severely affected Spain and France. This powerful storm represents an excellent study case to analyze the capabilities of the different ocean model systems available in the Spanish Mediterranean coasts to simulate extreme events, as well as to assess their suitability to enhance preparedness in maritime disasters with high impacts on coastal areas. Five different operational ocean forecasting services able to predict the storm-induced ocean circulation are evaluated. Three of the systems are delivered by the Copernicus Marine Service (hereafter CMEMS): the CMEMS global scale solution (GLO-1/12°), the specific Mediterranean basin scale one (MED-1/24°), and the regional solution for the Atlantic façade (IBI-1/36°), which includes also part of the western Mediterranean. These CMEMS core products are complemented with two higher resolution models focused on more limited areas, which provide operational forecasts for coastal applications: the WMOP system developed at the Balearic Islands Coastal Observing and Forecasting System (SOCIB) with a horizontal resolution of roughly 2 km and the Puertos del Estado (PdE) SAMOA systems with a 350-m resolution that cover the coastal domains of the Spanish Port Authorities of Barcelona, Tarragona, Castellón and Almeria. Both the WMOP and SAMOA models are nested in CMEMS regional systems (MED and IBI, respectively) and constitute good examples of coastal-scale-oriented CMEMS downstream services. The skill of these five ocean models in reproducing the surface dynamics in the area during Gloria is evaluated using met-ocean in situ measurements from numerous buoys (moored in coastal and open waters) and coastal meteorological stations as a reference to track the effects of the storm in essential ocean variables such as surface current, water temperature, and salinity throughout January 2020. Furthermore, modeled surface dynamics are validated against hourly surface current fields from the two high-frequency radar systems available in the zone (the SOCIB HF-Radar system covering the eastern part of the Ibiza Channel and the PdE one at Tarragona, which covers the Ebro Delta, one of the coastal areas most impacted by Gloria). The results assess the performance of the dynamical downscaling at two different levels: first, within the own CMEMS service (with their regional products, as enhanced solutions with respect to the global one) and second in the coastal down-streaming service side (with very high-resolution models reaching coastal scales). This multi-model study case focused on Storm Gloria has allowed to identify some strengths and limitations of the systems currently in operations, and it can help outlining future model service upgrades aimed at better forecasting extreme coastal events.

Comparison of direct measurements of sea ice thickness and snow height, CryoSat-2 observations and PIOMAS numerical estimates

The paper presents the results of comparison of contact measurements of ice thicknesses and snow heights performed at the points of the hydrological stations of the “Transarktika-2019” expedition in April 2019 north of Franz-Josef Land archipelago, with altimetry observations of the Cryosat-2 satellite and numerical estimates of the PIOMAS (Pan-Arctic numerical Ice and Ocean Modeling system and data Assimilation). A significantly better correspondence is predictably shown between the variability of the ice thickness directly measured and observed using the CryoSat-2 satellite than that for the numerical PIOMAS system estimates. A trial correction of the algorithm for calculating the ice thickness by replacing the climatic values of the ice density, snow density and height with data from direct measurements also predictably improves the quality of calculating the ice thickness from satellite observations. The mean / route mean square differences obtained for ice thicknesses (+44/+96 cm for uncorrected and +30/+95 cm for corrected CryoSat-2 satellite, –14/+81 cm for PIOMAS system) and snow height (–4/+12 cm for CryoSat-2 satellite, –15/+12 cm for PIOMAS system) show the scale of uncertainty in estimating sea ice thickness and snow height for areas dominated by medium and thick first-year ice. An anomaly of the ice thicknesses observed during the expedition is given in comparison with the background characteristics based on historical ice charting data for 1970s — 1990s, earlier High-Arctic aircraft “Sever” expeditions during 1950s — 1970s and the stated remote observations and numerical estimates for 2000s — 2019. Comparison shows that the AARI expedition was actually carried out in one of the most favorable years for ice research in the last decade for this region — the average ice thickness in April 2019 was 15 — 28 cm higher than that for the interval 2011 — 2019 with a slightly lower (1 — 2 cm) height of the snow cover. In the earlier period of the 1970 — 1990s this area was characterized by significantly thicker old ice with characteristic thicknesses ~ 60 cm more than in April 2019.

A Model Study of the Discharges Effects of Kaidu River on the Salinity Structure of Bosten Lake

The salinization of Bosten Lake, which is the largest lake in the arid or semi-arid region of Xinjiang, has increased. To study the effects of the inflow change of Kaidu River, the main recharge, on the salinity structure of Bosten Lake, the Estuarine, Coastal, and Ocean Modeling System with Sediments (ECOMSED), a basic three-dimensional numerical model, was used. The model is forced by realistic atmospheric forcing and river inflows, and verified by observational data. The model simulations can map the lake water movement processes and offer an understanding of the likely role of river runoff on the Bosten Lake salinity structure. The water mainly flows eastward at the surface and westward at the bottom. The river runoff of Kaidu River significantly affects the salinity structure of the southwestern part of the lake. The Kaidu River discharge mostly flowed northeastward along the west bank of the lake, so with decreasing Kaidu River discharge, the salinity of the region from the inlet of the river to its right (looking in the direction of the flow) subsequently increased. This study helps to the mastering of the dynamic change of salinity and provides some quantity information for controlling the salinization of Bosten Lake.

Impact of Model Physics on Seasonal Forecasts of Surface Air Temperature in the Arctic

The impacts of model physics and initial sea ice thickness on seasonal forecasts of surface energy budget and air temperature in the Arctic during summer were investigated based on Climate Forecast System, version 2 (CFSv2), simulations. The model physics changes include the enabling of a marine stratus cloud scheme and the removal of the artificial upper limit on the bottom heat flux from ocean to sea ice. The impact of initial sea ice thickness was examined by initializing the model with relatively realistic sea ice thickness generated by the Pan-Arctic Ice Ocean Modeling and Assimilation System (PIOMAS). Model outputs were compared to that from a control run that did not impose physics changes and used Climate Forecast System Reanalysis (CFSR) sea ice thickness. After applying the physics modification to either sea ice thickness initialization, the simulated total cloud cover more closely resembled the observed monthly variations of total cloud cover except for the midsummer reduction. Over the Chukchi–Bering Seas, the model physics modification reduced the seasonal forecast bias in surface air temperature by 24%. However, the use of initial PIOMAS sea ice thickness alone worsened the surface air temperature predictions. The experiment with physics modifications and initial PIOMAS sea ice thickness achieves the best surface air temperature improvement over the Chukchi–Bering Seas where the area-weighted forecast bias was reduced by 71% from 1.05 K down to −0.3 K compared with the control run. This study supports other results that surface temperatures and sea ice characteristics are highly sensitive to the Arctic cloud and radiation formulations in models and need priority in model formulation and validation.

Validation of newly designed regional earth system model (RegESM) for Mediterranean Basin

We present a validation analysis of a regional earth system model system (RegESM) for the Mediterranean Basin. The used configuration of the modeling system includes two active components: a regional climate model (RegCM4) and an ocean modeling system (ROMS). To assess the performance of the coupled modeling system in representing the climate of the basin, the results of the coupled simulation (C50E) are compared to the results obtained by a standalone atmospheric simulation (R50E) as well as several observation datasets. Although there is persistent cold bias in fall and winter, which is also seen in previous studies, the model reproduces the inter-annual variability and the seasonal cycles of sea surface temperature (SST) in a general good agreement with the available observations. The analysis of the near-surface wind distribution and the main circulation of the sea indicates that the coupled model can reproduce the main characteristics of the Mediterranean Sea surface and intermediate layer circulation as well as the seasonal variability of wind speed and direction when it is compared with the available observational datasets. The results also reveal that the simulated near-surface wind speed and direction have poor performance in the Gulf of Lion and surrounding regions that also affects the large positive SST bias in the region due to the insufficient horizontal resolution of the atmospheric component of the coupled modeling system. The simulated seasonal climatologies of the surface heat flux components are also consistent with the CORE.2 and NOCS datasets along with the overestimation in net long-wave radiation and latent heat flux (or evaporation, E), although a large observational uncertainty is found in these variables. Also, the coupled model tends to improve the latent heat flux by providing a better representation of the air–sea interaction as well as total heat flux budget over the sea. Both models are also able to reproduce the temporal evolution of the inter-annual anomaly of surface air temperature and precipitation (P) over defined sub-regions. The Mediterranean water budget (E, P and E–P) estimates also show that the coupled model has high skill in the representation of water budget of the Mediterranean Sea. To conclude, the coupled model reproduces climatological land surface fields and the sea surface variables in the range of observation uncertainty and allow studying air–sea interaction and main regional climate characteristics of the basin.

A fully coupled atmosphere–ocean wave modeling system for the Mediterranean Sea: interactions and sensitivity to the resolved scales and mechanisms

Abstract. It is commonly accepted that there is a need for a better understanding of the factors that contribute to air–sea interactions and their feedbacks. In this context it is important to develop advanced numerical prediction systems that treat the atmosphere and the ocean as a unified system. The realistic description and understanding of the exchange processes near the ocean surface requires knowledge of the sea state and its evolution. This can be achieved by considering the sea surface and the atmosphere as a continuously cross-talking dynamic system. Following and adapting concepts already developed and implemented in large-scale numerical weather models and in hurricane simulations, this study aims to present the effort towards developing a new, high-resolution, two-way fully coupled atmosphere–ocean wave model in order to support both operational and research activities. A specific issue that is emphasized is the determination and parameterization of the air–sea momentum fluxes in conditions of extremely high and time-varying winds. Software considerations, data exchange as well as computational and scientific performance of the coupled system, the so-called WEW (worketa-wam), are also discussed. In a case study of a high-impact weather and sea-state event, the wind–wave parameterization scheme reduces the resulted wind speed and the significant wave height as a response to the increased aerodynamic drag over rough sea surfaces. Overall, WEW offers a more realistic representation of the momentum exchanges in the ocean wind–wave system and includes the effects of the resolved wave spectrum on the drag coefficient and its feedback on the momentum flux.

Performance of One-Way and Two-Way Nesting Techniques Using the Shelf Circulation Modelling System for the Eastern Canadian Shelf

The performance of one-way and two-way nesting techniques is assessed in this study using model results produced by a regional ocean circulation modelling system for the eastern Canadian shelf. The assessment is made in terms of dynamical consistency between the parent model (PM) and child model (CM), representation of general circulation features, and reduction of numerical noise generated during the interaction of the PM and CM. It is demonstrated that the feedback from the CM to the PM in numerical experiments using two-way nesting ensures that the large-scale circulation produced by the PM and CM be dynamically consistent over the region where the two model domains overlap. In comparison with one-way nesting, two-way nesting leads to a better representation of coastal currents over the Gulf of St. Lawrence and the Scotian Shelf and improves the large-scale circulation in the results produced by the PM. This study also examines an alternative two-way nesting technique based on the semi-prognostic method in which differences between the PM and CM densities are used to adjust the horizontal momentum balance in the PM. Model results demonstrate the advantage of the semi-prognostic method in eliminating numerical noise during the feedback from the CM to PM while ensuring dynamical consistency between the two model components.

Wave modelling for the German Bight coastal-ocean predicting system

Ocean wave modeling has shown impressive developments, both in theoretical aspects as in the quality of the results available to users. The recent advances in development of the WAM wave model for forecasts applications at operational services and climate assessments for the North Sea and the German Bight are presented here. Ocean waves control the exchange of energy, momentum, heat, moisture, gas, etc. between the ocean and the atmosphere in the earth system. The impact of waves on currents and water levels in coastal areas is demonstrated. This is an important step towards a fully coupled atmosphere-wave- ocean modelling system. The synergy between wave observations and models for the North Sea and German Bight is increased on the road to improving the ocean state estimate and predictions in the coastal areas and generating up-to-date information, products and knowledge.

A three-dimensional water quality model to evaluate the environmental capacity of nitrogen and phosphorus in Jiaozhou Bay, China

Jiaozhou Bay has recently suffered from serious problems with pollution and eutrophication. Thus, land-based pollutant load must be reduced through a national control program. In this study, we developed a 3D water quality model to determine the environmental capacity of nitrogen and phosphorus in Jiaozhou Bay. A 3D hydrodynamic model (the estuarine, coastal, and ocean modeling system with sediments) was coupled with a water quality model, which was adapted from the dynamic model of nitrogen and phosphorus for a mesocosm near Jiaozhou Bay. The water quality model is divided into seven components: dissolved inorganic nitrogen, phosphate, phytoplankton, zooplankton, detritus, dissolved organic nitrogen, and dissolved organic phosphorus. Furthermore, it was calibrated based on data collected from Jiaozhou Bay in 2003. The proposed model effectively reproduced the spatiotemporal variability in nutrient concentration, thus suggesting that a reasonable numerical representation of the prototype system must be developed for further evaluation of environmental capacity.

Modeling and forecasting the distribution of Vibrio vulnificus in Chesapeake Bay.

To construct statistical models to predict the presence, abundance and potential virulence of Vibrio vulnificus in surface waters of Chesapeake Bay for implementation in ecological forecasting systems. We evaluated and applied previously published qPCR assays to water samples (n = 1636) collected from Chesapeake Bay from 2007-2010 in conjunction with State water quality monitoring programmes. A variety of statistical techniques were used in concert to identify water quality parameters associated with V. vulnificus presence, abundance and virulence markers in the interest of developing strong predictive models for use in regional oceanographic modeling systems. A suite of models are provided to represent the best model fit and alternatives using environmental variables that allow them to be put to immediate use in current ecological forecasting efforts. Environmental parameters such as temperature, salinity and turbidity are capable of accurately predicting abundance and distribution of V. vulnificus in Chesapeake Bay. Forcing these empirical models with output from ocean modeling systems allows for spatially explicit forecasts for up to 48 h in the future. This study uses one of the largest data sets compiled to model Vibrio in an estuary, enhances our understanding of environmental correlates with abundance, distribution and presence of potentially virulent strains and offers a method to forecast these pathogens that may be replicated in other regions.

Development and Role of Flux Coupler in a Coupled Atmosphere - Ocean Modeling System

ABSTRACT The exchange of momentum and radiation fluxes plays a vital role at the interface of atmosphere and ocean. In the present study an attempt has been made to develop a flux coupling algorithm for coupling atmosphere and ocean models as a single model system. The developed coupler has been used for the Kalpakkam region located in the east coast of India as a case study to understand the coastal processes and effect of thermal water pollution on coastal processes due to the establishment of nuclear power plants. The annual and seasonal variation of air temperature, wind speed, sea surface temperature and radiation fluxes; latent and sensible heat fluxes over the Kalpakkam region for a period of five years (2007 - July 2011) has been observed with available OAFlux data. During this period increasing trend of dynamical and thermodynamical variables has been observed. The Weather Research and Forecast (WRF) model has been configured and used to simulate atmospheric parameters over the Kalpakkam region. Thereafter the developed coupler has been used to extract the wind components and radiation parameters from the WRF model output, wherein the wind and radiation stresses were computed using bulk aerodynamic formulations. The heat fluxes show a positive correlation with the temperature gradient and wind speed. The robustness of developed algorithm was assessed through numerical experiments and statistical methods. The model simulated output and coupler estimated variable of wind speed and radiation fluxes have been validated with the available observations. The model results showed the overestimation of radiation fluxes than the observational data. The overestimated fluxes have been added as a correction factor in the coupler to maintain the constant variation in the fluxes. The comparison ascertains the fact that model and coupler performance is satisfactory in comparison with observational data.

A model study of the effects of river discharges and interannual variation of winds on the plume front in winter in Pearl River Estuary

A three-dimensional numerical model, Estuarine, Coastal and Ocean Modeling System with Sediments (ECOMSED), is employed to study the mechanism of plume front in the Pearl River Estuary (PRE) in detail. The model is forced by winds, tides and river discharges. The modeled results of tidal elevation, current velocity and salinity are in reasonable agreement with observational data in the PRE. By analyzing momentum and saltwater transport balance equations, it is found that the wind stress term, the pressure gradient term and the local time derivative term of velocity are dominant in the momentum equation, while the local time derivative term, the horizontal advective term and the vertical mixing term of salinity are dominant in the salinity transport equation. The residual current at surface along the plume front is seaward and stronger, whilst that in the bottom layer is mainly landward. A series of sensitive experiments is also run to examine the responses of plume front to changes of river discharges at different inlets in Lingdingyang Bay and interannual variation of northeast winds in winter. The location of plume front responds differently to the change of river discharge at different inlets. An increase in the river discharge at Dahu inlet seems to affect the location of plume front most among the four river inlets, it makes the plume front move eastward and southward wholly; the variation of river discharge at Nansha or Fengmamiao inlet on the location of plume front is more local and weaker; whilst the variation of river discharge at Hengmen inlet has little effect on the plume front. The location of plume front also changes in response to the interannual variation of northeast winds in winter, the stronger or the more eastward the winds are, the more westward the plume front moves, and only in the northern PRE, the response of plume front to the variation of wind speeds is largely different from that to the variation of wind directions.

A model tool for bathing water quality management: A case study on Salmonella occurrence at the southern Baltic coast

Despite Salmonella bacteria being important and wide-spread human pathogenic micro-organisms in coastal waters, they are not used as indicator organisms in the new EU bathing water directive 2006/7/EC. However, this does not exclude Salmonella as a potential threat in European bathing waters. For this reason we are analysing a Salmonella pollution event at the Polish seaside resort Miedzyzdroje in summer 2008 and present here a generic workflow for bathing water quality management developed through our results. We combine water quality measurements and a 3-D ocean model system with the goal to locate possible sources via backtracking simulations. Forward scenario simulations lead to spatial Salmonella pollution maps. A continuous pollution source in town is identified here as the most likely origin of bacterial occurrence. There is a high probability that the contamination impacted beaches within a minimum radius of 10 km from the source. Varied wind scenario results provide transport ranges and therefore the potential to improve the managerial response during a future Salmonella pollution in Miedzyzdroje. The overall impact of hygienic water pollution shows the necessity of new management instruments. Supra-regional recommendations are summarised in the generic workflow. It presents first action recommendations, aimed at safeguarding the health of bathers, in the event of a contamination without implementing a model study. If the pollution is continuous, an extended model study is recommended. Furthermore, the results are able to identify concrete sampling locations, based on regions most affected by the pollution event as shown within the case study.

Observation-based evaluation of surface wave effects on currents and trajectory forecasts

Knowledge of upper ocean currents is needed for trajectory forecasts and is essential for search and rescue operations and oil spill mitigation. This paper addresses effects of surface waves on ocean currents and drifter trajectories using in situ observations. The data set includes colocated measurements of directional wave spectra from a wave rider buoy, ocean currents measured by acoustic Doppler current profilers (ADCPs), as well as data from two types of tracking buoys that sample the currents at two different depths. The ADCP measures the Eulerian current at one point, as modelled by an ocean general circulation model, while the tracking buoys are advected by the Lagrangian current that includes the wave-induced Stokes drift. Based on our observations, we assess the importance of two different wave effects: (a) forcing of the ocean current by wave-induced surface fluxes and the Coriolis–Stokes force, and (b) advection of surface drifters by wave motion, that is the Stokes drift. Recent theoretical developments provide a framework for including these wave effects in ocean model systems. The order of magnitude of the Stokes drift is the same as the Eulerian current judging from the available data. The wave-induced momentum and turbulent kinetic energy fluxes are estimated and shown to be significant. Similarly, the wave-induced Coriolis–Stokes force is significant over time scales related to the inertial period. Surface drifter trajectories were analysed and could be reproduced using the observations of currents, waves and wind. Waves were found to have a significant contribution to the trajectories, and we conclude that adding wave effects in ocean model systems is likely to increase predictability of surface drifter trajectories. The relative importance of the Stokes drift was twice as large as the direct wind drag for the used surface drifter.

Suspended sediment transport in the offshore near Yangtze Estuary

Based on the Estuarine, Coastal and Ocean Modeling System with Sediments (ECOMSED) model, a 3-D hydrodynamic-transport numerical model was established for the offshore area near the Yangtze Estuary in the East China Sea. The hydrodynamic module was driven by tide and wind. Sediment module included sediment resuspension, transport and deposition of cohesive and non-cohesive sediment. The settling of cohesive sediment in the water column was modeled as a function of aggregation (flocculation) and deposition. The numerical results were compared with observation data for August, 2006. It shows that the sediment concentration reduces gradually from the seashore to the offshore area. Numerical results of concentration time series in the observation stations show two peaks and two valleys, according with the observation data. It is mainly affected by tidal current. The suspended sediment concentration is related to the tidal current during a tidal cycle, and the maximum concentration appears 1 h-4 h after the current maximum velocity has reached.

Towards an Ocean Forecasting System for the Aegean Sea

A first approach towards an ocean forecasting system for the Aegean Sea is presented in this paper. The ocean forecasting component is part of the POSEIDON system which includes high resolution atmospheric and wave modeling of the Aegean Sea together with an extended network of observational buoys for the continuous monitoring of physical, biological and chemical parameters. The ocean forecast model system consists of a high-resolution (1/20° times; 1/20°) implementation of the Princeton Ocean Model (POM) forced by forecast surface fluxes of momentum, heat and freshwater provided by the regional ETA atmospheric model. The climatological runs of the model reproduce the basic features of the general circulation in the Aegean Sea. In the forecasting mode the model shows skill in forecasting the sea surface temperature (SST) phenomenology when compared with satellite-derived SST. Sensitivity tests (such as the inclusion of sea surface atmospheric pressure) have been performed in order to improve the model's n...

Assimilation of Satellite Altimetry into a Western North Pacific Operational Model

An ocean data assimilation system, COMPASS-K (the Comprehensive Ocean Modeling, Prediction, Analysis and Synthesis System in the Kuroshio-region), has been developed at the Meteorological Research Institute (MRI). The purposes of the development are understanding ocean variability in the Kuroshio region as a local response to a global climate change with assimilated four-dimensional data sets, development of an operational system in the Japan Meteorological Agency, and for the GODAE (Global Ocean Data Assimilation Experiment) project. The model is an eddy permitting version of an MRI-OGCM. Space-time decorrelation scales of ocean variability are estimated with TOPEX I POSEIDON (TIP) altimeter data. Subsurface temperature and salinity fields are projected from the TIP altimeter data with a statistical correlation method and are assimilated into the model with a time-retrospective nudging scheme. Seasonal variation in the western North Pacific is investigated. Realistic space-time distribution of the physical quantities, the path of Kuroshio and its separation from Honshu are captured well. The Kuroshio volume transport is well reproduced in a reanalysis experiment of 1993. Preliminary predictability experiments are done in February and March, 1994. Predictability diagram shows the time scale of the predictability for temperature field is about 17 days in the Kuroshio south of Japan. This time scale is smaller than that in the North Atlantic.