Geostrophic Wind Research Articles

Weak geostrophic wind driven ventilation in street canyons with trees and green walls: Cooperating or opposing dispersions of airborne pollutants?

Urban high-rise buildings and dense street canyons significantly disrupt synchronised wind flows, potentially exacerbating urban heat island effect. Fortunately, urban heat island effect can be mitigated by the green infrastructure in cities, particularly through the presence of roadside trees and green walls. However, the combined mechanisms of tree resistance, shading effects, and cooling effect, along with sensible thermal buoyancy flows within canyon ventilation and pollution dispersion, remain undisclosed. This study employs refined computational fluid dynamics numerical simulations to analyse airflow and pollutant dispersion within urban street canyons. Air exchange rate and pollutant retention time are employed to evaluate ventilation and pollutant dispersion, respectively, inside the street canyons.In scenarios with leeward green wall layout, higher tree canopy spread and trunk height lead to a shift in high pollutant concentrations from the region near the windward side to that near the leeward side of the typical canyon (H/W = 1). Conversely, in deep canyons (H/W = 5), most pollutants are retained near the bottom leeward side. Moreover, differences in pollutant dispersion between typical and deep canyons increase with with rising tree canopy spread. Additionally, higher tree canopy spread and trunk height values result in longer pollutant retention time, which are unfavourable for pollutant removal, and particularly limiting pollutant dispersion as street canyon depth increases. Similarly, increasing tree canopy spread reduces the difference in air exchange rate between windward and leeward green wall layout patterns. Regarding Richardson number variations —indicating thermal buoyancy effects—lower Richardson number values lead to reduced differences in pollutant retention time between typical and deep canyons. Furthermore, ventilation performance of the windward green pattern surpasses that of the leeward pattern. This research provides valuable insights for implementing green infrastructure to locally mitigate urban air pollution.

Impact of surface current and temperature feedback on kinetic energy over the North-East Atlantic from a coupled ocean / atmospheric boundary layer model

A one-dimensional Atmospheric Boundary Layer (ABL1D) model is coupled with the NEMO ocean model and implemented over the Iberian–Biscay–Ireland (IBI) area at 1/36° resolution to investigate the damping effect of the current and the thermal feedback on the kinetic energy (KE) at the mesoscale. This type of coupling between an ocean model and an ABL1D is a newly proposed approach as an alternative of intermediate complexity between bulk forcing and full coupling with an atmosphere model. In ABL1D, the prognostic tracers are nudged toward large-scale variables and the wind is guided by a low-frequency geostrophic wind provided from the ERA-Interim reanalyses. First, the ABL1D is successfully validated against satellite observations regarding the wind, and the dynamic coupling coefficient (linking the near surface wind and wind-stress to the of the surface currents) are consistent with the literature, over the period 2016–2017. Our results show that the thermal feedback has a negligible impact on kinetic energy (KE) and does not influence the strength of the current feedback in the region. Given the ABL1D physics, this further indicates that the changes in the vertical wind structure caused by CFB are primarily governed by local mechanical mechanisms associated with surface wind-stress condition, rather than by thermodynamic or non-local processes within the planetary boundary layer. The induced KE reduction by the current feedback amounts to 14% at the surface and propagates down to 2000 m, indicating that it can modify the vertical distribution of KE throughout the water column. KE reductions in the surface boundary layer (0 – 300 m) and in the interior (300 – 2000 m) are attributed to a reduction of the surface wind work by 4%, and of the pressure work by 7%, respectively. The Ekman pumping anomalies induced by the current feedback tend to attenuate eddy activity and horizontal pressure gradients at depth, illustrating the potential of the current feedback to induce a geostrophic adjustment on the water column. These results illustrate the relevance of the proposed ABL1D coupling approach for reproducing the wind-current coupling (a.k.a. current feedback effect) which cannot be taken into account straightforwardly with simple bulk forcing.

The dynamics of the transpolar drift current

The governing equations for a viscous, incompressible fluid, in the thin-shell approximation, written in a coordinate frame which ensures validity at the North Pole, is simplified by invoking the f-plane approximation. Together with suitable stress conditions at the surface, describing the interaction of wind over ice and ice over water, a solution is developed for the Transpolar Drift current. This involves prescribing the near-surface geostrophic ocean flow and also, for ease of calculation, a simple version of the Ekman flow. Then, via the stress conditions at the surface, the near-surface wind consistent with these flows is obtained. An example, which exhibits the reduction in ice thickness as the flow moves over the North Pole towards the Fram Strait, and also describes the relative directions and speeds of the wind, ice and water, is presented in graphical form; this solution recovers what is observed. The considerable freedom available in this approach, allowing choices to be made for the geostrophic flow and the wind (both variable), and the ice profile at some initial position, is explained, opening the way to further, extensive investigations.

Improving seasonal predictions of German Bight storm activity

Abstract. Extratropical storms are one of the major coastal hazards along the coastline of the German Bight, the southeastern part of the North Sea, and a major driver of coastal protection efforts. However, the predictability of these regional extreme events on a seasonal scale is still limited. We therefore improve the seasonal prediction skill of the Max Planck Institute Earth System Model (MPI-ESM) large-ensemble decadal hindcast system for German Bight storm activity (GBSA) in winter. We define GBSA as the 95th percentiles of three-hourly geostrophic wind speeds in winter, which we derive from mean sea-level pressure (MSLP) data. The hindcast system consists of an ensemble of 64 members, which are initialized annually in November and cover the winters of 1960/61–2017/18. We consider both deterministic and probabilistic predictions of GBSA, for both of which the full ensemble produces poor predictions in the first winter. To improve the skill, we observe the state of two physical predictors of GBSA, namely 70 hPa temperature anomalies in September, as well as 500 hPa geopotential height anomalies in November, in areas where these two predictors are correlated with winter GBSA. We translate the state of these predictors into a first guess of GBSA and remove ensemble members with a GBSA prediction too far away from this first guess. The resulting subselected ensemble exhibits a significantly improved skill in both deterministic and probabilistic predictions of winter GBSA. We also show how this skill increase is associated with better predictability of large-scale atmospheric patterns.

Using Shelf‐Edge Transport Composition and Sensitivity Experiments to Understand Processes Driving Sea Level on the Northwest European Shelf

AbstractVariability in ocean currents, temperature and salinity drive dynamic sea level (DSL) variability on the Northwest European Shelf (NWES). It is dominated by mass variations, with steric signals relatively small. A mechanistic explanation of how ocean dynamics relates to the mass component of NWES sea level variability is required. We use regional ocean model experiments to isolate sources of variability and then investigate the effect on monthly to‐interannual DSL variability together with the simulated momentum budgets along the shelfbreak. Regional (local) wind and non‐regional (remote) forcing are important on the NWES. For the local wind forcing, the net mass flux onto the shelf, which drives a shelf‐mean mode of DSL variability, results from a combination of surface Ekman, bottom Ekman and geostrophic flows and explains 73% of the variance in transport across the shelf‐edge. The geostrophic flow is closely related to wind stress with a flow about half that of surface Ekman transport but in the opposite direction. For the remotely forced mass‐flux across the shelf‐edge, the geostrophic component explains 62% of the variance and bottom friction plays an important indirect role. The remotely forced variability, while relatively spatially uniform, is more important for explaining DSL variance over the western NWES. This mode of variability is sensitive to signals propagating northward via a thin strip of the southern boundary near the Portuguese coast, consistent with coastal trapped wave signals. It also appears to drive steric height in the Bay of Biscay, which is related to DSL on the shelf.

Stepwise Subduction Observed at a Front in the Marginal Ice Zone in Fram Strait

AbstractAt high latitudes, submesoscale dynamics act on scales of (100 m–1 km) and are associated with the breakdown of geostrophic balance, vertical velocities, and energy cascading to small scales. Submesoscale features such as fronts, filaments, and eddies are ubiquitous in marginal ice zones forced by the large horizontal density gradients. In July 2020, we identified multiple fronts and filaments using a towed undulating vehicle near the sea ice edge in central Fram Strait, the oceanic gateway to the Arctic Ocean between Greenland and Svalbard. Sea ice covered the entire study region 1–2 weeks earlier, and a stratified meltwater layer was present. We observed a front between warm and saline Atlantic Water (AW) and cold and fresh Polar Water (PW) at 30–85 m depth, where we identified a subsurface maximum in chlorophyll fluorescence and other biogeochemical properties extending along the tilted isopycnals down to 75 m, indicating subduction of AW (mixed with meltwater) that had previously occurred. The meltwater layer also featured multiple shallow fronts, one of which exhibited high velocities and a subsurface maximum in chlorophyll fluorescence, possibly indicating subduction of PW below the meltwater layer. The fronts at different depth levels suggest a stepwise subduction process near the ice edge, where water subducts from the surface below the meltwater and then further down along subsurface fronts. The submesoscale features were part of a larger‐scale mesoscale pattern in the marginal ice zone. As sea ice continuously retreats, such features may become more common in the Arctic Ocean.

Late summer northwestward Amazon plume pathway under the action of the North Brazil Current rings

The North Brazil Current (NBC) flows offshore of the mouth of the Amazon River and seasonally sheds anticyclonic rings (NBC rings) that propagate northwestward and interact with the Amazon River plume (ARP). Mesoscale features have a high temporal variability that is hard to monitor from current weekly and monthly sea surface salinity (SSS) satellite fields. Novel SSS fields with a higher temporal resolution analyzed together with satellite geostrophic currents, chlorophyll-a, and wind speed and in-situ data from the “Microbiomes cruise” on the SV Tara in August–September 2021 revealed a late summer freshwater pathway, which was not well documented in earlier studies. By combining these datasets, we improved the characterization of summer ARP pathways. In 2021, the ARP was a succession of freshwater patches cut off from the main plume by the NBC rings. A patch of about 200.000 km2 with salinity <33.5 pss was observed in September 2021, bringing 0.5 Sv of Amazon water northwestward in a period where the mean ocean currents lead to eastward transport. This patch was shallow, very stratified, and it created a surface steric-height anomaly that was identified as an anticyclonic feature in altimetric sea level products. Once separated from the NBC retroflection, it was mainly driven by Ekman currents. Other similar patches were observed during the 2021 summer, leading to a strong intermittency of the ARP transport. They strongly contributed to make 2021 the year with the largest northwestward freshwater transport in late summer within the 2010–2021 time-period investigated. This freshwater transport pathway is important for all plume-related phenomena, and show the ability of combined SMOS and SMAP data to accurately represent the day-to-day SSS variability.

The high-frequency and rare events barriers to neural closures of atmospheric dynamics

Recent years have seen a surge in interest for leveraging neural networks to parameterize small-scale or fast processes in climate and turbulence models. In this short paper, we point out two fundamental issues in this endeavor. The first concerns the difficulties neural networks may experience in capturing rare events due to limitations in how data is sampled. The second arises from the inherent multiscale nature of these systems. They combine high-frequency components (like inertia-gravity waves) with slower, evolving processes (geostrophic motion). This multiscale nature creates a significant hurdle for neural network closures. To illustrate these challenges, we focus on the atmospheric 1980 Lorenz model, a simplified version of the Primitive Equations that drive climate models. This model serves as a compelling example because it captures the essence of these difficulties.

Variability and trends of near-surface wind speed over the Tibetan Plateau: The role played by the westerly and Asian monsoon

Near-surface wind speed exerts profound impacts on many environmental issues, while the long-term (≥60 years) trend and multidecadal variability in the wind speed and its underlying causes in global high-elevation and mountainous areas (e.g., Tibetan Plateau) remain largely unknown. Here, by examining homogenized wind speed data from 104 meteorological stations over the Tibetan Plateau for 1961–2020 and ERA5 reanalysis datasets, we investigated the variability and long-term trend in the near-surface wind speed and revealed the role played by the westerly and Asian monsoon. The results show that the homogenized annual wind speed displays a decreasing trend (−0.091 m s−1 per decade, p < 0.05), with the strongest in spring (−0.131 m s−1 per decade, p < 0.05), and the weakest in autumn (−0.071 m s−1 per decade, p < 0.05). There is a distinct multidecadal variability of wind speed, which manifested in an prominent increase in 1961–1970, a sustained decrease in 1970–2002, and a consistent increase in 2002–2020. The observed decadal variations are likely linked to large-scale atmospheric circulation, and the correlation analysis unveiled a more important role of westerly and East Asian winter monsoon in modulating near-surface wind changes over the Tibetan Plateau. The potential physical processes associated with westerly and Asian monsoon changes are in concordance with wind speed change, in terms of overall weakened horizontal air flow (i.e., geostrophic wind speed), declined vertical thermal and dynamic momentum transfer (i.e., atmospheric stratification thermal instability and vertical wind shear), and varied Tibetan Plateau vortices. This indicates that to varying degrees these processes may have contributed to the changes in near-surface wind speed over the Tibetan Plateau. This study has implications for wind power production and soil wind erosion prevention in the Tibetan Plateau.

A Deep Learning Approach to Extract Balanced Motions From Sea Surface Height Snapshot

AbstractExtracting balanced geostrophic motions (BM) from sea surface height (SSH) observations obtained by wide‐swath altimetry holds great significance in enhancing our understanding of oceanic dynamic processes at submesoscale wavelength. However, SSH observations derived from wide‐swath altimetry are characterized by high spatial resolution while relatively low temporal resolution, thereby posing challenges to extract the BM from a single SSH snapshot. To address this issue, this paper proposes a deep learning model called the BM‐UBM Network, which takes an instantaneous SSH snapshot as input and outputs the projection corresponding to the BM. Training experiments are conducted both in the Gulf Stream and South China Sea, and three metrics are considered to diagnose model's outputs. The favorable results highlight the potential capability of the BM‐UBM Network to process SSH measurements obtained by wide‐swath altimetry.

The seasonal characteristics of English Channel storminess have changed since the 19th Century

Information from a variety of sources has suggested that increased storminess was experienced across the British Isles in the late eighteenth/early nineteenth century. However, it is not clear how stormy that period was relative to current conditions. Using newly recovered barometric pressure data that extend back to 1748 we have constructed a measure of geostrophic wind speed for the English Channel region using a pressure-triangle approach. We show that the 1790−1820s was a period of increased storminess across the region. This storminess extended throughout the year, which is different to comparable increases observed since the 1990s, which were confined to the winter season. While a strengthened North Atlantic jet stream is implicated in both periods, in the earlier period it is likely that the storm track shifted slightly to a more southerly location. We discuss the potential forcing mechanisms responsible for the changes in storminess over this multi-century timeframe.

A new perspective on Saturn’s polar vortices: Coupling the thermosphere and stratosphere

We argue that the observed temperature and altitude structure at the base of Saturn’s polar thermosphere implies the presence of a deep, warm-cored, anticyclonic vortex that extends deep into the mesosphere, and possibly into the stratosphere. This would be a unique structure in the solar system: a permanent, large-scale middle atmospheric feature that is a direct consequence of magnetosphere-atmosphere coupling. We investigate the plausibility of the suggestion that this deep vortex is linked to the observed non-seasonal polar hotspots in Saturn’s stratosphere. We develop a simple model of the dynamics and thermal structure of the middle atmosphere polar vortex, assuming vertical viscous transport is dominant and that the zonal winds are close to geostrophic balance. The model requires the existence of a small poleward ageostrophic wind which causes compressional heating, producing the required pressure gradients. This scenario is supported by a re-analysis of previously published thermospheric modelling results. Using this model we estimate the extent to which the thermospheric vortex penetrates into the deeper atmosphere. The results are very sensitive to the distribution of eddy viscosity. The models which best match the thermospheric temperatures produce a polar hotspot in the stratosphere that is hotter and broader than observed, indicating that our model, while a useful starting point, is incomplete. We also point out a puzzle concerning the half-widths of the polar hotspot and the polar altitude bulge at the base of the thermosphere. We are able to explain the greater width of the hotspot provided that there is enhanced downwards vertical transport of angular momentum closer to the pole, a situation that is expected given the strong convergence and downwelling predicted in this region. Finally, we speculate that the tropospheric polar cyclone could be connected to the anticyclonic vortex in the thermosphere.

Spatiotemporal variability of the coastal circulation in the northern Gulf of Cadiz from Copernicus Sentinel-3A satellite radar altimetry measurements

This study presents a generalised characterisation of the surface circulation over the northern shelf of the Gulf of Cadiz, based on 4 years of high-resolution satellite altimetry data from Sentinel-3A and wind model data. The altimetry-based surface zonal currents, adjusted for bottom-drag and wind effects, are compared with a generic CMEMS product and validated against in-situ ADCP measurements. The proposed altimetry product demonstrates superior performance than the CMEMS product, accurately reflecting surface circulation direction compared to in-situ measurements (r = 0.77, RMSE = 0.10 m/s, bias = 0.01 m/s). The use of the bottom-drag and wind-corrected/uncorrected altimetry product for spatiotemporal analysis of the shelf circulation revealed the distinct contributions of wind-driven and geostrophic components in different basin sectors. The results show that over the western basin, positive (eastward) surface currents were predominantly driven by westerly winds, while only occasionally, westward flows coincided with easterly winds, suggesting a higher control of the geostrophic component over the westward flows. In contrast, over the eastern basin, both eastward and westward flows were found to be primarily driven by favourable winds. Additionally, the analysis of Absolute Dynamic Topography (ADT) values along the whole basin showed the presence of ADT gradients both along-shore and cross-shore over the shelf, contributing to geostrophic flows. Finally, the seasonal analysis showed that eastward circulation tends to dominate during the spring and summer months, related to the upwelling season in the Gulf of Cadiz and associated westerly winds. Westward flows prevail during the winter months, related to easterly winds and the rebalancing of the along-shore sea level gradient during relaxed upwelling conditions. The findings demonstrate a significant improvement in the use of satellite altimetry data to study complex oceanographic dynamics in coastal areas, where both spatial and temporal variability are high. Moreover, the similarity of our results to those obtained from in-situ systems supports the use of altimetry data and publicly available wind models to support oceanographic studies in remote or resource-limited areas.

The North Equatorial Current and rapid intensification of super typhoons

Super Typhoon Mangkhut, which traversed the North Equatorial Current (NEC; 8–17 °N) in the western North Pacific in 2018, was the most intense Category-5 tropical cyclone (TC) with the longest duration in history—3.5 days. Here we show that the combination of two factors—high ocean heat content (OHC) and increased stratification — makes the NEC region the most favored area for a rapid intensification (RI) of super typhoons, instead of the Eddy Rich Zone (17–25 °N), which was considered the most relevant for RI occurrence. The high OHC results from a northward deepening thermocline in geostrophic balance with the westward-flowing NEC. The stratification is derived from precipitation associated with the Inter-Tropical Convergence Zone in the summer peak typhoon season. These factors, which are increasingly significant over the past four decades, impede the TC-induced sea surface cooling, thus enhancing RI of TCs and simultaneously maintaining super typhoons over the NEC region.

Spline Model: A Hydrostatic/Non-Hydrostatic Dynamic Core with Space-Time Second-Order Precision and Its Exact Tests

We present a new explicit quasi-Lagrangian integration scheme with the three-dimensional cubic spline function transform (transform = fitting + interpolation, referred to as the “spline format”) on a spherical quasi-uniform longitude–latitude grid. It is a consistent longitude–latitude grid, and to verify the feasibility, accuracy, convergence, and stability of the spline format interpolation scheme for the upstream point on the longitude–latitude grid, which may map a quasi-uniform longitude–latitude grid, a set of ideal, exact test schemes is adopted, which are recognized and proven to be effective internationally. The equilibrium flow test, cross-polar flow test, and Rossby–Haurwitz wave test are used to illustrate the spline scheme uniformity to the linear scheme and to overcome the over-dense grid in the polar region and the non-singularity of the poles. The cross-polar flow test demonstrates that the geostrophic wind crosses the polar area correctly, including the South Pole and North Pole. A non-hydrostatic, fully compressible dynamic core is used to complete the density flow test, demonstrating the existence of a time-varying reference atmosphere and that the spline format can simulate highly nonlinear fine-scale transient flows. It can be compared for the two results of the density flow test between the solution with the spline format and the benchmark reference solution with the linear format. Based on the findings, the non-hydrostatic dynamic core with the spline format is recommended for adoption. When simulated for the flow over an ideal mountain, through the “topographic gravity wave test”, the bicubic surface terrain and terrain-following height coordinates, time-split integration, and vector discrete decomposition can be derived successfully. These may serve as the foundations for a global, unified spline-format numerical model in the future.

Growth and Breakdown of Kelvin–Helmholtz Billows in the Stable Atmospheric Boundary Layer

The development and breakdown of Kelvin–Helmholtz (KH) waves (billows) in the stable atmospheric boundary layer (SABL) and their impact on vertical transport of momentum and scalars have been examined utilizing large eddy simulations. These simulations are initialized with a vertically uniform geostrophic wind and a constant potential temperature lapse rate. An Ekman type of boundary layer develops, and an inflection point forms in the SABL, which triggers the KH instability (KHI). KHI develops with the kinetic energy (KE) in the KH billows growing exponentially with time. The subsequent onset of secondary shear instability along S-shaped braids leads to the turbulent breakdown of the KH billow cores and braids. The frictional ground surface tends to slow down the growth of KE near the surface, reduce the KH billow core depth, and likely suppress other types of secondary instability. KH billows induce substantial down-gradient transport of momentum and sensible heat, which can be further enhanced by the onset of secondary shear instability. Although the KHI-induced strong transport only lasts for around 10–20 min, it reduces vertical shear and stratification in the SABL, enhances surface winds, and results in a 2–3-fold increase in the SABL depth.

A Diagnostic Study on a Southward Southwest Vortex Track Forecasted by CMA-GFS: The Role of Initial Field

Abstract It is known that the southwest vortex (SWV) is an important weather system that may induce severe weather. The southward deviation of an SWV track forecasted by the Global Assimilation and Prediction System of the China Meteorological Administration (CMA-GFS) is systematically diagnosed in this study. The southward shift of the SWV is directly attributed to the deviation of the steering flow caused by the weak forecast of the upper-level trough. According to the diagnosis of potential tendency, the underestimation of the initial vorticity advection forecasted by CMA-GFS dominates the weak development of the upper-level trough. The underestimation of the vorticity advection is eventually sourced to the weak geostrophic wind caused by the weak initial meridional and zonal gradients of the midlevel height in front of the trough. The assimilation process on the initial field of the CMA-GFS acts a negative effect on forecasting this SWV track. It weakens the π field at midmodel level, resulting in the weak midlevel height gradient in front of the trough. A verified numerical experiment initialized by a more reasonable field is carried out and the southward shift of the SWV is obviously modified. This study suggests that a reasonable analysis field is crucial for the accurate forecast of the SWV track. Significance Statement The important impact of initial field deviation in key regions on the forecast in the late period is highlighted. A systematic diagnosis process for identifying and addressing forecast issues on SWV track is proposed. This research provides a comprehensive approach for diagnosing the forecast deviation associated with SWV track.

Characterization of physical properties of a coastal upwelling filament with evidence of enhanced submesoscale activity and transition from balanced to unbalanced motions in the Benguela upwelling region

Abstract. We combine high-resolution in situ data (acoustic Doppler current profiler (ADCP), Scanfish, and surface drifters) and remote sensing to investigate the physical characteristics of a major filament observed in the Benguela upwelling region. The 30–50 km wide and about 400 km long filament persisted for at least 40 d. Mixed-layer depths were less than 40 m in the filament and over 60 m outside of it. Observations of the Rossby number Ro from the various platforms provide the spatial distribution of Ro for different resolutions. Remote sensing focuses on geostrophic motions of the region related to the mesoscale eddies that drive the filament formation and thereby reveals |Ro|&lt;0.1. Ship-based measurements in the surface mixed layer reveal 0.5&lt;|Ro|&lt;1, indicating the presence of unbalanced, ageostrophic motions. Time series of Ro from triplets of surface drifters trapped within the filament confirm these relatively large Ro values and show a high variability along the filament. A scale-dependent analysis of Ro, which relies on the second-order velocity structure function, was applied to the latter drifter group and to another drifter group released in the upwelling zone. The two releases explored the area nearly distinctly and simultaneously and reveal that at small scales (&lt;15 km) Ro values are twice as large in the filament in comparison to its environment with Ro&gt;1 for scales smaller than ∼500 m. This suggests that filaments are hotspots of ageostrophic dynamics, pointing to the presence of a forward energy cascade. The different dynamics indicated by our Ro analysis are confirmed by horizontal kinetic energy wavenumber spectra, which exhibit a power law k−α with α∼5/3 for wavelengths 2π/k smaller than a transition scale of 15 km, supporting significant submesoscale energy at scales smaller than the first baroclinic Rossby radius (Ro1∼30 km). The detected transition scale is smaller than those found in regions with less mesoscale eddy energy, consistent with previous studies. We found evidence for the processes which drive the energy transfer to turbulent scales. Positive Rossby numbers 𝒪(1) associated with cyclonic motion inhibit the occurrence of positive Ertel potential vorticity (EPV) and stabilize the water column. However, where the baroclinic component of EPV dominates, submesoscale instability analysis suggests that mostly gravitational instabilities occur and that symmetric instabilities may be important at the filament edges.

Spatial and temporal variations in the micronutrient Fe across the Peruvian shelf from 1984 to 2017

High dissolved iron (dFe) concentrations of the order of 10–100 nmol L−1 are a feature of waters influenced by sedimentary inputs in oxygen minimum zones (OMZ). However, the temporal development of dFe concentrations is poorly defined due to a general reliance on snapshot cross-shelf sections to study marine trace metal dynamics. Multiple cruise campaigns since the 1980s have investigated Fe dynamics over the Peruvian shelf, particularly between 9° S and 17°S where the shelf is broad, extremely productive and known to feature benthic dFe effluxes which are amongst the highest measured globally. This extensive long-term dataset uniquely allows us to study the interannual variability in dFe concentrations and their response to El Niño–Southern Oscillation (ENSO) events. By combining data from 11 cruises during the period 1984–2017 we are able to evaluate dFe dynamics on interannual timescales in a major OMZ. The region where average dFe concentrations are sensitive to variations in ENSO is confined to a subsurface layer at depths between 50 and 150 m, particularly in the narrow coastal region within 50 km of the coastline. Subsurface dFe concentrations were generally low during El Niño events (0.7–15.4 nmol L−1) and relatively high with a wider range of variability during the cold ENSO phase (1.1–52.1 nmol L−1). Inverse relationships between wind speed and surface/subsurface dFe were evident. In the subsurface layer, this may be attributable to enhanced dFe offshore transport along isopycnals when upwelling-favorable winds relax in accordance with previously outlined theories. Surface layer (<40 m) dFe variability was likely associated with a dilution and/or oxidation effect depending on the strength of wind driven water column mixing. Upwelling brings macronutrient-rich water into the euphotic zone, but its intensity had a limited impact on upper layer dFe concentrations possibly due to the influence of an onshore geostrophic flow. Interannual variability in surface chlorophyll-a (Chl-a) was found to correlate with dFe concentration in the offshore zone of northern Peru. This is consistent with bioassay experiments and climatological residual nitrate concentrations which both indicate proximal Fe limitation of phytoplankton growth over and beyond the northern Peruvian shelf. Overall, our work highlights the importance of physical factors driving short-term variations in Fe availability in one of the world’s most economically important fishery regions and suggests that, despite pronounced spatial and temporal variability in dFe concentrations, the ENSO phase has an impact on dFe availability.

Conservation and stability in a discontinuous Galerkin method for the vector invariant spherical shallow water equations

We develop a novel and efficient discontinuous Galerkin spectral element method (DG-SEM) for the spherical rotating shallow water equations in vector invariant form. We prove that the DG-SEM is energy stable, and discretely conserves mass, absolute vorticity, and linear geostrophic balance on general curvlinear meshes. These theoretical results are possible due to our novel entropy stable numerical DG fluxes for the shallow water equations in vector invariant form. We experimentally verify these results on a cubed sphere mesh. Additionally, we show that our method is robust, that is it can be run stably without any dissipation. The entropy stable fluxes are sufficient to control the grid scale noise generated by geostrophic turbulence without the need for artificial stabilisation.