Browse
Abstract
The equatorial Pacific zonal circulation is composed of westward surface currents, the eastward equatorial undercurrent (EUC) along the thermocline, and upwelling in the eastern cold tongue. Part of this upwelling arises from water flowing along isotherms sloping up to the east, but it also includes water mass transformation and consequent diabatic (cross-isothermal) flow (w ci) that is a key element of surface-to-thermocline communication. In this study we investigate the mean seasonal cycle and subseasonal variability of cross-isothermal flow in the cold tongue using heat budget output from a high-resolution forced ocean model. Diabatic upwelling is present throughout the year with surface-layer solar-penetration-driven diabatic upwelling strongest in boreal spring and vertical mixing in the thermocline dominating during the rest of the year. The former constitutes warming of the surface layer by solar radiation rather than exchange of thermal energy between water parcels. The mixing-driven regime allows heat to be transferred to the core of the EUC by warming parcels at depth. On subseasonal time scales the passage of tropical instability waves (TIWs) enhances diabatic upwelling on and north of the equator. On the equator the TIWs enhance vertical shear and induce vertical-mixing-driven diabatic upwelling, while off the equator TIWs enhance the sub-5-daily eddy heat flux which enhances diabatic upwelling. Comparing the magnitudes of TIW, seasonal, and interannual w ci variability, we conclude that each time scale is associated with sizeable variance. Variability across all of these time scales needs to be taken into account when modeling or diagnosing the effects of mixing on equatorial upwelling.
Abstract
The equatorial Pacific zonal circulation is composed of westward surface currents, the eastward equatorial undercurrent (EUC) along the thermocline, and upwelling in the eastern cold tongue. Part of this upwelling arises from water flowing along isotherms sloping up to the east, but it also includes water mass transformation and consequent diabatic (cross-isothermal) flow (w ci) that is a key element of surface-to-thermocline communication. In this study we investigate the mean seasonal cycle and subseasonal variability of cross-isothermal flow in the cold tongue using heat budget output from a high-resolution forced ocean model. Diabatic upwelling is present throughout the year with surface-layer solar-penetration-driven diabatic upwelling strongest in boreal spring and vertical mixing in the thermocline dominating during the rest of the year. The former constitutes warming of the surface layer by solar radiation rather than exchange of thermal energy between water parcels. The mixing-driven regime allows heat to be transferred to the core of the EUC by warming parcels at depth. On subseasonal time scales the passage of tropical instability waves (TIWs) enhances diabatic upwelling on and north of the equator. On the equator the TIWs enhance vertical shear and induce vertical-mixing-driven diabatic upwelling, while off the equator TIWs enhance the sub-5-daily eddy heat flux which enhances diabatic upwelling. Comparing the magnitudes of TIW, seasonal, and interannual w ci variability, we conclude that each time scale is associated with sizeable variance. Variability across all of these time scales needs to be taken into account when modeling or diagnosing the effects of mixing on equatorial upwelling.
Abstract
The sea surface expressions of Mediterranean Water eddies, known as “meddies,” are observed in satellite data, and their main characteristics are measured. Satellite altimeter observations of surface expressions are detected over the meddies observed in situ using the MEDTRANS meddy dataset (1950–2013). In this study 209 observed meddy cores in the North Atlantic Ocean, selected over the period of the 22 years of sea surface height measurements with satellite altimetry (1993–2013), were analyzed. Results show relatively good agreement between the theoretical estimates of the meddy surface signals as reported by Bashmachnikov and Carton and the measured surface expressions. It was found that, on average, the theoretical results underestimate the measured sea surface elevations of the meddy surface expressions by a factor of 2. Although the variability of the measured expressions is reasonably well described by the combination of meddy core and the ocean background parameters of the theoretical expression, we cannot define a single individual parameter of the meddy core, which chiefly shapes the magnitude of the meddy surface signal. Interestingly, the overall distribution of characteristics of meddy surface expressions in the Atlantic shows that the sea level anomalies formed by meddies intensify westward, growing both in magnitude and radius. This opposes the expected theoretical decrease of meddy surface signals due to a known progressive decay of the meddy cores with distance from their generation region at the Iberian continental slope. This observed tendency is attributed to meddy interaction with the upper-ocean currents and other eddies (in particular in the region of the North Atlantic Current and Azores Current) that are not considered by the theory.
Abstract
The sea surface expressions of Mediterranean Water eddies, known as “meddies,” are observed in satellite data, and their main characteristics are measured. Satellite altimeter observations of surface expressions are detected over the meddies observed in situ using the MEDTRANS meddy dataset (1950–2013). In this study 209 observed meddy cores in the North Atlantic Ocean, selected over the period of the 22 years of sea surface height measurements with satellite altimetry (1993–2013), were analyzed. Results show relatively good agreement between the theoretical estimates of the meddy surface signals as reported by Bashmachnikov and Carton and the measured surface expressions. It was found that, on average, the theoretical results underestimate the measured sea surface elevations of the meddy surface expressions by a factor of 2. Although the variability of the measured expressions is reasonably well described by the combination of meddy core and the ocean background parameters of the theoretical expression, we cannot define a single individual parameter of the meddy core, which chiefly shapes the magnitude of the meddy surface signal. Interestingly, the overall distribution of characteristics of meddy surface expressions in the Atlantic shows that the sea level anomalies formed by meddies intensify westward, growing both in magnitude and radius. This opposes the expected theoretical decrease of meddy surface signals due to a known progressive decay of the meddy cores with distance from their generation region at the Iberian continental slope. This observed tendency is attributed to meddy interaction with the upper-ocean currents and other eddies (in particular in the region of the North Atlantic Current and Azores Current) that are not considered by the theory.
Abstract
Changes in dynamic manometric sea level ζm
represent mass-related sea level changes associated with ocean circulation and climate. We use twin model experiments to quantify magnitudes and spatiotemporal scales of ζm
variability caused by barometric pressure pa
loading at long periods (
Abstract
Changes in dynamic manometric sea level ζm
represent mass-related sea level changes associated with ocean circulation and climate. We use twin model experiments to quantify magnitudes and spatiotemporal scales of ζm
variability caused by barometric pressure pa
loading at long periods (
Abstract
Oceanic flows are energetically dominated by low vertical modes. However, disturbances in the form of atmospheric storms, eddy interactions with various forms of boundaries, or spontaneous emission by coherent structures can generate weak high-baroclinic modes. The feedback of the low-energy high-baroclinic modes on large-scale energetically dominant low modes may be weak or strong depending on the flow Rossby number. In this paper we study this interaction using an idealized setup by constraining the flow dynamics to a high-energy barotropic mode and a single low-energy high-baroclinic mode. Our investigation points out that at low Rossby numbers the barotropic flow organizes into large-scale coherent vortices via an inverse energy flux while the baroclinic flow accumulates predominantly in anticyclonic barotropic vortices. In contrast, with increasing Rossby number, the baroclinic flow catalyzes a forward flux of barotropic energy. The barotropic coherent vortices decrease in size and number, with a strong preference for cyclonic coherent vortices at higher Rossby numbers. On partitioning the flow domain into strain-dominant and vorticity-dominant regions based on the barotropic flow, we find that at higher Rossby numbers baroclinic flow accumulates in strain-dominant regions, away from vortex cores. Additionally, a major fraction of the forward energy flux of the flow takes place in strain-dominant regions. Overall, one of the key outcomes of this study is the finding that even a low-energy high-baroclinic flow can deplete and dissipate large-scale coherent structures at O(1) Rossby numbers.
Abstract
Oceanic flows are energetically dominated by low vertical modes. However, disturbances in the form of atmospheric storms, eddy interactions with various forms of boundaries, or spontaneous emission by coherent structures can generate weak high-baroclinic modes. The feedback of the low-energy high-baroclinic modes on large-scale energetically dominant low modes may be weak or strong depending on the flow Rossby number. In this paper we study this interaction using an idealized setup by constraining the flow dynamics to a high-energy barotropic mode and a single low-energy high-baroclinic mode. Our investigation points out that at low Rossby numbers the barotropic flow organizes into large-scale coherent vortices via an inverse energy flux while the baroclinic flow accumulates predominantly in anticyclonic barotropic vortices. In contrast, with increasing Rossby number, the baroclinic flow catalyzes a forward flux of barotropic energy. The barotropic coherent vortices decrease in size and number, with a strong preference for cyclonic coherent vortices at higher Rossby numbers. On partitioning the flow domain into strain-dominant and vorticity-dominant regions based on the barotropic flow, we find that at higher Rossby numbers baroclinic flow accumulates in strain-dominant regions, away from vortex cores. Additionally, a major fraction of the forward energy flux of the flow takes place in strain-dominant regions. Overall, one of the key outcomes of this study is the finding that even a low-energy high-baroclinic flow can deplete and dissipate large-scale coherent structures at O(1) Rossby numbers.
Abstract
The Antarctic Ice Sheet is losing mass as a result of increased ocean-driven melting of its fringing ice shelves. Efforts to represent the effects of basal melting in sea level projections are undermined by poor understanding of the turbulent ice shelf–ocean boundary layer (ISOBL), a meters-thick layer of ocean that regulates heat and salt transfer between the ocean and ice. To address this shortcoming, we perform large-eddy simulations of the ISOBL formed by a steady, geostrophic flow beneath horizontal ice. We investigate melting and ISOBL structure and properties over a range of free-stream velocities and ocean temperatures. We find that the melting response to changes in thermal and current forcing is highly nonlinear due to the effects of meltwater on ISOBL turbulence. Three distinct ISOBL regimes emerge depending on the relative strength of current shear and buoyancy forcing: “well-mixed,” “stratified,” or “diffusive-convective.” We present expressions for mixing-layer depth for each regime and show that the transitions between regimes can be predicted with simple nondimensional parameters. We use these results to develop a novel regime diagram for the ISOBL which provides insight into the varied melting responses expected around Antarctica and highlights the need to include stratified and diffusive-convective dynamics in future basal melting parameterizations. We emphasize that melting in the diffusive-convective regime is time dependent and is therefore inherently difficult to parameterize.
Significance Statement
The purpose of this study is to investigate the processes that control ocean-driven melting of Antarctic ice shelves (100–1000-m-thick floating extensions of the Antarctic ice sheet). Currently, these processes are poorly understood due to the difficulty of accessing the ocean beneath ice shelves. Using an ocean model, we determine the melting response to different ocean conditions, including feedbacks whereby cold, fresh meltwater can enhance or suppress turbulent eddies beneath the ice, depending on the ocean state. Our results point the way to improvements in the representation of ocean-driven melting in ocean/climate models, which will allow more accurate predictions of future climate and sea level.
Abstract
The Antarctic Ice Sheet is losing mass as a result of increased ocean-driven melting of its fringing ice shelves. Efforts to represent the effects of basal melting in sea level projections are undermined by poor understanding of the turbulent ice shelf–ocean boundary layer (ISOBL), a meters-thick layer of ocean that regulates heat and salt transfer between the ocean and ice. To address this shortcoming, we perform large-eddy simulations of the ISOBL formed by a steady, geostrophic flow beneath horizontal ice. We investigate melting and ISOBL structure and properties over a range of free-stream velocities and ocean temperatures. We find that the melting response to changes in thermal and current forcing is highly nonlinear due to the effects of meltwater on ISOBL turbulence. Three distinct ISOBL regimes emerge depending on the relative strength of current shear and buoyancy forcing: “well-mixed,” “stratified,” or “diffusive-convective.” We present expressions for mixing-layer depth for each regime and show that the transitions between regimes can be predicted with simple nondimensional parameters. We use these results to develop a novel regime diagram for the ISOBL which provides insight into the varied melting responses expected around Antarctica and highlights the need to include stratified and diffusive-convective dynamics in future basal melting parameterizations. We emphasize that melting in the diffusive-convective regime is time dependent and is therefore inherently difficult to parameterize.
Significance Statement
The purpose of this study is to investigate the processes that control ocean-driven melting of Antarctic ice shelves (100–1000-m-thick floating extensions of the Antarctic ice sheet). Currently, these processes are poorly understood due to the difficulty of accessing the ocean beneath ice shelves. Using an ocean model, we determine the melting response to different ocean conditions, including feedbacks whereby cold, fresh meltwater can enhance or suppress turbulent eddies beneath the ice, depending on the ocean state. Our results point the way to improvements in the representation of ocean-driven melting in ocean/climate models, which will allow more accurate predictions of future climate and sea level.
Abstract
The impact of El Niño–Southern Oscillation (ENSO) on the southern African climate is well documented and provides skill in the seasonal forecast of rainfall, but less is known about the impact of ENSO on the Benguela Current west of southern Africa. There is a significant weak correlation between ENSO and the Benguela Current upwelling sea surface temperature (SST) in austral summer. Correlation is positive for southern Benguela and negative for northern Benguela. A significant correlation exists with up to 8 months lag when ENSO leads. The impact of ENSO is due to weaker-than-normal upwelling favorable southeasterly winds during El Niño in southern Benguela, leading to warmer-than-normal coastal SST. In contrast, during La Niña, stronger-than-normal southeasterly winds lead to cooler-than-normal SST. The opposite effect applies to northern Benguela. The coastal wind change is part of an ENSO large-scale basinwide perturbation in the tropical and South Atlantic. However, non-ENSO-related SST variation in the Benguela upwelling can be as important as ENSO-related SST perturbation, and some ENSO events do not lead to the expected changes. Changes in the Benguela upwelling are linked to changes in the intensity of the trade winds associated with a change of the South Atlantic anticyclone intensity and position. In southern Benguela, changes are also associated with variations in midlatitude low pressure systems and associated upwelling unfavorable westerly winds. La Niñas favor the development of Benguela Niños in Angola and Namibia. This study shows the potential for SST seasonal predictability in the Benguela upwelling due to the leading lag correlation between ENSO and the Benguela upwelling SST.
Abstract
The impact of El Niño–Southern Oscillation (ENSO) on the southern African climate is well documented and provides skill in the seasonal forecast of rainfall, but less is known about the impact of ENSO on the Benguela Current west of southern Africa. There is a significant weak correlation between ENSO and the Benguela Current upwelling sea surface temperature (SST) in austral summer. Correlation is positive for southern Benguela and negative for northern Benguela. A significant correlation exists with up to 8 months lag when ENSO leads. The impact of ENSO is due to weaker-than-normal upwelling favorable southeasterly winds during El Niño in southern Benguela, leading to warmer-than-normal coastal SST. In contrast, during La Niña, stronger-than-normal southeasterly winds lead to cooler-than-normal SST. The opposite effect applies to northern Benguela. The coastal wind change is part of an ENSO large-scale basinwide perturbation in the tropical and South Atlantic. However, non-ENSO-related SST variation in the Benguela upwelling can be as important as ENSO-related SST perturbation, and some ENSO events do not lead to the expected changes. Changes in the Benguela upwelling are linked to changes in the intensity of the trade winds associated with a change of the South Atlantic anticyclone intensity and position. In southern Benguela, changes are also associated with variations in midlatitude low pressure systems and associated upwelling unfavorable westerly winds. La Niñas favor the development of Benguela Niños in Angola and Namibia. This study shows the potential for SST seasonal predictability in the Benguela upwelling due to the leading lag correlation between ENSO and the Benguela upwelling SST.
Abstract
Large-amplitude internal solitary wave (ISW) shoaling, breaking, and run-up was tracked continuously by a dense and rapidly sampling array spanning depths from 500 m to shore near Dongsha Atoll in the South China Sea. Incident ISW amplitudes ranged between 78 and 146 m with propagation speeds between 1.40 and 2.38 m s−1. The ratio between wave amplitude and a critical amplitude A 0 controlled breaking type and was related to wave speed cp and depth. Fissioning ISWs generated larger trailing elevation waves when the thermocline was deep and evolved into onshore propagating bores in depths near 100 m. Collapsing ISWs contained significant mixing and little upslope bore propagation. Bores contained significant onshore near-bottom kinetic and potential energy flux and significant offshore rundown and relaxation phases before and after the bore front passage, respectively. Bores on the shallow forereef drove bottom temperature variation in excess of 10°C and near-bottom cross-shore currents in excess of 0.4 m s−1. Bores decelerated upslope, consistent with upslope two-layer gravity current theory, though run-up extent Xr was offshore of the predicted gravity current location. Background stratification affected the bore run-up, with Xr farther offshore when the Korteweg–de Vries nonlinearity coefficient α was negative. Fronts associated with the shoaling local internal tide, but equal in magnitude to the soliton-generated bores, were observed onshore of 20-m depth.
Abstract
Large-amplitude internal solitary wave (ISW) shoaling, breaking, and run-up was tracked continuously by a dense and rapidly sampling array spanning depths from 500 m to shore near Dongsha Atoll in the South China Sea. Incident ISW amplitudes ranged between 78 and 146 m with propagation speeds between 1.40 and 2.38 m s−1. The ratio between wave amplitude and a critical amplitude A 0 controlled breaking type and was related to wave speed cp and depth. Fissioning ISWs generated larger trailing elevation waves when the thermocline was deep and evolved into onshore propagating bores in depths near 100 m. Collapsing ISWs contained significant mixing and little upslope bore propagation. Bores contained significant onshore near-bottom kinetic and potential energy flux and significant offshore rundown and relaxation phases before and after the bore front passage, respectively. Bores on the shallow forereef drove bottom temperature variation in excess of 10°C and near-bottom cross-shore currents in excess of 0.4 m s−1. Bores decelerated upslope, consistent with upslope two-layer gravity current theory, though run-up extent Xr was offshore of the predicted gravity current location. Background stratification affected the bore run-up, with Xr farther offshore when the Korteweg–de Vries nonlinearity coefficient α was negative. Fronts associated with the shoaling local internal tide, but equal in magnitude to the soliton-generated bores, were observed onshore of 20-m depth.
Abstract
Cross-shelf transport plays an important role in the heat, salt, and nutrient budgets of the continental shelf. In this study, we quantify cross-shelf volume transport and explore its dynamics within a high-resolution (2.5–6 km) regional ocean model of the East Australian Current (EAC) System, a western boundary current with a high level of mesoscale eddy activity. We find that the largest time-mean cross-shelf flows (>4 Sv per 100 km; 1 Sv ≡ 106 m3 s−1) occur inshore of the coherent western boundary current, between 26° and 30°S, while the strongest time-varying flows occur in the EAC southern extension, poleward of 32°S, associated with mesoscale eddies. Using a novel diagnostic equation derived from the momentum budget we show that the cross-shelf transport is dominated by the baroclinic and geostrophic component of the velocities, as the EAC jet is relatively free to flow over the variable shelfbreak topography. However, topographic interactions are also important and act through the bottom pressure torque term as a secondary driver of cross-shelf transport. The importance of topographic interaction also increases in shallower water inshore of the coherent jet. Downstream of separation, cross-shelf transport is more time-varying and associated with the interaction of mesoscale eddies with the shelf. The identification of the change in nature and drivers of cross-shelf transport in eddy versus jet dominated regimes may be applicable to understanding cross-shelf transport dynamics in other boundary current systems.
Significance Statement
Cross-shelf transport, i.e., the movement of water from the open ocean on or off the continental shelf, is not reported often as it is difficult to measure and model. We demonstrate a simple but effective method to do this and, using an ocean model, apply it to the East Australian Current System and show what drives it. The results show two distinct regimes, which differ depending on which part of the current system you are in. Our results help to place observations of cross-shelf transport in better context and provide a framework within which to consider the transport of other things such as heat and carbon from the open ocean to the continental shelf.
Abstract
Cross-shelf transport plays an important role in the heat, salt, and nutrient budgets of the continental shelf. In this study, we quantify cross-shelf volume transport and explore its dynamics within a high-resolution (2.5–6 km) regional ocean model of the East Australian Current (EAC) System, a western boundary current with a high level of mesoscale eddy activity. We find that the largest time-mean cross-shelf flows (>4 Sv per 100 km; 1 Sv ≡ 106 m3 s−1) occur inshore of the coherent western boundary current, between 26° and 30°S, while the strongest time-varying flows occur in the EAC southern extension, poleward of 32°S, associated with mesoscale eddies. Using a novel diagnostic equation derived from the momentum budget we show that the cross-shelf transport is dominated by the baroclinic and geostrophic component of the velocities, as the EAC jet is relatively free to flow over the variable shelfbreak topography. However, topographic interactions are also important and act through the bottom pressure torque term as a secondary driver of cross-shelf transport. The importance of topographic interaction also increases in shallower water inshore of the coherent jet. Downstream of separation, cross-shelf transport is more time-varying and associated with the interaction of mesoscale eddies with the shelf. The identification of the change in nature and drivers of cross-shelf transport in eddy versus jet dominated regimes may be applicable to understanding cross-shelf transport dynamics in other boundary current systems.
Significance Statement
Cross-shelf transport, i.e., the movement of water from the open ocean on or off the continental shelf, is not reported often as it is difficult to measure and model. We demonstrate a simple but effective method to do this and, using an ocean model, apply it to the East Australian Current System and show what drives it. The results show two distinct regimes, which differ depending on which part of the current system you are in. Our results help to place observations of cross-shelf transport in better context and provide a framework within which to consider the transport of other things such as heat and carbon from the open ocean to the continental shelf.
Abstract
The paper presents a combined numerical–deep learning (DL) approach for improving wind and wave forecasting. First, a DL model is trained to improve wind velocity forecasts by using past reanalysis data. The improved wind forecasts are used as forcing in a numerical wave forecasting model. This novel approach, used to combine physics-based and data-driven models, was tested over the Mediterranean. The correction to the wind forecast resulted in ∼10% RMSE improvement in both wind velocity and wave height over reanalysis data. This significant improvement is even more substantial at the Aegean Sea when Etesian winds are dominant, improving wave height forecasts by over 35%. The additional computational costs of the DL model are negligible compared to the costs of either the atmospheric or wave numerical model by itself. This work has the potential to greatly improve the wind and wave forecasting models used nowadays by tailoring models to localized seasonal conditions, at negligible additional computational costs.
Significance Statement
Wind and wave forecasting models solve a set of complicated physical equations. Improving forecasting accuracy is usually achieved by using a higher-resolution, empirical coefficients calibration or better physical formulations. However, measurements are rarely used directly to achieve better forecasts, as their assimilation can prove difficult. The presented work bridges this gap by using a data-driven deep learning model to improve wind forecasting accuracy, and the resulting wave forecasting. Testing over the Mediterranean Sea resulted in ∼10% RMSE improvement. Inspecting the Aegean Sea when the Etesian wind is dominant shows an outstanding 35% improvement. This approach has the potential to improve the operational atmospheric and wave forecasting models used nowadays by tailoring models to localized seasonal conditions, at negligible computational costs.
Abstract
The paper presents a combined numerical–deep learning (DL) approach for improving wind and wave forecasting. First, a DL model is trained to improve wind velocity forecasts by using past reanalysis data. The improved wind forecasts are used as forcing in a numerical wave forecasting model. This novel approach, used to combine physics-based and data-driven models, was tested over the Mediterranean. The correction to the wind forecast resulted in ∼10% RMSE improvement in both wind velocity and wave height over reanalysis data. This significant improvement is even more substantial at the Aegean Sea when Etesian winds are dominant, improving wave height forecasts by over 35%. The additional computational costs of the DL model are negligible compared to the costs of either the atmospheric or wave numerical model by itself. This work has the potential to greatly improve the wind and wave forecasting models used nowadays by tailoring models to localized seasonal conditions, at negligible additional computational costs.
Significance Statement
Wind and wave forecasting models solve a set of complicated physical equations. Improving forecasting accuracy is usually achieved by using a higher-resolution, empirical coefficients calibration or better physical formulations. However, measurements are rarely used directly to achieve better forecasts, as their assimilation can prove difficult. The presented work bridges this gap by using a data-driven deep learning model to improve wind forecasting accuracy, and the resulting wave forecasting. Testing over the Mediterranean Sea resulted in ∼10% RMSE improvement. Inspecting the Aegean Sea when the Etesian wind is dominant shows an outstanding 35% improvement. This approach has the potential to improve the operational atmospheric and wave forecasting models used nowadays by tailoring models to localized seasonal conditions, at negligible computational costs.
Abstract
Observational surface data are utilized to reconstruct the subsurface density and geostrophic velocity fields via the “interior + surface quasigeostrophic” (isQG) method in a subdomain of the Antarctic Circumpolar Current (ACC). The input variables include the satellite-derived sea surface height (SSH), satellite-derived sea surface temperature (SST), satellite-derived or Argo-based sea surface salinity (SSS), and a monthly estimate of the stratification. The density reconstruction is assessed against a newly released high-resolution in situ dataset that is collected by a southern elephant seal. The results show that the observed mesoscale structures are reasonably reconstructed. In the Argo-SSS-based experiment, pattern correlations between the reconstructed and observed density mostly exceed 0.8 in the upper 300 m. Uncertainties in the SSS products notably influence the isQG performance, and the Argo-SSS-based experiment yields better density reconstruction than the satellite-SSS-based one. Through the two-dimensional (2D) omega equation, we further employ the isQG reconstructions to diagnose the upper-ocean vertical velocities (denoted w isQG2D), which are then compared against the seal-data-based 2D diagnosis of w seal. Notable discrepancies are found between w isQG2D and w seal, primarily because the density reconstruction does not capture the seal-observed smaller-scale signals. Within several subtransects, the Argo-SSS-based w isQG2D reasonably reproduce the spatial structures of w seal, but present smaller magnitude. We also apply the isQG reconstructions to the 3D omega equation, and the 3D diagnosis of w isQG3D is very different from w isQG2D, indicating the limitations of the 2D diagnostic equation. With reduced uncertainties in satellite-derived products in the future, we expect the isQG framework to achieve better subsurface estimations.
Abstract
Observational surface data are utilized to reconstruct the subsurface density and geostrophic velocity fields via the “interior + surface quasigeostrophic” (isQG) method in a subdomain of the Antarctic Circumpolar Current (ACC). The input variables include the satellite-derived sea surface height (SSH), satellite-derived sea surface temperature (SST), satellite-derived or Argo-based sea surface salinity (SSS), and a monthly estimate of the stratification. The density reconstruction is assessed against a newly released high-resolution in situ dataset that is collected by a southern elephant seal. The results show that the observed mesoscale structures are reasonably reconstructed. In the Argo-SSS-based experiment, pattern correlations between the reconstructed and observed density mostly exceed 0.8 in the upper 300 m. Uncertainties in the SSS products notably influence the isQG performance, and the Argo-SSS-based experiment yields better density reconstruction than the satellite-SSS-based one. Through the two-dimensional (2D) omega equation, we further employ the isQG reconstructions to diagnose the upper-ocean vertical velocities (denoted w isQG2D), which are then compared against the seal-data-based 2D diagnosis of w seal. Notable discrepancies are found between w isQG2D and w seal, primarily because the density reconstruction does not capture the seal-observed smaller-scale signals. Within several subtransects, the Argo-SSS-based w isQG2D reasonably reproduce the spatial structures of w seal, but present smaller magnitude. We also apply the isQG reconstructions to the 3D omega equation, and the 3D diagnosis of w isQG3D is very different from w isQG2D, indicating the limitations of the 2D diagnostic equation. With reduced uncertainties in satellite-derived products in the future, we expect the isQG framework to achieve better subsurface estimations.