Influence of Pacific on Southern Indian Ocean Rossby Waves

The sea-surface height anomalies derived from Simple Ocean Data Assimilation (SODA) during 1958–2001, Topex/Poseidon satellite during 1993–2001 and the SODA heat content anomalies (125 m depth) during 1958–2001 are filtered into annual and biennial Rossby wave components using a two-dimensional Finite Impulse Response filter. The filtered Rossby wave components (both annual and biennial) in the southern Pacific and Indian Oceans have considerable strength and variability. The propagation of annual and biennial Rossby waves in the Indonesian through-flow region [12.5°S–7.5°S] of the Indian Ocean is in phase with the southern Pacific Ocean waves. So it is speculated that the Pacific Ocean influences the Indian Ocean, especially through the region 17.5°S to 7.5°S and thus the southern Pacific Rossby waves could be an unexplored contributor to the Indian Ocean Rossby waves. We also carried out Fast Fourier Transform (FFT) and wavelet analysis on the tide gauge sea-level data along the Australian coast to support our claim. Filtered annual and biennial components of SODA heat content anomalies (125 m depth) also support these findings.     

A finite element sea ice model of the Canadian Arctic Archipelago

The Canadian Arctic Archipelago (CAA) is a complex area formed by narrow straits and islands in the Arctic. It is an important pathway for freshwater and sea-ice transport from the Arctic Ocean to the Labrador Sea and ultimately to the Atlantic Ocean. The narrow straits are often crudely represented in coupled sea-ice–ocean models, leading to a misrepresentation of transports through these straits. Unstructured meshes are an alternative in modelling this complex region, since they are able to capture the complex geometry of the CAA. This provides higher resolution in the flow field and allows for more accurate transports (but not necessarily better modelling). In this paper, a finite element sea-ice model of the Arctic region is described and used to estimate the sea-ice fluxes through the CAA. The model is a dynamic–thermodynamic sea-ice model with elastic–viscous–plastic rheology and is coupled to a slab ocean, where the temperature and salinity are restored to climatology, with no velocities and surface elevation. The model is spun-up from 1973 to 1978 with NCEP/NARR reanalysis data. From 1979 to 2007, the model is forced by NCEP/DoE reanalysis data. The large scale sea-ice characteristics show good agreement with observations. The total sea-ice area agrees very well with observations and shows a sensitivity to the Arctic oscillation (AO). For 1998–2002, we find estimates for the sea-ice volume and area fluxes through Admunsen Gulf, McClure Strait and the Queen Elizabeth Islands that compare well with observation and are slightly better than estimates from other models. For Nares Strait, we find that the fluxes are much lower than observed, due to the missing effect of topographic steering on the atmospheric forcing fields. The 1979–2007 fluxes show large seasonal and interannual variability driven primarily by variability in the ice velocity field and a sensitivity to the AO and other large-scale atmospheric variability, which suggests that accurate atmospheric forcing might be crucial to modelling the CAA.     

Global upper ocean heat content and climate variability

Observational data from the Global Temperature and Salinity Profile Program were used to calculate the upper ocean heat content (OHC) anomaly. The thickness of the upper layer is taken as 300 m for the Pacific/Atlantic Ocean and 150 m for the Indian Ocean since the Indian Ocean has shallower thermoclines. First, the optimal spectral decomposition scheme was used to build up monthly synoptic temperature and salinity dataset for January 1990 to December 2009 on 1° × 1° grids and the same 33 vertical levels as the World Ocean Atlas. Then, the monthly varying upper layer OHC field (H) was obtained. Second, a composite analysis was conducted to obtain the total-time mean OHC field ([`([`(H)])] \bar{\bar{H}} ) and the monthly mean OHC variability ( [(\textH)\tilde] \widetilde{\text{H}} ), which is found an order of magnitude smaller than [^(\textH)] \widehat{\text{H}} . Third, an empirical orthogonal function (EOF) method is conducted on the residue data ( [^(\textH)] \widehat{\text{H}} ), deviating from [(\textH)\tilde] \widetilde{\text{H}}  +  [(\textH)\tilde] \widetilde{\text{H}} , in order to obtain interannual variations of the OHC fields for the three oceans. In the Pacific Ocean, the first two EOF modes account for 51.46% and 13.71% of the variance, representing canonical El Nino/La Nina (EOF-1) and pseudo-El Nino/La Nina (i.e., El Nino Modoki; EOF-2) events. In the Indian Ocean, the first two EOF modes account for 24.27% and 20.94% of the variance, representing basin-scale cooling/warming (EOF-1) and Indian Ocean Dipole (EOF-2) events. In the Atlantic Ocean, the first EOF mode accounts for 49.26% of the variance, representing a basin-scale cooling/warming (EOF-1) event. The second EOF mode accounts for 8.83% of the variance. Different from the Pacific and Indian Oceans, there is no zonal dipole mode in the tropical Atlantic Ocean. Fourth, evident lag correlation coefficients are found between the first principal component of the Pacific Ocean and the Southern Oscillation Index with a maximum correlation coefficient (0.68) at 1-month lead of the EOF-1 and between the second principal component of the Indian Ocean and the Dipole Mode Index with maximum values (around 0.53) at 1–2-month advance of the EOF-2. It implies that OHC anomaly contains climate variability signals.     

Decadal scale variability of sea surface temperature in the Mediterranean Sea in relation to atmospheric variability

Twenty-four years of AVHRR-derived sea surface temperature (SST) data (1985–2008) and 35 years of NOCS (V.2) in situ-based SST data (1973–2008) were used to investigate the decadal scale variability of this parameter in the Mediterranean Sea in relation to local air–sea interaction and large-scale atmospheric variability. Satellite and in situ-derived data indicate a strong eastward increasing sea surface warming trend from the early 1990s onwards. The satellite-derived mean annual warming rate is about 0.037°C year^–1 for the whole basin, about 0.026°C year^–1 for the western sub-basin and about 0.042°C year^–1 for the eastern sub-basin over 1985–2008. NOCS-derived data indicate similar variability but with lower warming trends for both sub-basins over the same period. The long-term Mediterranean SST spatiotemporal variability is mainly associated with horizontal heat advection variations and an increasing warming of the Atlantic inflow. Analysis of SST and net heat flux inter-annual variations indicates a negative correlation, with the long-term SST increase, driving a net air–sea heat flux decrease in the Mediterranean Sea through a large increase in the latent heat loss. Empirical orthogonal function (EOF) analysis of the monthly average anomaly satellite-derived time series showed that the first EOF mode is associated with a long-term warming trend throughout the whole Mediterranean surface and it is highly correlated with both the Eastern Atlantic (EA) pattern and the Atlantic Multidecadal Oscillation (AMO) index. On the other hand, SST basin-average yearly anomaly and NAO variations show low and not statistically significant correlations of opposite sign for the eastern (negative correlation) and western (positive correlation) sub-basins. However, there seems to be a link between NAO and SST decadal-scale variations that is particularly evidenced in the second EOF mode of SST anomalies. NOCS SST time series show a significant SST rise in the western basin from 1973 to the late 1980s following a large warming of the inflowing surface Atlantic waters and a long-term increase of the NAO index, whereas SST slowly increased in the eastern basin. In the early 1990s, there is an abrupt change from a very high positive to a low NAO phase which coincides with a large change in the SST spatiotemporal variability pattern. This pronounced variability shift is followed by an acceleration of the warming rate in the Mediterranean Sea and a change in the direction (from westward to eastward) of its spatial increasing tendency.     

Sensitivity of climate response to variations in freshwater hosing location

In a recent intercomparison of the response of general circulation models (GCMs) to high-latitude freshwater forcing (Stouffer et al., J Climate 19(8):1365–1387, 2006), a number of the GCMs investigated showed a localised warming response in the high-latitude North Atlantic, as opposed to the cooling that the other models showed. We investigated the causes for this warming by testing the sensitivity of the meridional overturning circulation (MOC) to variations in freshwater forcing location, and then analysing in detail the causes of the warming. By analysing results from experiments with HadCM3, we are able to show that the high-latitude warming is independent of the exact location of the additional freshwater in the North Atlantic or Arctic Ocean basin. Instead, the addition of freshwater changes the circulation in the sub-polar gyre, which leads to enhanced advection of warm, saline, sub-surface water into the Greenland–Iceland–Norwegian Sea despite the overall slowdown of the MOC. This sub-surface water is brought to the surface by convection, where it leads to a strong warming of the surface waters and the overlying atmosphere.     

Modification of water masses in the Barents Sea and its coupling to ice dynamics: a model study

The physical and biological environment of the Barents Sea is characterised by large variability on a wide range of scales. Results from a numerical ocean model, SINMOD, are presented showing that the physical variability is partly forced by changes in annual net ice import. The mean contribution from ice import in the simulation period (1979–2007) is about 40% of the total amount of ice melted each year. The annual ice import into the Barents Sea varies between 143 and 1,236 km³, and this causes a substantial variability in the amount of annual ice melt in the Barents Sea. This in turn impacts the freshwater content. The simulated freshwater contribution from ice is 0.02 Sv on average and 0.04 Sv at maximum. When mixed into a mean net Atlantic Water (AW) inflow of 1.1 Sv with a salinity of 35.1, this freshwater addition decreases the salinity of the modified AW to 34.4 and 33.9 for the mean and maximum freshwater fluxes, respectively. Ice import may thus be important for the Barents Sea production of Arctic Ocean halocline water which has salinity of about 34.5. The changes in the ice melt the following summer due to ice import also affect the formation of dense water in the Barents Sea by changing stratification, altering the vertical mixing rates and affecting heat loss from the warm AW. The model results thus indicate that ice import from the Arctic has a great impact on water mass modification in the Barents Sea which in turn impacts the ventilation of the Arctic Ocean.     

Modeling studies of the upper ocean response to a tropical cyclone

A coupled ocean and boundary layer flux numerical modeling system is used to study the upper ocean response to surface heat and momentum fluxes associated with a major hurricane, namely, Hurricane Dennis (July 2005) in the Gulf of Mexico. A suite of experiments is run using this modeling system, constructed by coupling a Navy Coastal Ocean Model simulation of the Gulf of Mexico to an atmospheric flux model. The modeling system is forced by wind fields produced from satellite scatterometer and atmospheric model wind data, and by numerical weather prediction air temperature data. The experiments are initialized from a data assimilative hindcast model run and then forced by surface fluxes with no assimilation for the time during which Hurricane Dennis impacted the region. Four experiments are run to aid in the analysis: one is forced by heat and momentum fluxes, one by only momentum fluxes, one by only heat fluxes, and one with no surface forcing. An equation describing the change in the upper ocean hurricane heat potential due to the storm is developed. Analysis of the model results show that surface heat fluxes are primarily responsible for widespread reduction (0.5°–1.5°C) of sea surface temperature over the inner West Florida Shelf 100–300 km away from the storm center. Momentum fluxes are responsible for stronger surface cooling (2°C) near the center of the storm. The upper ocean heat loss near the storm center of more than 200 MJ/m² is primarily due to the vertical flux of thermal energy between the surface layer and deep ocean. Heat loss to the atmosphere during the storm’s passage is approximately 100–150 MJ/m². The upper ocean cooling is enhanced where the preexisting mixed layer is shallow, e.g., within a cyclonic circulation feature, although the heat flux to the atmosphere in these locations is markedly reduced.     

Warming in Arctic intermediate and deep waters around Chukchi Plateau and its adjacent regions in 1999

Warming in the Arctic Intermediate Water (AIW) has been reported in a series of articles in recent years. Prior to 1992, the water temperature of AIW off the Chukchi Continental Slope had never exceeded 0.5℃. Aagaard et al.[1] and Carmack[2] were the first to report that the temperature in AIW in the area in 1993 was close to 1℃, based on the data observed in a Canadian expedition. In 1994, the temperature of AIW around the Chukchi Sea and Mendeleyev Ridge again rose by at least an…     

Regional patterns of observed sea level change: insights from a 1/4° global ocean/sea-ice hindcast

A global eddy-admitting ocean/sea-ice simulation driven over 1958–2004 by daily atmospheric forcing is used to evaluate spatial patterns of sea level change between 1993 and 2001. In the present study, no data assimilation is performed. The model is based on the Nucleus for European Models of the Ocean code at the 1/4° resolution, and the simulation was performed without data assimilation by the DRAKKAR project. We show that this simulation correctly reproduces the observed regional sea level trend patterns computed using satellite altimetry data over 1993–2001. Generally, we find that regional sea level change is best simulated in the tropical band and northern oceans, whereas the Southern Ocean is poorly simulated. We examine the respective contributions of steric and bottom pressure changes to the total regional sea level changes. For the steric component, we analyze separately the contributions of temperature and salinity changes as well as upper and lower ocean contributions. Generally, the model results show that most regional sea level changes arise from temperature changes in the upper 750 m of the ocean. However, contributions of salinity changes and deep steric changes can be locally important. We also propose a map of ocean bottom pressure changes. Finally, we assess the robustness of such a model by comparing this simulation with a second simulation performed by MERCATOR-Ocean based on the same core model, but differing by its short length of integration (1992–2001) and its surface forcing data set. The long simulation presents better performance over 1993–2001 than the short simulation, especially in the Southern Ocean where a long adjustment time seems to be needed. In memory of my little brother Jean-Eudes, whose thirst for science filled out the rich discussions we had about my investigations and his job as user-service provider for MERCATOR-Ocean.     

Oceanic environmental changes of subarctic Bering Sea in recent 100 years: Evidence from molecular fossils

The objectives of the past studies of global cli-matological changes were to obtain and interpret in-formation about various macroclimates and paleoenvi-ronments, and to understand changing processes andmechanisms of the global system environment, so asto provide basic materials for establishing a physicalmodel of forecasting climate environment. To reduceuncertainty of forecast, it is particularly important totake variation sequences in short dimension and highresolution in global climatologic…     

Interannual variability in the Biannual Rossby waves in the tropical Indian Ocean and its relation to Indian Ocean Dipole and El Nino forcing

The interannual variability of the tropical Indian Ocean is studied using Simple Ocean Data Assimilation (SODA) sea surface height anomalies (SSHA) and Hadley Centre Ice Sea Surface Temperature anomalies. Biannual Rossby waves (BRW) were observed along the 1.5° S and 10.5° S latitudes during the Indian Ocean Dipole (IOD) years. The SODA SSHA and its BRW components were comparable with those of Topex/Poseidon. The phase speed of BRW along 1.5° S is −28 cm/s, which is comparable with the theoretical speed of first mode baroclinic (equatorially trapped) Rossby waves. This is the first study to show that no such propagation is seen along 1.5° S during El Nino years in the absence of IOD. Thus the westward propagating downwelling BRW in the equatorial Indian Ocean is hypothesized as a potential predictor for IOD. These waves transport heat from the eastern equatorial Indian Ocean to west, long before the dipole formation. Along 10.5° S, the BRW formation mechanisms during the El Nino and IOD years were found to be different. The eastern boundary variations along 10.5° S, being localized, do not influence the ocean interior considerably. Major portion of the interannual variability of the thermocline, is caused by the Ekman pumping integrated along the characteristic lines of Rossby waves. The study provides evidence of internal dynamics in the IOD formation. The positive trend in the downwelling BRW (both in SODA and Topex/Poseidon) is of great concern, as it contributes to the Indian Ocean warming.     

Effects of Arctic Sea Ice Decline on Weather and Climate: A Review

The areal extent, concentration and thickness of sea ice in the Arctic Ocean and adjacent seas have strongly decreased during the recent decades, but cold, snow-rich winters have been common over mid-latitude land areas since 2005. A review is presented on studies addressing the local and remote effects of the sea ice decline on weather and climate. It is evident that the reduction in sea ice cover has increased the heat flux from the ocean to atmosphere in autumn and early winter. This has locally increased air temperature, moisture, and cloud cover and reduced the static stability in the lower troposphere. Several studies based on observations, atmospheric reanalyses, and model experiments suggest that the sea ice decline, together with increased snow cover in Eurasia, favours circulation patterns resembling the negative phase of the North Atlantic Oscillation and Arctic Oscillation. The suggested large-scale pressure patterns include a high over Eurasia, which favours cold winters in Europe and northeastern Eurasia. A high over the western and a low over the eastern North America have also been suggested, favouring advection of Arctic air masses to North America. Mid-latitude winter weather is, however, affected by several other factors, which generate a large inter-annual variability and often mask the effects of sea ice decline. In addition, the small sample of years with a large sea ice loss makes it difficult to distinguish the effects directly attributable to sea ice conditions. Several studies suggest that, with advancing global warming, cold winters in mid-latitude continents will no longer be common during the second half of the twenty-first century. Recent studies have also suggested causal links between the sea ice decline and summer precipitation in Europe, the Mediterranean, and East Asia.     

Structure of summer atmospheric boundary layer in the center of Arctic Ocean and its relation with sea ice extent change

The atmospheric vertical structure and changed characteristics of boundary layer parameters, as well as their relations with sea ice and temperature changes in the center of Arctic Ocean(80°–88°N) are presented by adopting GPS sounding data obtained by the 4th–6th Arctic expeditions of China and NCEP(National Centre for Environmental Prediction) reanalysis data. Obvious differences are observed regarding the tropopause, boundary layer height, temperature inversion, and vertical structure of wind speed and direction in the center Arctic Ocean in the summer of 2012, 2010, and 2014. These differences can be explained by the relations between temperature and changes in sea ice extent in September from 1979 to 2014. In September 2012, the Arctic sea ice extent decreased by 44% an with obvious warming process. In September 2010 and 2014, it decreased by 22.6% and 17% with an obvious cooling process, respectively. A comparison of the two processes shows that sea ice change has a significant influence on the structure of the atmospheric boundary layer. In the recent 30 years, the temperature changes of 1000 and 850 h Pa in the center of the Arctic Ocean have displayed an obvious warming trend and negative correlation with sea ice extent. These changes indicate that the continuous reduction of Arctic sea ice will continue the warming of the troposphere middle layer.     

Temporal Variability of Lower Stratosphere Temperature

Inter-monthly to inter-decadal global variability of lower stratosphere temperature (LST) is studied in order to improve current knowledge on its variability and trends, as well as natural and anthropogenic influences upon it. Principal Component Analysis (PCA) with S-mode Varimax rotated PCA were used. The first seven components, which explain 70% of variance make it possible to determine homogeneous LST behaviour zones with little overlap between areas, and practically no unclassified areas. Composite time series, referred to as reference series, in the core of the subregions defined by each of the PCs, were calculated in order to obtain the temporal patterns. The equatorial-tropical zone and the subtropical area display warmings caused by the eruptions of El Chichon and Mt. Pinatubo volcanoes as well as the strong influence of the Quasi-Biennial Oscillation (QBO) which leads to equatorial warming (cooling) in the west (east) phase and cooling (warming) in subtropical latitudes. Only low latitudes show some kind of global teleconnection between hemispheres. Significant correlation with several ocean/atmosphere index time-series like ENSO, Antarctic and Arctic Oscillations (AAO, AO), Arctic Circumpolar Vortex was detected over latitudinally separate regions. Antarctic and Arctic ozone hole values were contrasted with warming and cooling features registered in mid and high latitudes in both hemispheres. The LST reference series exhibit a negative trend, commonly attributed to the increase in greenhouse gases that lead to a warming of the troposphere and a cooling of the stratosphere, in all sub regions. The highest cooling rate of − 0.65 °C/ decade is detected in the Gobi desert, and the lowest values of −0.1 °C/ decade over the NE of Canada and Greenland which indicates the great longitudinal variability that the LST trends may present. The difference with other authors is mainly due to the fact that results are based either on latitudinal averages or radiosonde data.     

Basin patterns of upper ocean warming for 1993–2009

A previous study (Lyman et al., Nature 465:334–337, 2010) showed a robust warming signal of the global upper ocean (0–700 m). They examined several sources of uncertainty that contribute to differences among heat content estimations. However, their focus was limited to globally averaged estimation. This study presents the spatial pattern of the global heat content change based on observed gridded datasets (Levitus et al., Geophys Res Lett 36:L07608, 2009). The western Pacific, Atlantic, and Indian Oceans showed significant warming trends, whereas eastern Pacific and some areas of the Gulf Stream experienced negative trends during 1993–2009. Steady warming trend was obtained from the first EOF mode when El Nino and Southern Oscillation (ENSO)-related signals were removed. This result implies that the rapid increase in heat content of the upper ocean around 2000–2005 is not related to a sampling transition from XBT to Argo observations but is associated with a natural variability dominated by strong ENSO-related signals.     

Long-term ocean simulations in FESOM: evaluation and application in studying the impact of Greenland Ice Sheet melting

The Finite Element Sea-ice Ocean Model (FESOM) is formulated on unstructured meshes and offers geometrical flexibility which is difficult to achieve on traditional structured grids. In this work, the performance of FESOM in the North Atlantic and Arctic Ocean on large time scales is evaluated in a hindcast experiment. A water-hosing experiment is also conducted to study the model sensitivity to increased freshwater input from Greenland Ice Sheet (GrIS) melting in a 0.1-Sv discharge rate scenario. The variability of the Atlantic Meridional Overturning Circulation (AMOC) in the hindcast experiment can be explained by the variability of the thermohaline forcing over deep convection sites. The model also reproduces realistic freshwater content variability and sea ice extent in the Arctic Ocean. The anomalous freshwater in the water-hosing experiment leads to significant changes in the ocean circulation and local dynamical sea level (DSL). The most pronounced DSL rise is in the northwest North Atlantic as shown in previous studies, and also in the Arctic Ocean. The released GrIS freshwater mainly remains in the North Atlantic, Arctic Ocean and the west South Atlantic after 120 model years. The pattern of ocean freshening is similar to that of the GrIS water distribution, but changes in ocean circulation also contribute to the ocean salinity change. The changes in Arctic and sub-Arctic sea level modify exchanges between the Arctic Ocean and subpolar seas, and hence the role of the Arctic Ocean in the global climate. Not only the strength of the AMOC, but also the strength of its decadal variability is notably reduced by the anomalous freshwater input. A comparison of FESOM with results from previous studies shows that FESOM can simulate past ocean state and the impact of increased GrIS melting well.

     

Observation of atmospheric methane in the Arctic Ocean up to 87° north

Although the Arctic methane reservoir is large,the emission of methane from the Arctic Ocean into the atmosphere remains poorly constrained.Continuous ship-borne measurements of atmospheric methane near the surface ocean were carried out during two cruises to investigate methane emission from the Arctic Ocean up to the latitude of 87°N.Three-day air mass back trajectories along the cruise tracks indicated that the surface Arctic Ocean could be a potentially important source of methane to the atmosphere.Rapid bursts in methane concentration occurred mainly along the ocean frontal area,suggesting that frontal upwelling in the upper layer of the Arctic Ocean might contribute to methane emissions into the atmosphere.     

A one-dimensional modeling study of the diurnal cycle in the equatorial Atlantic at the PIRATA buoys during the EGEE-3 campaign

A one-dimensional model is used to analyze, at the local scale, the response of the equatorial Atlantic Ocean under different meteorological conditions. The study was performed at the location of three moored buoys of the Pilot Research Moored Array in the Tropical Atlantic located at 10° W, 0° N; 10° W, 6° S; and 10° W, 10° S. During the EGEE-3 (Etude de la circulation océanique et de sa variabilité dans le Golfe de Guinee) campaign of May–June 2006, each buoy was visited for maintenance during 2 days. On board the ship, high-resolution atmospheric parameters were collected, as were profiles of temperature, salinity, and current. These data are used here to initialize, force, and validate a one-dimensional model in order to study the diurnal oceanic mixed-layer variability. It is shown that the diurnal variability of the sea surface temperatures is mainly driven by the solar heat flux. The diurnal response of the near-surface temperatures to daytime heating and nighttime cooling has an amplitude of a few tenths of degree. The computed diurnal heat budget experiences a net warming tendency of 31 and 27 W m⁻² at 0° N and 10° S, respectively, and a cooling tendency of 122 W m⁻² at 6° S. Both observed and simulated mixed-layer depths experience a jump between the nighttime convection phase and the well-stabilized diurnal water column. Its amplitude changes dramatically depending on the meteorological conditions occurring at the stations and reaches its maximum amplitude (~50 m) at 10° S. At 6° and 10° S, the presence of barrier layers is observed, a feature that is clearer at 10° S. Simulated turbulent kinetic energy (TKE) dissipation rates, compared to independent microstructure measurements, show that the model tracks their diurnal evolution reasonably well. It is also shown that the shear and buoyancy productions and the vertical diffusion of TKE all contribute to the supply of TKE, but the buoyancy production is the main source of TKE during the period of the simulation.     

Mechanisms controlling the seasonal mixed-layer temperature and salinity of the Indonesian seas

We examine the seasonal mixed-layer temperature (MLT) and salinity (MLS) budgets in the Banda–Arafura Seas region (120–138° E, 8–3° S) using an ECCO ocean-state estimation product. MLT in these seas is relatively high during November–May (austral spring through fall) and relatively low during June–September (austral winter and the period associated with the Asian summer monsoon). Surface heat flux makes the largest contribution to the seasonal MLT tendency, with significant reinforcement by subsurface processes, especially turbulent vertical mixing. Temperature declines (the MLT tendency is negative) in May–August when seasonal insolation is smallest and local winds are strong due to the southeast monsoon, which causes surface heat loss and cooling by vertical processes. In particular, Ekman suction induced by local wind stress curl raises the thermocline in the Arafura Sea, bringing cooler subsurface water closer to the base of the mixed layer where it is subsequently incorporated into the mixed layer through turbulent vertical mixing; this has a cooling effect. The MLT budget also has a small, but non-negligible, semi-annual component since insolation increases and winds weaken during the spring and fall monsoon transitions near the equator. This causes warming via solar heating, reduced surface heat loss, and weakened turbulent mixing compared to austral winter and, to a lesser extent, compared to austral summer. Seasonal MLS is dominated by ocean processes rather than by local freshwater flux. The contributions by horizontal advection and subsurface processes have comparable magnitudes. The results suggest that ocean dynamics play a significant part in determining both seasonal MLT and MLS in the region, such that coupled model studies of the region should use a full ocean model rather than a slab ocean mixed-layer model.     

A three-dimensional surface wave–ocean circulation coupled model and its initial testing

A theoretical framework to include the influences of nonbreaking surface waves in ocean general circulation models is established based on Reynolds stresses and fluxes terms derived from surface wave-induced fluctuation. An expression for the wave-induced viscosity and diffusivity as a function of the wave number spectrum is derived for infinite and finite water depths; this derivation allows the coupling of ocean circulation models with a wave number spectrum numerical model. In the case of monochromatic surface wave, the wave-induced viscosity and diffusivity are functions of the Stokes drift. The influence of the wave-induced mixing scheme on global ocean circulation models was tested with the Princeton Ocean Model, indicating significant improvement in upper ocean thermal structure and mixed layer depth compared with mixing obtained by the Mellor–Yamada scheme without the wave influence. For example, the model–observation correlation coefficient of the upper 100-m temperature along 35° N increases from 0.68 without wave influence to 0.93 with wave influence. The wave-induced Reynolds stress can reach up to about 5% of the wind stress in high latitudes, and drive 2–3 Sv transport in the global ocean in the form of mesoscale eddies with diameter of 500–1,000 km. The surface wave-induced mixing is more pronounced in middle and high latitudes during the summer in the Northern Hemisphere and in middle latitudes in the Southern Hemisphere.