Publications with links

Export 129 results:
Sort by: Author Title Type [ Year  (Desc)]
Llanillo, PJ, Pelegri JL, Talley LD, Pena-Izquierdo J, Cordero RR.  2018.  Oxygen pathways and budget for the Eastern South Pacific Oxygen Minimum Zone. Journal of Geophysical Research-Oceans. 123:1722-1744.   10.1002/2017jc013509   AbstractWebsite

Ventilation of the eastern South Pacific Oxygen Minimum Zone (ESP-OMZ) is quantified using climatological Argo and dissolved oxygen data, combined with reanalysis wind stress data. We (1) estimate all oxygen fluxes (advection and turbulent diffusion) ventilating this OMZ, (2) quantify for the first time the oxygen contribution from the subtropical versus the traditionally studied tropical-equatorial pathway, and (3) derive a refined annual-mean oxygen budget for the ESP-OMZ. In the upper OMZ layer, net oxygen supply is dominated by tropical-equatorial advection, with more than one-third of this supply upwelling into the Ekman layer through previously unevaluated vertical advection, within the overturning component of the regional Subtropical Cell (STC). Below the STC, at the OMZ's core, advection is weak and turbulent diffusion (isoneutral and dianeutral) accounts for 89% of the net oxygen supply, most of it coming from the oxygen-rich subtropical gyre. In the deep OMZ layer, net oxygen supply occurs only through turbulent diffusion and is dominated by the tropical-equatorial pathway. Considering the entire OMZ, net oxygen supply (3.8 +/- 0.42 mu mol kg(-1) yr(-1)) is dominated by isoneutral turbulent diffusion (56.5%, split into 32.3% of tropical-equatorial origin and 24.2% of subtropical origin), followed by isoneutral advection (32.0%, split into 27.6% of tropical-equatorial origin and 4.4% of subtropical origin) and dianeutral diffusion (11.5%). One-quarter (25.8%) of the net oxygen input escapes through dianeutral advection (most of it upwelling) and, assuming steady state, biological consumption is responsible for most of the oxygen loss (74.2%).

Tamsitt, V, Abernathey RP, Mazloff MR, Wang J, Talley LD.  2018.  Transformation of deep water masses along Lagrangian upwelling pathways in the Southern Ocean. Journal of Geophysical Research: Oceans.   10.1002/2017JC013409   AbstractWebsite

Upwelling of northern deep waters in the Southern Ocean is fundamentally important for the closure of the global meridional overturning circulation and delivers carbon and nutrient‐rich deep waters to the sea surface. We quantify water mass transformation along upwelling pathways originating in the Atlantic, Indian, and Pacific and ending at the surface of the Southern Ocean using Lagrangian trajectories in an eddy‐permitting ocean state estimate. Recent related work shows that upwelling in the interior below about 400 m depth is localized at hot spots associated with major topographic features in the path of the Antarctic Circumpolar Current, while upwelling through the surface layer is more broadly distributed. In the ocean interior upwelling is largely isopycnal; Atlantic and to a lesser extent Indian Deep Waters cool and freshen while Pacific deep waters are more stable, leading to a homogenization of water mass properties. As upwelling water approaches the mixed layer, there is net strong transformation toward lighter densities due to mixing of freshwater, but there is a divergence in the density distribution as Upper Circumpolar Deep Water tends become lighter and dense Lower Circumpolar Deep Water tends to become denser. The spatial distribution of transformation shows more rapid transformation at eddy hot spots associated with major topography where density gradients are enhanced; however, the majority of cumulative density change along trajectories is achieved by background mixing. We compare the Lagrangian analysis to diagnosed Eulerian water mass transformation to attribute the mechanisms leading to the observed transformation.

Rosso, I, Mazloff MR, Verdy A, Talley LD.  2017.  Space and time variability of the Southern Ocean carbon budget. Journal of Geophysical Research-Oceans. 122:7407-7432.   10.1002/2016jc012646   AbstractWebsite

The upper ocean dissolved inorganic carbon (DIC) concentration is regulated by advective and diffusive transport divergence, biological processes, freshwater, and air-sea CO2 fluxes. The relative importance of these mechanisms in the Southern Ocean is uncertain, as year-round observations in this area have been limited. We use a novel physical-biogeochemical state estimate of the Southern Ocean to construct a closed DIC budget of the top 650 m and investigate the spatial and temporal variability of the different components of the carbon system. The dominant mechanisms of variability in upper ocean DIC depend on location and time and space scales considered. Advective transport is the most influential mechanism and governs the local DIC budget across the 10 day-5 year timescales analyzed. Diffusive effects are nearly negligible. The large-scale transport structure is primarily set by upwelling and downwelling, though both the lateral ageostrophic and geostrophic transports are significant. In the Antarctic Circumpolar Current, the carbon budget components are also influenced by the presence of topography and biological hot spots. In the subtropics, evaporation and air-sea CO2 flux primarily balances the sink due to biological production and advective transport. Finally, in the subpolar region sea ice processes, which change the seawater volume and thus the DIC concentration, compensate the large impact of the advective transport and modulate the timing of biological activity and air-sea CO2 flux.

Johnson, KS, Plant JN, Dunne JP, Talley LD, Sarmiento JL.  2017.  Annual nitrate drawdown observed by SOCCOM profiling floats and the relationship to annual net community production. Journal of Geophysical Research-Oceans. 122:6668-6683.   10.1002/2017jc012839   AbstractWebsite

Annual nitrate cycles have been measured throughout the pelagic waters of the Southern Ocean, including regions with seasonal ice cover and southern hemisphere subtropical zones. Vertically resolved nitrate measurements were made using in situ ultraviolet spectrophotometer (ISUS) and submersible ultraviolet nitrate analyzer (SUNA) optical nitrate sensors deployed on profiling floats. Thirty-one floats returned 40 complete annual cycles. The mean nitrate profile from the month with the highest winter nitrate minus the mean profile from the month with the lowest nitrate yields the annual nitrate drawdown. This quantity was integrated to 200 m depth and converted to carbon using the Redfield ratio to estimate annual net community production (ANCP) throughout the Southern Ocean south of 30 degrees S. A well-defined, zonal mean distribution is found with highest values (3-4 mol C m(-2) yr(-1)) from 40 to 50 degrees S. Lowest values are found in the subtropics and in the seasonal ice zone. The area weighted mean was 2.9 mol C m(-2) yr(-1) for all regions south of 40 degrees S. Cumulative ANCP south of 50 degrees S is 1.3 Pg C yr(-1). This represents about 13% of global ANCP in about 14% of the global ocean area. Plain Language Summary This manuscript reports on 40 annual cycles of nitrate observed by chemical sensors on SOCCOM profiling floats. The annual drawdown in nitrate concentration by phytoplankton is used to assess the spatial variability of annual net community production in the Southern Ocean. This ANCP is a key component of the global carbon cycle and it exerts an important control on atmospheric carbon dioxide. We show that the results are consistent with our prior understanding of Southern Ocean ANCP, which has required decades of observations to accumulate. The profiling floats now enable annual resolution of this key process. The results also highlight spatial variability in ANCP in the Southern Ocean.

Johnson, KS, Plant JN, Coletti LJ, Jannasch HW, Sakamoto CM, Riser SC, Swift DD, Williams NL, Boss E, Haentjens N, Talley LD, Sarmiento JL.  2017.  Biogeochemical sensor performance in the SOCCOM profiling float array. Journal of Geophysical Research-Oceans. 122:6416-6436.   10.1002/2017jc012838   AbstractWebsite

The Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) program has begun deploying a large array of biogeochemical sensors on profiling floats in the Southern Ocean. As of February 2016, 86 floats have been deployed. Here the focus is on 56 floats with quality-controlled and adjusted data that have been in the water at least 6 months. The floats carry oxygen, nitrate, pH, chlorophyll fluorescence, and optical backscatter sensors. The raw data generated by these sensors can suffer from inaccurate initial calibrations and from sensor drift over time. Procedures to correct the data are defined. The initial accuracy of the adjusted concentrations is assessed by comparing the corrected data to laboratory measurements made on samples collected by a hydrographic cast with a rosette sampler at the float deployment station. The long-term accuracy of the corrected data is compared to the GLODAPv2 data set whenever a float made a profile within 20 km of a GLODAPv2 station. Based on these assessments, the fleet average oxygen data are accurate to 1 +/- 1%, nitrate to within 0.5 +/- 0.5 mu mol kg(-1), and pH to 0.005 +/- 0.007, where the error limit is 1 standard deviation of the fleet data. The bio-optical measurements of chlorophyll fluorescence and optical backscatter are used to estimate chlorophyll a and particulate organic carbon concentration. The particulate organic carbon concentrations inferred from optical backscatter appear accurate to with 35 mg C m(-3) or 20%, whichever is larger. Factors affecting the accuracy of the estimated chlorophyll a concentrations are evaluated.

Haentjens, N, Boss E, Talley LD.  2017.  Revisiting Ocean Color algorithms for chlorophyll a and particulate organic carbon in the Southern Ocean using biogeochemical floats. Journal of Geophysical Research-Oceans. 122:6583-6593.   10.1002/2017jc012844   AbstractWebsite

The Southern Ocean (SO) ecosystem plays a key role in the carbon cycle by sinking a major part (43%) of the ocean uptake of anthropogenic CO2, and being an important source of nutrients for primary producers. However, undersampling of SO biogeochemical properties limits our understanding of the mechanisms taking place in this remote area. The Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) project has been deploying a large number of autonomous biogeochemical floats to study the SO (as of December 2016, 74 floats out of 200 have been deployed). SOCCOM floats measurements can be used to extend remote sensing chlorophyll a (chl a) and particulate organic carbon (POC) products under clouds or during the polar night as well as adding the depth dimension to the satellite-based view of the SO. Chlorophyll a concentrations measured by a sensor embedded on the floats and POC concentrations derived from backscattering coefficients were calibrated with samples collected during the floats' deployment cruise. Float chl a and POC were compared with products derived from observations of MODIS and VIIRS sensors. We find the Ocean Color Index (OCI) global algorithm to agree well with the matchups (within 9%, on average, for the Visible Infrared Imaging Radiometer Suite (VIIRS) and 12%, on average, for the Moderate Resolution Imaging Spectroradiometer Aqua (MODIS)). SO-specific algorithms estimating chl a are offset by similar to 45% south of the Sea Ice Extent Front (similar to 60 degrees S). In addition, POC estimates based on floats agree well with NASA's POC algorithm.

Tamsitt, V, Drake HF, Morrison AK, Talley LD, Dufour CO, Gray AR, Griffies SM, Mazloff MR, Sarmiento JL, Wang J, Weijer W.  2017.  Spiraling pathways of global deep waters to the surface of the Southern Ocean. Nature Communications. 8:172.   10.1038/s41467-017-00197-0   Abstract

Upwelling of global deep waters to the sea surface in the Southern Ocean closes the global overturning circulation and is fundamentally important for oceanic uptake of carbon and heat, nutrient resupply for sustaining oceanic biological production, and the melt rate of ice shelves. However, the exact pathways and role of topography in Southern Ocean upwelling remain largely unknown. Here we show detailed upwelling pathways in three dimensions, using hydrographic observations and particle tracking in high-resolution models. The analysis reveals that the northern-sourced deep waters enter the Antarctic Circumpolar Current via southward flow along the boundaries of the three ocean basins, before spiraling southeastward and upward through the Antarctic Circumpolar Current. Upwelling is greatly enhanced at five major topographic features, associated with vigorous mesoscale eddy activity. Deep water reaches the upper ocean predominantly south of the Antarctic Circumpolar Current, with a spatially nonuniform distribution. The timescale for half of the deep water to upwell from 30° S to the mixed layer is ~60–90 years.

Holte, J, Talley LD, Gilson J, Roemmich D.  2017.  An Argo mixed layer climatology and database. Geophysical Research Letters. 44:5618-5626.   10.1002/2017gl073426   AbstractWebsite

A global climatology and database of mixed layer properties are computed from nearly 1,250,000 Argo profiles. The climatology is calculated with both a hybrid algorithm for detecting the mixed layer depth (MLD) and a standard threshold method. The climatology provides accurate information about the depth, properties, extent, and seasonal patterns of global mixed layers. The individual profile results in the database can be used to construct time series of mixed layer properties in specific regions of interest. The climatology and database are available online at . The MLDs calculated by the hybrid algorithm are shallower and generally more accurate than those of the threshold method, particularly in regions of deep winter mixed layers; the new climatology differs the most from existing mixed layer climatologies in these regions. Examples are presented from the Labrador and Irminger Seas, the Southern Ocean, and the North Atlantic Ocean near the Gulf Stream. In these regions the threshold method tends to overestimate winter MLDs by approximately 10% compared to the algorithm.

Centurioni, LR, Hormann V, Talley LD, Arzeno I, Beal L, Caruso M, Conry P, Echols R, Fernando HJS, Giddings SN, Gordon A, Graber H, Harcourt RR, Jayne SR, Jensen TG, Lee CM, Lermusiaux PFJ, L'Hegaret P, Lucas AJ, Mahadevan A, McClean JL, Pawlak G, Rainville L, Riser SC, Seo H, Shcherbina AY, Skyllingstad E, Sprintall J, Subrahmanyam B, Terrill E, Todd RE, Trott C, Ulloa HN, Wang H.  2017.  Northern Arabian Sea Circulation Autonomous Research (NASCar): A research initiative based on autonomous sensors. Oceanography. 30:74-87.   10.5670/oceanog.2017.224   AbstractWebsite

The Arabian Sea circulation is forced by strong monsoonal winds and is characterized by vigorous seasonally reversing currents, extreme differences in sea surface salinity, localized substantial upwelling, and widespread submesoscale thermohaline structures. Its complicated sea surface temperature patterns are important for the onset and evolution of the Asian monsoon. This article describes a program that aims to elucidate the role of upper-ocean processes and atmospheric feedbacks in setting the sea surface temperature properties of the region. The wide range of spatial and temporal scales and the difficulty of accessing much of the region with ships due to piracy motivated a novel approach based on state-of-the-art autonomous ocean sensors and platforms. The extensive data set that is being collected, combined with numerical models and remote sensing data, confirms the role of planetary waves in the reversal of the Somali Current system. These data also document the fast response of the upper equatorial ocean to monsoon winds through changes in temperature and salinity and the connectivity of the surface currents across the northern Indian Ocean. New observations of thermohaline interleaving structures and mixing in setting the surface temperature properties of the northern Arabian Sea are also discussed.

Williams, NL, Juranek LW, Feely RA, Johnson KS, Sarmiento JL, Talley LD, Dickson AG, Gray AR, Wanninkhof R, Russell JL, Riser SC, Takeshita Y.  2017.  Calculating surface ocean pCO(2) from biogeochemical Argo floats equipped with pH: An uncertainty analysis. Global Biogeochemical Cycles. 31:591-604.   10.1002/2016gb005541   AbstractWebsite

More than 74 biogeochemical profiling floats that measure water column pH, oxygen, nitrate, fluorescence, and backscattering at 10 day intervals have been deployed throughout the Southern Ocean. Calculating the surface ocean partial pressure of carbon dioxide (pCO(2sw)) from float pH has uncertainty contributions from the pH sensor, the alkalinity estimate, and carbonate system equilibrium constants, resulting in a relative standard uncertainty in pCO(2sw) of 2.7% (or 11 mu atm at pCO(2sw) of 400 mu atm). The calculated pCO(2sw) from several floats spanning a range of oceanographic regimes are compared to existing climatologies. In some locations, such as the subantarctic zone, the float data closely match the climatologies, but in the polar Antarctic zone significantly higher pCO(2sw) are calculated in the wintertime implying a greater air-sea CO2 efflux estimate. Our results based on four representative floats suggest that despite their uncertainty relative to direct measurements, the float data can be used to improve estimates for air-sea carbon flux, as well as to increase knowledge of spatial, seasonal, and interannual variability in this flux. Plain Language Summary The Southern Ocean is a key player in the global flow of carbon, yet it is hard to reach, and there are relatively few measurements there, especially in winter. Measuring the amount of carbon dioxide gas in seawater is key to advancing our understanding of the Southern Ocean. More than 74 robotic floats that use sensors to measure seawater properties have been deployed throughout the Southern Ocean, and each has a lifetime of around 5 years. It is currently not possible to directly measure carbon dioxide gas from these floats; however, it is possible to estimate carbon dioxide from things that the float can measure, like pH, a measure of ocean acidity. Here surface ocean carbon dioxide is estimated from several floats and compared to two ship-based estimates. In some locations, the floats closely match the existing estimates, but in other locations the floats see significantly higher surface ocean carbon dioxide in the wintertime, reinforcing the idea that the Southern Ocean's role in the global carbon cycle needs a closer look. Our results show that despite not measuring carbon dioxide directly, these floats will help scientists learn a lot about the Southern Ocean's part in the global flow of carbon.

Carter, BR, Feely RA, Mecking S, Cross JN, Macdonald AM, Siedlecki SA, Talley LD, Sabine CL, Millero FJ, Swift JH, Dickson AG, Rodgers KB.  2017.  Two decades of Pacific anthropogenic carbon storage and ocean acidification along Global Ocean Ship-lebased Hydrographic Investigations Program sections P16 and P02. Global Biogeochemical Cycles. 31:306-327.   10.1002/2016gb005485   AbstractWebsite

A modified version of the extended multiple linear regression (eMLR) method is used to estimate anthropogenic carbon concentration (C-anth) changes along the Pacific P02 and P16 hydrographic sections over the past two decades. P02 is a zonal section crossing the North Pacific at 30 degrees N, and P16 is a meridional section crossing the North and South Pacific at similar to 150 degrees W. The eMLR modifications allow the uncertainties associated with choices of regression parameters to be both resolved and reduced. Canth is found to have increased throughout the water column from the surface to similar to 1000 m depth along both lines in both decades. Mean column Canth inventory increased consistently during the earlier (1990s-2000s) and recent (2000s-2010s) decades along P02, at rates of 0.53 +/- 0.11 and 0.46 +/- 0.11 mol Cm-2 a(-1), respectively. By contrast, Canth storage accelerated from 0.29 +/- 0.10 to 0.45 +/- 0.11 mol Cm-2 a(-1) along P16. Shifts in water mass distributions are ruled out as a potential cause of this increase, which is instead attributed to recent increases in the ventilation of the South Pacific Subtropical Cell. Decadal changes along P16 are extrapolated across the gyre to estimate a Pacific Basin average storage between 60 degrees S and 60 degrees N of 6.1 +/- 1.5 PgC decade(-1) in the earlier decade and 8.8 +/- 2.2 PgC decade(-1) in the recent decade. This storage estimate is large despite the shallow Pacific Canth penetration due to the large volume of the Pacific Ocean. By 2014, Canth storage had changed Pacific surface seawater pH by -0.08 to -0.14 and aragonite saturation state by -0.57 to -0.82.

Billheimer, S, Talley LD.  2016.  Annual cycle and destruction of Eighteen Degree Water. Journal of Geophysical Research-Oceans. 121:6604-6617.   10.1002/2016jc011799   AbstractWebsite

Eighteen Degree Water (EDW), the subtropical mode water of the western North Atlantic, is a voluminous, weakly stratified upper ocean water mass that acts as a subsurface reservoir of heat, nutrients, and CO2. This thick layer persists throughout the year, but nearly half of its volume is dispersed or mixed away, diffusing its properties into the thermocline, from the time it outcrops in winter until it is renewed the following year. CTD observations from Argo profiling floats and acoustically tracked, isothermally bound profiling floats are used to quantify EDW destruction rates and investigate the relevant processes responsible for the large annual cycle of EDW. EDW destruction occurs primarily at the top of the EDW layer, with the highest EDW destruction rates occurring during early summer. Slower, steadier EDW destruction is observed in early winter. EDW destruction is dominated by 1-D vertical diffusion, while mesoscale, along-isopycnal stirring is also significant, explaining approximately 1/3 of the total annual EDW destruction. Destruction via along-isopycnal processes is more prevalent near the Gulf Stream than in the southern Sargasso Sea, due to higher potential vorticity gradients and enhanced mesoscale activity.

Delman, AS, Sprintall J, McClean JL, Talley LD.  2016.  Anomalous Java cooling at the initiation of positive Indian Ocean Dipole events. Journal of Geophysical Research: Oceans.   10.1002/2016JC011635   AbstractWebsite

Anomalous sea surface temperature (SST) cooling south of Java, initiated during May–July, is an important precursor to positive Indian Ocean Dipole (pIOD) events. As shown previously, the Java SST anomalies are spatially and temporally coincident with seasonal upwelling induced locally by southeasterly trade winds. However, we confirm earlier findings that interannual variability of the Java cooling is primarily driven by remote wind forcing from coastal Sumatra and the equatorial Indian Ocean (EqIO); we also find an influence from winds along the Indonesian Throughflow. The wind forcing in the EqIO and along coastal Sumatra does not initiate SST cooling locally due to a deep thermocline and thick barrier layer, but can force upwelling Kelvin waves that induce substantial surface cooling once they reach the seasonally shallower thermocline near the coast of Java. Satellite altimetry is used to obtain a Kelvin wave coefficient that approximates Kelvin wave amplitude variations along the equator. All pIOD years in the satellite record have anomalous levels of upwelling Kelvin wave activity along the equator during April–June, suggesting that upwelling waves during this season are necessary for pIOD event development. However, a change to wind-forced downwelling Kelvin waves during July–August can abruptly terminate cool Java SST anomalies and weaken the pIOD event. Upwelling Kelvin wave activity along the equator and wind stress anomalies west of Sumatra are both robust predictors of the IOD index later in the calendar year, while values of the Kelvin wave coefficient are the most reliable predictor of pIOD events specifically.

Hernandez-Guerra, A, Talley LD.  2016.  Meridional overturning transports at 30 degrees S in the Indian and Pacific Oceans in 2002-2003 and 2009. Progress in Oceanography. 146:89-120.   10.1016/j.pocean.2016.06.005   AbstractWebsite

The meridional circulation and transports at 30 degrees S in the Pacific and Indian Oceans for the years 20022003 and 2009 are compared, using GO-SHIP hydrographic section data with an inverse box model and several choices of constraints. Southward heat transport across the combined Indian-Pacific sections, reflecting net heating north of these sections, doubled from -0.7 +/- 0.2 PW in 2002-2003 to -1.4 +/- 0.1 PW in 2009 (negative sign is southward), with the increase concentrated in the Indian Ocean (-0.6 PW compared with similar to 0.2 PW in the Pacific), and was insensitive to model choices for the Indonesian Throughflow. Diagnosed net evaporation also more than doubled in the Indian Ocean, from 0.21-0.27 Sv in 2002-2003 to 0.51-0.58 in 2009, with a smaller but significant increase in net evaporation in the Pacific, from 0.06-0.08 Sv to 0.16-0.32 Sv. These increased heat and freshwater exports coincided with Indian Ocean warming, a shift in the Indian's shallow gyre overturning transport to lower densities, and an increase in southward Agulhas Current transport from 75 Sv in 2002 to 92 Sv in 2009. The Indian's deep overturn weakened from about 11 Sv in 2002 to 7 Sv in 2009. In contrast, the Pacific Ocean overturning circulation was, nearly unchanged from 2003 to 2009, independent of model within the uncertainties. The East Australian Current transport decreased only slightly, from 52 Sv to 46 Sv. The southward Pacific Deep Water transport was at a higher density than the southward Indian Deep Water transport in both years and all models, similar to prior results. Estimated diapycnal diffusivity and velocity are strongly enhanced near the ocean bottom and are higher farther up in the water column in the Indian than in the Pacific, likely extending the reach of Indian Ocean overturning up to shallower depths than in the Pacific. The horizontal distribution of transports in the Pacific at all depths changed notably from 2003 to 2009, despite the stability of its meridional overturning structure. The 2009 horizontal structure resembles a "bowed gyre"; the hydrographic section data show that this disturbance extends to the abyss and disrupts the Deep Western Boundary Current structure in the Southwest Pacific Basin. Satellite altimetry suggests association with slow westward Rossby wave propagation generated in the eastern Pacific, with no apparent effect on the net overturning circulation. The Indian Ocean's upper ocean horizontal structure was stable between the two years even though its shallow gyre overturning transports changed significantly. On the other hand, northward abyssal transports concentrated' in the central Indian Ocean (Crozet Basin) in 2002 shifted westward to the Mozambique and Madagascar Basins in 2009, although the Crozet Basin's Deep Western Boundary Current existed in both years. (C) 2016 Elsevier Ltd. All rights reserved.

Abernathey, RP, Cerovecki I, Holland PR, Newsom E, Mazlo M, Talley LD.  2016.  Water-mass transformation by sea ice in the upper branch of the Southern Ocean overturning. Nature Geoscience. 9:596-+.   10.1038/ngeo2749   AbstractWebsite

Ocean overturning circulation requires a continuous thermodynamic transformation of the buoyancy of seawater. The steeply sloping isopycnals of the Southern Ocean provide a pathway for Circumpolar Deep Water to upwell from mid depth without strong diapycnal mixing(1-3), where it is transformed directly by surface fluxes of heat and freshwater and splits into an upper and lower branch(4-6). While brine rejection from sea ice is thought to contribute to the lower branch(7), the role of sea ice in the upper branch is less well understood, partly due to a paucity of observations of sea-ice thickness and transport(8,9). Here we quantify the sea-ice freshwater flux using the Southern Ocean State Estimate, a state-of-the-art data assimilation that incorporates millions of ocean and ice observations. We then use the water-mass transformation framework(10) to compare the relative roles of atmospheric, sea-ice, and glacial freshwater fluxes, heat fluxes, and upper-ocean mixing in transforming buoyancy within the upper branch. We find that sea ice is a dominant term, with differential brine rejection and ice melt transforming upwelled Circumpolar Deep Water at a rate of similar to 22 x 10(6) m(3) s(-1). These results imply a prominent role for Antarctic sea ice in the upper branch and suggest that residual overturning and wind-driven sea-ice transport are tightly coupled.

Tamsitt, V, Talley LD, Mazloff MR, Cerovecki I.  2016.  Zonal variations in the Southern Ocean heat budget. Journal of Climate. 29:6563-6579.   10.1175/JCLI-D-15-0630.1   AbstractWebsite

The spatial structure of the upper ocean heat budget in the Antarctic Circumpolar Current (ACC) is investigated using the ⅙°, data-assimilating Southern Ocean State Estimate (SOSE) for 2005–10. The ACC circumpolar integrated budget shows that 0.27 PW of ocean heat gain from the atmosphere and 0.38 PW heat gain from divergence of geostrophic heat transport are balanced by −0.58 PW cooling by divergence of Ekman heat transport and −0.09 PW divergence of vertical heat transport. However, this circumpolar integrated balance obscures important zonal variations in the heat budget. The air–sea heat flux shows a zonally asymmetric pattern of ocean heat gain in the Indian and Atlantic sectors and ocean heat loss in the Pacific sector of the ACC. In the Atlantic and Indian sectors of the ACC, the surface ocean heat gain is primarily balanced by divergence of equatorward Ekman heat transport that cools the upper ocean. In the Pacific sector, surface ocean heat loss and cooling due to divergence of Ekman heat transport are balanced by warming due to divergence of geostrophic heat advection, which is similar to the dominant heat balance in the subtropical Agulhas Return Current. The divergence of horizontal and vertical eddy advection of heat is important for warming the upper ocean close to major topographic features, while the divergence of mean vertical heat advection is a weak cooling term. The results herein show that topographic steering and zonal asymmetry in air–sea exchange lead to substantial zonal asymmetries in the heat budget, which is important for understanding the upper cell of the overturning circulation.

Williams, NL, Juranek LW, Johnson KS, Feely RA, Riser SC, Talley LD, Russell JL, Sarmiento JL, Wanninkhof R.  2016.  Empirical algorithms to estimate water column pH in the Southern Ocean. Geophysical Research Letters. 43:3415-3422.   10.1002/2016gl068539   AbstractWebsite

Empirical algorithms are developed using high-quality GO-SHIP hydrographic measurements of commonly measured parameters (temperature, salinity, pressure, nitrate, and oxygen) that estimate pH in the Pacific sector of the Southern Ocean. The coefficients of determination, R-2, are 0.98 for pH from nitrate (pH(N)) and 0.97 for pH from oxygen (pH(Ox)) with RMS errors of 0.010 and 0.008, respectively. These algorithms are applied to Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) biogeochemical profiling floats, which include novel sensors (pH, nitrate, oxygen, fluorescence, and backscatter). These algorithms are used to estimate pH on floats with no pH sensors and to validate and adjust pH sensor data from floats with pH sensors. The adjusted float data provide, for the first time, seasonal cycles in surface pH on weekly resolution that range from 0.05 to 0.08 on weekly resolution for the Pacific sector of the Southern Ocean.

Talley, LD, Feely RA, Sloyan BM, Wanninkhof R, Baringer MO, Bullister JL, Carlson CA, Doney SC, Fine RA, Firing E, Gruber N, Hansell DA, Ishii M, Johnson GC, Katsumata K, Key RM, Kramp M, Langdon C, Macdonald AM, Mathis JT, McDonagh EL, Mecking S, Millero FJ, Mordy CW, Nakano T, Sabine CL, Smethie WM, Swift JH, Tanhua T, Thurnherr AM, Warner MJ, Zhang J-Z.  2016.  Changes in Ocean Heat, Carbon Content, and Ventilation: A Review of the First Decade of GO-SHIP Global Repeat Hydrography. Annual Review of Marine Science. 8:185-215.   10.1146/annurev-marine-052915-100829   AbstractWebsite

Global ship-based programs, with highly accurate, full water column physical and biogeochemical observations repeated decadally since the 1970s, provide a crucial resource for documenting ocean change. The ocean, a central component of Earth's climate system, is taking up most of Earth's excess anthropogenic heat, with about 19% of this excess in the abyssal ocean beneath 2,000 m, dominated by Southern Ocean warming. The ocean also has taken up about 27% of anthropogenic carbon, resulting in acidification of the upper ocean. Increased stratification has resulted in a decline in oxygen and increase in nutrients in the Northern Hemisphere thermocline and an expansion of tropical oxygen minimum zones. Southern Hemisphere thermocline oxygen increased in the 2000s owing to stronger wind forcingand ventilation. The most recent decade of global hydrography has mapped dissolved organic carbon, a large, bioactive reservoir, for the first time and quantified its contribution to export production (∼20%) and deep-ocean oxygen utilization. Ship-based measurements also show that vertical diffusivity increases from a minimum in the thermocline to a maximum within the bottom 1,500 m, shifting our physical paradigm of the ocean's overturning circulation.

Williams, NL, Feely RA, Sabine CL, Dickson AG, Swift JH, Talley LD, Russell JL.  2015.  Quantifying anthropogenic carbon inventory changes in the Pacific sector of the Southern Ocean. Marine Chemistry. 174:147-160.   10.1016/j.marchem.2015.06.015   AbstractWebsite

The Southern Ocean plays a major role in mediating the uptake, transport, and long-term storage of anthropogenic carbon dioxide (CO2) into the deep ocean. Examining the magnitude and spatial distribution of this oceanic carbon uptake is critical to understanding how the earth's carbon system will react to continued increases in this greenhouse gas. Here, we use the extended multiple linear regression technique to quantify the total and anthropogenic change in dissolved inorganic carbon (DIC) along the S04P and P16S CLIVAR/U.S. Global Ocean Carbon and Repeat Hydrography Program lines south of 67 degrees S in the Pacific sector of the Southern Ocean between 1992 and 2011 using discrete bottle measurements from repeat occupations. Along the S04P section, which is located in the seasonal sea ice zone south of the Antarctic Circumpolar Current in the Pacific, the anthropogenic component of the DIC increase from 1992 to 2011 is mostly found in the Antarctic Surface Water (AASW, upper 100 m), while the increase in DIC below the mixed layer in the Circumpolar Deep Water can be primarily attributed to either a slowdown in circulation or decreased ventilation of deeper, high CO2 waters. In the AASW we calculate an anthropogenic increase in DIC of 12-18 mu mol kg(-1) and an average storage rate of anthropogenic CO2 of 0.10 +/- 0.02 mol m(-2) yr(-1) for this region compared to a global average of 0.5 +/- 0.2 mol m(-2) yr(-1). In surface waters this anthropogenic CO2 uptake results in an average pH decrease of 0.0022 +/- 0.0004 pH units yr(-1), a 0.47 +/- 0.10% yr(-1) decrease in the saturation state of aragonite (Omega(Aragonite)) and a 2.0 +/- 0.7 m yr(-1) shoaling of the aragonite saturation horizons (calculated for the Omega(Aragonite) = 1.3 contour). (C) 2015 Published by Elsevier B.V.

Delman, AS, McClean JL, Sprintall J, Talley LD, Yulaeva E, Jayne SR.  2015.  Effects of eddy vorticity forcing on the mean state of the Kuroshio Extension. Journal of Physical Oceanography. 45:1356-1375.   10.1175/jpo-d-13-0259.1   AbstractWebsite

Eddy-mean flow interactions along the Kuroshio Extension (KE) jet are investigated using a vorticity budget of a high-resolution ocean model simulation, averaged over a 13-yr period. The simulation explicitly resolves mesoscale eddies in the KE and is forced with air-sea fluxes representing the years 1995-2007. A mean-eddy decomposition in a jet-following coordinate system removes the variability of the jet path from the eddy components of velocity; thus, eddy kinetic energy in the jet reference frame is substantially lower than in geographic coordinates and exhibits a cross-jet asymmetry that is consistent with the baroclinic instability criterion of the long-term mean field. The vorticity budget is computed in both geographic (i. e., Eulerian) and jet reference frames; the jet frame budget reveals several patterns of eddy forcing that are largely attributed to varicose modes of variability. Eddies tend to diffuse the relative vorticity minima/maxima that flank the jet, removing momentum from the fast-moving jet core and reinforcing the quasi-permanent meridional meanders in the mean jet. A pattern associated with the vertical stretching of relative vorticity in eddies indicates a deceleration (acceleration) of the jet coincident with northward (southward) quasi-permanent meanders. Eddy relative vorticity advection outside of the eastward jet core is balanced mostly by vertical stretching of the mean flow, which through baroclinic adjustment helps to drive the flanking recirculation gyres. The jet frame vorticity budget presents a well-defined picture of eddy activity, illustrating along-jet variations in eddy-mean flow interaction that may have implications for the jet's dynamics and cross-frontal tracer fluxes.

Whalen, CB, MacKinnon JA, Talley LD, Waterhouse AF.  2015.  Estimating the mean diapycnal mixing using a finescale strain parameterization. Journal of Physical Oceanography. 45:1174-1188.   10.1175/jpo-d-14-0167.1   AbstractWebsite

Finescale methods are currently being applied to estimate the mean turbulent dissipation rate and diffusivity on regional and global scales. This study evaluates finescale estimates derived from isopycnal strain by comparing them with average microstructure profiles from six diverse environments including the equator, above ridges, near seamounts, and in strong currents. The finescale strain estimates are derived from at least 10 nearby Argo profiles (generally <60 km distant) with no temporal restrictions, including measurements separated by seasons or decades. The absence of temporal limits is reasonable in these cases, since the authors find the dissipation rate is steady over seasonal time scales at the latitudes being considered (0 degrees-30 degrees and 40 degrees-50 degrees). In contrast, a seasonal cycle of a factor of 2-5 in the upper 1000m is found under storm tracks (30 degrees-40 degrees) in both hemispheres. Agreement between the mean dissipation rate calculated using Argo profiles and mean from microstructure profiles is within a factor of 2-3 for 96% of the comparisons. This is both congruous with the physical scaling underlying the finescale parameterization and indicates that the method is effective for estimating the regional mean dissipation rates in the open ocean.

Oka, E, Uehara K, Nakano T, Suga T, Yanagimoto D, Kouketsu S, Itoh S, Katsura S, Talley LD.  2014.  Synoptic observation of Central Mode Water in its formation region in spring 2003. Journal of Oceanography. 70:521-534.   10.1007/s10872-014-0248-2   AbstractWebsite

Hydrographic data east of Japan from five research cruises and Argo profiling floats in spring 2003 have been analyzed to examine the relationship of the formation of Central Mode Water (CMW) and Transition Region Mode Water (TRMW) in late winter 2003 to thermohaline fronts and mesoscale eddies. TRMW and the denser variety of CMW (D-CMW) were formed continuously just south of the subarctic frontal zone between 155 degrees E and 165 degrees W with little relation to eddies, suggesting that the absence of the permanent thermocline and halocline in this area is essential for the formation. The lighter variety of CMW (L-CMW) was formed south of the Kuroshio bifurcation front and east of 165 degrees E, partly in an anticyclonic eddy associated with the Kuroshio Extension. Some portion of D-CMW and L-CMW likely had been subducted to the permanent pycnocline by crossing southward the Kuroshio bifurcation front and the Kuroshio Extension front, respectively. In contrast, the formation of these waters in the western regions was inactive and was significantly different from that described previously using multiyear Argo float data. West of 155 degrees E, TRMW and D-CMW were formed only in two anticyclonic eddies that had been detached from the Kuroshio Extension 1-2 years ago. L-CMW was hardly formed west of 165 degrees E, which might be related to the upstream Kuroshio Extension being in its stable state characterized by low regional eddy activity.

Waterhouse, AF, MacKinnon JA, Nash JD, Alford MH, Kunze E, Simmons HL, Polzin KL, St Laurent LC, Sun OM, Pinkel R, Talley LD, Whalen CB, Huussen TN, Carter GS, Fer I, Waterman S, Garabato ACN, Sanford TB, Lee CM.  2014.  Global patterns of diapycnal mixing from measurements of the turbulent dissipation rate. Journal of Physical Oceanography. 44:1854-1872.   10.1175/jpo-d-13-0104.1   AbstractWebsite

The authors present inferences of diapycnal diffusivity from a compilation of over 5200 microstructure profiles. As microstructure observations are sparse, these are supplemented with indirect measurements of mixing obtained from(i) Thorpe-scale overturns from moored profilers, a finescale parameterization applied to (ii) shipboard observations of upper-ocean shear, (iii) strain as measured by profiling floats, and (iv) shear and strain from full-depth lowered acoustic Doppler current profilers (LADCP) and CTD profiles. Vertical profiles of the turbulent dissipation rate are bottom enhanced over rough topography and abrupt, isolated ridges. The geography of depth-integrated dissipation rate shows spatial variability related to internal wave generation, suggesting one direct energy pathway to turbulence. The global-averaged diapycnal diffusivity below 1000-m depth is O(10(-4))m(2) s(-1) and above 1000-m depth is O(10(-5))m(2) s(-1). The compiled microstructure observations sample a wide range of internal wave power inputs and topographic roughness, providing a dataset with which to estimate a representative global-averaged dissipation rate and diffusivity. However, there is strong regional variability in the ratio between local internal wave generation and local dissipation. In some regions, the depth-integrated dissipation rate is comparable to the estimated power input into the local internal wave field. In a few cases, more internal wave power is dissipated than locally generated, suggesting remote internal wave sources. However, at most locations the total power lost through turbulent dissipation is less than the input into the local internal wave field. This suggests dissipation elsewhere, such as continental margins.

Carter, BR, Talley LD, Dickson AG.  2014.  Mixing and remineralization in waters detrained from the surface into Subantarctic Mode Water and Antarctic Intermediate Water in the southeastern Pacific. Journal of Geophysical Research-Oceans. 119:4001-4028.   10.1002/2013jc009355   AbstractWebsite

A hydrographic data set collected in the region and season of Subantarctic Mode Water and Antarctic Intermediate Water (SAMW and AAIW) formation in the southeastern Pacific allows us to estimate the preformed properties of surface water detrained into these water masses from deep mixed layers north of the Subantarctic Front and Antarctic Surface Water south of the front. Using 10 measured seawater properties, we estimate: the fractions of SAMW/AAIW that originate as surface source waters, as well as fractions that mix into these water masses from subtropical thermocline water above and Upper Circumpolar Deep Water below the subducted SAMW/AAIW; ages associated with the detrained surface water; and remineralization and dissolution rates and ratios. The mixing patterns imply that cabbeling can account for similar to 0.005-0.03 kg m(-3) of additional density in AAIW, and similar to 0-0.02 kg m(-3) in SAMW. We estimate a shallow depth (similar to 300-700 m, above the aragonite saturation horizon) calcium carbonate dissolution rate of 0.4 +/- 0.2 mmol CaCO3 kg(-1) yr(-1), a phosphate remineralization rate of 0.031 +/- 0.009 mu mol P kg(-1) yr(-1), and remineralization ratios of P:N:-O-2:C-org of 1:(15.5 +/- 0.6):(143 +/- 10):(104 +/- 22) for SAMW/AAIW. Our shallow depth calcium carbonate dissolution rate is comparable to previous estimates for our region. Our -O-2:P ratio is smaller than many global averages. Our model suggests neglecting diapycnal mixing of preformed phosphate has likely biased previous estimates of -O-2:P and C-org:P high, but that the C-org:P ratio bias may have been counteracted by a second bias in previous studies from neglecting anthropogenic carbon gradients.

Talley, LD.  2013.  Closure of the Global Overturning Circulation Through the Indian, Pacific, and Southern Oceans: Schematics and Transports. Oceanography. 26:80-97. AbstractWebsite

The overturning pathways for the surface-ventilated North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW) and the diffusively formed Indian Deep Water (IDW) and Pacific Deep Water (PDW) are intertwined. The global overturning circulation (GOC) includes both large wind-driven upwelling in the Southern Ocean and important internal diapycnal transformation in the deep Indian and Pacific Oceans. All three northern-source Deep Waters (NADW, IDW, PDW) move southward and upwell in the Southern Ocean. AABW is produced from the denser, salty NADW and a portion of the lighter, low oxygen IDW/PDW that upwells above and north of NADW. The remaining upwelled IDW/PDW stays near the surface, moving into the subtropical thermoclines, and ultimately sources about one-third of the NADW. Another third of the NADW comes from AABW upwelling in the Atlantic. The remaining third comes from AABW upwelling to the thermocline in the Indian-Pacific. Atlantic cooling associated with NADW formation (0.3 PW north of 32 degrees S; 1 PW = 1015 W) and Southern Ocean cooling associated with AABW formation (0.4 PW south of 32 degrees S) are balanced mostly by 0.6 PW of deep diffusive heating in the Indian and Pacific Oceans; only 0.1 PW is gained at the surface in the Southern Ocean. Thus, while an adiabatic model of NADW global overturning driven by winds in the Southern Ocean, with buoyancy added only at the surface in the Southern Ocean, is a useful dynamical idealization, the associated heat changes require full participation of the diffusive Indian and Pacific Oceans, with a basin-averaged diffusivity on the order of the Munk value of 10(-4) m(2) s(-1).