Thermocline Model for Estimating Argo Sea Surface Temperature

Argo has become an important constituent of the global ocean observation system. However, due to the lack of sea surface measurements from most Argo profiles, the application of Argo data is still limited. In this study, a thermocline model was constructed based on three key thermocline parameters, i.e, thermocline upper depth, the thermocline bottom depth, and thermocline temperature gradient. Following the model, we estimated the sea surface temperature of Argo profiles by providing the relationship between sea surface and subsurface temperature. We tested the effectiveness of our proposed model using statistical analysis and by comparing the sea surface temperature with the results obtained from traditional methods and in situ observations in the Pacific Ocean. The root mean square errors of results obtained from thermocline model were found to be significantly reduced compared to the extrapolation results and satellite retrieved temperature results. The correlation coefficient between the estimation result and in situ observation was 0.967. Argo surface temperature, estimated by the thermocline model, has been theoretically proved to be reliable. Thus, our model generates theoretically feasible data present the mesoscale phenomenon in more detail. Overall, this study compensates for the lack surface observation of Argo, and provides a new tool to establish complete Argo data sets.

Practical Dynamical-Statistical Reconstruction of Ocean’s Interior from Satellite Observations

The algorithms based on Surface Quasi-Geostrophic (SQG) dynamics have been developed and validated by many researchers through model products, however it is still doubtful whether these SQG-based algorithms are worth using in terms of observed data. This paper analyzes the factors impeding the practical application of SQG and makes amends by a simple “first-guess (FG) framework”. The proposed framework includes the correction of satellite salinity and the estimation of the FG background, making the SQG-based algorithms applicable in realistic circumstances. The dynamical-statistical method SQG-mEOF-R is thereafter applied to satellite data for the first time. The results are compared with two dynamical algorithms, SQG and isQG, and three empirical algorithms, multivariate linear regression (MLR), random forest (RF), and mEOF-R. The validation against Argo profiles showed that the SQG-mEOF-R presents a robust performance in mesoscale reconstruction and outperforms the other five algorithms in the upper layers. It is promising that the SQG-mEOF-R and the FG framework are applicable to operational reconstruction.

Can assimilation of satellite observations improve subsurface biological properties in a numerical model? A case study for the Gulf of Mexico

Given current threats to ocean ecosystem health, there is a growing demand for accurate biogeochemical hindcasts, nowcasts, and predictions. Provision of such products requires data assimilation, i.e., a comprehensive strategy for incorporating observations into biogeochemical models, but current data streams of biogeochemical observations are generally considered insufficient for the operational provision of such products. This study investigates to what degree the assimilation of satellite observations in combination with a priori model calibration by sparse BGC-Argo profiles can improve subsurface biogeochemical properties. The multivariate deterministic ensemble Kalman filter (DEnKF) has been implemented to assimilate physical and biological observations into a three-dimensional coupled physical–biogeochemical model, the biogeochemical component of which has been calibrated by BGC-Argo float data for the Gulf of Mexico. Specifically, observations of sea surface height, sea surface temperature, and surface chlorophyll were assimilated, and profiles of both physical and biological variables were updated based on the surface information. We assessed whether this leads to improved subsurface distributions, especially of biological properties, using observations from five BGC-Argo floats that were not assimilated. An alternative light parameterization that was tuned a priori using BGC-Argo observations was also applied to test the sensitivity of data assimilation impact on subsurface biological properties. Results show that assimilation of the satellite data improves model representation of major circulation features, which translate into improved three-dimensional distributions of temperature and salinity. The multivariate assimilation also improves the agreement of subsurface nitrate through its tight correlation with temperature, but the improvements in subsurface chlorophyll were modest initially due to suboptimal choices of the model's optical module. Repeating the assimilation run by using the alternative light parameterization greatly improved the subsurface distribution of chlorophyll. Therefore, even sparse BGC-Argo observations can provide substantial benefits for biogeochemical prediction by enabling a priori model tuning. Given that, so far, the abundance of BGC-Argo profiles in the Gulf of Mexico and elsewhere has been insufficient for sequential assimilation, updating 3D biological properties in a model that has been well calibrated is an intermediate step toward full assimilation of the new data types.

Combining Argo and Satellite Data Using Model-Derived Covariances: Blue Maps

Blue Maps aims to exploit the versatility of an ensemble data assimilation system to deliver gridded estimates of ocean temperature, salinity, and sea-level with the accuracy of an observation-based product. Weekly maps of ocean properties are produced on a 1/10°, near-global grid by combining Argo profiles and satellite observations using ensemble optimal interpolation (EnOI). EnOI is traditionally applied to ocean models for ocean forecasting or reanalysis, and usually uses an ensemble comprised of anomalies for only one spatiotemporal scale (e.g., mesoscale). Here, we implement EnOI using an ensemble that includes anomalies for multiple space- and time-scales: mesoscale, intraseasonal, seasonal, and interannual. The system produces high-quality analyses that produce mis-fits to observations that compare well to other observation-based products and ocean reanalyses. The accuracy of Blue Maps analyses is assessed by comparing background fields and analyses to observations, before and after each analysis is calculated. Blue Maps produces analyses of sea-level with accuracy of about 4 cm; and analyses of upper-ocean (deep) temperature and salinity with accuracy of about 0.45 (0.15) degrees and 0.1 (0.015) practical salinity units, respectively. We show that the system benefits from a diversity of ensemble members with multiple scales, with different types of ensemble members weighted accordingly in different dynamical regions.

Global Oceanic Eddy Identification: A Deep Learning Method From Argo Profiles and Altimetry Data

The inadequate spatial resolution of altimeter results in low identification efficiency of oceanic eddies, especially for small-scale eddies. It is well known that eddies can not only induce sea surface signal but more importantly have typical vertical structure characteristics. However, although the vertical structure characteristics are usually used for statistical analysis, they are seldom considered in the process of eddy recognition. This study is devoted to identifying eddies from the perspective of their vertical signal derived from the 18-year Argo data. Due to the irregular and noisy profile pattern, the direct identification of eddy core from Argo profile is deemed to be a challenge. With the popularity of artificial intelligence, a new hybrid method that combines the advantages of convolutional neural network (CNN) with extreme gradient boosting (XGBoost) is proposed to extract the representative vertical feature and identify eddy from a profile. First, CNN is employed as a feature extractor to automatically obtain vertical features from the input profile at the bottom of the network. Second, the obtained high-dimensional feature vectors are inputted into the XGBoost model, combined with other profile features for classifying profiles that are outside altimeter-identified eddies (Alt eddy). Finally, extensive experiments are implemented to demonstrate the efficiency of the proposed method. The results show that the classification accuracy of CNN-XGBoost model can reach 98%, and about 36% eddies are recaptured. These eddies, dubbed CNN-XGB eddies, are benchmarked against Alt eddies for the vertical structure and geographical distribution, demonstrating a similar or even stronger vertical signal and a prominent eddy belt in the tropical ocean. Within the proposed theory framework, there are various potentials to obtain a better outlook for eddy identification and in situ float observations.

Can assimilation of satellite observations improve subsurface biological properties in a numerical model? A case study for the Gulf of Mexico

Given current threats to ocean ecosystem health, there is a growing demand for accurate biogeochemical hindcasts, nowcasts, and predictions. Provision of such products requires data assimilation, i.e., a comprehensive strategy for incorporating observations into biogeochemical models, but current data streams of biogeochemical observations are generally considered insufficient for the operational provision of such products. This study investigates to what degree the satellite observations in combination with sparse BGC Argo profiles can improve subsurface biogeochemical properties. The multivariate Deterministic Ensemble Kalman Filter (DEnKF) has been implemented to assimilate physical and biological observations into a biogeochemical model of the Gulf of Mexico. First, the biogeochemical model component was tuned using BGC-Argo observations. Then, observations of sea surface height, sea surface temperature, and surface chlorophyll were assimilated, and profiles of both physical and biological variables were updated based on the surface information. We assessed whether this leads to improved subsurface distributions, especially of biological properties, using observations from five BGC-Argo floats that were not assimilated, but used in the a priori tuning. Results show that assimilation of the satellite data improves model representation of major circulation features, which translate into improved three-dimensional distributions of temperature and salinity. The multivariate assimilation also improves the agreement of subsurface nitrate through its tight correlation with temperature, but the improvements in subsurface chlorophyll were modest initially due to suboptimal choices of the model’s optical module. Repeating the assimilation run after adjusting light attenuation parameterization through further a priori tuning greatly improved the subsurface distribution of chlorophyll. Therefore, even sparse BGC-Argo observations can provide substantial benefits to biogeochemical prediction by enabling a priori model tuning. Given that, so far, the abundance of BGC-Argo profiles in the Gulf of Mexico and elsewhere is insufficient for sequential assimilation, updating 3D biological properties in a model that has been well calibrated is an intermediate step toward full assimilation of the new data types.

Variability of the hydrological characteristics of the Eastern Mediterranean over the last decade

<p><span>A hydrodynamic simulation is carried out over the entire Mediterranean basin at a resolution of 3 to 4 km and a duration of about 10 years (2011-2020). The results are systematically evaluated using Argo profiles focusing on the spatial distribution of water mass properties along their path, the main mesoscale structures, the mean vertical temperature and salinity profiles by sub-basins as well as their "pseudo temporal evolution" biased by the variability of the spatial and temporal distribution of Argo observations.</span></p><p><span>The simulation has generally very low mean biases (of the order of 0.01 for salinity) and correlations on the monthly time series reconstructed from the observations, of the order of 0.9 at the scale of the eastern basin, both in surface waters and at 200 m in intermediate waters. </span></p><p><span>The evolution of salinity over the decade is then analyzed from the simulation. Particular attention is paid to the main basins of water mass formation, the Adriatic, the Levantine basin and the South Aegean Sea. The factors driving this evolution are analyzed in each of these basins. The propagation of the changes from these formation areas to the entire eastern basin is then examined, with a particular focus on the intermediate waters. </span></p>

Sea Surface Temperature and Ocean Heat Content during Tropical Cyclones Pam (2015) and Winston (2016) in the Southwest Pacific Region

The sea surface temperature (SST) and upper ocean heat content (OHC) have been explored along the track of two tropical cyclones (TCs), TC Pam (2015) and TC Winston (2016). These TCs severely affected the islands of Vanuatu and Fiji, in the South Pacific Region (8°–30°S, 140°E– 170°W). The SST decreased by as much as 5.4°C along the tracks of the TCs with most cooling occurring to the left of the TCs tracks relative to TCs motion. SST cooling of 1-5°C has generally been observed during both the forced and relaxation stages of TC passage. The Argo profiles near the TC revealed observable mixed layer deepening. Subsurface warming was also observed post-TC passage from the temperature profile of one of the floats after the passage of both TCs. The OHC and heat fluxes are seen to play an important part in TC intensification as both these TCs intensified after passing over the regions of high OHC and enhanced heat fluxes. Apart from the traditionally used OHC obtained up to the depth of the 26°C isotherm (QH), the OHC was also determined up to the depth of the 20°C isotherm (QH,20). The QH and QH,20 values decreased in the majority of cases post TC passage while QH,20 increased in one instance post-TC passage for both the TCs. QH,20 has also been used to identify heat energy changes at deeper levels and correlated well with the traditionally used OHC during the weaker stages of the TCs.