Observations and operational model simulations reveal the impact of Hurricane Matthew (2016) on the Gulf Stream and coastal sea level

In October 7–9, 2016, Hurricane Matthew moved along the southeastern coast of the U.S., causing major flooding and significant damage, even to locations farther north well away from the storm’s winds. Various observations, such as tide gauge data, cable measurements of the Florida Current (FC) transport, satellite altimeter data and high-frequency radar data, were analyzed to evaluate the impact of the storm. The data show a dramatic decline in the FC flow and increased coastal sea level along the U.S. coast. Weakening of the Gulf Stream (GS) downstream from the storm’s area contributed to high coastal sea levels farther north. Analyses of simulations of an operational hurricane-ocean coupled model reveal the disruption that the hurricane caused to the GS flow, including a decline in transport of ∼20 Sv (1 Sv = 10⁶ m³ s⁻¹). In comparison, the observed FC reached a maximum transport of ∼40 Sv before the storm on September 10 and a minimum of ∼20 Sv after the storm on October 12. The hurricane impacts both the geostrophic part of the GS and the wind-driven currents, generating inertial oscillations with velocities of up to ±1 m s⁻¹. Analysis of the observed FC transport since 1982 indicated that the magnitude of the current weakening in October 2016 was quite rare (outside 3 standard deviations from the mean). Such a large FC weakening in the past occurred more often in October and November, but is extremely rare in June-August. Similar impacts on the FC from past tropical storms and hurricanes suggest that storms may contribute to seasonal and interannual variations in the FC. The results also demonstrated the extended range of coastal impacts that remote storms can cause through their influence on ocean currents.     

Circulation analysis in the northwest Indian Ocean using ARGO floats and surface drifter observations,and SODA reanalysis output

This study incorporates observations from Array of Real-time Geostrophic Oceanography (ARGO) floats and surface drifters to identify seasonal circulation patterns at the surface, 1000 m, 1500 m, and 2000 m in the northwest Indian Ocean, and quantify velocities associated with them. A skill comparison of the Simple Ocean Data Assimilation (SODA) reanalysis output was also performed to contribute to the understanding of the circulation dynamics in this region.Subsurface currents were quantified and validated using the ARGO float data. Surface currents were identified using surface drifter data and compared to the subsurface observations to enhance our previous understanding of surface circulations. Quantified Southwest Monsoon surface currents include the Somali Current (v_max = 179.5 cm/s), the East Arabian Current (v_max = 52.3 cm/s), and the Southwest Monsoon Current (v_max = 51.2 cm/s). Northeastward flow along the Somali coast is also observed at 1000 m (v_max = 26.1 cm/s) and 1500 m (v_max = 12.7 cm/s). Currents associated with the Great Whirl are observed at the surface (v_max = 161.4 cm/s) and at 1000 m (v_max = 16.2 cm/s). In contrast to previous studies, both ARGO and surface drifter data show the Great Whirl can form as early as the boreal Spring intermonsoon, lasting until the boreal Fall intermonsoon. The Arabian Sea exhibits eastward/southeastward flow at the surface, 1000 m, 1500 m, and 2000 m. Quantified Northeast Monsoon surface currents include the Somali Current (v_max = 97.3 cm/s), Northeast Monsoon Current (v_max = 30.0 cm/s), and the North Equatorial Current (v_max = 28.5 cm/s). Southwestward flow along the Somali coast extends as deep as 1500 m.Point-by-point vector and scalar correlations of SODA output to ARGO and surface drifter data showed that surface SODA output and surface drifter data generally produced a strong correlation attributed to surface currents strongly controlled by the monsoons, while subsurface correlations of SODA output and ARGO were mostly insignificant due to variability associated with intermonsoonal transitions. SODA output produced overall smaller velocities than both observational datasets. Assimilating ARGO velocities into the SODA reanalysis could improve subsurface velocity assimilation, especially during the boreal fall and spring when ARGO observations suggest that flow is highly variable.     

Kinematic metrics of the Loop Current in the Gulf of Mexico from satellite altimetry

We analyzed a 20-year time series (January 1st, 1993 through December 31st, 2012) of Loop Current (LC) surface area derived from satellite altimetry in the eastern Gulf of Mexico to estimate kinematical metrics of this potent flow. On average the LC intrudes to its maximum northward position about 216 ± 126 days after the previous eddy separation; and ∼30 ± 31 days later sheds a large anticyclonic eddy. When the northern extent of the LC intrusion following the previous eddy separation is greater than 27°N, the current retreats very quickly until it sheds another eddy with the entire separation process occurring on the order of 30 days. To first order the change in areal extent of the LC during intrusion into the Gulf occurs at an average rate of 225 km² day⁻¹, which corresponds to an intrusion velocity of 1.7 cm s⁻¹ of the LC front, and adds Caribbean water to the Gulf at a rate of 2.6 ± 0.7 Sv.     

Variability of the Somali Current and eddies during the southwest monsoon regimes

The meso-scale eddies and currents in the Arabian Sea are analyzed using different satellite observations, Simple Oceanic Data Assimilation (SODA) reanalysis, and Ocean Reanalysis System 4 (ORAS4) from 1993 to 2016 to investigate the impacts of Southwest (SW) Monsoon strength on Somali Current (SC) mesoscale circulations such as the Great Whirl (GW), the Socotra Eddy (SE), the Southern Gyre (SG), and smaller eddies. Increased Ekman pumping during stronger SW monsoons strengthens coastal upwelling along the Somali coast. The Arabian Sea basin-wide anticyclonic circulation and presence of the GW form mesoscale circulation patterns favourable to advection of upwelled waters eastward into the central Arabian Sea. In September, after the SW monsoon winds reach peak strength in July and August, a higher number of discrete anticyclonic eddies with higher ( > 20 cm) sea surface height anomalies develop in strong and normal intensity SW monsoon seasons than weaker SW monsoon seasons.     

The origin of sub-surface source waters define the sea–air flux of methane in the Mauritanian Upwelling,NW Africa

Concentrations and flux densities of methane were determined during a Lagrangian study of an advective filament in the permanent upwelling region off western Mauritania. Newly upwelled waters were dominated by the presence of North Atlantic Central Water and surface CH₄ concentrations of 2.2 ± 0.3 nmol L⁻¹ were largely in equilibrium with atmospheric values, with surface saturations of 101.7 ± 14%. As the upwelling filament aged and was advected offshore, CH₄ enriched South Atlantic Central Water from intermediate depths of 100–350 m was entrained into the surface mixed layer of the filament following intense mixing associated with the shelf break. Surface saturations increased to 198.9 ± 15% and flux densities increased from a mean value over the shelf of 2.0 ± 1.1 μmol m⁻² d⁻¹ to a maximum of 22.6 μmol m⁻² d⁻¹. Annual CH₄ emissions for this persistent filament were estimated at 0.77 ± 0.64 Gg which equates to a maximum of 0.35% of the global oceanic budget. This raises the known outgassing intensity of this area and highlights the importance of advecting filaments from upwelling waters as efficient vehicles for air-sea exchange.     

A Loop Current experiment: Field and remote measurements

An overview of a new comprehensive observational study of the Loop Current (LC) in the eastern Gulf of Mexico that encompassed full-depth and near-bottom moorings, pressure-equipped inverted echo sounders (PIES) and remote sensing is presented. The study array was designed to encompass the LC from the Campeche Bank to the west Florida escarpment. This overview centers about principal findings as they pertain to mesoscale dynamics. Two companion papers provide in-depth analyses. Three LC anticyclonic eddy separation events were observed with good 3D spatial coverage over the 2½ year extent of the field study; the three separations exhibited similar processes after the LC had extended into the eastern Gulf. Large scale (∼300 km wavelength, 40–60 day periods) southward propagating meanders developed on the eastern side of the LC over deep (∼3000 m) water that were the result of baroclinic instability between the upper layer meandering jet and lower layer cyclones and anticyclones. The lower layer was only highly energetic during relatively short (∼2–3 months) intervals just prior to or during eddy detachments because of baroclinic instability. The steepening of the meanders lead to a pinch-off of LC eddies. The deep lower-layer eddies, constrained by the closed topography of the southeastern Gulf, propagated westward across the detachment zone and appear to assist in achieving separation. Small scale (∼50–100 km, periods ∼10 days) frontal eddies, observed on the western side of the LC along the Campeche Bank slope, decay over the deep water of the northern part of an extended LC, and have little influence on lower layer eddies, the east side meanders and the eddy detachment processes.     

Using the mean pressure gradient and NCEP/NCAR reanalysis to estimate the strength of the South Atlantic Anticyclone

A new methodology is proposed to estimate the strength of the South Atlantic Anticyclone (SAA), using the gridded sea level pressure (SLP) of the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) Reanalysis data. The top quartile (1017.3 hPa) of the SLP data was found a reasonable criterion to delimit the SAA area. Consequently, we defined the SAA area as the quadrangle containing 80% of the observations with pressure >1017.3 hPa. In this quadrangle, an area weighted pressure gradient (AWPG) was computed for the whole area and for the north–south and west–east halves. When compared with maximum pressure, the AWPG showed a better correlation with the significant wave height (SWH) and wind speed (WS) derived from altimetry. The mean value of the AWPG was 8 × 10⁻⁴ Pa/m, with representative values of 9.1 × 10⁻⁴ Pa/m and 7.4 × 10⁻⁴ Pa/m for austral winter and summer, respectively. The phase difference between the monthly AWPG in the north and south sub-quadrangles accounts for the evolution of the spatial pattern of the anticyclone throughout a year. This quantitative approach proved to be a useful estimate of the strength of South Atlantic Anticyclone. Further improvements of this approach are discussed.     

Distribution of CO2 parameters in the Western Tropical Atlantic Ocean

The variability of sea surface Total Alkalinity (TA) and sea surface Total Inorganic Carbon (C_T) is examined using all available data in the western tropical Atlantic (WTA: 20°S-20°N, 60°W-20°W). Lowest TA and C_T are observed for the region located between 0°N-15°N/60°W-50°W and are explained by the influence of the Amazon plume during boreal summer. In the southern part of the area, 20°S-10°S/40°W-60°W, the highest values of TA and C_T are linked to the CO₂–rich waters due to the equatorial upwelling, which are transported by the South Equatorial Current (SEC) flowing from the African coast to the Brazilian shore. An increase of C_T of 0.9 ± 0.3 μmol kg⁻¹yr⁻¹ has been observed in the SEC region and is consistent with previous published estimates. A revised C_T-Sea Surface Salinity (SSS) relationship is proposed for the WTA to take into account the variability of C_T at low salinities. This new C_T-SSS relationship together with a published TA-SSS relationship allow to calculate pCO₂ values that compare well with observed pCO₂ (R² = 0.90).     

Increasing flood risk and wetland losses due to global sea-level rise: regional and global analyses

To develop improved estimates of (1) flooding due to storm surges, and (2) wetland losses due to accelerated sea-level rise, the work of Hoozemans et al. (1993) is extended to a dynamic analysis. It considers the effects of several simultaneously changing factors, including: (1) global sea-level rise and subsidence; (2) increasing coastal population; and (3) improving standards of flood defence (using GNP/capita as an “ability-to-pay” parameter). The global sea-level rise scenarios are derived from two General Circulation Model (GCM) experiments of the Hadley Centre: (1) the HadCM2 greenhouse gas only ensemble experiment and (2) the more recent HadCM3 greenhouse gas only experiment. In all cases there is a global rise in sea level of about 38 cm from 1990 to the 2080s. No other climate change is considered. Relative to an evolving reference scenario without sea-level rise, this analysis suggests that the number of people flooded by storm surge in a typical year will be more than five times higher due to sea-level rise by the 2080s. Many of these people will experience annual or more frequent flooding, suggesting that the increase in flood frequency will be more than nuisance level and some response (increased protection, migration, etc.) will be required. In absolute terms, the areas most vulnerable to flooding are the southern Mediterranean, Africa, and most particularly, South and South-east Asia where there is a concentration of low-lying populated deltas. However, the Caribbean, the Indian Ocean islands and the Pacific Ocean small islands may experience the largest relative increase in flood risk. By the 2080s, sea-level rise could cause the loss of up to 22% of the world's coastal wetlands. When combined with other losses due to direct human action, up to 70% of the world's coastal wetlands could be lost by the 2080s, although there is considerable uncertainty. Therefore, sea-level rise would reinforce other adverse trends of wetland loss. The largest losses due to sea-level rise will be around the Mediterranean and Baltic and to a lesser extent on the Atlantic coast of Central and North America and the smaller islands of the Caribbean. Collectively, these results show that a relatively small global rise in sea level could have significant adverse impacts if there is no adaptive response. Given the “commitment to sea-level rise” irrespective of any realistic future emissions policy, there is a need to start strategic planning of appropriate responses now. Given that coastal flooding and wetland loss are already important problems, such planning could have immediate benefits.     

Variability of CO2 fugacity at the western edge of the tropical Atlantic Ocean from the 8°N to 38°W PIRATA buoy

Hourly data of CO₂ fugacity (fCO₂) at 8°N–38°W were analyzed from 2008 to 2011. Analyses of wind, rainfall, temperature and salinity data from the buoy indicated two distinct seasonal periods. The first period (January to July) had a mean fCO₂ of 378.9 μatm (n = 7512). During this period, in which the study area was characterized by small salinity variations, the fCO₂ is mainly controlled by sea surface temperature (SST) variations (fCO₂ = 24.4*SST-281.1, r² = 0.8). During the second period (August–December), the mean fCO₂ was 421.9 μatm (n = 11571). During these months, the region is subjected to the simultaneous action of (a) rainfall induced by the presence of the Intertropical Convergence Zone (ITCZ); (b) arrival of fresh water from the Amazon River plume that is transported to the east by the North Equatorial Countercurrent (NECC) after the retroflection of the North Brazil Current (NBC); and (c) vertical input of CO₂-rich water due to Ekman pumping. The data indicated the existence of high-frequency fCO₂ variability (periods less than 24 h). This high variability is related to two different mechanisms. In the first mechanism, fCO₂ increases are associated to rapid increases in SST and are attributed to the diurnal cycle of solar radiation. In addition, low wind speed contributes to SST rising by inhibiting vertical mixing. In the second mechanism, fCO₂ decreases are associated to SSS decreases caused by heavy rainfall.     

Measurements and simulation of the M-component current and simultaneous electromagnetic fields at 60 m and 550 m

Simultaneous measurements of the M-component current (surges superimposed on lightning continuing currents) and the corresponding electromagnetic fields at 60 m and 550 m from the lightning channel are analyzed and simulated with a two-wave model. The measured results reveal that the M-component current at the bottom of the channel exhibits a V-shape character with a leading edge of 78 μs and a trailing edge of 194 μs, while the electric field pulses at 60 m and 550 m have trailing edges faster than leading edges. The peak of the M-component current lags behind the electric field peak by tens of microseconds, when the distance increases to 550 m, the disparity of the time shift increases as well. However, the waveshape of the M-component current is similar to that of the magnetic field pulse. The M-component electric fields at 60 m and 550 m are 1.16 kV/m and 0.17 kV/m, respectively, and exhibit a logarithmic distance dependence which implies that the M-component charge density increases with height. Additionally, a two-wave model is used to examine the sensitivity of the predicted electric and magnetic fields to the speed and current reflection coefficient variations of the M-component. The simulated results show that the effects are different for the electric and magnetic fields. The M-component speed essentially controls the electric field, but has little effect on the magnetic field. Larger reflection coefficient results in a larger magnetic field, but a smaller electric field.     

The seasonal variation of the upper layers of the South China Sea (SCS) circulation and the Indonesian through flow (ITF): An ocean model study

The upper layer, wind-driven circulation of the South China Sea (SCS), its through-flow (SCSTF) and the Indonesian through flow (ITF) are simulated using a high resolution model, FVCOM (finite volume coastal ocean model) in a regional domain comprising the Maritime Continent. The regional model is embedded in the MIT global ocean general circulation model (ogcm) which provides surface forcing and boundary conditions of all the oceanographic variables at the lateral open boundaries in the Pacific and Indian oceans. A five decade long simulation is available from the MITgcm and we choose to investigate and compare the climatologies of two decades, 1960–1969 and 1990–1999.The seasonal variability of the wind-driven circulation produced by the monsoon system is realistically simulated. In the SCS the dominant driving force is the monsoon wind and the surface circulation reverses accordingly, with a net cyclonic tendency in winter and anticyclonic in summer. The SCS circulation in the 90s is weaker than in the 60s because of the weaker monsoon system in the 90s. In the upper 50 m the interaction between the SCSTF and ITF is very important. The southward ITF can be blocked by the SCSTF at the Makassar Strait during winter. In summer, part of the ITF feeds the SCSTF flowing into the SCS through the Karimata Strait. Differently from the SCS, the ITF is primarily controlled by the sea level difference between the western Pacific and eastern Indian Ocean. The ITF flow, consistently southwestward below the surface layer, is stronger in the 90s.The volume transports for winter, summer and yearly are estimated from the simulation through all the interocean straits. On the annual average, there is a ∼5.6 Sv of western Pacific water entering the SCS through the Luzon Strait and ∼1.4 Sv exiting through the Karimata Strait into the Java Sea. Also, ∼2 Sv of SCS water enters the Sulu Sea through the Mindoro Strait, while ∼2.9 Sv flow southwards through the Sibutu Strait merging into the ITF. The ITF inflow occurs through the Makassar Strait (up to ∼62%) and the Lifamatola Strait (∼38%). The annual average volume transport of the ITF inflow from the simulation is ∼15 Sv in the 60s and ∼16.6 Sv in the 90s, very close to the long term observations. The ITF outflow through the Lombok, Ombai and Timor straits is ∼16.8 Sv in the 60s and 18.9 Sv in the 90s, with the outflow greater by 1.7 Sv and 2.3 Sv respectively. The transport estimates of the simulation at all the straits are in rather good agreement with the observational estimates.We analyze the thermal structure of the domain in the 60s and 90s and assess the simulated temperature patterns against the SODA reanalysis product, with special focus on the shallow region of the SCS. The SODA dataset clearly shows that the yearly averaged temperatures of the 90s are overall warmer than those of the 60s in the surface, intermediate and some of the deep layers and the decadal differences (90s − 60s) indicate that the overall warming of the SCS interior is a local effect. In the simulation the warm trend from the 60s to the 90s in well reproduced in the surface layer. In particular, the simulated temperature profiles at two shallow sites at midway in the SCSTF agree rather well with the SODA profiles. However, the warming trend in the intermediate (deep) layers is not reproduced in the simulation. We find that this deficiency is mostly due to a deficiency in the initial temperature fields provide by the MITgcm.     

The Gulf Stream pathway and the impacts of the eddy-driven abyssal circulation and the Deep Western Boundary Current

A hydrodynamic model of the subtropical Atlantic basin and the Intra-Americas Sea (9–47°N) is used to investigate the dynamics of Gulf Stream separation from the western boundary at Cape Hatteras and its mean pathway to the Grand Banks. The model has five isopycnal Lagrangian layers in the vertical and allows realistic boundary geometry, bathymetry, wind forcing, and a meridional overturning circulation (MOC), the latter specified via ports in the northern and southern boundaries. The northward upper ocean branch of the MOC (14 Sv) was always included but the southward Deep Western Boundary Current (DWBC) was excluded in some simulations, allowing investigation of the impacts of the DWBC and the eddy-driven mean abyssal circulation on Gulf Stream separation from the western boundary. The result is resolution dependent with the DWBC playing a crucial role in Gulf Stream separation at 1/16° resolution but with the eddy-driven abyssal circulation alone sufficient to obtain accurate separation at 1/32° resolution and a realistic pathway from Cape Hatteras to the Grand Banks with minimal DWBC impact except southeast of the Grand Banks. The separation from the western boundary is particularly sensitive to the strength of the eddy-driven abyssal circulation. Farther to the east, between 68°W and the Grand Banks, all of the 1/16° and 1/32° simulations with realistic topography (with or without a DWBC) gave similar generally realistic mean pathways with clear impacts of the topographically constrained eddy-driven abyssal circulation versus very unrealistic Gulf Stream pathways between Cape Hatteras and the Grand Banks from otherwise identical simulations run with a flat bottom, in reduced-gravity mode, or with 1/8° resolution and realistic topography. The model is realistic enough to allow detailed model-data comparisons and a detailed investigation of Gulf Stream dynamics. The corresponding linear solution with a Sverdrup interior and Munk viscous western boundary layers, including one from the northward branch of the MOC, yielded two unrealistic Gulf Stream pathways, a broad eastward pathway centered at the latitude of Cape Hatteras and a second wind plus MOC-driven pathway hugging the western boundary to the north. Thus, a high resolution model capable of simulating an inertial jet is required to obtain a single nonlinear Gulf Stream pathway as it separates from the coast. None of the simulations were sufficiently inertial to overcome the linear solution need for a boundary current north of Cape Hatteras without assistance from pathway advection by the abyssal circulation, even though the core speeds of the simulated currents were consistent with observations near separation. In the 1/16° simulation with no DWBC and a 1/32° simulation with high bottom friction and no DWBC the model Gulf Stream overshot the observed separation latitude. With abyssal current assistance the simulated (and the observed) mean Gulf Stream pathway between separation from the western boundary and 70°W agreed closely with a constant absolute vorticity (CAV) trajectory influenced by the angle of the coastline prior to separation. The key abyssal current crosses under the Gulf Stream at 68.5–69°W and advects the Gulf Stream pathway southward to the terminus of an escarpment in the continental slope. There the abyssal current crosses to deeper depths to conserve potential vorticity while passing under the downward-sloping thermocline of the stream and then immediately retroflects eastward onto the abyssal plain, preventing further southward pathway advection. Thus specific topographic features and feedback from the impact of the Gulf Stream on the abyssal current pathway determined the latitude of the stream at 68.5–69°W, a latitude verified by observations. The associated abyssal current was also verified by observations.     

Analysis of rainfall characteristics of the Madden–Julian oscillation using TRMM satellite data

Rainfall characteristics of the Madden–Julian oscillation (MJO) are analyzed primarily using tropical rainfall measuring mission (TRMM) precipitation radar (PR), TRMM microwave imager (TMI) and lighting imaging sensor (LIS) data. Latent heating structure is also examined using latent heating data estimated with the spectral latent heating (SLH) algorithm.The zonal structure, time evolution, and characteristic stages of the MJO precipitation system are described. Stratiform rain fraction increases with the cloud activity, and the amplitude of stratiform rain variation associated with the MJO is larger than that of convective rain by a factor of 1.7. Maximum peaks of both convective rain and stratiform rain precede the minimum peak of the outgoing longwave radiation (OLR) anomaly which is often used as a proxy for the MJO convection. Stratiform rain remains longer than convective rain until ∼4000 km behind the peak of the mature phase. The stratiform rain contribution results in the top-heavy heating profile of the MJO.Associated with the MJO, there are tri-pole convective rain top heights (RTH) at 10–11, ∼7 and ∼3 km, corresponding to the dominance of afternoon showers, organized systems, and shallow convections, respectively. The stratiform rain is basically organized with convective rain, having similar but slightly lower RTH and slightly lags the convective rain maximum. It is notable that relatively moderate (∼7 km) RTH is dominant in the mature phase of the MJO, while very tall rainfall with RTH over 10 km and lightning frequency increase in the suppressed phase. The rain-yield-per flash (RPF) varies about 20–100% of the mean value of ∼2–10 × 10⁹ kg fl⁻¹ over the tropical warm ocean and that of ∼2–5 × 10⁹ kg fl⁻¹ over the equatorial Islands, between the convectively suppressed phase and the active phase of MJO, in the manner that RPF is smaller in the suppressed phase and larger in the active phase.     

Simulated variability of the Atlantic meridional overturning circulation

To examine the multi-annual to decadal scale variability of the Atlantic Meridional Overturning Circulation (AMOC) we conducted a four-member ensemble with a daily reanalysis forced, medium-resolution global version of the isopycnic coordinate ocean model MICOM, and a 300-years integration with the fully coupled Bergen Climate Model (BCM). The simulations of the AMOC with both model systems yield a long-term mean value of 18 Sv and decadal variability with an amplitude of 1–3 Sv. The power spectrum of the inter-annual to decadal scale variability of the AMOC in BCM generally follows the theoretical red noise spectrum, with indications of increased power near the 20-years period. Comparison with observational proxy indices for the AMOC, e.g. the thickness of the Labrador Sea Water, the strength of the baroclinic gyre circulation in the North Atlantic Ocean, and the surface temperature anomalies along the mean path of the Gulf Stream, shows similar trends and phasing of the variability, indicating that the simulated AMOC variability is robust and real. Mixing indices have been constructed for the Labrador, the Irminger and the Greenland-Iceland-Norwegian (GIN) seas. While convective mixing in the Labrador and the GIN seas are in opposite phase, and linked to the NAO as observations suggest, the convective mixing in the Irminger Sea is in phase with or leads the Labrador Sea. Newly formed deep water is seen as a slow, anomalous cold and fresh, plume flowing southward along the western continental slope of the Atlantic Ocean, with a return flow of warm and saline water on the surface. In addition, fast-travelling topographically trapped waves propagate southward along the continental slope towards equator, where they go east and continue along the eastern rim of the Atlantic. For both types of experiments, the Northern Hemisphere sea level pressure and 2 m temperature anomaly patterns computed based on the difference between climate states with strong and weak AMOC yields a NAO-like pattern with intensified Icelandic low and Azores high, and a warming of 0.25–0.5 °C of the central North Atlantic sea-surface temperature (SST). The reanalysis forced simulations indicate a coupling between the Labrador Sea Water production rate and an equatorial Atlantic SST index in accordance with observations. This coupling is not identified in the coupled simulation.     

Westerly wind events,diurnal cycle and central Pacific El Niño warming

In the past three decades, the strongest central Pacific (CP) El Niño event was observed in 2009–2010 by satellites. When intensity of this CP El Niño reached its maximum, large diurnal variations of sea surface temperature (SST) were also observed from tropical atmosphere ocean moorings in the central equatorial Pacific. Solar radiation in the equatorial central Pacific is larger than 140 W/m², which leads to the amplitude of diurnal cycle of SST primarily determined by large-scale wind patterns. Intraseasonal westerly wind events (WWEs) can lead to an eastward displacement of the warm pool and also can weaken the trade winds in central Pacific. When the occurrence of equatorial WWEs is more than 20 days in a month, monthly mean wind speed in central equatorial Pacific has high possibility of wind speed less than 3 m/s, thus has pronounced diurnal cycle of SST. The diurnal cycle of SST will rectify daily mean SST. Reduced mixing at the base of the mixed layer and suppression of entrainment due to the accumulated effect of diurnal cycle may lead to warmer SST in the following month. This study suggests the occurrence of more diurnal SST events may contribute to the increasing intensity of the CP El Niño events.     

Impacts of the South China Sea Throughflow on seasonal and interannual variations of the Indonesian Throughflow

Impacts of the South China Sea Throughflow (SCST) on seasonal and interannual variations of the Indonesian Throughflow are studied by comparing outputs from ocean general circulation model (OGCM) experiments with and without the SCST. The observed subsurface maximum in the southward flow through the Makassar Strait is simulated only when the SCST, which is driven by the large-scale wind, is allowed in the model. The mean volume and heat transport by the Makassar Strait Throughflow are reduced by 1.7 Sv and 0.19 PW, respectively, by the existence of the SCST in the model. The difference is particularly remarkable during boreal winter when the SCST reaches its seasonal maximum. Furthermore, the SCST is strengthened during El Niño, leading to the weakening in the southward volume and heat transport through the Makassar Strait by 0.37 Sv and 0.05 PW, respectively. These findings from the OGCM experiments suggest that the SCST may play an important role in climate variability of the Indo-Pacific Ocean.     

Climate change impacts on ecosystems and the terrestrial carbon sink: a new assessment

Climate output from the UK Hadley Centre's HadCM2 and HadCM3 experiments for the period 1860 to 2100, with IS92a greenhouse gas forcing, together with predicted patterns of N deposition and increasing CO₂, were input (offline) to the dynamic vegetation model, Hybrid v4.1 (Friend et al., 1997; Friend and White, 1999). This model represents biogeochemical, biophysical and biogeographical processes, coupling the carbon, nitrogen and water cycles on a sub-daily timestep, simulating potential vegetation and transient changes in annual growth and competition between eight generalized plant types in response to climate.Global vegetation carbon was predicted to rise from about 600 to 800 PgC (or to 650 PgC for HadCM3) while the soil carbon pool of about 1100 PgC decreased by about 8%. By the 2080s, climate change caused a partial loss of Amazonian rainforest, C₄ grasslands and temperate forest in areas of southern Europe and eastern USA, but an expansion in the boreal forest area. These changes were accompanied by a decrease in net primary productivity (NPP) of vegetation in many tropical areas, southern Europe and eastern USA (in response to warming and a decrease in rainfall), but an increase in NPP of boreal forests. Global NPP increased from 45 to 50 PgC y⁻¹ in the 1990s to about 65 PgC y⁻¹ in the 2080s (about 58 PgC y⁻¹ for HadCM3). Global net ecosystem productivity (NEP) increased from about 1.3 PgC y⁻¹ in the 1990s to about 3.6 PgC y⁻¹ in the 2030s and then declined to zero by 2100 owing to a loss of carbon from declining forests in the tropics and at warm temperate latitudes — despite strengthening of the carbon sink at northern high latitudes. HadCM3 gave a more erratic temporal evolution of NEP than HadCM2, with a dramatic collapse in NEP in the 2050s.     

Sensitivity studies with the North Atlantic sigma coordinate Princeton Ocean Model

The sigma coordinate, Princeton Ocean Model (POM) has been configured for the North Atlantic Ocean between 5°N and 50°N as part of data assimilation, model predictability and intercomparison studies. The model uses a curvilinear orthogonal grid with higher resolution in the western North Atlantic and lower resolution in the eastern North Atlantic. A series of experiments, each one of a 10-year duration, are performed to evaluate the sensitivity of the ocean mean state and variability to model parameters and model configuration; these experiments include open vs. closed boundary conditions, low vs. high resolution grids, and different choices of diffusion and viscosity. The results show that the use of closed boundaries together with near-boundary buffer zones where temperature and salinity are relaxed towards the observed values give less realistic flows, weaker recirculation gyres and less realistic Gulf Stream separation than do open boundary conditions. The experiments show that the sensitivity of the ocean variability in the model to the choice of the Smagorinsky diffusion and viscosity coefficients significantly differs from one region to another and largely depends on other attributes such as the mean position of the Gulf Stream in each simulation. A 50% change in model resolution in the Gulf Stream region has a larger effect on ocean variability than a change of diffusivity by a factor of 10. In areas where either the high or the low resolution models have sufficient resolution, as in the Gulf of Mexico, they are able to produce variability comparable to that observed from altimeter data; elsewhere, model variability is underestimated.     

Decadal- to interannual-scale source water variations in the Caribbean Sea recorded by Puerto Rican coral radiocarbon

Water that forms the Florida Current, and eventually the Gulf Stream, coalesces in the Caribbean from both subtropical and equatorial sources. The equatorial sources are made up of, in part, South Atlantic water moving northward and compensating for southward flow at depth related to meridional overturning circulation. Subtropical surface water contains relatively high amounts of radiocarbon (¹⁴C), whereas equatorial waters are influenced by the upwelling of low ¹⁴C water and have relatively low concentrations of ¹⁴C. We use a 250 year record of Δ¹⁴C in a coral from southwestern Puerto Rico along with previously published coral Δ¹⁴C records as tracers of subtropical and equatorial water mixing in the northern Caribbean. Data generated in this study and from other studies indicate that the influence of either of the two water masses can change considerably on interannual to interdecadal time scales. Variability due to ocean dynamics in this region is large relative to variability caused by atmospheric ¹⁴C changes, thus masking the Suess effect at this site. A mixing model produced using coral Δ¹⁴C illustrates the time varying proportion of equatorial versus subtropical waters in the northern Caribbean between 1963 and 1983. The results of the model are consistent with linkages between multidecadal thermal variability in the North Atlantic and meridional overturning circulation. Ekman transport changes related to tradewind variability are proposed as a possible mechanism to explain the observed switches between relatively low and high Δ¹⁴C values in the coral radiocarbon records.