Temperature-salinity Relationships Research Articles

Technical note: An autonomous flow-through salinity and temperature perturbation mesocosm system for multi-stressor experiments

Abstract. The rapid environmental changes in aquatic systems as a result of anthropogenic forcings are creating a multitude of challenging conditions for organisms and communities. The need to better understand the interaction of environmental stressors now, and in the future, is fundamental to determining the response of ecosystems to these perturbations. This work describes an automated ex situ mesocosm perturbation system that can manipulate several variables of aquatic media in a controlled setting. This perturbation system was deployed in Kongsfjorden (Svalbard); within this system, ambient water from the fjord was heated and mixed with freshwater in a multifactorial design to investigate the response of mixed-kelp communities in mesocosms to projected future Arctic conditions. The system employed an automated dynamic offset scenario in which a nominal temperature increase was programmed as a set value above real-time ambient conditions in order to simulate future warming. A freshening component was applied in a similar manner: a decrease in salinity was coupled to track the temperature offset based on a temperature–salinity relationship in the fjord. The system functioned as an automated mixing manifold that adjusted flow rates of warmed and chilled ambient seawater, with unmanipulated ambient seawater and freshwater delivered as a single source of mixed media to individual mesocosms. These conditions were maintained via continuously measured temperature and salinity in 12 mesocosms (1 control and 3 treatments, all in triplicate) for 54 d. System regulation was robust, as median deviations from nominal conditions were < 0.15 for both temperature (∘C) and salinity across the three replicates per treatment. Regulation further improved during a second deployment that mimicked three marine heat wave scenarios in which a dynamic temperature regulation held median deviations to < 0.036 ∘C from the nominal value for all treatment conditions and replicates. This perturbation system has the potential to be implemented across a wide range of conditions to test single or multi-stressor drivers (e.g., increased temperature, freshening, and high CO2) while maintaining natural variability. The automated and independent control for each experimental unit (if desired) provides a large breadth of versatility with respect to experimental design.

Biogeochemical structure of the Laptev Sea in 2015-2020 associated with the River Lena plume

The discharge of rivers and the subsequent dispersion of their plumes play a pivotal role in the biogeochemical cycling of the Arctic Ocean. Based on the data collected during annual transects conducted in the autumn period (September-October) from 2015-2020, this study explores the effect of River Lena plume dispersion on the seasonal and interannual changes in the hydrophysical and biogeochemical structure of the southeastern Laptev Sea. The temperature-salinity relationship (T-S), Redfield ratio and multiparameter cluster analysis were used to investigate variations in the water mass structure along the transect. The results revealed that the plume’s interannual and seasonal spreading patterns play a crucial role in regulating the local physical, biogeochemical, and biological processes in the southern Laptev Sea. During September-October, the hydrochemical water mass structure along the transects shifted from highly stratified to unstratified as the plume’s mixing intensity increased. Anomalous hydrochemical distributions were observed due to coastal upwelling, which was primarily characterized by high total alkalinity and nitrate levels, and low organic phosphorus, nitrite, and ammonia levels in the seawater. Wind and cold weather conditions drive deep vertical mixing of seawater, causing the resuspension of bottom sediment and the subsequent enrichment of bottom water by nutrients. Multi-parameter cluster analysis is used to describe the details of water mass structures in the highly dynamic southern Laptev Sea, with water mass structures typically undergoing significant changes within two weeks between September and October. The migration and transformation of water masses throughout the seasons are influenced by the volume of river discharge, fall-winter cooling, and atmospheric circulation patterns. Furthermore, the general atmospheric circulation is confirmed to be the primary cause of the interannual variation in the spread of the Lena River plume over the southeast Laptev Sea.

Cyclonic eddies and upper thermocline fine-scale structures in the Antarctic Circumpolar Current

Mesoscale eddies in the open ocean are mostly formed by baroclinic instability, in which the available potential energy from the large-scale slope of the isopycnals is converted into the kinetic energy of the flow around the eddy. As a permissible form of motion within a rapidly rotating and stratified fluid eddies driven by baroclinic instability are important for the poleward and vertical transport, not only of physical properties, but also biogeochemical ones. In this paper, we present observations from four cyclonic eddies in the Antarctic Circumpolar Current. We have sorted them by apparent age, based on altimeter data and consideration of the degree of homogenisation of the potential temperature-salinity(\U0001d703S) relationship, and then looked at the spatial distribution of measures of fine-scale variability in the upper thermocline. The youngest eddy shows isopycnals which are domed upwards and it contains a variety of waters with differing temperature-salinity characteristics. The fine-scale variability is higher in the core of the eddy. The older eddies show a core which is more homogeneous in potential temperature and salinity. The isopycnals are flatter in the centre of the eddy, and in cross-section, they can be M-shaped, so that the steepest gradients are concentrated around the edge. The fine-scale variability is more concentrated around the edges where the density gradients are stronger. We hypothesise that lateral stirring and mixing processes within the eddy homogenise the water so that the temperature-salinity relationship becomes tighter. When the eddy eventually collapses, this modified water can be released back into the flow. Thus, we see how the interplay of mesoscale and small-scale processes are modifying water mass properties and, potentially, regulate biogeochemical processes.

Biological and physical controls on the oxygen cycle in the Kuroshio Extension from an array of profiling floats

The Kuroshio Extension (KE) is a current associated with the largest CO2 flux into the Pacific Ocean, a broad region of uptake that extends across the Pacific basin between the subarctic and subtropical regions. The relative importance of the biological and physical processes controlling this sink is uncertain. Because oxygen is stoichiometrically linked to changes in dissolved inorganic carbon due to photosynthesis and respiration and subject to many of the same physical drivers as the CO2 flux, in situ oxygen measurements help to determine the processes driving this large CO2 flux. We analyzed data from eight Argo profiling floats equipped with oxygen sensors to estimate oxygen fluxes in the upper ocean of the KE region (approximate bounds: 25°N to 45°N, 135°E to 155°E). In situ air calibrations of these sensors allowed us to accurately measure air-sea oxygen differences, which largely control the flux of oxygen to and from the atmosphere. To characterize distinct biogeographical regions in the Kuroshio Extension and to accommodate seasonal north-south shifts in the location of the regional boundaries, we averaged oxygen measurements from different floats along isopycnal surfaces into 3 regions based on temperature-salinity relationships: North KE, Central KE, and South KE. Using these regional concentration time series, we determined the physical fluxes using an upper ocean layered model and calculated the residual oxygen flux, a combination of column-integrated net community production and physical processes unexplained by the model. The annual oxygen budget is largely a balance of air-sea exchange and the residual oxygen term. Residual oxygen fluxes are -5.4 ± 1.1, -5.9 ± 0.1, and 2.1 ± 2.2 mol O2 m-2 yr-1 (where negative is a loss from the upper ocean) for the North, Central, and South KE regions, respectively. The North and Central KE are regions of mode water formation, which balances the large air-sea fluxes into the ocean. The South KE oxygen residual indicates a biologically produced flux to the atmosphere in two out of three years that agrees with previous estimates of subtropical annual net community production (ANCP) but exhibits high interannual variability. This study suggests that physical processes are the primary drivers for the annual uptake of the gases oxygen and carbon dioxide in the Central KE region where the annual CO2 uptake is strongest.

An Updated Estimate of Salinity for the Atlantic Ocean Sector Using Temperature–Salinity Relationships

AbstractSimultaneous temperature and salinity profile measurements are of extreme importance for research; operational oceanography; research and applications that compute content and transport of mass, heat, and freshwater in the ocean; and for determining water mass stratification and mixing rates. Historically, temperature profiles are much more abundant than simultaneous temperature and salinity profiles. Given the importance of concurrent temperature and salinity profiles, several methods have been developed to derive salinity solely based on temperature profile observations, such as expendable bathythermograph (XBT) temperature measurements, for which concurrent salinity observations are typically not available. These empirical methods used to date contain uncertainties as a result of temporal changes in salinity and seasonality in the mixed layer, and are typically regionally based. In this study, a new methodology is proposed to infer salinity in the Atlantic Ocean from the water surface to 2000-m depth, which addresses the seasonality in the upper ocean and makes inferences about longer-term changes in salinity. Our results show that when seasonality is accounted for, the variance of the residuals is reduced in the upper 150 m of the ocean and the dynamic height errors are smaller than 4 cm in the whole study domain. The sensitivity of the meridional heat and freshwater transport to different empirical methods of salinity estimation is studied using the high-density XBT transect across 34.5°S in the South Atlantic Ocean. Results show that accurate salinity estimates are more important on the boundaries, suggesting that temperature–salinity compensation may be also important in those regions.

Submesoscale Dynamics in the Northern Gulf of Mexico. Part II: Temperature–Salinity Relations and Cross-Shelf Transport Processes

AbstractThis paper, the second of three, investigates submesoscale dynamics in the northern Gulf of Mexico under the influence of the Mississippi–Atchafalaya River system, using numerical simulations at 500-m horizontal resolution with climatological atmospheric forcing. The Turner angle Tu, a measure of the relative effect of temperature and salinity on density, is examined with respect to submesoscale current generation in runs with and without riverine forcing. Surface Tu probability density functions in solutions including rivers show a temperature-dominated signal offshore, associated with Loop Current water, and a nearshore salinity-dominated signal, associated with fresh river water, without a clear compensating signal, as often found instead in the ocean’s mixed layer. The corresponding probability distribution functions in the absence of rivers differ, illustrating the key role played by the freshwater output in determining temperature–salinity distributions in the northern Gulf of Mexico during both winter and summer. A quantity referred to as temperature–salinity covariance is proposed to determine what fraction of the available potential energy that is released during the generation of submesoscale circulations leads to the destruction of density gradients while leaving spice gradients untouched, thereby leading to compensation. It is shown that the fresh river fronts to the east of the Bird’s Foot can evolve toward compensation in concert with a gradual release of available potential energy. It is further demonstrated that, during winter, the cross-shelf freshwater transport mechanism to the west of the Bird’s Foot is well approximated by a diffusive process, whereas to the east is better represented by a ballistic process associated with the Mississippi water that converges in a jetlike pattern.

A Microscale View of Mixing and Overturning across the Antarctic Circumpolar Current

AbstractThe relative roles of isoneutral stirring by mesoscale eddies and dianeutral stirring by small-scale turbulence in setting the large-scale temperature–salinity relation of the Southern Ocean against the action of the overturning circulation are assessed by analyzing a set of shear and temperature microstructure measurements across Drake Passage in a “triple decomposition” framework. It is shown that a picture of mixing and overturning across a region of the Antarctic Circumpolar Current (ACC) may be constructed from a relatively modest number of microstructure profiles. The rates of isoneutral and dianeutral stirring are found to exhibit distinct, characteristic, and abrupt variations: most notably, a one to two orders of magnitude suppression of isoneutral stirring in the upper kilometer of the ACC frontal jets and an order of magnitude intensification of dianeutral stirring in the subpycnocline and deepest layers of the ACC. These variations balance an overturning circulation with meridional flows of O(1) mm s−1 across the ACC’s mean thermohaline structure. Isoneutral and dianeutral stirring play complementary roles in balancing the overturning, with isoneutral processes dominating in intermediate waters and the Upper Circumpolar Deep Water and dianeutral processes prevailing in lighter and denser layers.

Chemical evidence for the origin of the cold water belt along the northeastern coast of Hokkaido

In the southwestern Okhotsk Sea, the cold water belt (CWB) is frequently observed on satellite images offshore of the Soya Warm Current flowing along the northeastern coast of Hokkaido, Japan, during summertime. It has been speculated that the CWB is upwelling cold water that originates from either subsurface water of the Japan Sea off Sakhalin or bottom water of the Okhotsk Sea. Hydrographic and chemical observations (nutrients, humic-type fluorescence intensity, and iron) were conducted in the northern Japan Sea and southwestern Okhotsk Sea in early summer 2011 to clarify the origin of the CWB. Temperature–salinity relationships, vertical distributions of chemical components, profiles of chemical components against density, and the (NO3 + NO2)/PO4 relationship confirm that water in the CWB predominantly originates from Japan Sea subsurface water.

Thermodynamic modeling of ternary CH4[sbnd]H2O[sbnd]NaCl fluid inclusions

This paper reports the application of thermodynamic models, including equations of state, to ternary CH4H2ONaCl fluid inclusions. A simple equation describing pressure–temperature–salinity relations on the CH4 hydrate-liquid-vapor surface has been developed to calculate the NaCl contents (salinities) of inclusions, where the dissociation pressure of CH4 hydrate coexisting with vapor and liquid at a given temperature is calculated with a pressure equation of pure CH4. The pressure equation is a function of temperature and CH4 Raman peak position shift corrected by Ne lamp. With these relations and the latest CH4 solubility and PVTx models, a new iterative approach is presented to calculate the CH4 contents of CH4H2ONaCl inclusions on the assumption that the bulk molar volume of an inclusion at the melting temperature of CH4 hydrate and at the vapor bubble disappearance (homogenization) temperature are identical. A prominent merit of this method is that the compositions, molar volumes and homogenization pressures of CH4H2ONaCl inclusions can be simultaneously obtained without having to use volume fractions of vapor bubbles at the dissociation temperatures of CH4 hydrates determined based on optical observations or measurements. The homogenization pressures and isochores of CH4H2ONaCl fluid inclusions from updated models are briefly discussed. The code to estimate PVTx properties of inclusions in the ternary system CH4H2ONaCl, based on microthermometric and Raman data, can be obtained from Chemical Geology or the corresponding author (maoshide@163.com).

Assessing the Application of Argo Profiling Float Data to the Study of the Seasonal Variation of the Hydrological Parameters and the Current Field East of Luzon Strait

The seasonal characteristics of water masses, mixed-layer depth and the current field east of Luzon Strait are studied using Argo profiling float data. Based on data from March 2006 to November 2010, the temperature-salinity relation indicates that the seasonal variation of water masses is not obvious, except that the surface temperature is lower and the salinity is higher in spring and winter; the water masses in other layers change little. The seasonal variation of the mixed-layer depth, which is deepest at over 170 m in winter and shallowest at about 25 m in summer, is also discussed. Based on the P-vector method, multi-year seasonal mean Argo data and Levitus data are used to calculate the respective current fields. The surface geostrophic current derived from altimetric data is compared with those computed from Argo data and Levitus data. It is shown that, the current field computed from Argo data, which shows the Kuroshio and eddies clearly, is similar to that obtained from altimetric data, except that the current speed is less in the former case than in the latter one; this may be related to the scarcity and uneven distribution of Argo floats, the subsequent interpolation and extrapolation errors, the shortcomings of the P-vector method, the lack of a barotropic component in the calculation or the altimetry-derived current error associated with the marine geoid, but the resultant current field is much better than that obtained from Levitus data. Moreover, Argo data have an advantage in that they can be used to derive a three-dimensional current field, whereas the altimetric data can only be used to derive the surface current field. The seasonal variations in the vertical current structure and the Kuroshio transport east of Taiwan Island are also discussed; it is found that the thickness of the northward Kuroshio is about 700 m, and the seasonal zonal variation in the main axis of the Kuroshio is obvious. RÉSUMÉ [Traduit par la rédaction] Nous étudions les caractéristiques saisonnières des masses d'eau, la profondeur de la couche de mélange et le champ de courant à l'est du détroit de Luçon à l'aide des données des flotteurs profileurs Argo. D'après les données recueillies entre mars 2006 et novembre 2010, la relation température-salinité indique que la variation saisonnière des masses d'eau n'est pas évidente, sauf que la température de surface est plus basse et la salinité plus élevée au printemps et en hiver; les masses d'eau dans les autres couches changent peu. Nous discutons aussi de la variation saisonnière de la profondeur de la couche de mélange, qui est maximale à plus de 170 m en hiver et minimale à environ 25 m en été. En nous basant sur la méthode des vecteurs P, nous utilisons les données Argo et les données Levitus pluriannuelles de moyennes saisonnières pour calculer les champs de courant respectifs. Nous comparons le courant géostrophique de surface déduit des données altimétriques avec ceux calculés à l'aide des données Argo et des données Levitus. Il apparait que le champ de courant calculé d'après les données Argo, qui montre clairement le Kuroshio et les remous, est semblable à celui obtenu des données altimétriques, sauf que la vitesse du courant est moindre dans le premier cas que dans le dernier; cela peut être lié à la rareté et à la distribution inégale des flotteurs Argo, aux erreurs subséquentes d'interpolation et d'extrapolation, aux faiblesses de la méthode des vecteurs P, au manque d'une composante barotropique dans les calculs ou à l'erreur dans le courant déduit des données altimétriques associée au géoïde marin, mais le champ de courant résultant est bien meilleur que celui obtenu à l'aide des données Levitus. De plus, les données Argo présentent l'avantage de pouvoir être utilisées pour déduire un champ de courant tridimensionnel, alors que les données altimétriques ne permettent que de déduire le champ de courant de surface. Nous discutons aussi des variations saisonnières dans la structure verticale du courant et dans le transport par le Kuroshio à l'est de l’île de Taïwan; nous trouvons que la profondeur du Kuroshio circulant vers le nord est d'environ 700 m et la variation saisonnière zonale dans l'axe principal du Kuroshio est évidente.

Variational assimilation of HF radar surface currents in a coastal ocean model off Oregon

The impact of assimilation of sea surface velocity fields observed by a set of high-frequency (HF) radars is studied using a three-dimensional ocean circulation model configured along the Oregon coast. The study period is June–July 2008 featuring upwelling and separation of the coastal currents into the adjacent interior ocean. The nonlinear model is based on the Regional Ocean Modeling System (ROMS) and the data assimilation (DA) component on the AVRORA system utilizing the representer-based variational algorithm. Assimilation proceeds in a series of 3-day windows, providing an analysis solution in each window and a 3-day forecast into the next window. Experiments with two different initial condition error covariances are compared (one is dynamically balanced, based on the linearized equation of state, temperature-salinity relation, and geostrophic and thermal wind balance relations and the other is multivariate unbalanced). While the assimilation impact is statistically better in the case of the balanced covariance, the case with the unbalanced covariance also provides sensible improvement in terms of surface velocity and sea surface temperature (SST) model-data forecast statistics. The analysis of representer functions shows that even if the initial condition error covariance is unbalanced, the correction fields at the model initial time are partially balanced after each dynamical field is smoothed independently, due to inherent dynamical properties of the adjoint model. Assimilation of the HF radar surface currents improves not only surface velocity forecasts, but also geometry of the upwelling SST front and the sea surface height (SSH) slope near the coast, as verified against unassimilated satellite SSH and SST data. The assimilation also alters the latitudinal distribution of the time-averaged offshore transport. Combined HF radar velocity and other observations, e.g., altimetry, is needed to better constrain surface geostrophic currents in the entire model domain, including the area not covered by the HF radar data.

Several Important Issues in Salinity Quality Control of Argo Float

Several floats that deployed by different countries during the last years over the global ocean, which served as a contribution to the Array for real time geostrophic oceanography (Argo), are studied in this paper with respect to their uniqueness in the Delayed Mode Quality Control (DMQC) process. Of these floats, issues were found in the temperature-salinity relationships during the DMQC, which would easily be misinterpreted or overlooked. This study lists and summarizes the possible situations of the floats that might confront in their observation life, together with the reasonable explanations. The authors present some thoughts on evaluating the performance of the conductivity sensors and calibrating the float salinity data.

Fifty-Year Trends in Global Ocean Salinities and Their Relationship to Broad-Scale Warming

Abstract Using over 1.6 million profiles of salinity, potential temperature, and neutral density from historical archives and the international Argo Program, this study develops the three-dimensional field of multidecadal linear change for ocean-state properties. The period of analysis extends from 1950 to 2008, taking care to minimize the aliasing associated with the seasonal and major global El Niño–Southern Oscillation modes. Large, robust, and spatially coherent multidecadal linear trends in salinity to 2000-dbar depth are found. Salinity increases at the sea surface are found in evaporation-dominated regions and freshening in precipitation-dominated regions, with the spatial pattern of change strongly resembling that of the mean salinity field, consistent with an amplification of the global hydrological cycle. Subsurface salinity changes on pressure surfaces are attributable to both isopycnal heave and real water-mass modification of the temperature–salinity relationship. Subduction and circulation by the ocean’s mean flow of surface salinity and temperature anomalies appear to account for most regional subsurface salinity changes on isopycnals. Broad-scale surface warming and the associated poleward migration of isopycnal outcrops drive a clear and repeating pattern of subsurface isopycnal salinity change in each independent ocean basin. Qualitatively, the observed global multidecadal salinity changes are thus consonant with both broad-scale surface warming and the amplification of the global hydrological cycle.

Subduction on the Northern and Southern Flanks of the Gulf Stream

Abstract Sections of temperature, salinity, dissolved oxygen, and velocity were made crossing the Gulf Stream in late January 2006 to investigate the role of frontal processes in the formation of Eighteen Degree Water (EDW), the Subtropical Mode Water of the North Atlantic. The sections were nominally perpendicular to the stream and measured in a Lagrangian frame by following a floating spar buoy drifting in the Gulf Stream’s warm core. During the survey, EDW was isolated from the mixed layer by the stratified seasonal pycnocline, suggesting that EDW was not yet actively being formed at this time in the season and at the longitudes over which the survey was conducted (64°–70°W). However, in two of the sections, the seasonal pycnocline in the core of the Gulf Stream was broken by an intrusion of cold, fresh, weakly stratified water, nearly saturated in oxygen, that appears to have been subducted from the surface mixed layer north of the stream. The intrusion was identified in three of the sections in profiles with a nearly identical temperature–salinity relation. From the western-to-easternmost sections, where the intrusion was observed, the depth of the intrusion’s salinity minimum descended by ∼90 m in the 71 h it took to complete this part of the survey. This apparent subduction occurred primarily on the upstream side of a meander trough, where the cross-stream velocity was confluent and frontogenetic. Using a variant of the omega equation, the vertical velocity driven by the confluent flow was inferred and yielded downwelling in the vicinity of the intrusion spanning 10–40 m day−1, a range of values consistent with the intrusion’s observed descent, suggesting that frontal subduction was responsible for the formation of the intrusion. In the easternmost section located downstream of the meander trough, the flow was diffluent, driving an inferred vertical circulation that was of the opposite sense to that in the section upstream of the trough. In transiting the two sides of the trough, the intrusion was observed to move toward the center of the stream between the downwelling branches of the opposing vertical circulations, resulting in a downward Lagrangian mean vertical velocity and net subduction. Hydrographic evidence of the subduction of weakly stratified surface waters was seen in the southern flank of the Gulf Stream as well. The solution of the omega equation suggests that this subduction was associated with a relatively shallow vertical circulation confined to the upper 200 m of the water column in the proximity of the front marking the southern edge of the warm core.

Freshwater export from the Labrador Current to the North Atlantic Current at the Tail of the Grand Banks of Newfoundland

Historical hydrographic data, spanning the period 1896–2006, are used to examine the annual mean and seasonal variations in the distribution of freshwater along and across the shelf/slope boundary along the Labrador and Newfoundland Shelves and the Grand Banks of Newfoundland. Particular attention is paid to the export of freshwater along the eastern Grand Banks, between Flemish Cap and the Tail of the Grand Banks, as this has long been identified as a preferential region for the loss of mass and freshwater from the boundary. The data are combined into isopycnally averaged long-term annual and monthly mean gridded property fields and the evolving distribution of fresh arctic-origin water is analyzed in fields of salinity anomaly, expressed as departures from the “central water” temperature–salinity relation of the Gulf Stream. The climatology confirms that cold/fresh northern-source waters are advected offshore within the retroflecting Labrador Current along the full length of the boundary between Flemish Cap and the Tail of the Grand Banks. In fact, it is estimated that most of the equatorward baroclinic transport at the boundary must retroflect back toward the north in order to explain the annual mean distribution of salinity in the climatology. While the retroflection of the Labrador Current appears seasonally robust, the freshwater distribution within the retroflection region varies in response to (1) the freshness of the water available for export which is set by the arrival and rapid flushing of the seasonal freshwater pulse at the boundary, (2) seasonal buoyancy forcing at the surface which alters the vertical stratification across the retroflection region, restricting certain isopycnal export pathways, and (3) the density structure along the eastern Grand Banks, which defines the progressive retroflection of the Labrador Current.

Evaluation of the S(T) assimilation method with the Argo dataset

AbstractRecent observations from the Argo dataset of temperature and salinity profiles are used to evaluate a series of 3‐year data assimilation experiments in a global ice–ocean general circulation model. The experiments are designed to evaluate a new data assimilation system whereby salinity is assimilated along isotherms, S(T). In addition, the role of a balancing salinity increment to maintain water mass properties is investigated. This balancing increment is found to effectively prevent spurious mixing in tropical regions induced by univariate temperature assimilation, allowing the correction of isotherm geometries without adversely influencing temperature–salinity relationships. In addition, the balancing increment is able to correct a fresh bias associated with a weak subtropical gyre in the North Atlantic using only temperature observations. The S(T) assimilation method is found to provide an important improvement over conventional depth level assimilation, with lower root‐mean‐squared forecast errors over the upper 500 m in the tropical Atlantic and Pacific Oceans. An additional set of experiments is performed whereby Argo data are withheld and used for independent evaluation. The most significant improvements from Argo assimilation are found in less well‐observed regions (Indian, South Atlantic and South Pacific Oceans). When Argo salinity data are assimilated in addition to temperature, improvements to modelled temperature fields are obtained due to corrections to model density gradients and the resulting circulation. It is found that observations from the Argo array provide an invaluable tool for both correcting modelled water mass properties through data assimilation and for evaluating the assimilation methods themselves. Copyright © 2009 Royal Meteorological Society

Warming and Freshening in the Abyssal Southeastern Indian Ocean*

Abstract Warming and freshening of abyssal waters in the eastern Indian Ocean between 1994/95 and 2007 are quantified using data from two closely sampled high-quality occupations of a hydrographic section extending from Antarctica northward to the equator. These changes are limited to abyssal waters in the Princess Elizabeth Trough and the Australian–Antarctic Basin, with little abyssal change evident north of the Southeast Indian Ridge. As in previous studies, significant cooling and freshening is observed in the bottom potential temperature–salinity relations in these two southern basins. In addition, analysis on pressure surfaces shows abyssal warming of about 0.05°C and freshening of about 0.01 Practical Salinity Scale 1978 (PSS-78) in the Princess Elizabeth Trough, and warming of 0.1°C with freshening of about 0.005 in the abyssal Australian–Antarctic Basin. These 12-yr differences are statistically significant from zero at 95% confidence intervals over the bottom few to several hundred decibars of the water column in both deep basins. Both warming and freshening reduce the density of seawater, contributing to the vertical expansion of the water column. The changes below 3000 dbar in these basins suggest local contributions approaching 1 and 4 cm of sea level rise, respectively. Transient tracer data from the 2007 occupation qualitatively suggest that the abyssal waters in the two southern basins exhibiting changes have significant components that have been exposed to the ocean surface within the last few decades, whereas north of the Southeast Indian Ridge, where changes are not found, the component of abyssal waters that have undergone such ventilation is much reduced.

In Situ Calibration of Moored CTDs Used for Monitoring Abyssal Water

Abstract To monitor changes in heat content and geostrophic volume transport of abyssal water accurately, 50 moored conductivity–temperature–depth (CTD) recorders used for density measurements were calibrated in situ by simultaneous observations with accurate shipboard CTDs. Comparisons of the data from the moored and shipboard CTDs showed pressure sensitivities of 0–3 mK at 6000 dbar for the temperature sensors of the moored CTDs. From the in situ calibrations, the uncertainties of the moored CTD data for the deep ocean (≥3000 dbar) were estimated to be 0.6 dbar, 0.6 mK, and 0.0026 for pressure, temperature, and salinity, respectively, relative to the shipboard CTD reference. Time drifts of the moored CTD data, estimated from the in situ calibrations before and after 17- or 14-month mooring deployments in the deep ocean, were considerably smaller than typical stabilities as specified by the manufacturer. However, time drifts of the pressure sensors tended to be negative and the result suggests that pressure data from most present Argo floats, which use the same type of pressure sensor, may have a systematic negative bias. Time series salinity data calculated from the in situ–calibrated CTDs were slightly biased (mean of +0.0014) with respect to the shipboard CTD salinity data, based on potential temperature–salinity relationships, possibly due to a disequilibrium of the moored CTD conductivity sensors during the in situ calibrations.

Living conditions, abundance and composition of the metazoan fauna in surface and sub-ice layers in pack ice of the western Weddell Sea during late spring

Abstract The surface and sub-ice layer habitats and their metazoan fauna were studied on a drifting pack-ice floe in the western Weddell Sea from 29 November 2004 to 1 January 2005 during the “Ice Station POLarstern” (ISPOL). Flooding of the floe occurred at some places, and the establishment of surface layers with a brownish colour due to growing algae was observed at several sampling sites. The average surface-layer temperature, brine salinity and brine volume were −1.4 °C, 25.3 and 54%, respectively. The temperature–salinity relationship in the surface layer was seldom at equilibrium conditions. Chlorophyll a (Chl a) concentrations in the brine varied between 1.0 and 53.5 μg L−1. Surface-layer thickness, salinity, Chl a concentration and copepod abundances were generally higher at the edge of the floe than in the inner part. The sympagic copepod species Drescheriella glacialis/racovitzai and Stephos longipes, with abundances ranging between 0 and 3830 ind. L−1 (median: 2 ind. L−1) and 0 and 1293 ind. L−1 (median: 4 ind. L−1), respectively, were the dominant members of the surface-layer meiofauna. Their populations consisted mainly of adults and early naupliar stages, which points to an active reproduction of these species within the surface layer. Other taxa found in the surface layer were undetermined turbellarians, the gastropod Tergipes antarcticus, and, for the first time, the ctenophore Callianira antarctica, and the amphipods Eusirus antarcticus and Eusirus tridentatus. During the course of our study, slight melting at the ice underside took place, releasing sympagic organisms to the water column. Chl a concentrations in the sub-ice water layer were very low (0.1–0.5 μg L−1), except for 25 December when the Chl a concentration at 0 m depth increased to 2.3 μg L−1. The most dominant sympagic copepod species found in the sub-ice layer was Ectinosoma sp., with abundances ranging between 1 and 599 ind. m−3 (median: 25 ind. m−3). Other sympagic copepod species occurring regularly in this habitat were D. glacialis/racovitzai, Diarthrodes cf. lilacinus, Idomene antarctica and S. longipes. All of these sympagic species were generally found in higher abundances at 0 m depth underneath the ice than at 5 m depth, in contrast to pelagic copepod species that occurred more frequently at 5 m depth. Niche separation and probable life-cycle strategies of dominant sympagic metazoans are discussed.

The evolution of volcano-hosted geothermal systems based on deep wells from Karaha-Telaga Bodas, Indonesia

Temperature and pressure surveys, fluid samples, and petrologic analyses of rock samples from deep drill holes at the Karaha - Telaga Bodas geothermal field on the volcanic ridge extending northward from Galunggung Volcano, West Java, have provided a unique opportunity to characterize the evolution of an active volcano-hosted geothermal system. Wells up to 3 km in depth have encountered temperatures as high as 353°C and a weakly altered granodiorite that intruded to within 2 to 3 km of the surface. The intrusion is shallowest beneath the southern end of the field where an acid lake overlies a nearly vertical low resistivity structure (<10 ohm-m) defined by magnetotelluric measurements. This structure is interpreted to represent a vapor-dominated chimney that provides a pathway to the surface for magmatic gases. Four distinct hydrothermal mineral assemblages document the evolution of the geothermal system and the transition from liquid- to vapor-dominated conditions. The earliest assemblage represents the initial liquid-dominated system generated during emplacement of the granodiorite between 5910 ± 76 and 4200 ± 150 y BP. Tourmaline, biotite, actinolite, epidote and clay minerals were deposited contemporaneously at progressively greater distances from the intrusive contact (assemblage 1). At 4200 ± 150 y BP, flank collapse and the formation of the volcano9s crater, Kawah Galunggung, resulted in catastrophic decompression and boiling of the hydrothermal fluids. This event initiated development of the modern vapor-dominated regime. Chalcedony and then quartz were deposited as the early low salinity liquids boiled (assemblage 2). Both vapor- and liquid-rich fluid inclusions were trapped in the quartz crystals. Liquid-rich fluid inclusions from the southern part of the field record salinities ranging from 0 to 26 weight percent NaCl-CaCl₂ equivalent and locally contain fluorite daughter crystals. We suggest, based on temperature-salinity relationships and evidence of boiling, that these fluids were progressively concentrated as steam was lost from the system. However, mixing with fluids derived from the underlying intrusion or generated during the formation of acid SO₄ water on the vapor-dominated chimney margins could have contributed to the observed salinities. As pressures declined, CO₂–and SO₄-rich steam-heated water drained downward, depositing anhydrite and calcite (assemblage 3) in the fractures, limiting further recharge. Fluid inclusions with salinities up to 31 weight percent NaCl equivalent were trapped in these minerals as the descending water vaporized. The final assemblage is represented by precipitates of NaCl, KCl and FeCl_x deposited on rock surfaces in portions of the vapor-dominated zone that boiled dry. Vapor-dominated conditions extend over a distance of at least 10 km and to depths below sea level. Deep wells drilled into the underlying liquid-dominated reservoir in the northern and central part of the volcanic ridge produce low salinity fluids representing recent recharge of meteoric and steam-heated water. The evolution of volcanic-hosted vapor-dominated geothermal systems can be described by a five stage model. Stage 1 involves the formation of an over-pressured liquid-dominated geothermal system soon after magmatic intrusion. In Stages 2 and 3, pressures progressively decrease, and a curtain of steam-heated water surrounding a magmatic vapor-dominated chimney at 350°C and 14 ± 2 MPa develops. The relatively low pressure near the base of the chimney causes liquid inflow adjacent to the intrusion and the development of a secondary marginal vapor-dominated zone. In Stage 4, the magmatic vapor discharge from the intrusion becomes small, vapor pressure declines, and the secondary vapor-dominated zone expands above the intrusion. In Stage 5, the vapor-dominated zone floods because heat from the intrusion is insufficient to boil all liquid inflow. A more common, liquid-dominated volcanic-hosted system then develops.