Modeling seasonal circulation, upwelling and tidal mixing in the Arafura and Timor Seas

The extensive shallow tropical seas off northern Australia, encompassing the Arafura and Timor Seas, have been identified as one of the most pristine marine environments on the planet. However, the remoteness and the absence of major industrial development that has contributed to this status have the additional consequence that relatively little is known about these systems. This study is the first to model oceanographic conditions across the tidally dominated Arafura and Timor Seas, and their seasonal variability. The results are based on a high-resolution (0.05°) ocean circulation model forced by realistic winds, waves and tides. The main focus of the study is on physical processes that influence the distributions of sediments and primary productivity across the system. Regions of high bottom stress and tidal mixing have been identified, including a large offshore area around Van Diemen Rise (Timor Sea). Lagrangian particle tracks have revealed a seasonal overturning cell that stretches across the Gulf of Carpentaria (Arafura Sea) with upwelling and downwelling on either side of the Gulf. The presence of coastal upwelling and downwelling is shown to provide a dynamically consistent explanation for the persistent turbid boundary layer observed around the shallow coastal waters of the Gulf.     

Existence of eddies at crossroad of the Indonesian seas

An eddy-resolving Indo-Pacific ocean circulation model was applied to highlight the behavior of eddies throughout the Indonesian seas. The complexity of the topography and coastline at the entrance of the Makassar Strait induce an eddy-type throughflow, instead of a straightforward flow. A sill and a narrow passage in the Makassar strait creates a barrier and impedes the continuation of eddies from the Pacific ocean, but the existence of a steep deep basin (>500 m depth) between the Java and Flores seas indicates a possible area for eddies. Based on our numerical results, we described the presence of a unique eddy structure north of Lombok Island, which we designated the “Lombok Eddy” and verified it by performing a drifter release field experiment and reviewing monthly mean climatology data from the World Ocean Atlas 2001 and XBT PX2 track data. NCEP/NCAR reanalysis, satellite observation data, and mixed layer depth analysis were also used to confirm these processes. By analyzing numerical simulation results and available temperature datasets, two additional eddies were found. All eddies form primarily due to eastward local winds correlated with seasonal monsoon winds during the austral summer. These eddies vary synchronously at an interannual time scale. Together, they are referred to as the Lombok Eddy Train (LET), which affects the surface layer down to a depth of 60 m, and the intensity of the eddy system is strongly affected by mixed layer depth variability from December to February.     

An assessment of the mixed layer salinity budget in the tropical Pacific Ocean. Observations and modelling (1990–2009)

This paper investigates mechanisms controlling the mixed-layer salinity (MLS) in the tropical Pacific during 1990–2009. We use monthly 1°?×?1° gridded observations of salinity, horizontal current and fresh water flux, and a validated ocean general circulation model with no direct MLS relaxation in both its full resolution (0.25° and 5 days) and re-sampled as the observation time/space grid resolution. The present study shows that the mean spatial distribution of MLS results from a subtle balance between surface forcing (E???P, evaporation minus precipitation), horizontal advection (at low and high frequencies) and subsurface forcing (entrainment and mixing), all terms being of analogous importance. Large-scale seasonal MLS variability is found mainly in the Intertropical and South Pacific Convergence Zones due to changes in their meridional location (and related heavy P), in the North Equatorial Counter Currents, and partly in the subsurface forcing. Maximum interannual variability is found in the western Pacific warm pool and in both convergence zones, in relation to El Niño Southern Oscillation (ENSO) events. In the equatorial band, this later variability is due chiefly to the horizontal advection of low salinity waters from the western to the central-eastern basin during El Niño (and vice versa during La Niña), with contrasted evolution for the Eastern and Central Pacific ENSO types. Our findings reveal that all terms of the MLS equation, including high-frequency (<1 month) salinity advection, have to be considered to close the salinity budget, ruling out the use of MLS (or sea surface salinity) only to directly infer the mean, seasonal and/or interannual fresh water fluxes.     

High-altitude data assimilation system experiments for the northern summer mesosphere season of 2007

A global numerical weather prediction system is extended to the mesosphere and lower thermosphere (MLT) and used to assimilate high-altitude satellite measurements of temperature, water vapor and ozone from MLS and SABER during May–July 2007. Assimilated temperature and humidity from 100 to 0.001 hPa show minimal biases compared to satellite data and existing analysis fields. Saturation ratios derived diagnostically from these assimilated temperature and water vapor fields at PMC altitudes and latitudes compare well with seasonal variations in PMC frequency measured from the aeronomy of ice in the mesosphere (AIM) satellite. Synoptic maps of these diagnostic saturation ratios correlate geographically with three independent transient mesospheric cloud events observed at midlatitudes by SHIMMER on STPSat-1 and by ground observers during June 2007. Assimilated temperatures and winds reveal broadly realistic amplitudes of the quasi 5-day wave and migrating tides as a function of latitude and height. For example, analyzed winds capture the dominant semidiurnal MLT wind patterns at 55°N in June 2007 measured independently by a meteor radar. The 5-day wave and migrating diurnal tide also modulate water vapor mixing ratios in the polar summer MLT. Possible origins of this variability are discussed.     

Advanced insights into sources of vertical velocity in the ocean

Estimating vertical velocity in the oceanic upper layers is a key issue for understanding ocean dynamics and the transport of biogeochemical elements. This paper aims to identify the physical sources of vertical velocity associated with sub-mesoscale dynamics (fronts, eddies) and mixed-layer depth (MLD) structures, using (a) an ocean adaptation of the generalized Q-vector form of the ω-equation deduced from a primitive equation system which takes into account the turbulent buoyancy and momentum fluxes and (b) an application of this diagnostic method for an ocean simulation of the Programme Océan Multidisciplinaire Méso Echelle (POMME) field experiment in the North-Eastern Atlantic. The approach indicates that w-sources can play a significant role in the ocean dynamics and strongly depend on the dynamical structure (anticyclonic eddy, front, MLD, etc.). Our results stress the important contribution of the ageostrophic forcing, even under quasi-geostrophic conditions. The turbulent w-forcing was split into two components associated with the spatial variability of (a) the buoyancy and momentum (Ekman pumping) surface fluxes and (b) the MLD. Process (b) represents the trapping of the buoyancy and momentum surface energy into the MLD structure and is identified as an atmosphere/oceanic mixed-layer coupling. The momentum-trapping process is 10 to 100 times stronger than the Ekman pumping and is at least 1,000 times stronger than the buoyancy w-sources. When this decomposition is applied to a filamentary mixed-layer structure simulated during the POMME experiment, we find that the associated vertical velocity is created by trapping the surface wind-stress energy into this structure and not by Ekman pumping.     

Seasonal circulations in the Malay Peninsula Eastern continental shelf from a wave–tide–circulation coupled model

A wave–tide–circulation coupled model based on Princeton Ocean Model is established to study the seasonal circulation in the Malay Peninsula Eastern Continental Shelf region. The model successfully reconstructs the observed seasonal variation of the circulation in the region, as well as the main currents. The simulated tidal harmonic constants, sea surface temperature, and sea surface height anomaly agree with the observations well. The model results show that the upper-layer circulation in the region is mainly controlled by the monsoon winds, while there are two transitions in spring and fall. An anti-cyclonic eddy is present off the Peninsular Malaysia’s east coast in summer, centered at 5°N and 105.5°E, both in the TOPEX/Poseidon data and in the model. Numerical experiments show that the wind stress curl and bathymetry steering are responsible for its formation.     

Impacts of boundary layer parameterization schemes and air-sea coupling on WRF simulation of the East Asian summer monsoon

The planetary boundary layer(PBL)scheme in the regional climate model(RCM)has a significant impact on the interactions and exchanges of moisture,momentum,and energy between land,ocean,and atmosphere;however,its uncertainty will cause large systematic biases of RCM.Based on the four different PBL schemes(YSU,ACM2,Boulac,and MYJ)in Weather Research and Forecasting(WRF)model,the impacts of these schemes on the simulation of circulation and precipitation during the East Asian summer monsoon(EASM)are investigated.The simulated results of the two local turbulent kinetic energy(TKE)schemes,Boulac and MYJ,are more consistent with the observations than those in the two nonlocal closure schemes,YSU and ACM2.The former simulate more reasonable low-level southwesterly flow over East China and west pacific subtropical high(WPSH)than the latter.As to the modeling of summer monsoon precipitation,both the spatial distributions and temporal evolutions from Boulac and MYJ are also better than those in YSU and ACM2 schemes.In addition,through the comparison between YSU and Boulac experiments,the differences from the results of EASM simulation are more obvious over the oceanic area.In the experiments with the nonlocal schemes YSU and ACM2,the boundary layer mixing processes are much stronger,which lead to produce more sea surface latent heat flux and enhanced convection,and finally induce the overestimated precipitation and corresponding deviation of monsoon circulation.With the further study,it is found that the absence of air-sea interaction in WRF may amplify the biases caused by PBL scheme over the ocean.Consequently,there is a reduced latent heat flux over the sea surface and even more reasonable EASM simulation,if an ocean model coupled into WRF.     

Significance of the normalized seasonality of wind field and its rationality for characterizing the monsoon

The significance of the normalized seasonality of wind field is theoretically studied and the intrinsic relationship between its critical value and the definition of the monsoon region is revealed. As a result, the rationality which the monsoon region is characterized by the significant normalized seasonality is proved. Based on this, by use of the NECP/NCAR reanalysis data during 1958–1997, the spatial structure of the significant normalized seasonality of wind field is investigated, and the significant baroclinity of the seasonal variation of the atmospheric general circulation in the vertical direction is verified. Moreover, in the tropics there exists an anti-symmetric distribution between the significant seasonality in the eastern hemisphere and that in the western hemisphere, and the monsoon is linked closely to some important systems such as the subtropical highs, the night jet and the west wind channels.     

A climatic dataset of ocean vertical turbulent mixing coefficient based on real energy sources

Using data on wind stress, significant height of combined wind waves and swell, potential temperature, salinity and seawater velocity, as well as objectively-analyzed in situ temperature and salinity, we established a global ocean dataset of calculated wind- and tide-induced vertical turbulent mixing coefficients. We then examined energy conservation of ocean vertical mixing from the point of view of ocean wind energy inputs, gravitational potential energy change due to mixing (with and without artificially limiting themixing coefficient), and K-theory vertical turbulent parameterization schemes regardless of energy inputs. Our research showed that calculating the mixing coefficient with average data and artificial limiting the mixing coefficient can cause a remarkable lack of energy conservation, with energy losses of up to 90% and changes in the energy oscillation period. The data also show that wind can introduce a huge amount of energy into the upper layers of the Southern Ocean, and that tidesdo so in regions around underwater mountains. We argue that it is necessary to take wind and tidal energy inputs into account forlong-term ocean climate numerical simulations. We believe that using this ocean vertical turbulent mixing coefficient climatic dataset is a fast and efficient method to maintain the ocean energy balance in ocean modeling research.     

Internal-wave-induced and double-diffusive nutrient fluxes to the nutrient-consuming surface layer in the oligotrophic subtropical North Atlantic

In the literature, an inconsistency exists between estimates of biotically-effected carbon export inferred from large-scale geochemical studies (Jenkins 1982; 47 gC m^–2 a^–1) and local measurements of turbulent nutrient supply (Lewis et al. 1986; 4 gC m^–2 a^–1) in the eastern subtropical North Atlantic. Nutrient supply to the upper ocean by turbulent mixing is reexamined using local standard oceanographic measurements and high-resolution vertical profiles of nutrients averaged over a large region directly comparable to that investigated by Jenkins (1982).Responsible Editor: Franciscus Colijn     

Ocean circulation on the North Australian Shelf

The ocean circulation on Australia's Northern Shelf is dominated by the Monsoon and influenced by large-scale interannual variability. These driving forces exert an ocean circulation that influences the deep Timor Sea Passage of the Indonesian Throughflow, the circulation on the Timor and Arafura Shelves and, further downstream, the Leeuwin Current. Seasonal maxima of northeastward (southwestward) volume transports on the shelf are almost symmetric and exceed 10⁶ m³/s in February (June). The associated seasonal cycle of vertical upwelling from June to August south of 8.5°S and between 124°E and 137.5°E exceeds 1.5×10⁶ m³/s across 40 m depth. During El Niño events, combined anomalies from the seasonal means of high regional wind stresses and low inter-ocean pressure gradients double the northeastward volume transport on the North Australian Shelf to 1.5×10⁶ m³/s which accounts for 20% of the total depth-integrated transport across 124°E and reduce the total transport of the Indonesian Throughflow. Variability of heat content on the shelf is largely determined by Pacific and Indian Ocean equatorial wind stress anomalies with some contribution from local wind stress forcing.     

SST biases over the Northwest Pacific and possible causes in CMIP5 models

In this paper, the features and possible causes of sea surface temperature(SST) biases over the Northwest Pacific are investigated based on a mixed-layer heat budget analysis in 21 coupled general circulation models(CGCMs) from phase 5 of the Coupled Model Inter-comparison Project(CMIP5). Most CMIP5 models show cold SST biases throughout the year over the Northwest Pacific. The largest biases appear during summer, and the smallest biases occur during winter. These cold SST biases are seen at the basin scale and are mainly located in the inner region of the low and mid-latitudes. According to the mixed-layer heat budget analysis, overestimation of upward net sea surface heat fluxes associated with atmospheric processes are primarily responsible for the cold SST biases. Among the different components of surface heat fluxes, overestimated upward latent heat fluxes induced by the excessively strong surface winds contribute the most to the cold SST biases during the spring, autumn, and winter seasons. Conversely, during the summer, overestimated upward latent heat fluxes and underestimated downward solar radiations at the sea surface are equally important. Further analysis suggests that the overly strong surface winds over the Northwest Pacific during winter and spring are associated with excessive precipitation over the Maritime Continent region,whereas those occurring during summer and autumn are associated with the excessive northward extension of the intertropical convergence zone(ITCZ). The excessive precipitation over the Maritime Continent region and the biases in the simulated ITCZ induce anomalous northeasterlies, which are in favor of enhancing low-level winds over the North Pacific. The enhanced surface wind increases the sea surface evaporation, which contributes to the excessive upward latent heat fluxes. Thus, the SST over the Northwest Pacific cools.     

Seasonal upwelling on the Western and Southern Shelves of the Gulf of Mexico

An 8-year database of sea surface temperature (SST), 7 years of Sea-viewing Wide Field-of-view Sensor (SeaWiFS) ocean color images, wind fields, and numerical model results are analyzed to identify regions and periods of coastal upwelling on the western and southern shelves of the Gulf of Mexico. On the seasonal scale, it is found that on the Tamaulipas, Veracruz, and southwestern Texas–Louisiana shelves there are upwelling favorable winds from April to August, when southeasterly winds are dominant and cold SST anomalies associated with upwelling are observed along their coasts. However, during summer, values of chlorophyll-a concentration are lower than those in autumn and winter, which are high due to advection of old bloom biological material from upstream. During winter, there is a cold front on the Tamaulipas shelf produced by advection of cold water from the Texas–Louisiana shelf and not due to upwelling. On the eastern Campeche Bank, persistent upwelling is observed due to favorable winds throughout the year with cold SST and large chlorophyll-a content along the inner shelf from May to September. On the Tamaulipas shelf, the summer upwelling delays the annual SST peak until September, while in most of the Gulf SST peaks in August. This difference is due to the end of the upwelling favorable wind conditions and the September seasonal current reversal.     

A numerical study of the dynamics of the wave-driven circulation within a fringing reef system

The circulation driven by wave breaking, tides and winds within a fringing coral reef system (Ningaloo Reef) in Western Australia was investigated using the ocean circulation model ROMS two-way coupled to the wave model SWAN. Currents within the system were dominantly forced by wave breaking, with flow driven over the shallow reefs and towards the lagoon, which returned to the ocean through channels in the reef. Hindcast model simulations were compared against an extensive field dataset, revealing that the coupled wave–circulation model could accurately predict the waves and currents throughout this morphologically complex reef–lagoon system. A detailed momentum budget analysis showed that, over the reef, a dominant cross-shore balance was established between radiation stress gradients and a pressure (mean water level) gradient (similar to a beach). Within the lagoon, alongshore currents were primarily balanced by alongshore gradients in wave setup, which drove flow towards (and ultimately out) the channels. The importance of these wave-driven currents to Ningaloo Reef was quantified over a full seasonal cycle, during periods when wave and wind conditions significantly differed. These results showed that wave breaking still overwhelmingly dominated the circulation and flushing of Ningaloo Reef throughout the year, with winds playing an insignificant role.     

On the response of a turbulent coastal buoyant current to wind events: the case of the Western Adriatic Current

This numerical study focuses on the response of the Western Adriatic Current to wind forcing. The turbulent buoyant surface current is induced by the Po river outflow in the Adriatic Sea. Idealized and realistic wind conditions are considered by retaining the complex geomorphology of the middle Adriatic basin. In the absence of wind, the Adriatic Promontories force the current to separate from the coast and induce instabilities. Persistent 7-m s^− 1 downwelling favorable northwesterly winds thicken and narrow the current. Instabilities whose size is ~10 km develop but ultimately vanish, since there is not enough across-shore space to grow. On the contrary, 7-m s^− 1 upwelling favorable southeasterly winds thin and widen the current, and instabilities can grow to form mesoscale (~35 km) features. When realistic winds are considered, the same trends are observed, but the state of the sea set up by previous wind events also plays a crucial role. The turbulent regimes set up by different winds affect mixing and the WAC meridional transport. With downwelling winds, the transport is generally southward and mixing happens mostly between the fresher (S ≤ 38) salinity classes. With upwelling winds, the transport decreases and changes sign, and mixing mainly involves saltier (S > 38) waters. In all cases, mixing is enhanced when a finer 0.5-km horizontal resolution is employed.     

Atmospheric tides and mean winds in the meteor region over Santa Maria (29.7°S; 53.8°W)

We compare meteor radar measurements of the MLT region winds at Santa Maria, Brazil (29.7°S, 53.8°W) with the Horizontal Neutral Wind Model (HWM-93) and the Global Scale Wave Model (GSWM-00). The observed annual variation of the prevailing zonal wind disagrees in some respects with the HWM-93 model. Also, the zonal diurnal tide amplitude shows an annual variation, whereas that of the GSWM-00 is semiannual, and its vertical wavelength is smaller than that suggested by the model. The observed semidiurnal tide shows seasonal and inter-annual variations and the phase is evanescent during almost the whole year.     

Late Quaternary clay minerals off Middle Vietnam in the western South China Sea: Implications for source analysis and East Asian monsoon evolution

High-resolution clay mineral records combined with oxygen isotopic stratigraphy over the past 450 ka during late Quaternary from Core MD05-2901 off Middle Vietnam in the western South China Sea are re-ported to reconstruct a history of East Asian monsoon evolution. Variations in Illite, chlorite, and kaolinite contents indicate a strong glacial-interglacial cyclicity, while changes in smectite content present a higher frequency cyclicity. The provenance analysis indicates a mixture of individual clay minerals from various sources surrounding the South China Sea. Smectite derived mainly from the Sunda shelf and its major source area of the Indonesian islands. Illite and chlorite originated mainly from the Mekong and Red rivers. Kaolinite was provided mainly by the Pearl River. Spectral analysis of the kaolin-ite/(illite chlorite) ratio displays a strong eccentricity period of 100 ka, implying the ice sheet-forced win-ter monsoon evolution; whereas higher frequency changes in the smectite content show an ice sheet-forced obliquity period of 41 ka, and precession periods of 23 and 19 ka and a semi-precession period of 13 ka as well, implying the tropical-forced summer monsoon evolution. The winter monsoon evolution is generally in coherence with the glacial-interglacial cyclicity, with intensified winter monsoon winds during glacials and weakened winter monsoon winds during interglacials; whereas the summer monsoon evolution provides an almost linear response to the summer insolation of low latitude in the Northern Hemisphere, with strengthened summer monsoon during higher insolation and weakened summer monsoon during lower insolation. The result suggests that the high-latitude ice sheet and low-latitude tropical factor could drive the late Quaternary evolution of East Asian winter and summer monsoons, respectively, implying their diplex and self-contained forcing mechanism.     

Observations of the upper ocean response to storm forcing in the South Atlantic Roaring Forties

In the austral summer of 1992–1993 the passage of a storm system drove a strong upper ocean response at 45°S in the mid-South Atlantic. Good in situ observations were obtained. CTD casts revealed that the mixed layer deepened by \sim40 m over 4 days. Wind stirring dominated over buoyancy flux-driven mixing during the onset of high winds. Doppler shear currents further reveal this to be intimately related to inertial dynamics. The penetration depth of inertial currents, which are confined to the mixed layer, increases with time after a wind event, matched by a downward propagation of low values of the Richardson number. This suggests that inertial current shear is instrumental in producing turbulence at the base of the mixed layer. Evolution of inertial transport is simulated using a time series of ship-observed wind stress. Simulated transport is only 30-50% of the observed transport, suggesting that much of the observed inertial motion was forced by an earlier (possibly remote) storm. Close proximity of the subtropical front further complicates the upper ocean response to the storm. A simple heat balance for the upper 100 m reveals that surface cooling and mixing (during the storm) can account for only a small fraction of an apparent \sim1 °C mixed layer cooling.     

Long-Lead Prediction Skill of Indian Summer Monsoon Rainfall Using Outgoing Longwave Radiation (OLR): an Application of Canonical Correlation Analysis

There is a statistical linkage between tropical outgoing longwave radiation (OLR) and all-India summer monsoon rainfall (AISMR). A positive and significant correlation is observed over the surrounding areas of northeast Australia and the Arafura Sea during the months of January and February (J&F) which drops down as the lead month decreases. The OLR index as an area average over the surrounding areas of northeast Australia and the Arafura Sea is found to have 0.4 correlation with AISMR. The index is also found to be strongly correlated with the Indian monsoon index. In view of the teleconnection pattern, the OLR index is used for the development of statistical models using the concept of linear regression (LR) and canonical correlation analysis (CCA). Potential of CCA over LR is observed for the prediction of seasonal rainfall over the northwest, west central and over the whole country as well. The seasonal rainfall predictability basically comes from the months of June and September.