Search

Modern climate projections lack adequate spatial and temporal resolution due to computational constraints, leading to inaccuracies in representing critical processes like thunderstorms that occur on the sub-resolution scale. Hybrid methods combining physics with machine learning (ML) offer faster, higher fidelity climate simulations by outsourcing compute-hungry, high-resolution simulations to ML emulators. However, these hybrid ML-physics simulations require domain-specific data and workflows that have been inaccessible to many ML experts. As an extension of the ClimSim dataset (Yu et al., 2024), we present ClimSim-Online, which also includes an end-to-end workflow for developing hybrid ML-physics simulators. The ClimSim dataset includes 5.7 billion pairs of multivariate input/output vectors, capturing the influence of high-resolution, high-fidelity physics on a host climate simulator's macro-scale state. The dataset is global and spans ten years at a high sampling frequency. We provide a cross-platform, containerized pipeline to integrate ML models into operational climate simulators for hybrid testing. We also implement various ML baselines, alongside a hybrid baseline simulator, to highlight the ML challenges of building stable, skillful emulators. The data (https://huggingface.co/datasets/LEAP/ClimSim_high-res) and code (https://leap-stc.github.io/ClimSim and https://github.com/leap-stc/climsim-online) are publicly released to support the development of hybrid ML-physics and high-fidelity climate simulations., Comment: This manuscript is an expanded version of our paper that received the Outstanding Paper Award at the NeurIPS 2023 conference

The geographical variability, frequency content, and vertical structure of near-surface oceanic kinetic energy (KE) are important for air-sea interaction, marine ecosystems, operational oceanography, pollutant tracking, and interpreting remotely sensed velocity measurements. Here, KE in high-resolution global simulations (HYbrid Coordinate Ocean Model; HYCOM, and Massachusetts Institute of Technology general circulation model; MITgcm), at the sea surface (0 m) and 15 m, are respectively compared with KE from undrogued and drogued surface drifters. Global maps and zonal averages are computed for low-frequency ($<$ 0.5 cpd), near-inertial, diurnal, and semi-diurnal bands. Both models exhibit low-frequency equatorial KE that is low relative to drifter values. HYCOM near-inertial KE is higher than in MITgcm, and closer to drifter values, probably due to more frequently updated atmospheric forcing. HYCOM semi-diurnal KE is lower than in MITgcm, and closer to drifter values, likely due to inclusion of a parameterized topographic internal wave drag. A concurrent tidal harmonic analysis in the diurnal band demonstrates that much of the diurnal flow is non-tidal. We compute a simple proxy of near-surface vertical structure, the ratio of 0 m KE to 0 m KE plus 15 m KE in model outputs, and undrogued KE to undrogued KE plus drogued KE in drifter observations. Over most latitudes and frequency bands, model ratios track the drifter ratios to within error bars. Values of this ratio demonstrate significant vertical structure in all frequency bands except the semidiurnal band. Latitudinal dependence in the ratio is greatest in diurnal and low-frequency bands., Comment: revised for AGU JGR: Oceans

We describe the application of a new technique from nonlinear dynamical systems to infer the Lagrangian connectivity of the deep global ocean. We approximate the dynamic Laplacian using Argo trajectories from January 2011 to January 2017 and extract the eight dominant coherent (or dynamically self-connected) regions at 1500m depth. Our approach overcomes issues such as sparsity of observed data, and floats continually leaving and entering the dataset; only 10\% of floats record for the full six years. The identified coherent regions maximally trap water within them over the six-year time frame, providing a distinct analysis of the deep global ocean, and relevant information for planning future float deployment. While our study is concerned with ocean circulation at a multi-year, global scale, the dynamic Laplacian approach may be applied at any temporal or spatial scale to identify coherent structures in ocean flow from positional time series information arising from observations or models.

The transport by materially coherent eddies is studied in a two-layer quasigeostrophic model of geophysical turbulence. The coherent eddies are identified by closed contours of the Lagrangian-averaged vorticity deviation obtained from Lagrangian particles advected by the flow. A series of flow regimes with different bottom friction strengths are considered---it is found that coherent eddies become more prevalent and longer-lasting as the strength of the bottom drag increases. These coherent eddies, with average core radius close to the deformation radius, propagate zonally with speeds close to the long baroclinic Rossby wave speed and meridionally with a preference for cyclones to propagate poleward and anticyclones to propagate equatorward. The meridional propagation preference of the coherent eddies gives rise to a systematic upgradient potential vorticity (PV) transport, which is in the opposite direction as the background PV transport and not captured by standard Lagrangian diffusivity estimates. The upgradient PV transport by coherent eddy cores is less than 10% of the total PV transport, but the PV transport by the periphery flow induced by the PV inside coherent eddies is significant and downgradient. This clarifies the distinct roles of the trapping and stirring effect of coherent eddies in PV transport in geophysical turbulence.

In large-eddy simulations, subgrid-scale (SGS) processes are parameterized as a function of filtered grid-scale variables. First-order, algebraic SGS models are based on the eddy-viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse-grained velocity from direct numerical simulations (DNSs) of the atmospheric boundary layer at friction Reynolds numbers Re_{\tau} up to 1243 without invoking the eddy-viscosity assumption. The DNN model was found to produce higher correlation of SGS stresses compared to the Smagorinsky model and the Smagorinsky-Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The additional information on potential temperature and pressure were found not to be useful for SGS modeling. Deep learning thus demonstrates great potential for LESs of geophysical turbulence., Comment: 33 pages, 11 figures, 3 tables

We should be in a golden age of scientific discovery, given that we have more data and more compute power available than ever before, plus a new generation of algorithms that can learn effectively from data. But paradoxically, in many data-driven fields, the eureka moments are becoming increasingly rare. Scientists are struggling to keep pace with the explosion in the volume and complexity of scientific data. We describe here a few simple architectural principles that we believe are essential in order to create effective, robust, and flexible platforms that make the best use of emerging technology to deal with the exponential growth of scientific data.

Rotationally coherent Lagrangian vortices (RCLVs) are identified from satellite-derived surface geostrophic velocities in the Eastern Pacific (180$^\circ$-130$^\circ$ W) using the objective (frame-invariant) finite-time Lagrangian-coherent-structure detection method of Haller et al. (2016) based on the Lagrangian-averaged vorticity deviation. RCLVs are identified for 30, 90, and 270 day intervals over the entire satellite dataset, beginning in 1993. In contrast to structures identified using Eulerian eddy-tracking methods, the RCLVs maintain material coherence over the specified time intervals, making them suitable for material transport estimates. Statistics of RCLVs are compared to statistics of eddies identified from sea-surface height (SSH) by Chelton et al. 2011. RCLVs and SSH eddies are found to propagate westward at similar speeds at each latitude, consistent with the Rossby wave dispersion relation. However, RCLVs are uniformly smaller and shorter-lived than SSH eddies. A coherent eddy diffusivity is derived to quantify the contribution of RCLVs to meridional transport; it is found that RCLVs contribute less than 1% to net meridional dispersion and diffusion in this sector, implying that eddy transport of tracers is mostly due to incoherent motions, such as swirling and filamentation outside of the eddy cores, rather than coherent meridional translation of eddies themselves. These findings call into question prior estimates of coherent eddy transport based on Eulerian eddy identification methods.