skip to main content


Title: An evaluation of the LLC4320 global-ocean simulation based on the submesoscale structure of modeled sea surface temperature fields
Abstract. We have assembled 2 851 702 nearly cloud-free cutout images (sized 144 km × 144 km) of sea surface temperature (SST) data from the entire 2012–2020 Level-2 Visible Infrared Imaging Radiometer Suite (VIIRS) dataset to perform a quantitative comparison to the ocean model output from the MIT General Circulation Model (MITgcm). Specifically, we evaluate outputs from the LLC4320 (LLC, latitude–longitude–polar cap) 148∘ global-ocean simulation for a 1-year period starting on 17 November 2011 but otherwise matched in geography and the day of the year to the VIIRS observations. In lieu of simple (e.g., mean, standard deviation) or complex (e.g., power spectrum) statistics, we analyze the cutouts of SST anomalies with an unsupervised probabilistic autoencoder (PAE) trained to learn the distribution of structures in SST anomaly (SSTa) on ∼ 10–80 km scales (i.e., submesoscale to mesoscale). A principal finding is that the LLC4320 simulation reproduces, over a large fraction of the ocean, the observed distribution of SSTa patterns well, both globally and regionally. Globally, the medians of the structure distributions match to within 2σ for 65 % of the ocean, despite a modest, latitude-dependent offset. Regionally, the model outputs reproduce mesoscale variations in SSTa patterns revealed by the PAE in the VIIRS data, including subtle features imprinted by variations in bathymetry. We also identify significant differences in the distribution of SSTa patterns in several regions: (1) in an equatorial band equatorward of 15∘; (2) in the Antarctic Circumpolar Current (ACC), especially in the eastern half of the Indian Ocean; and (3) in the vicinity of the point at which western boundary currents separate from the continental margin. It is clear that region 3 is a result of premature separation in the simulated western boundary currents. The model output in region 2, the southern Indian Ocean, tends to predict more structure than observed, perhaps arising from a misrepresentation of the mixed layer or of energy dissipation and stirring in the simulation. The differences in region 1, the equatorial band, are also likely due to model errors, perhaps arising from the shortness of the simulation or from the lack of high-frequency and/or wavenumber atmospheric forcing. Although we do not yet know the exact causes for these model–data SSTa differences, we expect that this type of comparison will help guide future developments of high-resolution global-ocean simulations.  more » « less
Award ID(s):
2150228
NSF-PAR ID:
10484222
Author(s) / Creator(s):
; ; ; ;
Editor(s):
Ham, David; Añel, Juan Antonio; Kerkweg, Astrid; Lo, Min-Hui; Neale, Richard; Sander, Rolf; Ullrich, Paul
Publisher / Repository:
European Scientific Union
Date Published:
Journal Name:
Geoscientific Model Development
Edition / Version:
1
Volume:
16
Issue:
23
ISSN:
1991-9603
Page Range / eLocation ID:
7143 to 7170
Subject(s) / Keyword(s):
["sea surface temperature","remote sensing"]
Format(s):
Medium: X Size: 26.5MB Other: pdf
Size(s):
["26.5MB"]
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Terrestrial‐marine dust fluxes, pedogenic carbonate δ13C values, and various paleovegetation proxies suggest that Africa experienced gradual cooling and drying across the Pliocene‐Pleistocene (Plio‐Pleistocene) boundary (2.58 million years ago [Ma]). However, the timing, magnitude, resolution, and relative influences of orbitally‐driven changes in high latitude glaciations and low latitude insolation differ by region and proxy. To disentangle these forcings and investigate equatorial eastern African climate across the Plio‐Pleistocene boundary, we generated a high‐resolution (∼3,000‐year) data set of compound‐specificn‐alkane leaf wax δ2H values—a robust proxy for atmospheric circulation and precipitation amount—from the HSPDP‐BTB13‐1A core, which spans a ∼3.3–2.6 Ma sequence in the Baringo‐Tugen Hills‐Barsemoi Basin of central Kenya. In combination with the physical sedimentology, our data indicate that precipitation varied strongly with orbital obliquity, not precession, during the late Pliocene, perhaps imparted by variations in the cross‐equatorial insolation gradient. We also observe a marked shift toward wetter conditions beginning ∼3 Ma that corresponds with global cooling, drying in western Australia, and a steepening of the west‐east zonal Indian Ocean (IO) sea surface temperature (SST) gradient. We propose that northward migration of the Subtropical Front reduced Agulhas current leakage, warming the western IO and causing changes in the IO zonal SST gradient at 3 Ma, a process that has been observed in the latest Pleistocene‐Holocene but not over longer timescales. Thus, the late Cenozoic moisture history of eastern Africa is driven by a complex mixture of low‐latitude insolation, the IO SST gradient, and teleconnections to distal high‐latitude cooling.

     
    more » « less
  2. Abstract

    Thermospheric mass density perturbations are commonly observed during geomagnetic storms and fundamental to upper atmosphere dynamics, but the sources of these perturbations are not well understood. Large neutral density perturbations during storms greatly affect the drag experienced by low Earth orbit. We investigated the thermospheric density perturbations at all latitudes observed along the CHAMP and GRACE satellite trajectories during the August 24–25, 2005 geomagnetic storm. Observations show that large neutral density enhancements occurred not only at high latitudes, but also globally. Large density perturbations were seen in the equatorial regions away from the high‐latitude, magnetospheric energy sources. We used the high‐resolution Multiscale Atmosphere Geospace Environment (MAGE) model to simulate consecutive neutral density changes observed by satellites during the storm. The MAGE simulation, which resolved mesoscale high‐latitude convection electric fields and field‐aligned currents, and included physics‐based specification of auroral precipitation, was contrasted with a standalone ionosphere‐thermosphere simulation driven by a high‐latitude electrodynamics empirical model. The comparison demonstrates that first‐principles representations of highly dynamic and localized Joule heating events in a fully coupled whole geospace model is critical to accurately capture both generation and propagation of traveling atmospheric disturbances (TADs) that produce neutral density perturbations globally. The MAGE simulation shows that larger density peaks in the equatorial region observed by CHAMP and GRACE are the result of TADs generated at high‐latitudes in both hemispheres, and intersect at low‐latitudes. This study reveals the importance of investigating thermospheric density variations at all latitudes in a fully coupled geospace model with sufficiently high resolving power.

     
    more » « less
  3. Abstract. Wind work at the air-sea interface is the transfer of kinetic energy between the ocean and the atmosphere and, as such, is an important part of the ocean-atmosphere coupled system. Wind work is defined as the scalar product of ocean wind stress and surface current, with each of these two variables spanning, in this study, a broad range of spatial and temporal scales, from 10 km to more than 3000 km and hours to months. These characteristics emphasize wind work's multiscale nature. In the absence of appropriate global observations, our study makes use of a new global, coupled ocean-atmosphere simulation, with horizontal grid spacing of 2–5 km for the ocean and 7 km for the atmosphere, analyzed for 12 months.We develop a methodology, both in physical and spectral spaces, to diagnose three different components of wind work that force distinct classes of ocean motions, including high-frequency internal gravity waves, such as near-inertial oscillations, low-frequency currents such as those associated with eddies, and seasonally averaged currents, such as zonal tropical and equatorial jets.The total wind work, integrated globally, has a magnitude close to 5 TW, a value that matches recent estimates. Each of the first two components that force high-frequency and low-frequency currents, accounts for ∼ 28 % of the total wind work and the third one that forces seasonally averaged currents, ∼ 44 %. These three components, when integrated globally, weakly vary with seasons but their spatial distribution over the oceans has strong seasonal and latitudinal variations. In addition, the high-frequency component that forces internal gravity waves, is highly sensitive to the collocation in space and time (at scales of a few hours) of wind stresses and ocean currents. Furthermore, the low-frequency wind work component acts to dampen currents with a size smaller than 250 km and strengthen currents with larger sizes. This emphasizes the need to perform a full kinetic budget involving the wind work and nonlinear advection terms as small and larger-scale low-frequency currents interact through these nonlinear terms.The complex interplay of surface wind stresses and currents revealed by the numerical simulation motivates the need for winds and currents satellite missions to directly observe wind work.

     
    more » « less
  4. null (Ed.)
    Abstract The Indian Ocean has received increasing attention for its large impacts on regional and global climate. However, sea surface temperature (SST) variability arising from Indian Ocean internal processes has not been well understood particularly on decadal and longer timescales, and the external influence from the Tropical Pacific has not been quantified. This paper analyzes the interannual-to-decadal SST variability in the Tropical Indian Ocean in observations and explores the external influence from the Pacific versus internal processes within the Indian Ocean using a Linear Inverse Model (LIM). Coupling between Indian Ocean and tropical Pacific SST anomalies (SSTAs) is assessed both within the LIM dynamical operator and the unpredictable stochastic noise that forces the system. Results show that the observed Indian Ocean Basin (IOB)-wide SSTA pattern is largely a response to the Pacific ENSO forcing, although it in turn has a damping effect on ENSO especially on annual and decadal timescales. On the other hand, the Indian Ocean Dipole (IOD) is an Indian Ocean internal mode that can actively affect ENSO; ENSO also has a returning effect on the IOD, which is rather weak on decadal timescale. The third mode is partly associated with the Subtropical Indian Ocean Dipole (SIOD), and it is primarily generated by Indian Ocean internal processes, although a small component of it is coupled with ENSO. Overall, the amplitude of Indian Ocean internally generated SST variability is comparable to that forced by ENSO, and the Indian Ocean tends to actively influence the tropical Pacific. These results suggest that the Indian-Pacific Ocean interaction is a two-way process. 
    more » « less
  5. Abstract The Sea Surface Temperature Anomaly (SSTA) in tropical Atlantic during boreal spring and summer shows two dominant modes: a basin-warming and a meridional dipole mode, respectively. Observational and coupled model simulations indicate that the former induces a Pacific La Niña in the succeeding winter whereas the latter cannot. The basin-warming forcing induces a La Niña through a Kelvin wave response and the associated wind-evaporation-SST-convection (WESC) feedback over the northern Indian Ocean (NIO) and Maritime Continent (MC). Anomalous Kelvin wave easterly interacts with the monsoonal westerly, leading to a warm SSTA and a northwest-southeast oriented heating anomaly in NIO/MC, which further induces easterly and cold SSTA over the equatorial Pacific. In contrast, the dipole forcing has little impact on the Indian and Pacific Oceans due to the offsetting of the Kelvin wave to the asymmetric Atlantic heating. Further observational and modeling studies towards the Tropical North Atlantic (TNA) and Equatorial Atlantic (EA) SSTA modes indicate that the TNA (EA) forcing induces a CP- (EP-) type ENSO. In both cases, the Kelvin wave response and the WESC feedback over the NIO/MC are important in conveying the Atlantic’s impact. The difference lies in distinctive Rossby wave responses – A marked westerly anomaly appears in the equatorial eastern Pacific (EEP) for the TNA forcing (due to its westward location) while no significant wind response is observed in EEP for the EA forcing. The westerly anomaly prevents a cooling tendency in EEP through anomalous zonal and vertical advection according to a mixed-layer heat budget analysis. 
    more » « less