Horizontal Momentum Advection Research Articles

Understanding Energy Pathways in the Gulf Stream

Abstract The Gulf Stream (GS) is one of the strongest ocean currents on the planet. Eddy-rich resolution models are needed to properly represent the dynamics of the GS; however, kinetic energy (KE) can be in excess in these models if not dissipated efficiently. The question of how and how much energy is dissipated and in particular how it flows through ocean scales thus remains an important and largely unanswered question. Using a high-resolution (∼2 km) ocean model [Coastal and Regional Ocean Community (CROCO)], we characterize the spatial and temporal distribution of turbulent cascades in the GS based on a coarse-grained method. We show that the balanced flow is associated with an inverse cascade while the forward cascade is explained by ageostrophic advection associated with frontogenesis. Downscale fluxes are dominant at scales smaller than about 20 km near the surface and most intense at the GS North Wall. There is also strong seasonal variability in KE flux, with the forward cascade intensifying in winter and the inverse cascade later in spring. The forward cascade, which represents an interior route to dissipation, is compared with both numerical and boundary dissipation processes. The contribution of interior dissipation is an order of magnitude smaller than that of the other energy sinks. We thus evaluate the sensitivity of horizontal momentum advection schemes on energy dissipation and show that the decrease in numerical dissipation in a high-order scheme leads to an increase in dissipation at the boundaries, not in the downscale flux.

Investigation of the ΣO/N2 Depletion With Latitudinally Tilted Equatorward Boundary Observed by GOLD During the Geomagnetic Storm on April 20, 2020

AbstractA long‐lived (>10 hr) O/N2 column density ratio (∑O/N2) depletion at middle latitudes was observed by the Global‐scale Observations of the Limb and Disk (GOLD) during the 20 April 2020 geomagnetic storm. The observed ∑O/N2 depletion tilts latitudinally at equatorward boundary with the lowest latitude of ∼20°N at ∼75°W, and changes orientation near 75°W, with a north‐westward tilt to the west of ∼75W and a north‐eastward tilt to the east of ∼75W. The National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model shows the similar ∑O/N2 depletion patterns and magnitudes, but the longitudes with the most equatorward were ∼100°W and moved westward from ∼100°W to ∼140°W during the storm recovery period. Horizontal winds play major role in forming this middle‐latitude ∑O/N2 depletion structure. The strongest equatorward winds in the longitude sector near the magnetic pole resulted in the most equatorward expansion of all longitudes, generating the depletion with latitudinally tilted equatorward boundary. In addition, the high‐latitude westward zonal wind and the middle‐latitude eastward zonal wind maintained the long existence of this depletion in the GOLD Field of View. The further analysis on neutral winds shows that storm‐time meridional winds at middle latitudes in the mid‐low thermosphere are mainly dominated by pressure gradient force, while zonal winds are determined primarily by horizontal momentum advection.

Differences in Tropical Rainfall in Aquaplanet Simulations With Resolved or Parameterized Deep Convection

AbstractThis study investigates the effects of resolved deep convection on tropical rainfall and its multi‐scale variability. A series of aquaplanet simulations are analyzed using the Model for Prediction Across Scales‐Atmosphere with horizontal cell spacings from 120 to 3 km. The 3‐km experiment uses a novel configuration with 3‐km cell spacing between 20°S and 20°N and 15‐km cell spacing poleward of 30°N/S. A comparison of those experiments shows that resolved deep convection yields a narrower, stronger, and more equatorward intertropical convergence zone, which is supported by stronger nonlinear horizontal momentum advection in the boundary layer. There is also twice as much tropical rainfall variance in the experiment with resolved deep convection than in the experiments with parameterized convection. All experiments show comparable precipitation variance associated with Kelvin waves; however, the experiment with resolved deep convection shows higher precipitation variance associated with westward propagating systems. Resolved deep convection also yields at least two orders of magnitude more frequent heavy rainfall rates (>2 mm hr−1) than the experiments with parameterized convection. A comparison of organized precipitation systems demonstrates that tropical convection organizes into linear systems that are associated with stronger and deeper cold pools and upgradient convective momentum fluxes when convection is resolved. In contrast, parameterized convection results in more circular systems, weaker cold pools, and downgradient convective momentum fluxes. These results suggest that simulations with parameterized convection are missing an important feedback loop between the mean state, convective organization, and meridional gradients of moisture and momentum.

Role of Advection on the Evolution of Near-Surface Temperature and Wind in Urban-Aware Simulations

AbstractThe role of advection of heat and momentum on the evolution of near-surface temperature and wind is evaluated in urban-aware simulations over Houston, Texas, under dry conditions on a light-wind day. Two sets of experiments, each consisting of four simulations using different planetary boundary layer (PBL) schemes, were conducted over 48 h using the default urban scheme (BULK) and the single-layer urban canopy model (SLUCM) available within the Weather Research and Forecasting Model. We focus on understanding and quantifying the role played by temperature and momentum advection, particularly on the windward and leeward sides of the city. Previous studies have largely ignored any quantitative analysis of impacts from the advection of momentum over an urban area. The horizontal advection of temperature was found to be more important in the BULK because of the larger surface temperature gradient caused by warmer surface temperatures over urban areas than in the SLUCM. An analysis of the momentum budget shows that horizontal advection of zonal and meridional momentum plays a prominent role during the period of peak near-surface winds and that this effect is more pronounced in the windward side of the city. The local tendency in peak winds in the leeward side lags that in the windward side by about 1–2 h, similar to the lag found in horizontal momentum advection. The sensitivity of the results to different urban and PBL schemes was explored. The results imply that representation and influence of land-use patterns via sophisticated urban parameterizations generate locally driven winds that best resemble observations.

Diurnal Variations of Low-Level Winds and Precipitation Response to Large-Scale Circulations during a Heavy Rainfall Event

AbstractA succession of MCSs developed during the last week of October 2016 and produced extreme heavy rainfall in central China. The event underwent an evident shift from a mei-yu-like warm scenario to an autumn cold scenario. Diurnal cycles of rainfall and low-level winds may be modulated by the shifting of large-scale atmospheric conditions. We conducted observational analyses and numerical experiments to examine how large-scale circulations influenced rainfall systems through diurnally varying processes. The results show that, in the first half (warm) period of the event, intense rainfall mostly occurred in eastern-central China with an early morning peak. It was closely related to a nocturnal southwesterly low-level jet (NLLJ) on the flank of the western Pacific subtropical high. The NLLJ formed near midnight in southern China where ageostrophic wind rotated clockwise due to Blackadar’s inertial oscillation. The NLLJ extended downstream to central China during the predawn hours due to the horizontal advection of momentum. Both the formation and extension of the NLLJ were supported by an enhanced subtropical high that provided relatively warm conditions with surface heating for boundary layer inertial oscillation and strong background southwesterly winds for momentum transport. The NLLJ induced MCSs at its northern terminus where the low-level ascent, moisture flux convergence, and convective instability were enhanced during the predawn hours. In the second half period with an intrusion of cold air, the diurnal amplitude of low-level winds became small under relatively cold and cloudy conditions. Moderate rainfall tended to occur in western-central China with a peak after midnight, most likely due to frontogenetic processes, upslope lifting, and nighttime cloud-top cooling.

The Role of Equatorial Rossby Waves in Tropical Cyclogenesis. Part I: Idealized Numerical Simulations in an Initially Quiescent Background Environment

Abstract This two-part series of papers examines the role of equatorial Rossby (ER) waves in tropical cyclone (TC) genesis. To do this, a unique initialization procedure is utilized to insert n = 1 ER waves into a numerical model that is able to faithfully produce TCs. In this first paper, experiments are carried out under the idealized condition of an initially quiescent background environment. Experiments are performed with varying initial wave amplitudes and with and without diabatic effects. This is done to both investigate how the properties of the simulated ER waves compare to the properties of observed ER waves and explore the role of the initial perturbation strength of the ER wave on genesis. In the dry, frictionless ER wave simulation the phase speed is slightly slower than the phase speed predicted from linear theory. Large-scale ascent develops in the region of low-level poleward flow, which is in good agreement with the theoretical structure of an n = 1 ER wave. The structures and phase speeds of the simulated full-physics ER waves are in good agreement with recent observational studies of ER waves that utilize wavenumber–frequency filtering techniques. Convection occurs primarily in the eastern half of the cyclonic gyre, as do the most favorable conditions for TC genesis. This region features sufficient midlevel moisture, anomalously strong low-level cyclonic vorticity, enhanced convection, and minimal vertical shear. Tropical cyclogenesis occurs only in the largest initial-amplitude ER wave simulation. The formation of the initial tropical disturbance that ultimately develops into a tropical cyclone is shown to be sensitive to the nonlinear horizontal momentum advection terms. When the largest initial-amplitude simulation is rerun with the nonlinear horizontal momentum advection terms turned off, tropical cyclogenesis does not occur, but the convectively coupled ER wave retains the properties of the ER wave observed in the smaller initial-amplitude simulations. It is shown that this isolated wave-only genesis process only occurs for strong ER waves in which the nonlinear advection is large. Part II will look at the more realistic case of ER wave–related genesis in which a sufficiently intense ER wave interacts with favorable large-scale flow features.

Evolution of the quasi 2‐day wave during January 2006

We examine the evolution of the quasi 2‐day wave in the middle atmosphere during the period from 5 January to 5 February 2006 using global synoptic meteorological fields from the high‐altitude Navy Operational Global Atmospheric Prediction System Advanced Level Physics, High Altitude (NOGAPS‐ALPHA) forecast‐assimilation system. This period is characterized by a high level of planetary wave activity in the Northern Hemisphere (winter) extratropical stratosphere prior to a sudden stratospheric warming (SSW) on 20 January 2006. Space‐time spectral analysis of 6‐hourly NOGAPS‐ALPHA fields finds the largest quasi 2‐day wave amplitudes in the Southern Hemisphere (summer) extratropical upper mesosphere. Eliassen‐Palm flux diagnostics indicate that this extratropical quasi 2‐day wave is related to baroclinic instability along the equatorward flank of the summer easterly jet. The quasi 2‐day wave is also evident in NOGAPS‐ALPHA water vapor fields near the tropical stratopause and is related to barotropic instability. We find that the strong planetary wave activity leading up to the SSW produced an enhanced northward component of the residual meridional circulation that influenced the background zonal winds and, by extension, the quasi 2‐day wave forcing in both the tropical and extratropical regions. In the tropical region, the combination of enhanced horizontal momentum advection by the residual meridional circulation and inertially unstable circulations related to planetary wave breaking in the subtropics produced conditions favoring barotropic instability. In the extratropical region, the enhanced residual meridional circulation altered the zonal wind tendency through increased Coriolis torque.

How momentum advection schemes influence current-topography interactions at eddy permitting resolution

Recent studies have shown that the use of an enstrophy-and-energy-conserving momentum advection scheme substantially reduces widespread biases of mean currents in the global 1/4° DRAKKAR model. This paper investigates the origin of these improvements. A series of sensitivity simulations with different momentum advection schemes is performed with the North Atlantic 1/4° DRAKKAR model. Three second order momentum advection schemes conserving, respectively, enstrophy ( ens), energy ( efx) and both quantities ( een) are tested and their impact on the model solution are compared. The mean kinetic energy vertical profile is found to change up to 10% depending on the chosen scheme. This sensitivity is maximum in bottom layers. The analysis of the vorticity tendency due to horizontal momentum advection reveals that the three schemes differ mostly in bottom layers as well. The average magnitude of this term is enhanced with the efx scheme and reduced with the een scheme. These differences are found to be consistent with the instantaneous tendency of each scheme. In addition, we show that the differences between the schemes are related to the grid-scale irregularity of the velocity field. Both the grid scale irregularity and the differences between the schemes are found to be enhanced in bottom layers. We conclude that the model solution depends crucially on the ability of the momentum advection scheme to handle under-resolved flows close to the bottom topography. This work emphasizes the critical influence of topography in eddy-active regions on mean circulation features such as the position of the North-Atlantic current or the Gulf Stream separation.

An analysis of the momentum forcing in the high‐latitude lower thermosphere

We analyze the dynamics of the high‐latitude thermsopheric wind system below 170 km for negative IMF Bz by using a fully nonlinear model with a realistic distribution of the forcing. A transition of the forcing patterns and their relative contribution to the high‐latitude lower thermospheric wind system occurs around 123 km under various conditions, weak or strong IMF, summer or winter. Winds around and above 123 km are sustained by the gradient‐wind balance among divergent/convergent pressure gradient, Coriolis, and horizontal momentum advection (mainly centrifugal) accelerations. Below 123 km winds are maintained by the approximate balance of divergent/convergent pressure gradient, Coriolis, and Hall ion drag accelerations through modified geostrophy. The dominant contribution to the wind tendency (time rate of change) is the rotational component of the ion drag acceleration. The wind tendency above 123 km tends to resemble rotational Pedersen ion drag acceleration well, which reflects a rotated pattern of the E × B velocity. Near and below 123 km the wind tendency is also affected by the rotational component of the Hall ion drag acceleration whose pattern no longer closely resembles the pattern of the E × B velocity, and the wind pattern can differ significantly from that well above 123 km. Simulations for different strengths of the IMF and different seasons indicate that largely divergent/convergent Coriolis and horizontal momentum advection accelerations tend approximately to balance with the horizontal pressure gradient (as well as with divergent/convergent ion drag at lower altitudes) under various conditions. As the forcing increases the radius of curvature of the strong winds also tends to increase, so that the centrifugal acceleration does not increase quadratically with the maximum wind speed, and the tendency for a rough balance between the Coriolis and horizontal momentum advection accelerations in the duskside vortex above 123 km is maintained.

Nonlinear interaction between wind‐forced currents and near‐inertial oscillations and internal waves in nearshore frontal regions

A cross‐sectional model with a coastal boundary and bottom front is used to examine the nonlinear interaction between wind‐forced currents and internal waves at the near‐inertial frequency (f) and those at super‐ and subinertial frequencies (ωf). In the frontal region, nonlinear effects associated with vorticity in the frontal jet and differences in eddy viscosity gives rise to Ekman pumping. This allows the wind's momentum at super‐ and subinertial frequencies to penetrate to a greater depth than that due to momentum diffusion. Internal waves produced in the frontal region are trapped on the western side of the front but can propagate on the eastern side, although some trapping occurs at depth. Away from the frontal region, nonlinear interaction gives rise to energy at the frequencies (ωf + f) and (∣ωf − f∣) primarily through vertical shear associated with the f frequency and vertical velocity at the ωf frequency. In the frontal region, where internal waves at the ωf and f frequency are trapped at depth by the sloping isotherms, the horizontal momentum advection term is comparable to the nonlinear term involving vertical velocity, and there is maximum nonlinear interaction. Internal waves at the (ωf + f) frequency are superinertial and can propagate away from their generation region, while those at (∣ωf − f∣) are subinertial. Calculations show appreciable spatial variability in the interaction processes, which will influence the spectra derived from observations. The spatial characteristics of the interaction change in the frontal region depending upon ωf being super‐ or subinertial.

The diurnal wind variation in a variable eddy viscosity semi-geostrophic Ekman boundary-layer model: Analytical study

¶A time-dependent semi-geostrophic Ekman boundary-layer model (SG), including slowly varying eddy diffusivity with height and inertial term effects, is developed to investigate the diurnal wind variation in the planetary boundary layer (PBL). An approximate analytical solution of this model is derived by using the WKB method, which extends the Tan and Farahani (1998)’s solution by including the vertical variable eddy viscosity. The features of the diurnal wind variation in the PBL mainly depend on three factors: the latitude, horizontal momentum advection and eddy viscosity. The vertical variable eddy viscosity has little influence on diurnal wind variation in the PBL at the low latitude, however its effect may be exacerbated in the mid- and high latitudes. In comparing with the constant eddy viscosity case, the decreasing (increasing) with height eddy viscosity produces a large (small) maximum wind speed (MWS) in the PBL, however, the eddy viscosity that has a mid-layer peak in the vertical gives rise to a higher height of occurrence of MWS. For the boundary-layer wind structure, there is a singular point when the modified SG inertial oscillation frequency η equals the forcing frequency ω. The isotachs of boundary-layer wind speed have almost no tilt to left or right relative to time evolution and the occurrence time of the MWS is the earliest at the singular point. The feature will be enhanced in the decreasing with height eddy viscosity case and weakened in the eddy viscosity initially increasing with height case.

Numerical simulation on mesoscale convective system along Mei-Yu front in Southern China

A mesoscale convective system (MCS) was generated over the South China Sea at around 115° E, 21°N on 0000 UTC 7 June 1998. Riding along the Mei-Yu front, the system moved through the Taiwan Strait and finally hit the southern part of the Taiwan Island, produced over 300 mm of rainfall over parts of Taiwan in the next 24 hours. It was found that the Penn State-NCAR Mesoscale Model Version 5 (MM5) did quite well in simulating the evolution of the MCS. Diagnostic studies on the mesoscale structure of MCS, moisture concentration process and the momentum budget were based on the model output. The concept model of the MCS in the mature phase can be concluded as the following: At the 850 hPa level, the mesoscale low level jet (mLLJ) was found to the southwest of the MCS, which was also associated with a mesoscale low generated by convection, the mesoscale upper level jet (mULJ) was found due east of the MCS (and an upper-level mesoscale high) at the 300 hPa level. Horizontal advection of momentum acted to accelerate the exit of the mLLJ but decelerate the entrance of the mLLJ where the pressure gradient force acted to compensate this lost and maintain its strength. For maintenance of the mULJ, vertical advection of momentum and the pressure gradient force were both found to be important.

Generalized Inverse of a Reduced Gravity Primitive Equation Ocean Model and Tropical Atmosphere–Ocean Data

A nonlinear 2‰-layer reduced gravity primitive equations (PE) ocean model is used to assimilate sea surface temperature (SST) data from the Tropical Atmosphere‐Ocean (TAO) moored buoys in the tropical Pacific. The aim of this project is to hindcast cool and warm events of this part of the ocean, on seasonal to interannual timescales. The work extends that of Bennett et al., who used a modified Zebiak‐Cane coupled model. They were able to fit a year of 30-day averaged TAO data to within measurement errors, albeit with significant initial and dynamical residuals. They assumed a 100-day decorrelation timescale for the dynamical residuals. This long timescale for the residuals reflects the neglect of resolvable processes in the intermediate coupled model, such as horizontal advection of momentum. However, the residuals in the nonlinear PE model should be relatively short timescale errors in parameterizations. The scales for these residuals are crudely estimated from the upper ocean turbulence studies of Peters et al. and Moum. The assimilation is performed by minimizing a weighted least squares functional expressing the misfits to the data and to the model throughout the tropical Pacific and for 18 months. It is known that the minimum lies in the ‘‘data subspace’’ of the state or solution space. The minimum is therefore sought in the data subspace, by using the representer method to solve the Euler‐Lagrange (EL) system. Although the vector space decomposition and solution method assume a linear EL system, the concept and technique are applied to the nonlinear EL system (resulting from the nonlinear PE model), by iterating with linear approximations to the nonlinear EL system. As a first step, the authors verify that sequences of solutions of linear iterates of the forward PE model do converge. The assimilation is also used as a significance test of the hypothesized means and covariances of the errors in the initial conditions, dynamics, and data. A ‘‘strong constraint’’ inverse solution is computed. However, it is outperformed by the ‘‘weak constraint’’ inverse. A cross validation by withheld data is presented, as well as an inversion with the model forced by the Florida State University winds, in place of a climatological wind forcing used in the former inversions.

Surface Winds at Landfall of Hurricane Andrew (1992)—A Reply

and ocean, separately, and then merging them along the coastline without including full dynamics and physics interaction,’’ as mentioned in LZY. In this reply, we will show that their converting the flight-level (at an altitude of 2.5‐3.2 km) data with a simple, one-dimensional planetary boundary layer (PBL) model tended to introduce large errors in surface winds, while their merging, without considering a transition zone produced by horizontal momentum advection, accounts for the unrealistic sharp discontinuity. These problems could be clear

A coastal ocean model intercomparison study for a three-dimensional idealised test case

Several coastal ocean models have been used to compute the circulation on the Northwest European Continental Shelf. Five of them, developed within the European Union, are compared in the scope of an idealised three-dimensional test case, dealing with the geostrophic adjustment of a freshwater cylinder. As the central eddy adjusts, unstable baroclinic vortices start to grow. All the models are able to produce such unstable vortices. However, two of them produce an order-two instability, which is in accordance with a previous laboratory experiment, while the others exhibit an order-four instability. Using a simple scaling analysis, it is seen that the azimuthal wavenumber depends on the ratio of the kinetic energy to the available potential energy. It appears that the discrepancy in the azimuthal wavenumber is mainly due to the effect of the discretisation of the horizontal advection of momentum which could produce significant decrease of the total kinetic energy.

An Analytical Study of the Diurnal Variations of Wind in a Semi-geostrophic Ekman Boundary Layer Model

A time-dependent semi-geostrophic Ekman boundary-layer model based on the geostrophic momentum approximation is used to study the diurnal wind variation in the planetary boundary layer (PBL) and the evolution of the low-level nocturnal jet (LLJ). The coefficient of eddy viscosity varies periodically with time, varies linearly with height in the surface layer and is constant above the surface layer. The influence of horizontal advection of momentum on the diurnal wind variation in the PBL, the development of inertial oscillations (IOs) and the formation of the LLJ are examined. In comparison with the Ekman solutions, the diurnal wind variation in semi-geostrophic Ekman boundary-layer dynamics has the following features: (1) the phase angle of the diurnal wind wave shifts with height, the rate of shifting is increased in anticyclonic regions and decreased in cyclonic regions, (2) the time of occurrence of the low-level maximum wind speed is later in anticyclonic regions and earlier in cyclonic regions, (3) the height of occurrence of the maximum wind speed is higher in the anticyclonic and lower in cyclonic regions, (4) the wind speed maximum and the amplitude of the diurnal wind variation are larger in anticyclonic and smaller in cyclonic regions, (5) the period of IOs is larger in anticyclonic regions and smaller in cyclonic regions, (6) anticyclonic vorticity is conducive to the generation of LLJ in the PBL. These features are interpreted by means of the physical properties of semi-geostrophic Ekman boundary-layer dynamics and inertial oscillation dynamics.

Mesoscale Momentum Budget in a Midlatitude Squall Line: A Numerical Case Study

Abstract A mesoβ-scale momentum budget and its effect on larger-scale mean flow in a midlatitude mesoscale convective system are investigated using a numerical simulation of an intense squall line that occurred during 10–11 June 1985 PRE-STORM. It is found that the momentum generation normal to the line associated with the latent heating and cooling contributes most significantly to the momentum budget and determines the mesoβ-scale internal structure and evolution of the squall line. The momentum generation along the line contributes to the initial development of a mesovortex but has little effect on the final vertical structure of the along-line flow. Both vertical and horizontal momentum advection have significant contributions, particularly to the vertical mixing of the along-line flow; and this component of the horizontal momentum is locally transported down-gradient. It is also found that for midlatitude convective systems, convectively generated downdrafts can play as prominent a role as updrafts i...