Ship-Mounted Real-Time Surface Observational System on board Indian Vessels for Validation and Refinement of Model Forcing Fields

A network of ship-mounted real-time Automatic Weather Stations integrated with Indian geosynchronous satellites Indian National Satellites (INSATs)] 3A and 3C, named Indian National Centre for Ocean Information Services Real-Time Automatic Weather Stations (I-RAWS), is established. The purpose of I-RAWS is to measure the surface meteorological-ocean parameters and transmit the data in real time in order to validate and refine the forcing parameters (obtained from different meteorological agencies) of the Indian Ocean Forecasting System (INDOFOS). Preliminary validation and intercomparison of analyzed products obtained from the National Centre for Medium Range Weather Forecasting and the European Centre for Medium-Range Weather Forecasts using the data collected from I-RAWS were carried out. This I-RAWS was mounted on board oceanographic research vessel Sagar Nidhi during a cruise across three oceanic regimes, namely, the tropical Indian Ocean, the extratropical Indian Ocean, and the Southern Ocean. The results obtained from such a validation and intercomparison, and its implications with special reference to the usage of atmospheric model data for forcing ocean model, are discussed in detail. It is noticed that the performance of analysis products from both atmospheric models is similar and good; however, European Centre for Medium-Range Weather Forecasts air temperature over the extratropical Indian Ocean and wind speed in the Southern Ocean are marginally better.

Yanai waves in the western equatorial Indian Ocean

Observations and models have shown the presence of intraseasonal fluctuations in 20-30-day and 10-20-day bands in the equatorial Indian Ocean west of 60 degrees E (WEIO). Their spatial and temporal structures characterize them as Yanai waves, which we label low-frequency (LFYW) and high-frequency (HFYW) Yanai waves, respectively. We explore the dynamics of these intraseasonal signals, using an ocean general circulation model (Modular Ocean Model) and a linear, continuously stratified model. Yanai waves are forced by the meridional wind tau(y) everywhere in the WEIO most strongly during the monsoon seasons. They are forced both directly in the interior ocean and by reflection of the interior response from the western boundary; interference between the interior and boundary responses results in a complex surface pattern that propagates eastward and has nodes. Yanai waves are also forced by instabilities primarily during June/July in a region offshore from the western boundary (52-55 degrees E). At that time, eddies, generated by barotropic instability of the Southern Gyre, are advected southward to the equator. There, they generate a westward-propagating, cross-equatorial flow field, v(eq), with a wave number/frequency spectrum that fits the dispersion relation of a number of Yanai waves, and these waves are efficiently excited. Typically, Yanai waves associated with several baroclinic modes are excited by both wind and eddy forcing; and typically, they superpose to create beams that carry energy vertically and eastward along ray paths. The same processes generate LFYWs and HFYWs, and hence, their responses are similar; differences are traceable to the property that HFYWs have longer wavelengths than LFYWs for each baroclinic mode.

Annual cycle of stratification and tidal fronts in the Bohai Sea: A model study

In general, competition between buoyancy mechanisms and mixing dynamics largely determines the water column structure in a shelf sea. A three dimensional baroclinic ocean model forced by surface heat fluxes and the 2.5 order Mellor-Yamada turbulence scheme is used to simulate the annual cycle of the temperature in the Bohai Sea. The difference between the sea surface temperature (SST) and sea bottom temperature (SBT) is used to examine the evolution of its vertical stratification. It is found that the water column is well-mixed from October to March and that the seasonal thermocline appears in April, peaks in July and then weakens afterwards, closely following the heat budget. In addition, the Loder parameter based on the topography and tidal current amplitude is also computed in order to examine tidal fronts in the BS, which are evident in summer months when the wind stirring mechanism is weak.

Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

IMPACTS OF THE TROPICAL PACIFIC COUPLED PROCESS ON THE INTERANNUAL VARIABILITY IN THE INDIAN OCEAN

The basic features of climatology and interannual variations of tropical Pacific and Indian Oceans were analyzed using a coupled general circulation model (CGCM), which was constituted with an intermediate 2.5-layer ocean model and atmosphere model ECHAM4. The CGCM well captures the spatial and temporal structure of the Pacific El Nino-Southern Oscillation (ENSO) and the variability features in the tropical Indian Ocean. The influence of Pacific air-sea coupled process on the Indian Ocean variability was investigated carefully by conducting numerical experiments. Results show that the occurrence frequency of positive/negative Indian Ocean Dipole (IOD) event will decrease/increase with the presence/absence of the coupled process in the Pacific Ocean. Further analysis demonstrated that the air-sea coupled process in the Pacific Ocean affects the IOD variability mainly by influencing the zonal gradient of thermocline via modulating the background sea surface wind.

Submesoscale activity over the shelf of the northern South China Sea in summer: simulation with an embedded model

We applied a primitive equation ocean model to simulate submesoscale activities and processes over the shelf of the northern South China Sea (NSCS) with a one-way nesting technology for downscaling. The temperature and density fields showed that submesoscale activities were ubiquitous in the NSCS shelf. The vertical velocity was considerably enhanced in submesoscale processes and could reach an average of 58 m per day in the subsurface. At this point, the mixed layer depth also was deepened along the front, and the surface kinetic energy also increased with the intense vertical movement induced by submesoscale activity. Thus, submesoscale stirring/mixing is important for tracers, such as temperature, salinity, nutrients, dissolved organic, and inorganic carbon. This result may have implication for climate and biogeochemical investigations.

Seasonal and interannual variability of primary and export production in the South China Sea: a three-dimensional physical-biogeochemical model study

To investigate the seasonal and interannual variations in biological productivity in the South China Sea (SCS), a Pacific basin-wide physical - biogeochemical model has been developed and used to estimate the biological productivity and export flux in the SCS. The Pacific circulation model, based on the Regional Ocean Model Systems (ROMS), is forced with daily air-sea fluxes derived from the NCEP (National Centers for Environmental Prediction) reanalysis between 1990 and 2004. The biogeochemical processes are simulated with a carbon, Si(OH)(4), and nitrogen ecosystem (CoSiNE) model consisting of silicate, nitrate, ammonium, two phytoplankton groups (small phytoplankton and large phytoplankton), two zooplankton grazers (small micrograzers and large mesozooplankton), and two detritus pools. The ROMS-CoSiNE model favourably reproduces many of the observed features, such as ChI a, nutrients, and primary production (PP) in the SCS. The modelled depth-integrated PP over the euphotic zone (0-125 m) varies seasonally, with the highest value of 386 mg C m (-2) d (-1) during winter and the lowest value of 156 mg C m (-2) d (-1) during early summer. The annual mean value is 196 mg C m (-2) d (-1). The model-integrated annual mean new production (uptake of nitrate), in carbon units, is 64.4 mg C m (-2) d (-1) which yields an f-ratio of 0.33 for the entire SCS. The modelled export ratio (e-ratio: the ratio of export to PP) is 0.24 for the basin-wide SCS. The year-to-year variation of biological productivity in the SCS is weaker than the seasonal variation. The large phytoplankton group tends to dominate over the smaller phytoplankton group, and likely plays an important role in determining the interannual variability of primary and new production.

An individual-based model of the early life history of mackerel (Scomber scombrus) in the eastern North Atlantic, simulating transport, growth and mortality

The main purpose of this paper is to provide the core description of the modelling exercise within the Shelf Edge Advection Mortality And Recruitment (SEAMAR) programme. An individual-based model (IBM) was developed for the prediction of year-to-year survival of the early life-history stages of mackerel (Scomber scombrus) in the eastern North Atlantic. The IBM is one of two components of the model system. The first component is a circulation model to provide physical input data for the IBM. The circulation model is a geographical variant of the HAMburg Shelf Ocean Model (HAMSOM). The second component is the IBM, which is an i-space configuration model in which large numbers of individuals are followed as discrete entities to simulate the transport, growth and mortality of mackerel eggs, larvae and post-larvae. Larval and post-larval growth is modelled as a function of length, temperature and food distribution; mortality is modelled as a function of length and absolute growth rate. Each particle is considered as a super-individual representing 10 super(6) eggs at the outset of the simulation, and then declining according to the mortality function. Simulations were carried out for the years 1998-2000. Results showed concentrations of particles at Porcupine Bank and the adjacent Irish shelf, along the Celtic Sea shelf-edge, and in the southern Bay of Biscay. High survival was observed only at Porcupine and the adjacent shelf areas, and, more patchily, around the coastal margin of Biscay. The low survival along the shelf-edge of the Celtic Sea was due to the consistently low estimates of food availability in that area.

On the nudging terms at open boundaries in regional ocean models

For open boundary conditions (OBCs) in regional models, a nudging term added to radiative and/or advective conditions during the wave or flow propagation outward from the model domain of interest is widely used, to prevent the predicted boundary values from evolving to become quite different from the external data, especially for a long-term integration. However, nudging time scales are basically unknown, leading to many empirical selections. In this paper, a method for objectively estimating nudging time scales during outward propagation is proposed, by using internal model dynamics near the boundary. We tested this method and other several commonly used OBCs for cases of both an idealized model domain and a realistic configuration, and model results demonstrated that the proposed method improves the model solutions. Many similarities are found between the nudging and mixing time scales, in magnitude, spatial and temporal variations, since the nudging mainly replaces the effect of the mixing terms in this study. However, the mixing time scale is not an intrinsic property of the nudging term because in other studies the nudging term might replace terms other than the mixing terms and, thus, should reflect other characteristic time scales.

Assessment of a zoomed global model for the North Sea by comparison with a conventional nested regional model

Assessment of marine downscaling of global model simulations to the regional scale is a prerequisite for understanding ocean feedback to the atmosphere in regional climate downscaling. Major difficulties arise from the coarse grid resolution of global models, which cannot provide sufficiently accurate boundary values for the regional model. In this study, we first setup a stretched global model (MPIOM) to focus on the North Sea by shifting poles. Second, a regional model (HAMSOM) was performed with higher resolution, while the open boundary values were provided by the stretched global model. In general, the sea surface temperatures (SSTs) in the two experiments are similar. Major SST differences are found in coastal regions (root mean square difference of SST is reaching up to 2°C). The higher sea surface salinity in coastal regions in the global model indicates the general limitation of this global model and its configuration (surface layer thickness is 16 m). By comparison, the advantage of the absence of open lateral boundaries in the global model can be demonstrated, in particular for the transition region between the North Sea and Baltic Sea. On long timescales, the North Atlantic Current (NAC) inflow through the northern boundary correlates well between both model simulations (R~0.9). After downscaling with HAMSOM, the NAC inflow through the northern boundary decreases by ~10%, but the circulation in the Skagerrak is stronger in HAMSOM. The circulation patterns of both models are similar in the northern North Sea. The comparison suggests that the stretched global model system is a suitable tool for long-term free climate model simulations, and the only limitations occur in coastal regions. Regarding the regional studies focusing on the coastal zone, nested regional model can be a helpful alternative.

Transports and budgets in a 1/4 ° global ocean reanalysis 1989–2010

Large-scale ocean transports of heat and freshwater have not been well monitored, and yet the regional budgets of these quantities are important to understanding the role of the oceans in climate and climate change. In contrast, atmospheric heat and freshwater transports are commonly assessed from atmospheric reanalysis products, despite the presence of non-conserving data assimilation based on the wealth of distributed atmospheric observations as constraints. The ability to carry out ocean reanalyses globally at eddy-permitting resolutions of 1/4 ° or better, along with new global ocean observation programs, now makes a similar approach viable for the ocean. In this paper we examine the budgets and transports within a global high resolution ocean model constrained by ocean data assimilation, and compare them with independent oceanic and atmospheric estimates.

Atlantic meridional heat transports in two ocean reanalyes evaluated against the RAPID array

The Atlantic meridional overturning circulation in two versions of the NEMO ¼° global ocean model has been compared with the RAPID transport array at 26oN. Both model versions reproduce the mean MOC strength well although the Florida Straits flows differ because the pathway of the Gulf Stream is not strongly constrained at this resolution. Both models however have a mean meridional heat transport of 1.07PW, much lower than the 1.35PW from RAPID observations in Apr04-Oct07. Much of the heat transport discrepancy is due to lower transports in summer across the MidOcean (Bahamas-Africa) section, due to stronger southward geostrophic flows in the top 100m where the water is warmest. Seasonal thermocline changes increase temperature differences across the basin driving stronger geostrophic shear, but this effect is much weaker in the top 100m of the RAPID velocity data. The effect accounts for a reduction of 1.1Sv in MOC and 0.1PW in heat transports. The rest of the discrepancy comes from lower Ekman transports from using ERAInterim winds instead of QuikSCAT, a smaller zonally-varying “Eddy” heat transport component, estimated from repeat XBT sections in the observations, and the southward throughflow in the model. Other differences in depth structure of the model MOC and RAPID observations are described but have much less impact on heat transports.

What processes drive the ocean heat transport?

The ocean contributes to regulating the Earth’s climate through its ability to transport heat from the equator to the poles. In this study we use long simulations of an ocean model to investigate whether the heat transport is carried primarily by wind-driven gyres or whether it is dominated by deep circulations associated with abyssal mixing and high latitude convection. The heat transport is computed as a function of temperature classes. In the Pacific and Indian ocean, the bulk of the heat transport is associated with wind-driven gyres confined to the thermocline. In the Atlantic, the thermocline gyres account for only 40% of the total heat transport. The remaining 60% is associated with a circulation reaching down to cold waters below the thermocline. Using a series of sensitivity experiments, we show that this deep heat transport is primarily set by the strength and patterns of surface winds and only secondarily by diabatic processes at high latitudes in the North Atlantic. Abyssal mixing below 2000 m has hardly any impact on ocean heat transport. A major implication is that the role of the ocean in regulating Earth’s climate strongly depends on how surface winds change across different climates in both hemispheres at low and high latitudes.

Quantifying the roles of ocean circulation and biogeochemistry in governing ocean carbon-13 and atmospheric carbon dioxide at the last glacial maximum

We use a state-of-the-art ocean general circulation and biogeochemistry model to examine the impact of changes in ocean circulation and biogeochemistry in governing the change in ocean carbon-13 and atmospheric CO2 at the last glacial maximum (LGM). We examine 5 different realisations of the ocean's overturning circulation produced by a fully coupled atmosphere-ocean model under LGM forcing and suggested changes in the atmospheric deposition of iron and phytoplankton physiology at the LGM. Measured changes in carbon-13 and carbon-14, as well as a qualitative reconstruction of the change in ocean carbon export are used to evaluate the results. Overall, we find that while a reduction in ocean ventilation at the LGM is necessary to reproduce carbon-13 and carbon-14 observations, this circulation results in a low net sink for atmospheric CO2. In contrast, while biogeochemical processes contribute little to carbon isotopes, we propose that most of the change in atmospheric CO2 was due to such factors. However, the lesser role for circulation means that when all plausible factors are accounted for, most of the necessary CO2 change remains to be explained. This presents a serious challenge to our understanding of the mechanisms behind changes in the global carbon cycle during the geologic past.

The role of air–sea coupling in the simulation of the Madden–Julian oscillation in the Hadley Centre model

The role of air–sea coupling in the simulation of the Madden–Julian oscillation (MJO) is explored using two configurations of the Hadley Centre atmospheric model (AGCM), GA3.0, which differ only in F, a parameter controlling convective entrainment and detrainment. Increasing F considerably improves deficient MJO-like variability in the Indian and Pacific Oceans, but variability in and propagation through the Maritime Continent remains weak. By coupling GA3.0 in the tropical Indo-Pacific to a boundary-layer ocean model, KPP, and employing climatological temperature corrections, well resolved air–sea interactions are simulated with limited alterations to the mean state. At default F, when GA3.0 has a poor MJO, coupling produces a stronger MJO with some eastward propagation, although both aspects remain deficient. These results agree with previous sensitivity studies using AGCMs with poor variability. At higher F, coupling does not affect MJO amplitude but enhances propagation through the Maritime Continent, resulting in an MJO that resembles observations. A sensitivity experiment with coupling in only the Indian Ocean reverses these improvements, suggesting coupling in the Maritime Continent and West Pacific is critical for propagation. We hypothesise that for AGCMs with a poor MJO, coupling provides a “crutch” to artificially augment MJO-like activity through high-frequency SST anomalies. In related experiments, we employ the KPP framework to analyse the impact of air–sea interactions in the fully coupled GA3.0, which at default F shows a similar MJO to uncoupled GA3.0. This is due to compensating effects: an improvement from coupling and a degradation from mean-state errors. Future studies on the role of coupling should carefully separate these effects.