Ocean Interaction Research Articles

Oceanic internal solitary wave interactions via the KP equation in a three-layer fluid with shear flow

Oceanic internal solitary wave interactions via the KP equation in a three-layer fluid with shear flow

Modeling seasonal-to-decadal ocean–cryosphere interactions along the Sabrina Coast, East Antarctica

Abstract. The Totten Ice Shelf (TIS) and Moscow University Ice Shelf (MUIS), along the Sabrina Coast of Wilkes Land, are the floating seaward terminuses of the second-largest freshwater reservoir in the East Antarctic Ice Sheet. Being a marine ice sheet, it is vulnerable to the surrounding ocean conditions. Recent comprehensive oceanographic observations, including bathymetric measurements off the Sabrina Coast, have shed light on the widespread intrusion of warm modified Circumpolar Deep Water (mCDW) onto the continental shelf and the intense ice–ocean interaction beneath the TIS. However, the spatiotemporal coverage of the observation is very limited. Here, we use an ocean–sea ice–ice shelf model with updated bathymetry to better understand the regional ocean circulations and ocean–cryosphere interactions. The model successfully captured the widespread intrusions of mCDW, local sea ice production and the ocean heat and volume transports into the TIS cavity, facilitating an examination of the overturning ocean circulation within the ice shelf cavities and the resultant basal melting. We found notable differences in the temporal variability in ice shelf basal melting across the two adjacent ice shelves of the TIS and the western part of the MUIS. Ocean heat transport by mCDW controls the low-frequency interannual-to-decadal variability in ice–ocean interactions, but the sea ice production in the Dalton Polynya strongly modifies the signals, explaining the regional difference between the two ice shelves. The formation of a summertime eastward-flowing undercurrent beneath the westward-flowing Antarctic Slope Current is found to play an important role in the seasonal delivery of ocean heat to the continental shelf.

Simulation of MJO with improved deep convection scheme in different resolutions of BCC-CSM2 models

AbstractThis study investigates the impacts of modifying the deep convection scheme on the ability to simulate the Madden–Julian Oscillation (MJO) in the Beijing Climate Center Climate System Model version 2 with a medium resolution (BCC-CSM2-T159) and a high resolution (BCC-CSM2-T382). On the basis of the original deep convection scheme, a modified scheme is suggested, which involves the transport processes of deep convective cloud water. The liquid cloud water that is detrained is transferred horizontally to its neighboring grids, and a portion of the cloud water that is horizontally transported is allowed to be transported downward into the lower troposphere. Both BCC-CSM2-T159 and BCC-CSM2-T382 with the modified deep convection scheme perform better than that used the original deep convection scheme in reproducing the major features of the MJO, such as its spectrum, period, intensity, eastward propagation and life cycle. Further analysis shows that those pronounced improvements in the MJO features in both BCC-CSM2-T159 and BCC-CSM2-T382 with the modified scheme are caused by transport processes of deep convective cloud water. The modified deep convection scheme enhances moisture and energy exchange from the lower troposphere to the upper troposphere around convective cloud, and promotes the convergence of moisture in the lower troposphere to the east of the MJO convection center, and then induces eastward propagation of the MJO. The comparisons between the coupled experiments and their corresponding experiments following Atmospheric Model Intercomparison Project (AMIP) simulations indicated that atmosphere–ocean interactions are also important to improve MJO simulations in the models.

Helicity: A Possible Indicator of Negative Feedback Initiation of Tropical Cyclone–Ocean Interaction

AbstractThe development of large‐scale vortex dynamo during tropical cyclogenesis through the mutual intensification of primary and secondary circulation in terms of helical evolution is well studied. However, the influence of atmospheric helicity on ocean surface heat and moisture fluxes associated with tropical cyclone (TC) evolution is yet to be understood. At its development stage, the heat and moisture flux from the ocean surface increases with the increase in intensity of the TC. This time, TC‐Ocean interaction works in a positive feedback mechanism. However, after a particular stage, the very intense wind of the TC causes strong turbulent mixing in the underlying ocean. Hence, the upper ocean layer temperature starts decreasing, leading to a decrease in heat and moisture flux exchange to the TC. We have analyzed the emergence of vortex helical structure in TC and the role of helicity in modulating surface heat and moisture fluxes. The result shows that helical development is associated with the axisymetrization of the vertical velocity and diabatic heating at the middle and upper troposphere. There is no definite TC intensity at which the initiation of negative feedback of TC‐upper ocean interaction occurs. However, the increase in helicity above 200 × 10−2 ms−2 reversed the TC‐ocean interaction trend. Therefore, TC helicity evolution and upper ocean interaction are essential to understand TC evolution.

The Effect of Salinity on Ocean Circulation and Ice–Ocean Interaction on Enceladus

Observational data suggest that the ice shell on Enceladus is thicker at the equator than at the pole, indicating an equator-to-pole ice flow. If the ice shell is in an equilibrium state, the mass transport of the ice flow must be balanced by the freezing and melting of the ice shell, which in turn is modulated by the ocean heat transport. Here we use a numerical ocean model to study the ice–ocean interaction and ocean circulation on Enceladus with different salinities. We find that salinity fundamentally determines the ocean stratification. A stratified layer forms in the low-salinity ocean, affecting the ocean circulation and heat transport. However, in the absence of tidal heating in the ice shell, the ocean heat transport is found to always be toward lower latitudes, resulting in freezing at the poles, which cannot maintain the ice shell geometry against the equator-to-pole ice flow. The simulation results suggest that either the ice shell on Enceladus is not in an equilibrium state or tidal dissipation in the ice shell is important in maintaining the ice shell geometry. The simulations also suggest that a positive feedback between cross-equatorial ocean heat transport and ice melting results in spontaneous symmetry breaking between the two hemispheres. This feedback may play a role in the observed interhemispheric asymmetry in the ice shell.

Some finite difference methods for solving linear fractional KdV equation

The time-fractional Korteweg de Vries equation can be viewed as a generalization of the classical KdV equation. The KdV equations can be applied in modeling tsunami propagation, coastal wave dynamics, and oceanic wave interactions. In this study, we construct two standard finite difference methods using finite difference methods with conformable and Caputo approximations to solve a time-fractional Korteweg-de Vries (KdV) equation. These two methods are named as FDMCA and FDMCO. FDMCA utilizes Caputo's derivative and a finite-forward difference approach for discretization, while FDMCO employs conformable discretization. To study the stability, we use the Von Neumann Stability Analysis for some fractional parameter values. We perform error analysis using L1 & L∞ norms and relative errors, and we present results through graphical representations and tables. Our obtained results demonstrate strong agreement between numerical and exact solutions when the fractional operator is close to 1.0 for both methods. Generally, this study enhances our comprehension of the capabilities and constraints of FDMCO and FDMCA when used to solve such types of partial differential equations laying some ground for further research.

Investigating Atmospheric Responses to and Mechanisms Governing North Atlantic Sea Surface Temperatures over 10-Year Periods

Abstract North Atlantic sea surface temperature (SST) variability plays a critical role in modulating the climate system. However, characterizing patterns of North Atlantic SST variability and diagnosing the associated mechanisms is challenging because they involve coupled atmosphere–ocean interactions with complex spatiotemporal relationships. Here we address these challenges by applying a time-evolving self-organizing map approach to a long preindustrial coupled control simulation and identify a variety of 10-yr spatiotemporal evolutions of winter SST anomalies, including but not limited to those associated with the North Atlantic Oscillation–Atlantic multidecadal variability (NAO–AMV)-like interactions. To assess mechanisms and atmospheric responses associated with various SST spatiotemporal evolutions, composites of atmospheric and oceanic variables associated with these evolutions are investigated. Results show that transient-eddy activities and atmospheric circulation responses exist in almost all the evolutions that are closely correlated to the details of the SST pattern. In terms of the mechanisms responsible for generating various SST evolutions, composites of ocean heat budget terms demonstrate that contributions to upper-ocean temperature tendency from resolved ocean advection and surface heat fluxes rarely oppose each other over 10-yr periods in the subpolar North Atlantic. We further explore the potential for predictability for some of these 10-yr SST evolutions that start with similar states but end with different states. However, we find that these are associated with abrupt changes in atmospheric variability and are unlikely to be predictable. In summary, this study broadly investigates the atmospheric responses to and the mechanisms governing the North Atlantic SST evolutions over 10-yr periods. Significance Statement Climate variability in the North Atlantic Ocean has wide-ranging impacts on global and regional climate. However, the processes involved include interactions between the ocean and atmosphere that vary across both space and time, making it challenging to characterize and predict. Using a novel machine learning approach, this study identifies various time evolutions of North Atlantic sea surface temperature patterns over 10-yr periods. This includes evolutions with similar start states but different trajectories, which have important implications for predictability. Furthermore, we investigate the mechanisms responsible for these evolutions and how different sea surface temperature patterns affect atmospheric circulation through small-scale atmospheric disturbances. These new insights into the complex ocean–atmosphere interactions over time are critical for improving decadal prediction skill.

Dissolved Inorganic Nutrient Biogeochemistry in an Urbanized Coastal Region: A Study of Dapeng Cove, Shenzhen

Dissolved inorganic nutrients are pivotal in maintaining the material and energy balance of marine ecosystems, impacting the survival and dynamic succession of marine organisms. To gain a deeper understanding of the source and sink characteristics of dissolved inorganic nutrients in bays affected by human activities and to elucidate the processes involving filter-feeding shellfish in relation to these nutrients, this study investigated the source and sink dynamics of dissolved inorganic nutrients in the Dapeng Cove sea area of Shenzhen. Over the past decade, a significant change in the N/P ratio within the survey area has been observed, suggesting a shift in nutrient limitation from nitrogen to phosphorus or phosphorus–silicon limitation. This induced change in the N/P ratio, along with Si/N and Si/P ratios, may facilitate the growth of cyanobacteria and, subsequently, alter the proportions of diatoms, dinoflagellates, and cyanobacteria. Seasonal fluctuations in human disturbance intensity and precipitation determine the seasonal and spatial distribution of nutrients in the bay, thereby influencing the bay ecosystem metabolism. The Land–Ocean Interactions in the Coastal Zone (LOICZ) model analysis revealed that the bay represents a major source of inorganic nitrogen and a source of phosphate in spring, summer, and autumn, while acting as a sink for phosphate in winter. Furthermore, rivers and groundwater represent the primary sources of phosphate and inorganic nitrogen in the bay. The bay exhibits an annual net ecosystem metabolic rate of 7.06 mmol C/m2/d, with denitrification dominating the nitrogen cycle at 12.65 mmol C/m2/d. Overall, the Dapeng Cove ecosystem displays net production exceeding respiration, classifying it as an autotrophic system. Additionally, the nitrogen cycle in the sea area is predominantly driven by denitrification. The analysis also revealed that the impact of oyster proliferation on the physical and chemical factors in the surveyed area is relatively weaker than that of surface runoff and groundwater inputs.

Linking Arctic stratospheric polar vortex weakening to rising CO2-induced intensification of the Indo-Pacific warm pool during the past five decades

Accompanying the global rise in greenhouse gas emissions, a warming trend in the Indo-Pacific warm pool (IPWP) has exerted a discernible influence on tropical atmosphere–ocean interactions. However, the impact of this intensification of the IPWP on the Arctic stratospheric polar vortex (ASPV) remains unclear. In this study, we revealed a link between the changes in the IPWP and ASPV during the early winter months, with nearly half of the weakening in the ASPV attributable to the intensification of the IPWP from 1968 to 2020. Wave trains triggered by the elevated SST in the IPWP region lead to enhanced eastward-propagating flux convergence in the northern high-latitude stratosphere, ultimately resulting in a diminishing ASPV. With increasing atmospheric greenhouse gas, the ASPV is poised to further weaken in the future, particularly in the context of a more intense IPWP. Our finding has significant implications for early winter ASPV strength and location prediction and seasonal weather forecasting.

Hydroclimate variability in the central Mediterranean during MIS 17 interglacial (Middle Pleistocene) highlights timing offset with monsoon activity

Mediterranean climates are characterized by warm, dry summers and mild, rainy winters. Previous studies suggest that over the last 1.36 Myr, Mediterranean winter rainfalls were in phase with the African monsoon. Here we present a high-resolution terrestrial and marine dataset for the Marine Isotope Stage 17 interglacial (Middle Pleistocene) from Southern Italy, showing that precipitation rates and regimes in the central Mediterranean varied independently of the monsoon system. Specifically, events of extreme summer precipitation were promoted by increased regional insolation rates and/or extratropical cyclones, and their magnitude was further enhanced by the advection of cool and humid North Atlantic air during stadials. Our findings provide new information on the short- to mid-term natural hydroclimatic variability of the Mediterranean basin, and offer new critical insights on land–ocean interactions at the regional scale by complementing previous analyses on the displacement of storm tracks toward southern Europe.

Automatic Extraction of the Calving Front of Pine Island Glacier Based on Neural Network

Calving front location plays a crucial role in studying ice–ocean interaction, mapping glacier area change, and constraining ice dynamic models. However, relying solely on visual interpretation to extract annual changes in the calving front of ice shelves is a time-consuming process. In this study, a comparative analysis was conducted on the segmentation obtained from fully convolutional networks (FCN), U-Net, and U2-Net models, revealing that U2-Net exhibited the most effective classification. Notably, U2-Net outperformed the other two models by more than 30 percent in terms of the F1 parameter. Therefore, this paper introduces an automated approach that utilizes the U2-Net model to extract the calving front of ice shelves based on a Landsat image, achieving an extraction accuracy of 58 m. To assess the model’s performance on additional ice shelves in the polar region, the calving front of the Totten and Filchner ice shelves were also extracted for the past decade. The findings demonstrated that the ice velocity of the Filchner ice shelf exceeded that of the Totten ice shelf. Between February 2014 and March 2015, the majority of the calving fronts along the Filchner Ice Shelf showed an advancing trend, with the fastest-moving front measuring 3532 ± 58 m/yr.

Unleashing the power of the Sun: the increasing impact of the solar cycle on off-season super typhoons since the 1990s

The occurrence of super typhoons outside the normal typhoon season can result in devastating loss of life and property damage. Our research reveals that the 11-year solar cycle can affect the incidence of these off-season typhoons (from November to April) in the western North Pacific by influencing sea surface temperature (SST) through a footprint mechanism. The solar cycle, once amplified by atmospheric and ocean interactions, generates a noticeable SST footprint in the subtropical North Pacific during winter and spring, which eventually intrudes into the tropical central Pacific and affects the atmospheric conditions, resulting in an increase or decrease in the occurrence of super typhoons during active or inactive solar periods. This mechanism has become more effective since the Atlantic Multi-decadal Oscillation (AMO) shifted to a warm phase in the 1990s, intensifying the subtropical Pacific couplings. An example of this type of off-season super typhoon during an active solar period is Typhoon Haiyan in 2013. By incorporating information about the solar cycle, we can anticipate the likelihood of super typhoon occurrences, thus improving decadal disaster preparation and planning.

Roles of the Indo-Pacific subsurface Kelvin waves and volume transport in prolonging the triple-dip 2020–2023 La Niña

The rare triple-dip 2020–2023 La Niña event has resulted in a series of extreme climate events across the globe. Here, we reveal the role of tropical Indo-Pacific oceanic interactions in driving the first triple-dip La Niña of the twenty-first century. Specifically, we found that the eastern Indian Ocean subsurface warming anomalies were associated with the re-intensification of the subsequent La Niña event. The subsurface warming anomaly signals were propagated eastward by equatorial and coastal subsurface Kelvin waves from the eastern Indian Ocean to the western Pacific Ocean through the Indo-Pacific oceanic pathway, which contributes to the accumulation of heat content and deepens the thermocline in the western tropical Pacific. The westward Indonesian Throughflow (ITF) transported more heat during multi-year La Niña events from the western Pacific Ocean to the eastern Indian Ocean than during single-year events, resulting in the injection of more warm water into the eastern Indian Ocean. The combination of subsurface Kelvin wave propagation and increased ITF volume transport in the Indo-Pacific region acted to prolong the heat content in the western Pacific during the decay phase of La Niña, ultimately leading to the rare triple-dip 2020–2023 La Niña event.

High‐Resolution Simulations of the Plume Dynamics in an Idealized 79°N Glacier Cavity Using Adaptive Vertical Coordinates

AbstractFor better projections of sea level rise, two things are needed: an improved understanding of the contributing processes and their accurate representation in climate models. A major process is basal melting of ice shelves and glacier tongues by the ocean, which reduces ice sheet stability and increases ice discharge into the ocean. We study marine melting of Greenland's largest floating ice tongue, the 79° North Glacier, using a high‐resolution, 2D‐vertical ocean model. While our fjord model is idealized, the results agree with observations of melt rate and overturning strength. Our setup is the first application of adaptive vertical coordinates to an ice cavity. Their stratification‐zooming allows a vertical resolution finer than 1 m in the entrainment layer of the meltwater plume, which is important for the plume development. We find that the plume development is dominated by entrainment only initially. In the stratified upper part of the cavity, the subglacial plume shows continuous detrainment. It reaches neutral buoyancy near 100 m depth, detaches from the ice, and transports meltwater out of the fjord. Melting almost stops there. In a sensitivity study, we show that the detachment depth depends primarily on stratification. Our results contribute to the understanding of ice–ocean interactions in glacier cavities. Furthermore, we suggest that our modeling approach with stratification‐zooming coordinates will improve the representation of these interactions in global ocean models. Finally, our idealized model topography and forcing are close to a real fjord and completely defined analytically, making the setup an interesting reference case for future model developments.

The Southern Ocean carbon and climate observations and modeling (SOCCOM) project: A review

The Southern Ocean serves as the primary gateway through which the intermediate, deep, and bottom waters of the ocean interact with the surface ocean (and thus the atmosphere), and it has a profound influence on the oceanic uptake of anthropogenic carbon and heat as well as nutrient resupply from the abyss to the surface. Yet it has been the least observed and understood region of the world ocean. The Southern Ocean Carbon and Climate Observations and Modeling (SOCCOM) project was implemented in 2014 with a goal to help remedy this deficit in observations and understanding. The SOCCOM project is based on two major advances that have the potential to transform understanding of the Southern Ocean. The first is the development of new biogeochemical sensors mounted on autonomous profiling floats that allow sampling of ocean biogeochemistry in 3-dimensional space. Floats may detect processes with a temporal resolution that ranges from hours to years. The second is that the climate modeling community finally has the computational resources and physical understanding to develop fully coupled climate models that can represent crucial, mesoscale processes in the Southern Ocean, as well as corresponding models that assimilate observations to produce a state estimate. The observational component, based on deployment of profiling floats with oxygen, nitrate, pH and bio-optical sensors, is generating vast amounts of new biogeochemical data that provide a year-round view of the Southern Ocean from the surface to 2000 m. The modeling effort is applying these observations and enhancing our understanding of the current ocean, and reducing uncertainty in projections of future carbon and nutrient cycles and climate. After nine years of operation, including a project renewal in the sixth year, the SOCCOM project has deployed more than 260 profiling floats. These floats have collected over 27,000 vertical profiles throughout the Southern Ocean. A data-assimilating biogeochemical state estimate model has been implemented. Here, the design of the SOCCOM project is reviewed and the scientific results that have been obtained are described. The project’s capability to help meet the observing system priorities outlined for a notional UN Decade for Ocean Sciences Southern Ocean observing system is assessed.

EMPOWERING OCEANIC RESEARCH: LEVERAGING GPT-4 ARCHITECTURE IN THE STUDY OF MARINE AEROSOLS

In recent years, artificial intelligence has surged to the forefront of scientific research, proving indispensable in various disciplines. This article delves into the innovative utilization of the GPT-4 architecture, a state-of-the-art AI model, in the study of marine aerosols. Marine aerosols play a pivotal role in climatic systems, influencing cloud formation, precipitation, and radiative balance. Traditional methods of studying these aerosols often require extensive time, resources, and face challenges in predicting complex oceanic interactions. The GPT-4 architecture transcends these limitations, offering rapid data processing, in-depth pattern recognition, and predictions with unprecedented accuracy. By integrating AI into oceanic aerosol research, we not only enhance our understanding of marine environments but also pave the way for broader applications of AI in atmospheric and marine sciences.

Consistency of Intra-Centennial Oscillations in Length of Day and Oceanic Characteristics

The paper presents analysis of intra-centennial (inter-decadal and multidecadal) variations of the length of day (LOD) and some oceanic parameters such as sea surface temperature (SST) and sea level (SL). Methods of multivariate regression analysis and correlation analysis are used. Results of the regression analysis show a spatially coherent response of SST to LOD variations on the multidecadal time scale. The earlier response is peculiar to the north and tropical Atlantic where the multidecadal SST variations are approximately opposite to the LOD variations. In the most remaining parts of the oceans, except especially in the Nino 3.4 region of the equatorial east Pacific, the multidecadal SST variations are generally lagged relative to the antiphase variations of the LOD. Smoothing of SST averaged over different areas and of the global mean SL shows that the intra-annual variations include inter-decadal, 20–30-year, multidecadal, 60–70-year, components that correspond to similar oscillation components in the LOD. The most striking correspondence of the two components is observed between the LOD and SST averaged over the Nino 3.4 region. Generally, there are significant correlations of the intra-centennial variations on the averaged and smoothed SST series and global mean SL with the LOD variations. We propose that angular momentum exchange processes involving oceanic circulation and interactions between the Earth’s core and the mantle play probably a part in the observed relationships of intra-centennial variations in oceanic parameters with variations in the LOD.

Contourite stratigraphic models linked to the light intermediate versus dense deep Mediterranean water flow regime variations (Alboran Sea, SW Mediterranean)

Several water masses are involved in the circulation of oceans, their bottom layers impacting on sedimentation through contourites. The majority of palaeoceanographic studies on regional contourites are performed for one water mass despite that their joint study would offer relevant clues to understand past ocean and climate interaction. This works presents for the first time a an analysis about the impact of the Light Intermediate Mediterranean (LMW) and Dense Deep Mediterranean (DMW) bottom currents on the sedimentation in the Alboran Sea (SW Mediterranean) and its paleoceanographic significance in response to climatic oscillations from the last glacial period to the Holocene. To do this, an integration of chronostratigraphical, sedimentological, and compositional data is carried out from contourites formed by those water masses. That integration enable us to define three distinct contourite stratigraphic models. (I) The contourite terrace model, characterized by coarse-grained contourites, which is an archive of the interplay between the high-energy Atlantic Water-LMW interface and glacioeustasy from the Younger Dryas (YD) to the Holocene. (II) The contourite drift models, which are archives of rapid ocean-climate coupled fluctuations since 29.5 kyr. They comprise coarse-grained contourites formed by a relatively fast LMW and fine-grained contourites formed by a relatively weak DMW, except for the Heinrich Stadials HS3 to HS1 and YD when coarse-grained contourites were deposited. (III) The contourite/turbidite mixed model represents another archive of DMW and glacioeustasy interplay from the end of the late Pleistocene to Holocene.That contourite stratigraphy allows us to infer for the first time the relative variability of the LMW versus DMW flow regimes, which records differences and similarities. The similarities indicate that the LMW and DMW fluctuations occur in parallel at millennial and centennial time scales. The differences refer to the overall higher velocity of LMW versus DMW; the magnitude changes in velocities that are lower for LMW and higher for DMW; the recognition of three short ventilation events (a, b, c) during HS1 and HS2 for only DMW; and the distinct LMW and DMW responses to the onset of glacial conditions and return to interglacial conditions during the HSs, YD and Holocene cold periods.The proposed contourite stratigraphic models can be applied for other areas in the Mediterranean margins to identify and correlate the LMW and DMW palaeoceanographic events throughout this sea. The findings suggest that the different water masses that make up the water column must be seriously considered to fully understand palaeoceanographic and palaeoclimatic studies based on contourites. This is because their distinct impact on sedimentation may provide new insights into their different palaeoceanographic responses to rapid climatic oscillations and their triggering mechanisms.

Stronger Tropical Cyclone–Induced Ocean Cooling in Near-Coastal Regions Compared to the Open Ocean

Abstract Tropical cyclones (TC) often induce strong mixing in the upper ocean that generates a trail of cooler sea surface temperature (Twake) in their wakes. The Twake can affect TC intensity, so its prediction is important, especially in coastal regions where TCs can make landfall. Coastal Twakes are often more complex than those in the open ocean due to the influences of coastline geometry, highly variable water depth, continental runoff, and shelf processes. Using observational data since 2002, here we show a significantly stronger global mean Twake in coastal regions compared to offshore regions. Temperature stratification is the main driver of stronger coastal Twakes in the North Atlantic and east Pacific. In the northwest Pacific and north Indian Ocean, the differences between coastal and offshore Twakes are smaller due to compensation between TC forcings and ocean stratification. The north Indian Ocean is unique in the Northern Hemisphere because salinity stratification plays a major role on the spatial distribution of Twake. In the South Pacific Ocean, TC intensity and translation speed are crucial for explaining coastal–offshore Twake differences, while ocean stratification and mixed layer depth are more important for the coastal–offshore Twake differences in the south Indian Ocean. These findings suggest that coastal–offshore differences in ocean stratification need to be properly represented in models in order to capture changes in TC-induced ocean cooling as storms approach landfall. Significance Statement Landfalling tropical cyclones (TCs) often cause considerable damage in coastal regions with dense human populations. Understanding TC–ocean interaction and how it differs between coastal and offshore regions can help predict TC intensity prior to landfall. Sea surface cooling after TC passage is an important proxy for TC–ocean interaction. A global evaluation of coastal TC-induced cooling has not been conducted. Using data covering two decades, we show significantly stronger TC-induced surface cooling in coastal regions compared to offshore regions at the global scale and in all basins except the northwest Pacific and north Indian Ocean. The difference is driven mainly by upper-ocean conditions in the North Atlantic, east Pacific, and south Indian Ocean, and by TC characteristics in the South Pacific.

Simulation of flexural-gravity wave propagation for elastic plates in shallow water using an energy-stable finite difference method with weakly enforced boundary and interface conditions

We introduce an energy stable, high-order-accurate finite difference approximation of the dynamic, pure bending Kirchhoff plate equations for complex geometries and spatially variable properties. We utilize the summation-by-parts (SBP) framework to discretize the biharmonic operator with variable coefficients, with attention given to free and clamped boundary conditions and corner conditions. Energy conservation is established by combining SBP boundary closures with weak enforcement of the boundary and interface conditions using a penalty (simultaneous approximation term, SAT) technique. Then we couple the plate equations to the shallow water equations to study flexural-gravity wave propagation, and prove that the semi-discrete system of equations is self-adjoint. We demonstrate the stability and accuracy properties of the method on curvilinear multiblock grids using the method of manufactured solutions. The method, which we provide in an open-source code, is then used to model ocean wave interactions with the Thwaites Glacier and Pine Island Ice Shelf in the Amundsen Sea off the coast of West Antarctica.