Ocean Model Research Articles

A regional ocean–atmosphere coupled model using CMA-TRAMS and LICOM: Preliminary results for tropical cyclone gale prediction over the northern South China Sea

This paper provides a comparative analysis of the performance of a high-resolution regional ocean–atmosphere coupled model in predicting tropical cyclone (TC) gales over the northern South China Sea. The atmosphere and ocean components of the coupled system are represented by the China Meteorological Administration's Tropical Regional Atmosphere Model for the South China Sea (CMA-TRAMS) and the LASG/IAP Climate system Ocean Model (LICOM), respectively. The Ocean Atmosphere Sea Ice Soil VersionH 3 (OASIS3) software has been utilized for the exchange of momentum, heat, and freshwater fluxes between these two components. An assessment of the coupled model's three-day predictions for five TCs’ gales was conducted. Preliminary findings indicate that the predicted TC tracks show less sensitivity to oceanic influences than the predicted TC intensities. Significant improvement in predicting the surface TC gales has been achieved through coupling the ocean model. This improvement is attributed to the impact of the warmer ocean's effect on TC intensification, counteracting the cooling effect of the cold wake. In summary, coupling has enhanced the model's predictive capabilities for TC gales. A detailed assessment of the coupled model's performance in predicting other tropical weather phenomena is forthcoming.摘要本文对一个具有我国自主知识产权的高分辨率区域海气耦合模式在南海北部台风大风预报性能进行对比分析. 该耦合模式的大气和海洋部分分别为中国气象局南海台风模式 (CMA-TRAMS) 和中国科学院大气物理研究所的海洋模式 (LICOM), 并用OASIS耦合器实现海气界面动量, 热量和淡水通量的交换. 本文对海气耦合模式预报5个台风的24, 48和72小时大风的性能进行了检验评估. 初步结果表明, 预报的台风路径对海洋影响的敏感性低于台风强度. 耦合模式使预报台风大风的准确性得到显著提高, 这一改进归结于海气耦合模式更为准确的海表温度对台风增强的影响, 从而抵消了冷尾流的降温效应. 总而言之, 耦合模式增强了模式对台风大风的预报能力. 关于海气耦合模式在其他热带天气中的预报能力评估将在未来工作中展示.

Detecting Instantaneous Tidal Signals in Ocean Models Utilizing Streaming Band‐Pass Filters

AbstractThrough the implementation of a streaming filter, output of numerical ocean simulations can be band‐pass filtered at tidal frequencies while the model is running, yielding time series of sinusoidal motions consisting of tidal signals in the filter's target frequency band. The filtering algorithm is developed from a system of two ordinary differential equations that represents the motion of a damped harmonic oscillator. The filter's response to a broadband input signal is unity at its target frequency but vanishes toward the low and high frequency limits. The decay of the filter response is controlled by a dimensionless parameter, which determines the filter's bandwidth. As a result, the filter allows signals within a small frequency band around its target frequency to pass through, while blocking signals outside of its target frequency band. In this work, the filtering algorithm is implemented into the barotropic solver of the Modular Ocean Model version 6 (MOM6) for determining the instantaneous tidal velocities of the semi‐diurnal and diurnal tides. Utilizing the filters, the frequency‐dependent internal wave drag is applied to the semi‐diurnal and diurnal frequency bands separately. The simulation results suggest that the performance of the algorithm is consistent with the filter transfer function in Fourier space. Potential applications of the algorithm also include de‐tiding the model output for nested regional ocean models, especially those for the purpose of operational forecasting.

Sensitivity of Simulated Arctic Ocean Salinity and Strait Transport to Interannually Variable Hydrologic Model Based Runoff

AbstractAs the Arctic warms at an increased rate compared to the rest of the globe, freshwater runoff has been shown to be increasing into the Arctic Ocean. The effects of this contemporary increase in riverine freshwater into the Arctic Ocean are estimated from ocean model simulations, using two runoff data sets. One runoff data set is based on older climatological data, which has no inter‐annual variability after 2007 and as such does not represent the observed increases in river runoff into the Arctic. The other data set comes from a hydrological model developed for the Arctic drainage basin, which includes contemporary changes in the climate. At the Pan‐Arctic scale this new data set represents an approximately 11% increase in runoff, compared with the older climatological data. Comparing two ocean model runs forced with the different runoff data sets, overall changes in different freshwater markers across the basin were found to be between 5% and 10%, depending on the regions. The strongest increases were seen from the Siberian rivers, which in turn caused the strongest freshening in the Eastern Arctic. As the surface waters of the Arctic Ocean are sensitive to runoff, incorporating hydrological model data can help to better understand current changes and potential future impacts from increased runoff with climate change.

Structure and dynamics of the near-surface water layer under conditions of natural convection

A series of laboratory and numerical experiments have been performed to study the process of the surface water cooling. The process of the submergence of an ensemble of jets (a thermal-like structure) or a negative buoyancy flux was recorded using digital cameras (∼2 and 10 Hz). Quantitative estimates of the structure and dynamical features of the initial phase of the submergence of a density front of the ensemble of jets were obtained. These estimates were compared with the previous results of the study of individual thermals. It was found that the deepening velocity is characterized by the nonlinear dependence from the submergence depth in time of both the thermals and the density front of the ensemble. The numerical calculations also showed nonlinearity in this dependence. The process of submergence of the density front was divided into three stages, which were described for the first time. The deepening velocity of the density front of the ensemble of jets was ∼0.15 cm/s. The results are qualitatively consistent with the theoretical estimates of the evolution of the free convection. The observed laboratory flows have the laminar nature of convective movements (Reynolds numbers ranging from 3 to 60). A geometric model was proposed to describe linear quantity of both the thermal and jet. It was shown that the effects of turbulence on the seawater measurements were caused externally and not related to the process of the surface water cooling. Thus, it is assumed that either thermal or saline turbulence dominates over free convection conditions. This conclusion has an important role in the study of the heat exchange process in ocean models.

A regional NEMO 4.0 configuration of the subpolar North Atlantic

The dynamics of the North Atlantic subpolar gyre (SPG) largely control the Atlantic Meridional Overturning Circulation (AMOC), thus affecting the long-term climate variability of the global ocean. Many key processes in the North Atlantic SPG, such as Gulf Stream detachment, eddy activity, deep convection and overflows, are frequently incorrectly presented or misrepresented in modern ocean general circulation models. Here, we present a new regional eddy-resolving model of the subpolar North Atlantic based on the NEMO4 configuration and analyze its ability to adequately represent the dynamics of the subpolar North Atlantic. The model performance is assessed using in situ data (cross-section hydrography, Argo and sea gliders), remote sensing data (satellite sea surface temperature and salinity) and the ANDRO and ARMOR3D products based on optimal interpolation products. In the next step, we analyze the effects of high-resolution forcing and several recent developments in ocean modeling to improve the model solution. The implementation of the self-diffusive momentum advection scheme with additional viscosity allows for the simulation of more realistic Deep Western Boundary Current and Irminger Rings lifecycle characteristics. The surface current feedback parameterization implemented in the forcing function and the modification of the kinetic energy budget in the upper ocean, especially in intensively eddying areas, simulate the behavior predicted by coupled models. While the model is sensitive to the implementation of cool skin–warm layer parameterization, this parameterization does not significantly improve the solution compared to the observational data. Regional implementation of the local-sigma vertical coordinate system in cascading areas allows overflow waters to spread over the ocean bottom more realistically. The remote effect of the implementation of the local-sigma coordinate is traced in the vertical thermohaline structure at the western slopes of the Irminger Sea and Iceland Basin as well as in the Deep Western Boundary Current in the Labrador Sea. Compared to z-partial step setting, local-sigma allows for correct representation of continuous densification of the water in the overturning part of the AMOC and thus North Atlantic Deep-Water formation. The potential of the new model configuration for analyzing the subpolar North Atlantic, including convective activity, is discussed.

Estimation of tidal energy potential in the Vietnam East Sea: A comprehensive analysis using semi-empirical tide models

This study assesses the tidal energy potential in the Vietnam East Sea (VES), a region with complex hydrodynamics and significant interactions with the Western Pacific Ocean. A rigorous validation of six prominent semi-empirical global ocean tide models (FES2014, DTU16, EOT20, GOT4.10c, HAMTIDE12, OSU12.V1) was conducted against an extensive network of 34 tide gauge stations in the VES region. The validation identified FES2014 as the superior model, exhibiting the closest agreement with observations due to its advanced data assimilation techniques, high-resolution grids, and robust hydrodynamic representation influenced by the complex VES bathymetry. Leveraging the FES2014 model, the influential role of the semi-diurnal component (M2) in dictating energy distribution is prominently evident. Evaluation results of tidal amplitudes revealed three emerging hotspots: the Batang Lupar estuary (Malaysia), Anpu Gang (China), and Bac Lieu (Vietnam), with tidal amplitudes up to 326.2 cm. The Anpu Gang exhibited the highest theoretical potential tidal energy of 14.90 Wh/m2, making it the most promising site for localized tidal power development. Although moderate compared to globally renowned sites, the identified VES hotspots merit consideration for small-to-medium scale projects tailored to local conditions. While limited to potential energy assessments, this study provides a crucial baseline for the VES region, highlighting opportunities for sustainable tidal energy exploitation.

GRACE Satellite Observations of Antarctic Bottom Water Transport Variability

AbstractAntarctic Bottom Water (AABW) formation and transport constitute a key component of the global ocean circulation. Direct observations suggest that AABW volumes and transport rates may be decreasing, but these observations are too temporally or spatially sparse to determine the cause. To address this problem, we develop a new method to reconstruct AABW transport variability using data from the GRACE (Gravity Recovery and Climate Experiment) satellite mission. We use an ocean general circulation model to investigate the relationship between ocean bottom pressure and AABW: we calculate both of these quantities in the model, and link them using a regularized linear regression. Our reconstruction from modeled ocean bottom pressure can capture 65%–90% of modeled AABW transport variability, depending on the ocean basin. When realistic observational uncertainty values are added to the modeled ocean bottom pressure, the reconstruction can still capture 30%–80% of AABW transport variability. Using the same regression values, the reconstruction skill is within the same range in a second, independent, general circulation model. We conclude that our reconstruction method is not unique to the model in which it was developed and can be applied to GRACE satellite observations of ocean bottom pressure. These advances allow us to create the first global reconstruction of AABW transport variability over the satellite era. Our reconstruction provides information on the interannual variability of AABW transport, but more accurate observations are needed to discern AABW transport trends.

A Global Ocean Opal Ballasting–Silicate Relationship

AbstractOpal and calcium carbonate are thought to regulate the biological pump's transfer of organic carbon to the deep ocean. A global sediment trap database exhibits large regional variations in the organic carbon flux associated with opal flux. These variations are well‐explained by upper ocean silicate concentrations, with high opal ‘ballasting’ in the silicate‐deplete tropical Atlantic Ocean, and low ballasting in the silicate‐rich Southern Ocean. A plausible, testable hypothesis is that opal ballasting varies because diatoms grow thicker frustules where silicate concentrations are higher, carrying less organic carbon per unit opal. The observed pattern does not fully emerge in an advanced ocean biogeochemical model when diatom silicification is represented using a single global parameterization as a function of silicate and iron. Our results suggest a need for improving understanding of currently modeled processes and/or considering additional parameterizations to capture the links between elemental cycles and future biological pump changes.

Horizontal and vertical dispersion in a wind-driven oceanic gyre model

This study addresses the horizontal and vertical dispersion of passive tracers in idealized wind-driven subtropical gyres. Synthetic particles within a closed basin are numerically advected to analyze their dispersion under different theoretical velocity fields. Horizontal dispersion simulations incorporate the classic wind-driven Stommel circulation along with (i) surface Ekman drift associated with the Stommel wind field and (ii) inertial effects due to particle size and buoyancy. Results reveal that the Ekman drift inhibits particle dispersion across the entire domain leading to tracer concentration in a quasi-stable distribution skewed toward the western side of the basin. Similar behavior is observed with inertial particles. The equilibrium state is quantified for different diffusivity values, particle sizes, and buoyancies. For vertical dispersion, simulations incorporate the three-dimensional Ekman velocity, which includes a negative vertical component, while ignoring inertial effects. Initially, surface particles accumulate around the gyre center while slowly sinking, but they disperse across the basin once they surpass the Ekman layer and are free from surface effects. Tracers sink more on the western side of the basin, regardless of horizontal diffusivity. On average, ignoring inertial effects, particles sink less with higher diffusivity and more with lower diffusivity, suggesting a potential for high horizontal distribution of sunken tracers in the ocean.

Seasonality of Submesoscale Vertical Heat Transport Modulated by Oceanic Mesoscale Eddies in the Kuroshio Extension

AbstractEnergetic mesoscale eddies are often accompanied by strong submesoscale variability, which plays a significant role in connecting mesoscale and turbulent motions in the ocean and leads to strong vertical motions. The product of a high‐resolution (1/48°) oceanic numerical model, the LLC4320, is employed to investigate the seasonal variations of vertical heat transport induced by submesoscale processes within multiple mesoscale eddies in the Kuroshio Extension (KE) region. In different seasons, the submesoscale vertical heat transport exhibits a consistent upward pattern, with notably higher magnitudes observed during winter. In winter, the maxima value of submesoscale vertical heat flux (SVHF) can account for approximately 60% of the total vertical heat flux (VHF). This is equivalent to the average net sea surface heat flux in a single eddy region. In summer and autumn, the maxima absolute value of submesoscale vertical heat flux can account for approximately 30% of the total VHF. Energy analysis reveals that baroclinic instability associated with vertical buoyancy flux has a crucial effect on generating submesoscale processes within the eddy region. The submesoscale motions are influenced by the mixed layer instability, strain‐induced frontogenesis, turbulent thermal wind and turbulent thermal wind‐induced frontogenesis within the upper mixed layer, while they are largely associated with the strain‐induced frontogenesis in the ocean interior. Furthermore, the upward low‐frequency submesoscale vertical heat transport is generated by submesoscale secondary circulation at eddy peripheries.

ALBATROSS: Advancing Southern Ocean tide modelling with high resolution and enhanced bathymetry

The knowledge of bathymetry and ocean tides plays a pivotal role at the crossroads of various scientific fields, especially in the Polar regions. Its significance extends to ocean circulation modeling and understanding the coupled dynamical response of the ocean, sea-ice and ice-sheet systems. In the Southern Ocean, conventional satellite altimetry measurements are rare below the 66° parallel. Hydrodynamic models are thus useful tools to provide spatially continuous information about ocean tides. However, the accuracy of ocean tide models around the Antarctic continent is currently limited by the quality of bathymetry. Recent reprocessing of decade-long CryoSat-2 data has facilitated a new computation of bathymetry around Antarctica, bringing innovative information on bathymetry gradients. This, combined with new compilations of bathymetry, ice draft, coastline, and grounding line datasets in ice-shelf regions, allows improving models and knowledge of ocean tides in the Southern Ocean. We developed a new high-resolution tidal model that implements the improved bathymetry data and includes data assimilation of satellite-altimetry tidal retrievals computed from CryoSat-2, filling the gap between the 66°S-limited coverage of the TOPEX-Jason suite missions and the Antarctic coast. Comparisons with tidal estimates derived from tide gauge measurements showed very good consistencies with an RMSE of 3 cm.

Ocean Alkalinity Enhancement in Deep Water Formation Regions Under Low and High Emission Pathways

AbstractOcean Alkalinity Enhancement (OAE) is an ocean‐based Carbon Dioxide Removal (CDR) method to mitigate climate change. Studies to characterize regional differences in OAE efficiencies and biogeochemical effects are still sparse. As subduction regions play a pivotal role for anthropogenic carbon uptake and centennial storage, we here evaluate OAE efficiencies in the subduction regions of the Southern Ocean, the Northwest Atlantic, and the Norwegian‐Barents Sea region. Using the ocean biogeochemistry model FESOM2.1‐REcoM3, we simulate continuous OAE globally and in the subduction regions under high (SSP3‐7.0) and low (SSP1‐2.6) emission scenarios. The OAE efficiency calculated by two different metrics is higher (by 8%–30%) for SSP3‐7.0 than for SSP1‐2.6 due to a lower buffer factor in a high‐ world. All subduction regions show a CDR potential (0.23–0.31; PgC uptake per Pg alkaline material) consistent with global OAE for both emission scenarios. Calculating the efficiency as the ratio of excess dissolved inorganic carbon (DIC) to excess alkalinity shows that the Southern Ocean and the Northwest Atlantic are as efficient as the global ocean (0.79–0.85), while the Norwegian‐Barents Sea region has a lower efficiency (0.65–0.75). The subduction regions store a fraction of excess carbon below 1 km that is 1.9 times higher than the global ocean. The excess surface alkalinity and thus uptake and storage follow the mixed‐layer depth seasonality, with the majority of the excess flux occurring in summer at shallow mixed layer depths. This study therefore highlights that subduction regions can be efficient for OAE if optimal deployment strategies are developed.

The effect of Southern Ocean topography on the global MOC and abyssal water mass distribution

Abstract We investigate the role of Southern Ocean topography and wind stress in the deep and abyssal ocean overturning and water mass composition using a suite of idealized global ocean circulation models. Specifically, we address how the presence of a meridional ridge in the vicinity of Drake Passage and the formation of an associated Southern Ocean gyre influences the water mass composition of the abyssal cell. Our experiments are carried out using a numerical representation of the global ocean circulation in an idealized two-basin geometry under varying wind-stress and Drake Passage ridge height. In the presence of a low Drake Passage ridge the overall strength of the meridional overturning circulation is primarily influenced by wind stress, with a topographically-induced weakening of the mid-depth cell and concurrent strengthening of the abyssal cell occurring only after ridge height passes 2500m. Passive tracer experiments show that a strengthening mid-depth cell leads to increased abyssal ventilation by North Atlantic water masses, as more North Atlantic Deep Water (NADW) enters the Southern Ocean and then spreads into the Indo-Pacific. We repeat our tracer experiments without restoring in the high-latitude Southern Ocean in order to identify the origin of water masses that circulate through the Southern Ocean before sinking into the abyss as Antarctic Bottom Water. Our results from these “exchange” tracer experiments show that an increasing ridge height in Drake Passage and the concurrent gyre spin-up lead to substantially decreased NADW-origin waters in the abyssal ocean, as more surface waters from north of the ACC are transferred into the Antarctic Bottom Water formation region.

Long-term variation of the eddy kinetic energy in the Northeastern South China sea

The seasonal to interannual variability of eddy kinetic energy (EKE) in the Northeastern South China Sea (NE-SCS) has been widely studied and it is recognized that they are strongly related to the state of the Kuroshio pathway in the Luzon Strait. While, due to the lack of long-term observations and high-resolution simulations, the decadal change of EKE in NE-SCS remains unexplored. In this study, we show the EKE trend in the past ∼ 30 years in the NE-SCS by using satellite observation and global HYbrid Coordinate Ocean Model reanalysis with the Navy Coupled Ocean Data Assimilation. It is found that due to the weakening of the Kuroshio in the Luzon Strait since 1990 s, the Kuroshio shows an enhanced looping path in the NE-SCS, inducing stronger EKE in this region. Further analysis confirms that the energy transfer by baroclinic instability is dominant for the increasing of EKE, when the Kuroshio intrudes into the NE-SCS and brings more potential energy inside the circulation. The Kuroshio state along the Luzon Strait is the key for modulating the EKE in the NE-SCS. Furthermore, the long-term weakening of Kuroshio current along the Luzon strait during 1993–2020 is determined by the decreasing of subtropical mode water, corresponding to the positive phase of the Atlantic Multidecadal Oscillation. This study provides insight into the interaction between marginal sea (i.e., the SCS) and the open ocean (i.e., the western Pacific Ocean), finally linking to the global climate change.

The Role of Benthic Fluxes in Acidifying the Bottom Waters in the Northern Gulf of Mexico Hypoxic Zone Based on an Updated Water Column Biogeochemical‐Seabed Diagenetic and Sediment Transport Model

AbstractThe seabed and the water column are tightly coupled in shallow coastal environments. Numerical models of seabed‐water interaction provide an alternative to observational studies that require concurrent measurements in both compartments, which are hard to obtain and rarely available. Here, we present a coupled model that includes water column biogeochemistry, seabed diagenesis, sediment transport and hydrodynamics. Our model includes realistic representations of biogeochemical reactions in both seabed and water column, and fluxes at their interface. The model was built on algorithms for seabed‐water exchange in the Regional Ocean Modeling System and expanded to include carbonate chemistry in seabed. The updated model was tested for two sites where benthic flux and porewater concentration measurements were available in the northern Gulf of Mexico hypoxic zone. The calibrated model reproduced the porewater concentration‐depth profiles and benthic fluxes of O2, dissolved inorganic carbon (DIC), TAlk, NO3 and NH4. We used the calibrated model to explore the role of benthic fluxes in acidifying bottom water during fair weather and resuspension periods. Under fair weather conditions, model results indicated that bio‐diffusion in sediment, labile material input and sediment porosity have a large control on the importance of benthic flux to bottom water acidification. During resuspension, the model indicated that bottom water acidification would be enhanced due to the sharp increase of the DIC/TAlk ratio of benthic fluxes. To conclude, our model reproduced the seabed‐water column exchange of biologically important solutes and can be used for quantifying the role of benthic fluxes in driving bottom water acidification over continental shelves.

Effects of urban eutrophication on pelagic habitat capacity in the Southern California Bight

Land-based nutrient inputs to the ocean have been linked to increased coastal productivity, subsurface acidification and O2 loss, even in upwelling systems like the Southern California Bight. However, whether eutrophication alters the [environment’s] capacity to support key taxa has yet to be evaluated for this region. Here, we assess the impact of land-based nutrient inputs on the availability of aerobic and calcifying habitat for key pelagic taxa using ocean model simulations. We find that acute, lethal conditions are not commonly induced in epipelagic surface waters, but that sublethal, ecologically relevant changes are pervasive. Land-based nutrient inputs reduce the potential aerobic and calcifier habitat during late summer, when viable habitat is at its seasonal minimum. A region of annually recurring habitat compression is predicted 30 – 90 km from the mainland, southeast of Santa Catalina Island. Here, both aerobic and calcifier habitat is vertically compressed by, on average, 25%, but can be as much as 60%. This effect can be traced to enhanced remineralization of organic matter that originates from the coast. These findings suggest that effects of land-based nutrients are not restricted to chemistry but extend to habitat capacity for multiple taxa of ecological and economic importance. Considerable uncertainty exists, however, in how this habitat compression translates to population-level effects.

A Stable Implementation of a Data‐Driven Scale‐Aware Mesoscale Parameterization

AbstractOcean mesoscale eddies are often poorly represented in climate models, and therefore, their effects on the large scale circulation must be parameterized. Traditional parameterizations, which represent the bulk effect of the unresolved eddies, can be improved with new subgrid models learned directly from data. Zanna and Bolton (2020), https://doi.org/10.1029/2020gl088376 (ZB20) applied an equation‐discovery algorithm to reveal an interpretable expression parameterizing the subgrid momentum fluxes by mesoscale eddies through the components of the velocity‐gradient tensor. In this work, we implement the ZB20 parameterization into the primitive‐equation GFDL MOM6 ocean model and test it in two idealized configurations with significantly different dynamical regimes and topography. The original parameterization was found to generate excessive numerical noise near the grid scale. We propose two filtering approaches to avoid the numerical issues and additionally enhance the strength of large‐scale energy backscatter. The filtered ZB20 parameterizations led to improved climatological mean state and energy distributions, compared to the current state‐of‐the‐art energy backscatter parameterizations. The filtered ZB20 parameterizations are scale‐aware and, consequently, can be used with a single value of the non‐dimensional scaling coefficient for a range of resolutions. The successful application of the filtered ZB20 parameterizations to parameterize mesoscale eddies in two idealized configurations offers a promising opportunity to reduce long‐standing biases in global ocean simulations in future studies.

How Landmass Distribution Influences the Atmospheric Dynamics of Tidally Locked Terrestrial Exoplanets

Interpretation of the ongoing efforts to simulate the atmospheres of potentially habitable terrestrial exoplanets requires that we understand the underlying dynamics and chemistry of such objects to a much greater degree than 1D or even simple 3D models enable. Here, for the tidally locked habitable-zone planet TRAPPIST-1e, we explore one effect which can shape the dynamics and chemistry of terrestrial planets: the inclusion of an Earth-like land–ocean distribution with orography. To do this we use the Earth-system model WACCM6/CESM2 to run a pair of TRAPPIST-1e models with N2–O2 atmospheres and with the substellar point fixed over either land or ocean. The presence of orography shapes atmospheric transport, and in the case of Earth-like orography, breaks the symmetry between the Northern and Southern Hemispheres which was previously found in slab ocean models. For example, peak zonal jet speeds in the Southern Hemisphere are 50%–100% faster than similar jets in the Northern Hemisphere. This also affects the meridional circulation, transporting equatorial material toward the south pole. As a result we also find significant changes in the atmospheric chemistry, including the accumulation of potentially lethal quantities of ozone at both the south pole and the surface. Future studies which investigate the effects of landmass distribution on the dynamics of exoplanetary atmospheres should pay close attention to both the dayside land fraction as well as the orography of the land. Simply modeling a flat landmass will not give a complete picture of its dynamical impact.

Impact of the Fujiwhara effect from tropical cyclones Seroja and Odette on ocean dynamic in Southern Indonesia: Insights from argo data and model analysis

Tropical cyclone (TC) disturbances can cause changes in ocean dynamics before, during, and afterwards. We employed an Argo data and ocean model product of Copernicus Marine Environment Monitoring Service (CMEMS) to address the gap in observations of the ocean response from the severe TC Seroja and its interactions with the tropical low TC Odette, known as the rare Fujiwhara effect phenomenon. Our study in southern Indonesia seawater was conducted between March and April 2021, beginning with reconciling the deviation of the CMEMS ocean model product by assessing the model accuracy to the Argo float data to build confidence in using the model for further analysis. Comparison revealed that the temperature model datasets have a higher correlation coefficient (∼0.95) with the Argo observation than the salinity model datasets. Further, we investigated the temporal evolution of vertical profiles and upper ocean responses both inside and beyond the TCs route zone. We discovered that the upper ocean response results in a colder and saltier surface with a recovery time of more than 15 days. Furthermore, as a result of the strong surface cooling, mixed layers inside the Seroja track were rapidly thickening following the passage of TC. Our findings provide insights into how the TCs Seroja and Odette affect the ocean layers, especially the primary ocean features and behaviours (temperature and salinity), using the Argo and model datasets.

Distinct Impacts of Topographic versus Planetary PV Gradients on Baroclinic Turbulence

Abstract The impacts of topographic versus planetary potential vorticity (PV) gradients on fully developed geostrophic turbulence are often treated as dynamically equivalent in theoretical understanding and parameterizations of ocean mesoscale turbulence. Using a suite of homogeneous, two-layer quasigeostrophic (QG) turbulence simulations, we identify similarities and distinctions of topographic versus planetary PV gradients in modulating turbulent eddy energy and fluxes. We show that, while an elevated background PV gradient can suppress geostrophic turbulence, a positive (negative) bottom slope with respect to the orientation of isopycnal slope barotropizes (debarotropizes) the turbulent energy at large scales, which contrasts with a dynamically inert planetary PV gradient in the modewise energy transfer. Then, in the presence of weak bottom drag, a positive slope energizes large-scale along-slope jets and limits small-scale barotropic eddies, both of which yield stronger eddy suppression effects than from a planetary PV gradient; by contrast, a negative slope hinders along-slope jet formation by enhancing the dual-energy cascade cycling, which alleviates the topographic suppression on eddies. In the presence of strong bottom drag, a positive slope elevates barotropic eddy energy, which further enhances the eddy-driven fluxes; by contrast, a negative slope confines turbulent energy to the baroclinic mode, which is strongly damped, causing further weakened turbulent energy and eddy fluxes. A flow regime captured by linear QG theories also emerges as the turbulent energy cascade is jointly suppressed by negative slopes and strong bottom drag. This study provides insights into parameterizing mesoscale eddy effects across sloping bathymetry in predictive ocean models.