The Argo array provides nearly 4000 temperature and salinity profiles of the top 2000 m of the ocean every 10 days. Still, Argo floats will never be able to measure the ocean at all times, everywhere. Optimized Argo float distributions should match the spatial and temporal variability of the many societally important ocean features that they observe. Determining these distributions is challenging because float advection is difficult to predict. Using no external models, transition matrices based on existing Argo trajectories provide statistical inferences about Argo float motion. We use the 24 years of Argo locations to construct an optimal transition matrix that minimizes estimation bias and uncertainty. The optimal array is determined to have a 2° × 2° spatial resolution with a 90-day time step. We then use the transition matrix to predict the probability of future float locations of the core Argo array, the Global Biogeochemical Array, and the Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) array. A comparison of transition matrices derived from floats using Argos system and Iridium communication methods shows the impact of surface displacements, which is most apparent near the equator. Additionally, we demonstrate the utility of transition matrices for validating models by comparing the matrix derived from Argo floats with that derived from a particle release experiment in the Southern Ocean State Estimate (SOSE).
- Award ID(s):
- 1936222
- NSF-PAR ID:
- 10436780
- Date Published:
- Journal Name:
- Biogeosciences
- Volume:
- 20
- Issue:
- 8
- ISSN:
- 1726-4189
- Page Range / eLocation ID:
- 1671 to 1690
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract This paper describes the new Regional Arctic Ocean/sea ice Reanalysis (RARE) with a domain that spans a subpolar/polar cap poleward of 45°N. Sequential data assimilation constrains temperature and salinity using World Ocean Database profiles as well as in situ and satellite SST, and PIOMAS sea ice thickness estimates. The 41-yr (1980–2020) RARE1.15.2 reanalysis with resolution varying between 2 and 5 km horizontally and 1–10 m vertically in the upper 100 m is examined. To explore the impact of resolution RARE1.15.2 is compared to a coarser-resolution SODA3.15.2, which uses the same modeling and data assimilation system. Improving resolution in the reanalysis system improves agreement with observations. It produces stronger more compact currents, enhances eddy kinetic energy, and strengthens along-isopycnal heat and salt transports, but reduces vertical exchanges and thus strengthens upper ocean haline stratification. RARE1.15.2 and SODA3.15.2 are also compared to the Hadley Center EN4.2.2 statistical objective analysis. In regions of reasonable data coverage such as the Nordic seas the three products produce similar time-mean distributions of temperature and salinity. But in regions of poor coverage and in regions where the coverage changes in time EN4.2.2 suffers more from those inhomogeneities. Finally, the impact on the Arctic of interannual temperature fluctuations in the subpolar gyres on the Arctic Ocean is compared. The influence of the subpolar North Pacific is limited to a region surrounding Bering Strait. The influence of the subpolar North Atlantic, in contrast, spreads throughout the Nordic seas and Barents Sea in all three products within two years.
Significance Statement The Arctic Ocean/sea ice system plays crucial roles in climate variability and change by controlling the northern end of the oceanic overturning circulation, the equator to pole air pressure gradient, and Earth’s energy balance. Yet the historical ocean observation set is sparse and inhomogeneous, while ocean dynamics has challengingly fine horizontal and vertical scales. This paper introduces a new Regional Arctic Ocean/sea ice Reanalysis (RARE) whose goal is to use the combined constraints of mesoscale ocean dynamics, historical observations, surface meteorology, and continental runoff in a data assimilation framework to reconstruct historical variability. RARE is used to produce a 41-yr ocean/sea ice reanalysis 1980–2020 whose results are described here.
-
more » « less
Abstract Argo‐type profiling floats do not receive satellite positioning while under sea ice. Common practice is to approximate unknown positions by linearly interpolating latitude‐longitude between known positions before and after ice cover, although it has been suggested that some improvement may be obtained by interpolating along contours of planetary‐geostrophic potential vorticity. Profiles with linearly interpolated positions represent 16% of the Southern Ocean Argo data set; consequences arising from this approximation have not been quantified. Using three distinct data sets from the Weddell Gyre—10‐day satellite‐tracked Argo floats, daily‐tracked RAFOS‐enabled floats, and a particle release simulation in the Southern Ocean State Estimate—we perform a data withholding experiment to assess position uncertainty in latitude‐longitude and potential vorticity coordinates as a function of time since last fix. A spatial correlation analysis using the float data provides temperature and salinity uncertainty estimates as a function of distance error. Combining the spatial correlation scales and the position uncertainty, we estimate uncertainty in temperature and salinity as a function of duration of position loss. Maximum position uncertainty for interpolation during 8 months without position data is 116 ± 148 km for latitude‐longitude and 92 ± 121 km for potential vorticity coordinates. The estimated maximum uncertainty in local temperature and salinity over the entire 2,000‐m profiles during 8 months without position data is 0.66 ∘C and 0.15 psu in the upper 300 m and 0.16 ∘C and 0.01 psu below 300 m.
-
more » « less
Abstract From seasonal cruises in the NE Pacific Ocean during 2017, we (1) determined dissolved organic carbon concentrations; (2) calculated net community production (NCP) from nitrate drawdown; and (3) established relationships between NCP and seasonal dissolved organic carbon (DOC) accumulation in the upper 75 m. The fraction of NCP that accumulated as DOC, hereafter referred to as the net dissolved production ratio, was calculated for several stations during spring and summer. The net dissolved production ratio was about 0.26 at the oceanic station Ocean Station Papa during different seasons and years. Using nitrate concentration profiles obtained from Bio‐Argo floats during 2009–2018 operating near Ocean Station Papa, we calculated NCP at high temporal resolution and then applied the 0.26 constant in order to (4) estimate DOC variability for the 9‐year period. We found strong seasonality near Ocean Station Papa, with NCP maxima during summers ranging from 0.3 to 2.9 mol C/m2and surface DOC concentrations estimated from 56 μmol/kg in winters to 73 μmol/kg in summers. There was a 10‐fold interannual variability in the seasonally accumulated inventory of DOC, ranging from 0.078 to 0.75 mol C/m2. This study reinforces the value of deploying floats equipped with chemical sensors in order to better understand marine biogeochemical cycles, especially when high resolution data cannot be obtained otherwise. Given that ~26% of NCP accumulates as DOC in the central Gulf of Alaska, the remaining balance of ~74% is available for export as sinking biogenic particles.
-
more » « less
Abstract The ocean mixed layer plays an important role in the coupling between the upper ocean and atmosphere across a wide range of time scales. Estimation of the variability of the ocean mixed layer is therefore important for atmosphere‐ocean prediction and analysis. The increasing coverage of in situ Argo profile data allows for an increasingly accurate analysis of the mixed layer depth (MLD) variability associated with deviations from the seasonal climatology. However, sampling rates are not sufficient to fully resolve subseasonal (
day) MLD variability. Yet, many multivariate observations‐based analyses include implicit modeled subseasonal MLD variability. One analysis method is optimal interpolation of in situ data, but the interior analysis can be improved by leveraging surface data with regression or variational approaches. Here, we demonstrate how machine learning methods and satellite sea surface temperature, salinity, and height facilitate MLD estimation in a pilot study of two regions: the mid‐latitude southern Indian and the eastern equatorial Pacific Oceans. We construct multiple machine learning architectures to produce weekly 1/2° gridded MLD anomaly fields (relative to a monthly climatology) with uncertainty estimates. We test multiple traditional and probabilistic machine learning techniques to compare both accuracy and probabilistic calibration. We validate our methodology by applying it to ocean model simulations. We find that incorporating sea surface data through a machine learning model improves the performance of spatiotemporal MLD variability estimation compared to optimal interpolation of Argo observations alone. These preliminary results are a promising first step for the application of machine learning to MLD prediction.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.5194/bg-20-1671-2023
