Search

Based on the ERA5 reanalysis datasets during 1980–2019, a total of eleven zonal shear lines (ZSLs) that caused heavy precipitation and lasted more than 60 hours over the Tibetan Plateau in summer are selected for composite analysis. By decomposing the kinetic energy (K) near the ZSL into divergent and rotational kinetic energies (KD and KR) and the kinetic energy of interaction between the divergent wind and the rotational wind (KRD), the influence of the rotational and divergent winds on the evolution of the ZSL intensity is investigated from the perspective of KD and KR. The main results are as follows. The ZSL is a comprehensive reflection of rotation and convergence. The intensity evolution of ZSL is essentially synchronized with those of K, KR, and KRD but lags behind KD by about three hours. The enhancement of K is mainly contributed by KR, which is governed by the conversion from KD to KR. Furthermore, the increase in the conversion from KD to KR is controlled by the geostrophic effect term Af, which is determined by the joint enhancement of the zonal rotational and meridional divergent wind components (uR and vD). Therefore, the joint enhancement of uR and vD controls the increase of the ZSL intensity, leading to increased precipitation.

This second part of a numerical study on shallow-cumulus dilution focuses on the sensitivity of cloud dilution to changes in the vertical wind profile. Insights are obtained through large-eddy simulations of maritime and continental cloud fields. In these simulations, the speed of the initially uniform geostrophic wind and the strength of geostrophic vertical wind shear in the cloud and subcloud layer are varied. Increases in the cloud-layer vertical wind shear (up to 9 ms-1km-1) lead to 40 %–50 % larger cloud-core dilution rates compared to their respective unsheared counterparts. When the background wind speed, on the other hand, is enhanced by up to 10 m s−1 and subcloud-layer vertical wind shear develops or is initially prescribed, the dilution rate decreases by up to 25 %. The sensitivities of the dilution rate are linked to the updraft strength and the properties of the entrained air. Increases in the wind speed or vertical wind shear result in lower vertical velocities across all sets of experiments with stronger reductions in the cloud-layer wind shear simulation (27 %–47 %). Weaker updrafts are exposed to mixing with the drier surrounding air for a longer time period, allowing more entrainment to occur (i.e., the “core-exposure effect”). However, reduced vertical velocities, in concert with increased cloud-layer turbulence, also assist in widening the humid shell surrounding the cloud cores, leading to entrainment of more humid air (i.e., the “core–shell dilution effect”). In the experiments with cloud-layer vertical wind shear, the core-exposure effect dominates and the cloud-core dilution increases with increasing shear. Conversely, when the wind speed is increased and subcloud-layer vertical wind shear develops or is imposed, the core–shell dilution effect dominates to induce a buffering effect. The sensitivities are generally stronger in the maritime simulations, where weaker sensible heat fluxes lead to narrower, more tilted, and, therefore, more suppressed cumuli when cloud-layer shear is imposed. Moreover, in the experiments with subcloud wind shear, the weaker baseline turbulence in the maritime case allows for a larger turbulence enhancement, resulting in a widening of the transition zones between the cores and their environment, leading to the entrainment of more humid air.

Land surface processes are vital to the performance of regional climate models in dynamic downscaling application. In this study, we investigate the sensitivity of the simulation by using the weather research and forecasting (WRF) model at 10-km resolution to the land surface schemes over Central Asia. The WRF model was run for 19 summers from 2000 to 2018 configured with four different land surface schemes including CLM4, Noah-MP, Pleim-Xiu and SSiB, hereafter referred as Exp-CLM4, Exp-Noah-MP, Exp-PX and Exp-SSiB respectively. The initial and boundary conditions for the WRF model simulations were provided by the National Centers for Environmental Prediction Final (NCEP-FNL) Operational Global Analysis data. The ERA-Interim reanalysis (ERAI), the GHCN-CAMS and the CRU gridded data were used to comprehensively evaluate the WRF simulations. Compared with the reanalysis and observational data, the WRF model can reasonably reproduce the spatial patterns of summer mean 2-m temperature, precipitation, and large- scale atmospheric circulation. The simulations, however, are sensitive to the option of land surface scheme. The performance of Exp-CLM4 and Exp-SSiB are better than that of Exp-Noah-MP and Exp-PX assessed by Multivariable Integrated Evaluation (MVIE) method. To comprehensively understand the dynamic and physical mechanisms for the WRF model’s sensitivity to land surface schemes, the differences in the surface energy balance between Ave-CLM4-SSiB (the ensemble average of Exp-CLM4 and Exp-SSiB) and Ave-NoanMP-PX (the ensemble average of Exp-Noah-MP and Exp-PX) are analyzed in detail. The results demonstrate that the sensible and latent heat fluxes are respectively lower by 30.42 W·m−2 and higher by 14.86 W·m−2 in Ave-CLM4-SSiB than that in Ave-NoahMP-PX. As a result, large differences in geopotential height occur over the simulation domain. The simulated wind fields are subsequently influenced by the geostrophic adjustment process, thus the simulations of 2-m temperature, surface skin temperature and precipitation are respectively lower by about 2.08 ℃, 2.23 ℃ and 18.56 mm·month−1 in Ave-CLM4-SSiB than that in Ave-NoahMP-PX over Central Asia continent.

Knowledge of ocean surface geostrophic circulation has been greatly improved due to satellite observations. However, it is essential to use different observation sources (satellites and in situ data) in order to improve accuracy and resolution. This study combines along-track Sea Level Anomalies (SLA) with geostrophic velocity estimated from surface drifters to map SLA and associated geostrophic current anomalies in the Gulf of Mexico. Firstly, substantial pre-processing is needed on the drifter data to extract the geostrophic component of the signal in order to be consistent with the physical content provided by altimetry. This step includes estimating and removing the Ekman current, Stokes drift and wind slippage. Three kinds of drifters are used: – Drifters belonging to Woods Hole Group, a CLS Group company that launches their own drifters in the Gulf of Mexico. – Drifters launched in the framework of the Lagrangian Submesoscale ExpeRiment (LASER) campaign (January-April 2016). – Drifters from Atlantic Oceanographic and Meteorological Laboratory (AOML). Secondly, drifters and along-track SLA from Jason2, HY2, SARAL and Cryosat-2 are combined through multivariate objective analysis to map a daily time series of SLA and associated geostrophic current anomalies from 01/09/2015 to 30/04/2016. Finally, comparisons with independent data reveal the improved agreement of maps combining both altimetry and drifter data, especially for the meridional component of geostrophic current.

The global superrotation index S compares the integrated axial angular momentum of the atmosphere to that of a state of solid-body corotation with the underlying planet. The index S is similar to a zonal Rossby number, which suggests it may be a useful indicator of the circulation regime occupied by a planetary atmosphere. We investigate the utility of S for characterizing regimes of atmospheric circulation by running idealized Earthlike general circulation model experiments over a wide range of rotation rates Ω, 8ΩE to ΩE/512, where ΩE is Earth’s rotation rate, in both an axisymmetric and three-dimensional configuration. We compute S for each simulated circulation, and study the dependence of S on Ω. For all rotation rates considered, S is on the same order of magnitude in the 3D and axisymmetric experiments. For high rotation rates, S ≪ 1 and S ∝ Ω−2, while at low rotation rates S ≈ 1/2 = constant. By considering the limiting behavior of theoretical models for S, we show how the value of S and its local dependence on Ω can be related to the circulation regime occupied by a planetary atmosphere. Indices of S ≪ 1 and S ∝ Ω−2 define a regime dominated by geostrophic thermal wind balance, and S ≈ 1/2 = constant defines a regime where the dynamics are characterized by conservation of angular momentum within a planetary-scale Hadley circulation. Indices of S ≫ 1 and S ∝ Ω−2 define an additional regime dominated by cyclostrophic balance and strong equatorial superrotation that is not realized in our simulations.

The interannual variability of summertime eddy-induced heat transport (EHT) in the western South China Sea (WSCS) is investigated based on the downgradient eddy diffusivity method and explored its formation mechanism. Estimations of long-term mean EHT and its monthly evolution reveal that the largest EHT in the SCS occurs in the WSCS region during the summer. In the WSCS, enhanced EHT and eddy kinetic energy (EKE) levels are simultaneously observed in 1994, 1999, 2002, 2006, 2008, 2009, 2012, 2014 whereas the lower EHT and EKE levels are observed in 1995, 2000, 2001, 2003, 2004, 2007, 2010, 2015, 2017 during JAS (July, August and September) months. Analysis of the Simple Ocean Data Assimilation, version 3.3.1 (SODAv3) data along 110.75° E reveals a strong surface intensification of the summertime eastward jet (SEJ) in the EHT-strong years than the EHT-weak years. Linear stability analysis conducted by adopting a 2 ½-layer reduced gravity model shows that the increased EHT in EHT-strong years is due to the enhanced baroclinic instability caused by the strong vertical shear developed through the surface-intensification of SEJ. The cause for the interannually varying vertical shear can be sought in the interannually varying meridional temperature gradient which is influenced by the combined forcing of the meridional Ekman flux convergence, meridional geostrophic flux convergence, and convergence of the latitudinally dependent surface heat flux forcing. It is also found that the interannual variations of EHT in the SCS are partially influenced by the local wind stress curl and remote forcing from the eastern boundary.

Multidecadal variability of the Atlantic Meridional Overturning Circulation (AMOC) plays a vital role in Earth’s climate variability. Climate change has the potential to alter the causes and characteristics of AMOC multidecadal variability. Here we use a coupled climate model to simulate AMOC multidecadal variability under three distinct atmospheric CO2 concentrations: Last Glacial Maximum, preindustrial, and 4 × preindustrial levels. Firstly, we discover that AMOC multidecadal variability exhibits a shortened period and a reduced amplitude with increasing atmospheric CO2. We find that these changes in AMOC variability are largely related to enhanced ocean stratification in the subpolar North Atlantic with increasing CO2 which in turn changes the characteristics of westward propagating oceanic baroclinic Rossby waves. Our analysis indicates that the shortened period is primarily due to the increased speed of free oceanic Rossby waves, and the reduced amplitude is mainly due to the reduced magnitude of atmospherically-forced oceanic Rossby waves. Mean flow effects, in the form of eastward mean zonal advection and westward geostrophic self-advection, need to be considered as they largely increase the speed of Rossby waves and hence allow for a better estimate of the changes in the period and amplitude of AMOC variability. Secondly, to explore the possible linkage between atmospheric variability and AMOC fluctuations under each CO2 concentration in a qualitative manner, we analyze the relationship between the North Atlantic Oscillation (NAO) and the AMOC and find a significant negative correlation between the two only under the preindustrial levels where the NAO leads the AMOC by 3–11 years.

A linear two-layer model is used to elucidate the role of prognostic moisture on quasigeostrophic (QG) motions in the presence of a mean thermal wind (). Solutions to the basic equations reveal two instabilities that can explain the growth of moist QG systems. The well-documented baroclinic instability is characterized by growth at the synoptic scale (horizontal scale of ~1000 km) and systems that grow from this instability tilt against the shear. Moisture–vortex instability—an instability that occurs when moisture and lower-tropospheric vorticity exhibit an in-phase component—exists only when moisture is prognostic. The instability is also strongest at the synoptic scale, but systems that grow from it exhibit a vertically stacked structure. When moisture is prognostic and is easterly, baroclinic instability exhibits a pronounced weakening while moisture vortex instability is amplified. The strengthening of moisture–vortex instability at the expense of baroclinic instability is due to the baroclinic () component of the lower-tropospheric flow. In westward-propagating systems, lower-tropospheric westerlies associated with an easterly advect anomalous moisture and the associated convection toward the low-level vortex. The advected convection causes the vertical structure of the wave to shift away from one that favors baroclinic instability to one that favors moisture–vortex instability. On the other hand, a westerly reinforces the phasing between moisture and vorticity necessary for baroclinic instability to occur. Based on these results, it is hypothesized that moisture–vortex instability is an important instability in humid regions of easterly such as the South Asian and West African monsoons.

We examined the characteristics of nocturnal low-level wind maximum (NLWM) at a topographically complex site in Boseong, South Korea. The study site is located on the southern coast of the Korean peninsula. Although the area within a few kilometers of the site is flat, it is surrounded by mountains and faces the ocean on the southeast side. We defined NLWM as the lowest wind speed maximum within the lowest 300 m of the atmosphere that is at least 2 m s−1 greater than the next minimum above during nighttime. Climatology of NLWMs was derived from a 4-year half-hourly database of wind profiles (11 levels) obtained from a 300-m tower. The half-hourly data showed that NLWM occurred on approximately 18% of the nights. NLWMs were typically situated at 40 m above the ground. Half of the NLWMs had a speed of 3–5 m s−1, and its direction was predominantly west-northwesterly regardless of the height or season. The occurrence of NLWM showed a maximum in October and a minimum in July. The cause of NLWM at this site was discussed. It was found that the NLWM at this site occurred in the presence of valley wind forcing on nights with weak northerly geostrophic winds at 850 hPa.

In a recent study (Želi et al. in Bound Layer Meteorol 176:229–249, 2020), we have shown that the explicit algebraic Reynolds-stress (EARS) model, implemented in a single-column context, is able to capture the main features of a stable atmospheric boundary layer (ABL) for a range of stratification levels. We here extend the previous study and show that the same formulation and calibration of the EARS model also can be applied to a dry convective ABL. Five different simulations with moderate convective intensities are studied by prescribing surface heat flux and geostrophic forcing. The results of the EARS model are compared to large-eddy simulations of Salesky and Anderson (J Fluid Mech 856:135–168, 2018). It is shown that the EARS model performs well and is able to capture the counter-gradient heat flux in the upper part of the ABL due to the presence of the non-gradient term in the relation for vertical turbulent heat flux. The model predicts the full Reynolds-stress tensor and heat-flux vector and allows us to compare other important aspects of a convective ABL such as the profiles of vertical momentum variance. Together with the previous studies, we show that the EARS model is able to predict the essential features of the ABL. It also shows that the EARS model with the same model formulation and coefficients is applicable over a wide range of stable and moderately unstable stratifications.

A spectral problem of the Orr–Sommerfeld type is considered to describe unstable disturbances of oceanic geostrophic currents with a linear vertical velocity profile by taking into account the vertical diffusion of buoyancy and friction. Numerical solutions are obtained for different values of the dimensionless parameters of the problem. Calculations of the spectra of eigenvalues and growth rates are compared with those of a similar problem for an ideal fluid. It is shown that (a) dissipation expands the range of wave numbers of unstable disturbances, (b) dissipation can increase the growth rates of baroclinic disturbances, (c) disturbances driven by the instability of the critical layer can grow faster than baroclinic disturbances, and (d) currents with a width equal to or less than the Rossby radius can be unstable; nearly circular (axisymmetric) unstable disturbances (submesoscale eddies) can develop in narrow currents or frontal zones.

The predictability of high impact weather events over the North Atlantic is controlled by synoptic‐scale systems and the mesoscale structures embedded within them. Despite forecast uncertainty being greatest at small scales at the initial time, forecast error projects strongly onto synoptic and larger scales within days. Different stages of error growth have previously been identified including: convective instability, baroclinic instability and the influence of divergent outflow on the tropopause position, and interactions between disturbances at tropopause level.\ud \ud Evidence is presented for “predictability barriers” (PBs) identified with events on certain validation dates during the North Atlantic Waveguide and Downstream impact Experiment (NAWDEX) where ensemble spread grows more quickly than usual, but ensemble mean forecast error grows even faster. An advective mechanism for diabatic influence on the development of tropopause ridges is hypothesised to be linked to the PB events. A semi‐geostrophic balance tool is used to attribute the response of the 3‐D ageostrophic flow to geostrophic and diabatic forcing, enabling a novel diagnostic for Diabatically‐Induced Ageostrophic Advection of potential vorticity (DIAA).\ud \ud It is shown that predictability barriers are linked to events with strong diabatic influence on tropopause advection during the NAWDEX period. Error growth exceeds ensemble spread rate by approximately 4/3 during strong DIAA events, showing that predictive skill is considerably lower in these situations.

The zonal flow associated with cut-off lows (COLs) comprises two jet streaks of different spatial extents. The smaller scale jet streak, located north of the COLs, forms as a result of meridional divergence of vorticity advection and it is quasi-stationary, relative to the COLs. It dissipates as the COLs do the same. The larger scale jet streak gives rise to anticyclonic and equatorward Rossby wave breaking (RWB) as it propagates southeasterly to the base of the ridge, south of the COL and then northeasterly beyond that point. As the jet streak propagates it brings with it the anticyclonic barotropic shear that causes the Rossby waves to break. Its propagation is caused by zonal momentum advection by the zonal flow from jet streak entrance to its exit. As it propagates, its northwesterly/southeasterly orientation changes to one that is more zonal to become south-westerly/northeasterly at the end of the COL life cycle. This change in orientation is due to meridional advection of zonal momentum,where the meridional flow advects momentum southward (northward) at the jet streak entrance (exit). The jet streaks form a split jet structure and the winds between the streaks is decelerated by vorticity advection convergence. Because the flow and COL (and RWB) life cycle are coupled, understanding the dynamics that underlie the changes in the COL ambient flow contributes to resolving the outstanding RWB/COL causality problem.

Air–sea interactions play a critical role in the climate system. This study investigates wind-induced changes in the ocean surface temperature and sea ice cover feeding back onto the atmospheric circulation. This interaction was modeled in the Nordic seas, using a partial coupling method to constrain the ocean with prescribed wind forcing in an otherwise fully coupled Earth system model. This enabled the evaluation of not only the direct oceanic, but also the indirect atmospheric response to idealized forcing scenarios of perturbed winds over the Nordic seas. The results show that an anticyclonic wind anomaly forcing leads to significant surface cooling in the Greenland Sea mostly due to anomalous drift of sea ice. During winter, the cooling reduces the net surface heat flux to the atmosphere and increases sea level pressure. The pressure gradients result in anomalous geostrophic southerly winds, which locally are comparable both in direction and in velocity to the prescribed forcing anomalies, suggesting a positive feedback.

This study investigates the spatial structure and the seasonal occurrence of foehn wind to the east of the Andes using a flow blocking analysis in a 20-year climate simulation. The latter was performed by the Eta-CPTEC regional model at 50-km horizontal resolution. This version of the model includes a cut-cell scheme to represent topography and a finite-volume vertical advection scheme for dynamic variables. The results indicate that foehn wind more frequently blows during winter and spring on the eastern slopes of the Andes, except to the south of 37° S where it blows at all seasons. Higher mountains of the Central Andes (27° S–35° S) and the High Plateau (15° S–27° S) result in blocked foehn events, with a weak adjustment to the geostrophic balance. On the Central Andes, rain and snow on mountain tops may also contribute to generate foehn wind on the eastern slopes. The results show that a low pressure develops to the east of the Central Andes, and also to the east of the High Plateau when foehn blows. Lower mountains in Patagonia (to the south of 37° S) result in more frequent non-blocked foehn event, with better adjustment to the geostrophic balance., Facultad de Ciencias Astronómicas y Geofísicas

Using National Centers for Environmental Prediction reanalysis data for the period 28 June to 12 July during 2001 to 2013, the secondary circulation (SC) associated with the mei-yu front was quantitatively diagnosed by numerically solving a primitive version of the Sawyer-Eliassen equation. Results demonstrate that a direct SC exists near the mei-yu front zone during mid-summer and the synoptic-scale geostrophic deformations are the main factors determining SC structures. About 94% of the sinking strength and 61% of the ascending strength in the SC are induced by the geostrophic deformations. Other terms, such as diabatic heating, ageostrophic dynamical forcing, and frictional forcing, mainly influence the fine flow pattern of the SC. The forced SC produces a frontogenesis area tilting to the north with altitude. Further diagnosis clarifies the positive feedback involving the geostrophic shear forcing and vorticity frontogenesis in the upper-level mei-yu front zone. Furthermore, statistical results indicate that all 34 deep convection cases that occurred in the warm region of the meiyu front over the period 2004–2013 experienced high-level frontogenesis associated with along-jet cold advection. The cyclonic shear forcing “moved” the monsoon SC’s subsidence branch to the warm side of the mei-yu front and caused the subsidence branch to extend downwards to the lower troposphere, conducive to the initiation of deep convection in the warm region of the mei-yu front.

In the tropics, rainfall is coupled with waves in the form of, for example, convectively coupled equatorial waves (CCEWs) and the Madden–Julian oscillation (MJO). In perhaps the simplest viewpoint of CCEWs, the effects of moisture and convective adjustment can predict the basic aspects of their propagation and structure: reduced propagation speeds and reduced meridional length scales. Here, a similar simple viewpoint is investigated for the MJO’s propagation and structure. To do this investigation, budget analyses of a model MJO are first presented to illustrate and motivate the asymptotic scaling assumptions. Asymptotic models are then derived for the MJO. In brief, the structure of the asymptotic MJO is described by a tropical geostrophic balance, and the slow propagation arises from the dynamics of moist static energy. To be specific, if the moist static energy has a background vertical gradient that is asymptotically weak (i.e., a moist stability that is nearly neutral), then it supports a slowly propagating wave. Beyond these main aspects, other processes also have an influence, such as eddy diffusion of moisture. In comparing the simple viewpoints of CCEWs and the MJO, one main difference is in the propagation speeds: relative to a dry wave speed of 50 m s−1, the MJO has a speed of 5 m s−1, resulting from a reduction factor of 0.1 related to moist stability, whereas the basic CCEW speed is 15 m s−1, resulting from a reduction factor of the square root of 0.1, related to the square root of the moist stability.

Normal modes are used to investigate the contributions of geostrophic vortices and inertia–gravity waves to the energy spectrum of an idealized baroclinic wave simulation. The geostrophic and ageostrophic modal spectra (GE and AE, respectively) are compared to the rotational and divergent kinetic energy (RKE and DKE, respectively), which are often employed as proxies for vortex and wave energy. In our idealized f-plane framework, the horizontal modes are Fourier, and the vertical modes are found by solving an appropriate eigenvalue problem. For low vertical mode number n, both the GE and AE spectra are steep; however, for higher n, while both spectra are shallow, the AE is shallower than the GE and the spectra cross. The AE spectra are peaked at the Rossby deformation wavenumber knR, which increases with n. Analysis of the horizontal mode equations suggests that, for large wavenumbers k≫knR, the GE is approximated by the RKE, while the AE is approximated by the sum of the DKE and potential energy. These approximations are supported by the simulations. The vertically averaged RKE and DKE spectra are compared to the sum of the GE and AE spectra over all vertical modes; the spectral slopes of the GE and AE are close to those of the RKE and DKE, supporting the use of the Helmholtz decomposition to estimate vortices and waves in the midlatitudes. However, the AE is consistently larger than the DKE because of the contribution from the potential energy. Care must be taken when diagnosing the mesoscale transition from the intersection of the vortex and wave spectra; GE and AE will intersect at a different scale than RKE and DKE, despite their similar slopes.

International audience; An improved estimation of the surface currents in the Levantine Basin of the Mediterranean sea is crucial for a wide range of applications, including pollutants transport and nutrients distribution. This estimation remains challenging due to the scarcity or shortcomings of various data types used for this purpose. In this paper, we present an objective validation of a variational assimilation algorithm that blends geostrophic velocities derived from altimetry, wind-induced velocities, and drifter positions, to continuously obtain velocity corrections. The assessment of the validation impact was based on available independent in-situ data (current meters, gliders, and independent drifters) and satellite ocean color images. In all cases, the improvement was shown either qualitatively (position of the eddies) or quantitatively.

While wind is the primary driver of the oceanic general circulation, it kills the ocean’s most energetic motions., Wind is the primary driver of the oceanic general circulation, yet the length scales at which this energy transfer occurs are unknown. Using satellite data and a recent method to disentangle multiscale processes, we find that wind deposits kinetic energy into the geostrophic ocean flow only at scales larger than 260 km, on a global average. We show that wind removes energy from scales smaller than 260 km at an average rate of −50 GW, a process known as eddy killing. To our knowledge, this is the first objective determination of the global eddy killing scale. We find that eddy killing is taking place at almost all times but with seasonal variability, peaking in winter, and it removes a substantial fraction (up to 90%) of the wind power input in western boundary currents. This process, often overlooked in analyses and models, is a major dissipation pathway for mesoscales, the ocean’s most energetic scales.

Future changes in geostrophic winds over Europe and the North Atlantic region were studied utilizing output data from 21 CMIP5 global climate models (GCMs). Changes in temporal means, extremes, and the joint distribution of speed and direction were considered. In concordance with previous research, the time mean and extreme scalar wind speeds do not change pronouncedly in response to the projected climate change; some degree of weakening occurs in the majority of the domain. Nevertheless, substantial changes in high wind speeds are identified when studying the geostrophic winds from different directions separately. In particular, in northern Europe in autumn and in parts of northwestern Europe in winter, the frequency of strong westerly winds is projected to increase by up to 50%. Concurrently, easterly winds become less common. In addition, we evaluated the potential of the GCMs to simulate changes in the near-surface true wind speeds. In ocean areas, changes in the true and geostrophic winds are mainly consistent and the emerging differences can be explained (e.g., by the retreat of Arctic sea ice). Conversely, in several GCMs the continental wind speed response proved to be predominantly determined by fairly arbitrary changes in the surface properties rather than by changes in the atmospheric circulation. Accordingly, true wind projections derived directly from the model output should be treated with caution since they do not necessarily reflect the actual atmospheric response to global warming.

Estimates of the meridional mass transport of Antarctic Bottom Water, calculated using the coupled ocean-atmosphere Earth System Model on the basis of the original data assimilation method are presented. For assimilation, we use data of the latitudinal CTD sections of temperature and salinity of the WOCE international experiment in 1991–1995. Estimates of the current velocities of Antarctic Bottom Water with the assimilation of observational data are given. We used the author’s data-assimilation method, which was previously referred to as the generalized Kalman Filter (GKF) method. In this particular case, it coincides with the classical Kalman method (EnKF). We also present the estimates of mass transport based on a standard geostrophic dynamic scheme. It is shown that model calculations with data assimilation are qualitatively the same and are quantitatively close to the estimates of the geostrophic flow transport based on the dynamic method.

The increasing hub height of wind turbines aims at optimizing the wind energy yield at one location and offers the possibility to provide new areas for wind power, for example forests. Inhomogeneous environmental conditions of locations for wind turbines as well as the hub heights of more than 100 m cause challenges for flow models and their potential for wind power assessment. This includes special features of the wind field like low-level jets (LLJs), frequently observed local wind maxima in the nocturnal boundary layer. To characterize the dependencies of LLJs, the micro-scale model HIRVAC2D (HIgh Resolution Vegetation Atmosphere Coupler 2D) is applied in the study. The model HIRVAC2D is capable of modelling different vegetation types by explicitly considering the highly resolved structure of varying plant parameters. Beyond that, the model enables the resolution of temporally variable atmospheric circulation patterns during day- and night-time with typical thermal stratifications. In this way, HIRVAC2D is suitable to capture the nocturnal LLJ development and its characteristics. Results of several HIRVAC2D simulations are presented in order to deduce quantitatively the sensitivity of LLJs to vegetation and model parameters as well as meteorological quantities. It is shown that the geostrophic wind speed is an important criterion for the development of LLJs within a height range between 50 and 300 m. For a geostrophic wind speed of 4 m s−1, a nocturnal LLJ occurs remarkably more frequent as for a wind speed of 10 m s−1. To interpret and evaluate this result regarding possible wind power applications, a frequency distribution of the geostrophic wind speed was calculated over 30 years exemplarily at two locations using the meso-scale model COSMO in climate mode. Additionally, the type of land use has an impact on the height and intensity of LLJs. For a grassland site, the nocturnal LLJ is noticeably more frequent in the considered height range, but with a smaller wind speed and at a lower height above ground in comparison to deciduous or coniferous forests.

Atmospheric cyclones with strong winds significantly impact ocean circulation, regional sea surface temperature, and deep water formation across the global oceans. Thus they are expected to play a key role in a variety of energy transport mechanisms. Even though wind-generated internal gravity waves are thought to contribute significantly to the energy balance of the deep ocean, their excitation mechanisms are only partly understood. The present study investigates the generation of internal gravity waves during a geostrophic adjustment process in a Boussinesq model with axisymmetric geometry. The atmospheric disturbance is set by an idealized pulse of cyclonic wind stress with a Rankine vortex structure. Strength, radius and duration of the forcing are varied. The effect upon wave generation of stratification with variable mixed-layer depth is also examined. Results indicate that internal gravity waves are generated after approximately one inertial period. The outward radial energy flux is dominated by waves having structure close to vertical mode-1 and with frequency close to the inertial frequency. Less energetic higher mode waves are observed to be generated close to the sea floor underneath the storm. The total radiated energy corresponds to approximately 0.02% of the wind input. Deeper mixed-layer conditions as well as weaker stratification reduce this fraction. The low energy transfer rates suggest that other processes that drive vertical motion like surface heat fluxes, turbulent motion, mixed region collapse and storm translation are essential for significant energy extraction by internal gravity waves to occur.

—A spectral problem of the Orr–Sommerfeld type for describing stable and unstable perturbations of oceanic geostrophic flows with linear vertical velocity shear is considered (parabolic vertical velocity profile). Calculations of eigenvalues, increments of growth rate of unstable modes, and eigenfunctions of the fastest growing perturbations are presented. It is found that the instability of the geostrophic flow is observed over a wide range of horizontal scales: in addition to longwave perturbations with a phase velocity exceeding the maximum velocity of the geostrophic current and perturbations with the scales of the Rossby radius, shortwave modes exist with scales much smaller than the Rossby radius (submesoscale structures). The results of the model are used to describe intrusions in the Arctic basin, which are observed under conditions of absolutely stable stratification.

Based on daily precipitation data from the Chinese Meteorological Administration and reanalysis data from the National Centers for Environmental Prediction-Department of Energy, the character of low-frequency precipitation variability during the first rainy season (April–June) over South China and its corresponding atmospheric circulations in the mid-high latitudes are investigated. The results show that the precipitation anomalies during this period exhibit obvious quasi-biweekly oscillation (QBWO) features, with a period of 8–24 days. The influence of wave trains in the mid-high latitudes to low-frequency persistent heavy rain event (PHR-LF event, the 8–24-day filtered precipitation larger than one standard deviation of filtered time series and persisting at least three days over South China) is further discussed. During the first rainy season over South China, there are two low-frequency wave trains in the mid-high latitudes associated with the PHR-LF event—the wave train crossing the Eurasian continent and the wave train along the subtropical westerly jet. Analysis of wave activity flux indicates that the wave energy disperses toward eastern China along these two low-frequency wave trains from north to south and from west to east, and then propagates downward over South China. Accordingly, the disturbance of the relative vorticity of the cyclonic anomalies over eastern China is strengthened, which enhances the meridional gradient of relative vorticity. Owing to the transport of low-frequency relative vorticity and geostrophic vorticity by meridional wind, the ascending motion over South China intensifies and lasts for a long time, triggering a PHR-LF event. In addition, the tropical system is also a key factor to PHR-LF event. The QBWO of the convection over the South China Sea provide moisture for PHR-LF events, maintaining persistent rainfall and vertical ascending motion over South China.

Summary: This paper analyzes the mixing characteristics of the Subarctic Front (SAF) in the Kuroshio-Oyashio Confluence Region based on temperature, salinity, and current data obtained from surveys and remote sensing in June 2016. The frontal zone of the observed area is at 145°–151°E, 38°–41°N. The front is distributed between 25.5–26.7 σθ in a band pattern inclined from north to south and is deeper in the south. The region shallower than 200 m and distributed along the isopycnal of 25.9–26.1 σθ has the strongest horizontal temperature and salinity gradients, and the largest of the former can reach over 0.7°C/km. Diapycnal mixing of the SAF is mainly turbulent; it is stronger in the north than in the south. The region with stronger turbulence (Kρ > 10−3.5 m2/s) is distributed mainly in water layers within and under the front (26.1–26.7 σθ), showing that the SAF is shallower in the north and deeper in the south along the front. Symmetric instability may be the main factor causing strong turbulent mixing in the frontal zone. Double diffusion mixing is stronger in the south than in the north; the region with stronger double diffusion (Kθ > 10−4.5 m2/s) is distributed mainly in water layers within and above the front (25–26.5 σθ) on the southern side of the SAF. These water layers are dominated mainly by ``salt-fingering’’ double diffusion, with only a few water layers dominated by ``diffusive layering’’ double diffusion mixing in middle and lower waters deeper than 300 m. Keywords: Kuroshio–Oyashio confluence region, Subarctic front, Mixing, Turbulent eddy diffusivity, Thermal diffusivity

Horizontal velocity spirals with a clockwise rotation (downward looking) rate of 1.7° m−1, on average, were observed in the western and northern Yellow Sea from December 2006 to February 2007. With the observed thermal wind relation, the beta-spiral theory was used to explain the dynamics of spirals. It was found that the horizontal diffusion of geostrophic vortex stretching is likely to be a major mechanism for generating geostrophic spirals. Vertical advection associated with surface/bottom Ekman pumping and topography-induced upwelling is too weak to support these spirals. Strong wind stirring and large heat loss in wintertime lead to weak stratification and diminish the effects of vertical advection. The cooling effect and vertical diffusion are offset by an overwhelming contribution of horizontal diffusion in connection with vortex stretching. The Richardson number-dependent vertical eddy diffusivity reaches a magnitude of 10−4 m2 s−1 on average. An eddy diffusivity of 2870 m2 s−1 is required for dynamic balance by estimating the residual term. This obtained value of 10−4 m2 s−1 is in good agreement with the estimation in terms of observed eddy activities. The suppressed and unsuppressed diffusivities in the observation region are 2752 and 2881 m2 s−1, respectively, which supports a closed budget for velocity rotation.

Sea level measurements from an absolute dynamic typography dataset for the period of 1993–2015 were used to study the variability in geostrophic circulation in the Labrador and Irminger seas. Analysis of current intensity trends calculated by geostrophic formulas has shown a weakening in cyclonic circulation observed in these two convection centers. The intensity trends are negative in all the main currents making up the Subpolar Gyre: the Labrador, East Greenland, West Greenland, and Irminger currents. On average, the negative intensity trends in the indicated currents are –0.3 cm s–1 yr–1, in individual areas reaching –0.5 cm s–1 yr–1, which corresponds to a decrease in current speed of 6.9–11.5 cm/s for the 23-year interval. Weakening of cyclonic circulation is likely predetermined by weakening convection processes and a decrease in the dome-shaped curvature of isopycnal surfaces generated by upwelling of weakly stratified water closer to the surface. This is confirmed by analysis of satellite data: in periods when deep convection weakens, the intensity of cyclonic circulation decreases at the margin of the Labrador and Irminger seas. Simultaneously with weakening currents, an increase in intensity is observed from the outer (with respect to the convection centers) margin of the considered currents, including in coastal zones.

The turbulence generated by wind shear is described at grey-zone resolutions using a theoretical neutral boundary layer based on atmospheric conditions constructed from measurements from the CASES-99 field campaign. Six-metre-resolution large-eddy simulations (LES) are performed to access the “true” resolved turbulence for two cases, corresponding to a forcing of the boundary layer by zonal geostrophic wind speeds of $$10\,\text {m}\,\text {s}^{-1}$$ and $$20\,\text {m}\,\text {s}^{-1}$$ . The LES fields are subject to a coarse-graining procedure in order to compute turbulence diagnostics in the grey zone, with the robustness and weakness of various averaging procedures tested, for which simple top-hat averaging is found to be both suitable and accurate. In addition, the “true” resolved and subgrid-scale fluxes, variances, turbulent kinetic energy and production terms are quantified on various scales. The grey zone of turbulence is defined as the range of scales where 10–90% of turbulence is resolved, which here ranges from resolutions of 25– $$800\,\hbox {m}$$ ( $$0.031$$ ). The turbulence parametrizations, which are tested in the Meso-NH model by running simulations at resolutions from the LES scale to the mesoscale, fail to produce the correct turbulence regardless of resolution.

Characteristics of flow in the planetary boundary layer (PBL) strongly affect the design of tall structures. PBL modelling in building codes, based as it is on empirical data from the 1960s and 1970s, differs significantly from contemporary PBL models, which account for both “neutral” flows, and “conventionally neutral” flows. PBL heights estimated in these relatively sophisticated models are typically approximately half as large as those obtained using the classical asymptotic similarity approach, and are one order of magnitude larger than those specified in North American and Japanese building codes. A simple method is proposed for estimating the friction velocity and PBL height as functions of specified surface roughness and geostrophic wind speed. Based on published results, it is tentatively determined that, even at elevations as high as 800 m above the surface, the contribution to the resultant mean flow velocity of the component V normal to the surface stress is negligible and the veering angle is of the order of only 5\(^{\circ }\). This note aims to encourage dialogue between boundary-layer meteorologists and structural engineers.

The aim of this study is to analyse the influence of atmospheric circulation over Central Europe on the long-term variability of relative sunshine duration (RSD) and air temperature in Poland. The analysis of data from 23 stations located in different parts of the country has made it possible to indicate to what extent the trends of these two climate elements are related to atmospheric circulation and whether they occur throughout the year or only in certain months. To this end, the author's own classification of circulation types was used, which relies on a set of indices related to the velocity and direction of geostrophic flow, as well as total shear vorticity, calculated using a gridded data set of sea-level pressures. The results of the research showed that atmospheric circulation over Central Europe has a significant impact on RSD and air temperature in Poland, but not on their upward trends in recent decades. The increase in the average annual values of both these climate elements was marked in almost every circulation type. This contradicts the assumption of classical synoptic climatology, which assumes that circulation types themselves are usually considered to be constant, i.e. that changes in air temperature can be related only to an increase or decrease in the incidence of the circulation types.

The geostrophic drag law (GDL), which predicts the geostrophic drag coefficient and the cross-isobaric angle, is relevant for meteorological applications such as wind energy. For conventionally neutral atmospheric boundary layers (CNBLs) capped by an inversion, the GDL coefficients A and B are affected by the inversion strength and latitude, expressible via the ratio of the Brunt–Vaisala frequency (N) to the Coriolis parameter (f). We present large-eddy simulations (LES) covering a wider range of N/|f| than considered previously, and show that A and B obtained from carefully performed LES collapse to a single curve when plotted against N/|f|. This verifies the GDL for CNBLs over an extended range of N/|f| within LES. Additionally, in agreement with atmospheric observations, we show that using A = 1.9 and B = 4.4 accurately predicts the geostrophic drag coefficient in the limit of weak inversion strength or high latitude ((Formula presented.)). However, due to the strong dependence of B on N/|f|, corresponding predictions for the cross-isobaric angle are less accurate. As we find significant deviations between the LES results and the original parameterization of the GDL for CNBLs, we update the corresponding model coefficients.

Motivated by the abundance of low clouds in the subtropics, where the easterly trade winds prevail, we study the role of shallow convection in the momentum budget of the trades. To this end, we use ICON-LEM hindcasts run over the North Atlantic for 12 days corresponding to the NARVAL1 (winter) and NARVAL2 (summer) flight campaigns. The simulation protocol consists of several nested domains, and we focus on the inner domains (≈100 × 100 km2) which have been run at resolutions of 150–600 m and are forced by analysis data, thus exhibiting realistic conditions. Combined, the resolved advection and the subgrid stresses decelerate the easterly flow over a frictional layer that balances the prevailing geostrophic wind forcing. Irrespective of the horizontal resolution, this layer is about 2 km deep in the strong winter trades and 1 km in summer, as winds and geostrophic forcing weaken and cloudiness reduces. The unresolved processes are strongest near the surface and are well captured by traditional K-diffusion theory, but convective-scale motions which are not considered in K-diffusion theory contribute the most in the upper part of the mixed layer and are strongest just below cloud base. The results point out that convection in the mixed layer – the roots of trade-wind cumuli and subcloud-layer circulations – play an important role in slowing down easterly flow below cloud base (but little in the cloud layer itself), which helps make the zonal wind jet more distinct. Most of the friction within the clouds and near the wind jet stems from smaller-scale turbulence stresses.

This study investigates the evolution of German Bight (southeastern North Sea) storminess from 1897 to 2018 through analysing upper quantiles of geostrophic wind speeds, which act as a proxy for past storm activity. Here, geostrophic wind speeds are calculated from triplets of mean sea level pressure observations that form triangles over the German Bight. The data used in the manuscript are provided by the International Surface Pressure Databank and the national meteorological services of Denmark, Germany, and the Netherlands. The derivation of storm activity is achieved by enhancing the established triangle proxy method via combining and merging storminess time series from numerous partially overlapping triangles in an ensemble‐like manner. The utilized approach allows for the construction of robust, long‐term and subdaily German Bight storminess time series. Further, the method provides insights into the underlying uncertainty of the time series. The results show that storm activity over the German Bight is subject to multidecadal variability. The latest decades are characterized by an increase in activity from the 1960s to the 1990s, followed by a decline lasting into the 2000s and below‐average activity up until present. The results are backed through a comparison with reanalysis products from four datasets, which provide high‐resolution wind and pressure data starting in 1979 and offshore wind speed measurements taken from the FINO‐WIND project. This study also finds that German Bight storminess positively correlates with storminess in the NE Atlantic in general. In certain years, however, notably different levels of storm activity in the two regions can be found, which likely result from shifted large‐scale circulation patterns.

Given a flow domain D with subdomains D1 and D2, piecewise potential vorticity inversion (PPVI) inverts a potential vorticity (PV) anomaly in D2 and assumes vanishing PV in D1 where boundary conditions must be taken into account. It is a widely held view that the PV anomaly exerts a far-field influence on D1, which is revealed by PPVI. Tests of this assertion are conducted using a simple quasigeostrophic model where an upper layer D2 contains a PV anomaly and D1 is the layer underneath. This anomaly is inverted. Any downward physical impact of PV in D2 must also be represented in the results of a downward piecewise density inversion (PDI) based on the hydrostatic relation and the density in D2 as following from PPVI. There is no doubt about the impact of the mass in D2 on the flow in the lower layer D1. Thus results of PPVI and PDI have to agree closely. First, PPVI is applied to a locally confined PV anomaly in D2. There is no far-field “response” in D1 if stationarity is imposed. Modifications of boundary conditions lead to “induced” flows in D1 but the results of PPVI and PDI differ widely. This leads to a simple proof that there is no physical far-field influence of PV anomalies in D2. Wave patterns of the streamfunction restricted to D2 are prescribed in a second series of tests. The related PV anomalies are obtained by differentiation and are also confined to D2 in this case. This approach illustrates the basic procedure to derive PV fields from observations which excludes a far-field response.

During the past 25 years, altimetric observations of the ocean surface from space have been mapped to provide two dimensional sea surface height (SSH) fields that are crucial for scientific research and operational applications. The SSH fields can be reconstructed from conventional altimetric data using temporal and spatial interpolation. For instance, the standard Developing Use of Altimetry for Climate Studies (DUACS) products are created with an optimal interpolation method that is effective for both low temporal and low spatial resolution. However, the upcoming next-generation SWOT mission will provide very high spatial resolution but with low temporal resolution. The present paper makes the case that this temporal–spatial discrepancy induces the need for new advanced mapping techniques involving information on the ocean dynamics. An algorithm is introduced, dubbed the BFN-QG, that uses a simple data assimilation method, the back-and-forth nudging (BNF), to interpolate altimetric data while respecting quasigeostrophic (QG) dynamics. The BFN-QG is tested in an observing system simulation experiments and compared to the DUACS products. The experiments consider as reference the high-resolution numerical model simulation NATL60 from which are produced realistic data: four conventional altimetric nadirs and SWOT data. In a combined nadirs and SWOT scenario, the BFN-QG substantially improves the mapping by reducing the root-mean-square errors and increasing the spectral effective resolution by 40 km. Also, the BFN-QG method can be adapted to combine large-scale corrections from nadir data and small-scale corrections from SWOT data so as to reduce the impact of SWOT correlated noises and still provide accurate SSH maps.

Author(s): Vishnu, S; Boos, WR; Ullrich, PA; O'Brien, TA | Abstract: Cyclonic low-pressure systems (LPS) produce abundant rainfall in South Asia, where they are traditionally categorized as monsoon lows, monsoon depressions, and more intense cyclonic storms. The India Meteorological Department (IMD) has tracked monsoon depressions for over a century, finding a large decline in their number in recent decades, but their methods have changed over time and do not include monsoon lows. This study presents a fast, objective algorithm for identifying monsoon LPS and uses it to assess interannual variability and trends in reanalyses. Variables and thresholds used in the algorithm are selected to best match a subjectively analyzed LPS data set while minimizing disagreement between four reanalyses in a training period. The stream function of 850 hPa horizontal wind is found to be optimal in this sense; it is less noisy than vorticity and represents the complete nondivergent wind, even when flow is not geostrophic. Using this algorithm, LPS statistics are computed for five reanalyses, and none show a detectable trend in monsoon depression counts since 1979. Both the Japanese 55-year Reanalysis (JRA-55) and the IMD data set show a step-like reduction in depression counts when they began using geostationary satellite data, in 1979 and 1982, respectively; the 1958–2018 linear trend in JRA-55, however, is smaller than in the IMD data set, and its error bar includes 0. There are more LPS in seasons with above-average monsoon rainfall and in La Nina years, but few other large-scale modes of interannual variability are found to modulate LPS counts, lifetimes, or track length consistently across reanalyses.

In this study, we show the relationship between sea-level anomalies (SLA) and upper-ocean parameters in the Equatorial Indian Ocean (EIO). This work also focuses on the variability of SLA obtained from satellite altimeter data in different spatial and temporal scales and its relationship with computed ocean heat content (OHC), dynamic height (DH), and thermocline depth (20 °C isotherm: D20) during 1993–2015. SLA showed low Pearson’s correlation coefficient (CC) with upper-ocean parameters over central EIO resembling a “Monopole” pattern. The Array for Real-time Geostrophic Oceanography (ARGO) in situ profile data in the central EIO also confirmed this. SLA over this monopole showed low correlations with all parameters as compared with eastern and western EIO. These findings show a clear signature of a persisting sea-level monopole in the central EIO. Oscillating SLA over western and eastern EIO during summer and winter monsoon months is found to be responsible for locking this monopole in the central EIO. Both SLA and OHC increased in EIO during 2006–2015 compared with 1993–2005. The month of January showed different east–west trends at different times. This trend during 1993–2015 is neutral, but it shifted from negative during 1993–2005 to positive during 2006–2015.

The decay of the Convective Boundary Layer (CBL) is studied using large-eddy simulations of free and advective CBLs, in which surface heat supply is suddenly cut off. After the cutoff, coherent convective circulations last about one convective time scale and then fade away. In the mixed layer, the decay time scale increases with height, indicating that nonlocal eddies decay slower than near-surface local eddies. The slower decay of turbulence in the middle of CBL than near-surface turbulence is reconfirmed from the analysis of pattern correlations of perturbations of vertical velocity. Perturbations of potential temperature and scalar concentration decay faster and slower than vertical velocity perturbations, respectively. A downward propagation of negative heat flux and its oscillation are found and a quadrant analysis reveals that warmer air sinking events are responsible for the downward propagation. The fourth quadrant events seem to be induced by demixing of air parcels, entrained from above the CBL. The advective CBL simulation with geostrophic wind illustrates that near-surface eddies are mechanically generated and they decelerate flow from the bottom up in the CBL/residual layer. The two-dimensional spectra show the height- and scale-dependent characteristics of decaying convective turbulence again in the free and advective boundary layer simulations.

Inertia-gravity waves (IGWs) play an essential role in the terrestrial atmospheric dynamics as they can lead to energy and momentum flux when propagating upwards. An open question is to which extent IGWs contribute to the total energy and to flattening of the energy spectrum observed at the mesoscale. In this work, we present an experimental investigation of the energy distribution between the large-scale balanced flow and the small-scale imbalanced flow. Weakly nonlinear IGWs emitted from baroclinic jets are observed in the differentially heated rotating annulus experiment. Similar to the atmospheric spectra, the experimental kinetic energy spectra reveal the typical subdivision into two distinct regimes with slopes $k^{-3}$ for the large scales and $k^{-5/3}$ for the small scales. By separating the spectra into the vortex and the wave component, it emerges that at the large-scale end of the mesoscale the gravity waves observed in the experiment cause a flattening of the spectra and provide most of the energy. At smaller scales, our data analysis suggests a transition towards a turbulent regime with a forward energy cascade up to where dissipation by diffusive processes occurs., Comment: 18 pages, 6 figures

Surface wind speed (SWS) from meteorological observation, global atmospheric reanalysis, and geostrophic wind speed (GWS) calculated from surface pressure were used to study the stilling and recovery of SWS over China from 1960 to 2017. China experienced anemometer changes and automatic observation transitions in approximately 1969 and 2004, resulting in SWS inhomogeneity. Therefore, we divided the entire period into three sections to study the SWS trend, and found a near zero annual trend in the SWS in China from 1960 to 1969, a significant decrease of -0.24 m/s decade-1 from 1970 to 2004, and a weak recovery from 2005 to 2017. By defining the 95th and 5th percentiles of monthly mean wind speeds as strong and weak winds, respectively, we found that the SWS decrease was primarily caused by a strong wind decrease of -8 % decade-1 from 1960 to 2017, but weak wind showed an insignificant decreasing trend of -2 % decade-1. GWS decreased with a significant trend of -3 % decade-1 before the 1990s, during the 1990s, GWS increased with a trend of 3 % decade-1 whereas SWS continued to decrease with a trend of 10 % decade-1. Consistent with SWS, GWS demonstrated a weak increase after the 2000s. After detrended, both of SWS and GWS showed synchronous decadal variability, which is related to the intensity of Aleutian low pressure over the North Pacific. However, current reanalyses cannot reproduce the decadal variability, and can not capture the decreasing trend of SWS either.

Floating plastic debris is an increasing source of pollution in the world's oceans. Numerical simulations using models of ocean currents give insight into the transport and distribution of microplastics in the oceans, but most simulations do not account for the oscillating flow caused by global barotropic tides. Here, we investigate the influence of barotropic tidal currents on the transport and accumulation of floating microplastics, by numerically simulating the advection of virtual plastic particles released all over the world's oceans and tracking these for 13 years. We use geostrophic and surface Ekman currents from GlobCurrent and the currents caused by the four main tidal constituents (M 2, S 2, K 1, and O 1) from the FES model. We analyze the differences between the simulations with and without the barotropic tidal currents included, focusing on the open ocean. In each of the simulations, we see that microplastic accumulates in regions in the subtropical gyres, which is in agreement with observations. The formation and location of these accumulation regions remain unaffected by the barotropic tidal currents. However, there are a number of coastal regions where we see differences when the barotropic tidal currents are included. Due to uncertainties of the model in coastal regions, further investigation is required in order to draw conclusions in these areas. Our results suggest that, in the global open ocean, barotropic tidal currents have little impact on the transport and accumulation of floating microplastic and can thus be neglected in simulations aimed at studying microplastic transport in the open ocean., Key Points Barotropic tidal currents have little impact on transport and accumulation of floating microplastic in the global open oceanSimulations of plastic dispersion in the open ocean do not need to take barotropic tides into accountA computationally efficient implementation for advecting virtual particles by tidal currents is provided

Five planetary-boundary-layer parametrizations of the Weather Research and Forecasting model are compared with respect to their ability to simulate the very stable and the weakly stable regimes of the stable boundary layer. This is performed for single column models where the large-scale mechanical forcing is represented by geostrophic wind speeds ranging from 0.5 to 12 m $$\hbox {s}^{-1}$$ . The performance of the models is assessed by contrasting the relationships they produce between the turbulence velocity scale and the mean wind speed, between potential temperature gradient and the mean wind speed, and between the flux and gradient Richardson numbers. The level-2.5 Mellor–Yamada–Nakanishi–Niino parametrization simulates the very stable regime the best, mainly because its heat eddy diffusivity decreases with respect to the momentum eddy diffusivity as the stability increases, while the same is not true for the other parametrizations considered.

International audience; Using the two-layer moist-convective rotating shallow water model, we study the process of relaxation (adjustment) of localized large-scale pressure anomalies in the lower equatorial troposphere, and show that it engenders coherent structures strongly resembling the Madden Julian Oscillation (MJO) events, as seen in vorticity, pressure, and moisture fields. We demonstrate that baroclinicity and moist convection substantially change the scenario of the quasi-barotropic "dry" adjustment, which was established in the framework of one-layer shallow water model and consists, in the long-wave sector, in the emission of equatorial Rossby waves, with dipolar meridional structure, to the West, and of equatorial Kelvin waves, to the East. If moist convection is strong enough, a dipolar cyclonic structure, which appears in the process of adjustment as a Rossby-wave response to the perturbation, transforms into a coherent modon-like structure in the lower layer, which couples with a baroclinic Kelvin wave through a zone of enhanced convection and produces, at initial stages of the process, a self-sustained slowly eastward-propagating zonally-dissymmetrical quadrupolar vorticity pattern. At the same time, a weaker quadrupolar structure of opposite sign arises in the upper layer, the whole picture similar to the active phase of the MJO events. The baroclinic Kelvin wave then detaches from the dipole, which keeps slow eastward motion, and circumnavigates the Equator, catching up and interacting with the dipole.

Gap wind jets (Tehuano winds) trigger supersquirts of colder water and mesoscale asymmetric dipoles in the Gulf of Tehuantepec (GT). However, the effects of successive gap wind jets on dipoles and their effects inside eddies have not yet been studied. Based on the wind fields, geostrophic currents, and surface drifter dispersion, this research documented three dipoles triggered and modified by Tehuano winds. Once a dipole develops, successive gap wind jets strengthen the vortices, and the anticyclonic eddy migrates southwestward while the cyclonic eddy is maintained on the east side of the GT. During the wind relaxation stage, the cyclonic eddy may propagate westward, but due to the subsequent re-intensification of the Tehuano winds, the vortex could break down, as was suggested by surface drifter dispersion pattern and geostrophic field data. The effect of the Tehuano winds was evaluating via eddy-Ekman pumping. Under Tehuano wind conditions, Ekman downwelling (upwelling) inside the anticyclonic (cyclonic) eddies may reach ~ -2.0 (0.5) m d-1 and decrease as the wind weakens. In the absence of Tehuano winds, Ekman downwelling inside the anticyclonic eddy was ~ 0.1 (-0.1) m d-1. The asymmetry of downwelling and upwelling inside eddies during Tehuano wind events may be associated with Tehuano wind forcing.

The along-shelf momentum balance of the Mid-Atlantic Bight (MAB) coastal ocean includes a significant contribution from the along-shelf gradient in sea level. This sea level tilt, of order 10−7, and other features of the mean sea level are not captured well in global mean dynamic topography (MDT) derived from hydrographic observations or satellite altimetry and gravity data, and is poorly represented in global and basin scale dynamical models. This is problematic for applications that would use coastal satellite altimeter data to estimate total water level above datum. We have produced a MDT for the MAB using the Regional Ocean Modeling System (ROMS) with 4-dimensional variational (4D-Var) data assimilation configured to obtain climatological annual and monthly mean results. The observations assimilated were a regional hydrographic climatology of temperature and salinity, and long-term mean velocity from HF-radar, shipboard ADCP, and current-meters. Assimilation adjusts the 3-dimensional ocean state, boundary conditions, and air-sea fluxes to minimize the model-data misfit. The assimilation of mean velocity data is vital to obtaining a realistic circulation result. The MDT exhibits a strong across-shelf sea level slope in geostrophic balance with the southwestward mean flow. The subtle along-shelf tilt is recovered and is relatively uniform throughout the MAB inside the 50 m isobath, but on the southern outer shelf it reverses sign and drives significant across-isobath flow, partially draining the southward mean transport. In the north, across-shelf flow is offshore in the surface and bottom Ekman layers, but largely balanced locally by inflow in the interior depth range.

The meridional geostrophic transport (MGT) in the interior tropical North Pacific Ocean is estimated based on global ocean heat and salt content data. The decadal variations of the zonally and vertically integrated MGT in the tropical North Pacific Ocean are found to precede the Pacific decadal oscillation (PDO) by 1–3 years. The dynamics of the MGT are analyzed based on Sverdrup theory. It is found that the total meridional transport variability (MGT plus Ekman) is dominated by the MGT variability having positive correlations with the PDO index. The Sverdrup transports differ from the total meridional transport significantly and have insignificant correlations with PDO index, suggesting that the MGT variability is not controlled by the Sverdrup dynamics. In comparison, the simulated meridional transport variability in the models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) and the Ocean General Circulation Model for the Earth Simulator are dominated by the Sverdrup transports, having insignificant correlations with the simulated PDO indices. The comparison suggests that the non-Sverdrup component in the MGT is important for the predictability of PDO and that significant deficiencies exist in these models in simulating a realistic structure of the tropical ocean gyre variability and predicting the decadal climate variations associated with it.