The dependence of southern ocean meridional overturning on wind stress

An eddy-resolving numerical model of a zonal flow, meant to resemble the Antarctic Circumpolar Current, is described and analyzed using the framework of J. Marshall and T. Radko. In addition to wind and buoyancy forcing at the surface, the model contains a sponge layer at the northern boundary that permits a residual meridional overturning circulation (MOC) to exist at depth. The strength of the residual MOC is diagnosed for different strengths of surface wind stress. It is found that the eddy circulation largely compensates for the changes in Ekman circulation. The extent of the compensation and thus the sensitivity of the MOC to the winds depend on the surface boundary condition. A fixed-heat-flux surface boundary severely limits the ability of the MOC to change. An interactive heat flux leads to greater sensitivity. To explain the MOC sensitivity to the wind strength under the interactive heat flux, transformed Eulerian-mean theory is applied, in which the eddy diffusivity plays a central role in determining the eddy response. A scaling theory for the eddy diffusivity, based on the mechanical energy balance, is developed and tested; the average magnitude of the diffusivity is found to be proportional to the square root of the wind stress. The MOC sensitivity to the winds based on this scaling is compared with the true sensitivity diagnosed from the experiments.

Natural ventilation assessment in typical open an semi-open urban environments under various wind directions

Semi-open street roofs protect pedestrians from intense sunshine and rains. Their effects on natural ventilation of urban canopy layers (UCL) are less understood. This paper investigates two idealized urban models consisting of 4(2×2) or 16(4×4) buildings under a neutral atmospheric condition with parallel (0°) or non-parallel (15°,30°,45°) approaching wind. The aspect ratio (building height (H) / street width (W)) is 1 and building width is B=3H. Computational fluid dynamic (CFD) simulations were first validated by experimental data, confirming that standard k-ε model predicted airflow velocity better than RNG k-ε model, realizable k–ε model and Reynolds stress model. Three ventilation indices were numerically analyzed for ventilation assessment, including flow rates across street roofs and openings to show the mechanisms of air exchange, age of air to display how long external air reaches a place after entering UCL, and purging flow rate to quantify the net UCL ventilation capacity induced by mean flows and turbulence. Five semi-open roof types are studied: Walls being hung above street roofs (coverage ratio λa=100%) at z=1.5H, 1.2H, 1.1H ('Hung1.5H', 'Hung1.2H', 'Hung1.1H' types); Walls partly covering street roofs (λa=80%) at z=H ('Partly-covered' type); Walls fully covering street roofs (λa=100%) at z=H ('Fully-covered' type).They basically obtain worse UCL ventilation than open street roof type due to the decreased roof ventilation. 'Hung1.1H', 'Hung1.2H', 'Hung1.5H' types are better designs than 'Fully-covered' and 'Partly-covered' types. Greater urban size contains larger UCL volume and requires longer time to ventilate. The methodologies and ventilation indices are confirmed effective to quantify UCL ventilation.

How to parametrize urban-canopy drag to reproduce wind-direction effects within the canopy

The mean wind direction within an urban canopy changes with height when the incoming flow is not orthogonal to obstacle faces. This wind-turning effect is induced by complex processes and its modelling in urban-canopy (UC) parametrizations is difficult. Here we focus on the analysis of the spatially-averaged flow properties over an aligned array of cubes and their variation with incoming wind direction. For this purpose, Reynolds-averaged Navier–Stokes simulations previously compared, for a reduced number of incident wind directions, against direct numerical simulation results are used. The drag formulation of a UCparametrization ismodified and different drag coefficients are tested in order to reproduce the wind-turning effect within the canopy for oblique wind directions. The simulations carried out for a UC parametrization in one-dimensional mode indicate that a height-dependent drag coefficient is needed to capture this effect.

The influence of ocean flow on newly forming sea ice

The heat and mass balance of the Arctic Ocean is very sensitive to the growth and decay of sea ice and the interaction between the heat and salt fields in the oceanic boundary layer. The hydraulic roughness of sea ice controls the detailed nature of turbulent fluxes in the boundary layer and hence is an important ingredient in model parameterizations. We describe a novel mechanism for the generation of corrugations of the sea ice–ocean interface, present a mathematical analysis elucidating the mechanism, and present numerical calculations for geophysically relevant conditions. The mechanism relies on brine flows developing in the sea ice due to Bernoulli suction by flow of ocean past the interface. For oceanic shears at the ice interface of 0.2 s−1, we expect the corrugations to form with a wavelength dependent upon the permeability structure of the sea ice which is described herein. The mechanism should be particularly important during sea ice formation in wind-maintained coastal polynyas and in leads. This paper applies our earlier analyses of the fundamental instability to field conditions and extends it to take account of the anisotropic and heterogeneous permeability of sea ice.

Flow structure and near-field dispersion in arrays of building-like obstacles

Dispersion in the near-field region of localised releases in urban areas is difficult to predict because of the strong influence of individual buildings. Effects include upstream dispersion, trapping of material into building wakes and enhanced concentration fluctuations. As a result, concentration patterns are highly variable in time and mean profiles in the near field are strongly non-Gaussian. These aspects of near-field dispersion are documented by analysing data from direct numerical simulations in arrays of building-like obstacles and are related to the underlying flow structure. The mean flow structure around the buildings is found to exert a strong influence over the dispersion of material in the near field. Diverging streamlines around buildings enhance lateral dispersion. Entrainment of material into building wakes in the very near field gives rise to secondary sources, which then affect the subsequent dispersion pattern. High levels of concentration fluctuations are also found in this very near field; the fluctuation intensity is of order 2 to 5.

Solar cycle evolution of dipolar and pseudostreamer belts and their relation to the slow solar wind

Dipolar streamers are coronal structures formed by open solar flux converging from coronal holes of opposite polarity. Thus the dipolar streamer belt traces the coronal foot print of the heliospheric current sheet (HCS), and it is strongly associated with the origin of slow solar wind. Pseudostreamers, on the other hand, separate converging regions of open solar flux from coronal holes of the same polarity and do not contain current sheets. They have recently received a great deal of interest as a possible additional source of slow solar wind. Here we add to that growing body of work by using the potential-field source-surface model to determine the occurrence and location of dipolar and pseudostreamers over the last three solar cycles. In addition to providing new information about pseudostreamer morphology, the results help explain why the observations taken during the first Ulysses perihelion pass in 1995 showed noncoincidence between dipolar streamer belt and the locus of slowest flow. We find that Carrington rotation averages of the heliographic latitudes of dipolar and pseudostreamer belts are systematically shifted away from the equator, alternately in opposite directions, with a weak solar cycle periodicity, thus keeping slow wind from the web of combined streamer belts approximately symmetric about the equator. The largest separation of dipolar and pseudostreamer belts occurred close to the Ulysses pass, allowing a unique opportunity to see that slow wind from pseudostreamer belts north of the southward-displaced dipolar belt was responsible for the noncoincident pattern.

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Quantitative ventilation assessments of idealized urban canopy layers with various urban layouts and the same building packing density

This paper investigates urban canopy layers (UCL) ventilation under neutral atmospheric condition with the same building area density (λp=0.25) and frontal area density (λf=0.25) but various urban sizes, building height variations, overall urban forms and wind directions. Turbulent airflows are first predicted by CFD simulations with standard k-ε model evaluated by wind tunnel data. Then air change rates per hour (ACH) and canopy purging flow rate (PFR) are numerically analyzed to quantify the rate of air exchange and the net ventilation capacity induced by mean flows and turbulence. With a parallel approaching wind (θ=0o), the velocity ratio first decreases in the adjustment region, followed by the fully-developed region where the flow reaches a balance. Although the flow quantities macroscopically keep constant, however ACH decreases and overall UCL ventilation becomes worse if urban size rises from 390m to 5km. Theoretically if urban size is infinite, ACH may reach a minimum value depending on local roof ventilation, and it rises from 1.7 to 7.5 if the standard deviation of building height variations increases (0% to 83.3%). Overall UCL ventilation capacity (PFR) with a square overall urban form (Lx=Ly=390m) is better as θ=0o than oblique winds (θ=15o, 30o, 45o), and it exceeds that of a staggered urban form under all wind directions (θ=0o to 45o), but is less than that of a rectangular urban form (Lx=570m, Ly=270m) under most wind directions (θ=30o to 90o). Further investigations are still required to quantify the net ventilation efficiency induced by mean flows and turbulence.

Wind pressure on a single building immersed in a low-jet wind profile

This paper for the first time discuss the wind pressure distribution on the building surface immersed in wind profile of low-level jet rather than a logarithmic boundary-layer profile. Two types of building models are considered, low-rise and high-rise building, relative to the low-level jet height. CFD simulation is carried out. The simulation results show that the wind pressure distribution immersed in a low-jet wine profile is very different from the typical uniform and boundary-layer flow. For the low-rise building, the stagnation point is located at the upper level of windward façade for the low-level jet wind case, and the separation zone above the roof top is not as obvious as the uniform case. For the high-rise building model, the height of stagnation point is almost as high as the low-level jet height.

Impact of non-hydrostatic effects and trapped lee waves on mountain wave drag in directionally sheared flow

The orographic gravity wave drag produced in flow over an axisymmetric mountain when both vertical wind shear and non-hydrostatic effects are important was calculated using a semi-analytical two-layer linear model, including unidirectional or directional constant wind shear in a layer near the surface, above which the wind is constant. The drag behaviour is determined by partial wave reflection at the shear discontinuity, wave absorption at critical levels (both of which exist in hydrostatic flow), and total wave reflection at levels where the waves become evanescent (an intrinsically non-hydrostatic effect), which produces resonant trapped lee wave modes. As a result of constructive or destructive wave interference, the drag oscillates with the thickness of the constant-shear layer and the Richardson number within it (Ri), generally decreasing at low Ri and when the flow is strongly non-hydrostatic. Critical level absorption, which increases with the angle spanned by the wind velocity in the constant-shear layer, shields the surface from reflected waves, keeping the drag closer to its hydrostatic limit. While, for the parameter range considered here, the drag seldom exceeds this limit, a substantial drag fraction may be produced by trapped lee waves, particularly when the flow is strongly non-hydrostatic, the lower layer is thick and Ri is relatively high. In directionally sheared flows with Ri = O(1), the drag may be misaligned with the surface wind in a direction opposite to the shear, a behaviour which is totally due to non-trapped waves. The trapped lee wave drag, whose reaction force on the atmosphere is felt at low levels, may therefore have a distinctly different direction from the drag associated with vertically propagating waves, which acts on the atmosphere at higher levels.

Statistical-dynamical downscaling for wind energy potentials: evaluation and applications to decadal hindcasts and climate change projections

A statistical–dynamical downscaling (SDD) approach for the regionalization of wind energy output (Eout) over Europe with special focus on Germany is proposed. SDD uses an extended circulation weather type (CWT) analysis on global daily mean sea level pressure fields with the central point being located over Germany. Seventy-seven weather classes based on the associated CWT and the intensity of the geostrophic flow are identified. Representatives of these classes are dynamically downscaled with the regional climate model COSMO-CLM. By using weather class frequencies of different data sets, the simulated representatives are recombined to probability density functions (PDFs) of near-surface wind speed and finally to Eout of a sample wind turbine for present and future climate. This is performed for reanalysis, decadal hindcasts and long-term future projections. For evaluation purposes, results of SDD are compared to wind observations and to simulated Eout of purely dynamical downscaling (DD) methods. For the present climate, SDD is able to simulate realistic PDFs of 10-m wind speed for most stations in Germany. The resulting spatial Eout patterns are similar to DD-simulated Eout. In terms of decadal hindcasts, results of SDD are similar to DD-simulated Eout over Germany, Poland, Czech Republic, and Benelux, for which high correlations between annual Eout time series of SDD and DD are detected for selected hindcasts. Lower correlation is found for other European countries. It is demonstrated that SDD can be used to downscale the full ensemble of the Earth System Model of the Max Planck Institute (MPI-ESM) decadal prediction system. Long-term climate change projections in Special Report on Emission Scenarios of ECHAM5/MPI-OM as obtained by SDD agree well to the results of other studies using DD methods, with increasing Eout over northern Europe and a negative trend over southern Europe. Despite some biases, it is concluded that SDD is an adequate tool to assess regional wind energy changes in large model ensembles.

A wind tunnel model for quantifying fluxes in the urban boundary layer

Transport of pollution and heatout of streets into the boundary layer above is not currently understood and so fluxes cannot be quantified. Scalar concentration within the street is determined by the flux out of it and so quantifying fluxes for turbulent flow over a rough urban surface is essential. We have developed a naphthalene sublimation technique to measure transfer from a two-dimensional street canyon in a wind tunnel for the case of flow perpendicular to the street. The street was coated with naphthalene, which sublimes at room temperature, so that the vapour represented the scalar source. The transfer velocity wT relates the flux out of the canyon to the concentration within it and is shown to be linearly related to windspeed above the street. The dimensionless transfer coefficient wT/Uδ represents the ventilation efficiency of the canyon (here, wT is a transfer velocity,Uδ is the wind speed at the boundary-layer top). Observed values are between 1.5 and 2.7 ×10-3 and, for the case where H/W→0 (ratio of buildingheight to street width), values are in the same range as estimates of transfer from a flat plate, giving confidence that the technique yields accurate values for street canyon scalar transfer. wT/Uδ varies with aspect ratio (H/W), reaching a maximum in the wake interference regime (0.3 < H/W < 0.65). However, when upstream roughness is increased, the maximum in wT/Uδ reduces, suggesting that street ventilation is less sensitive to H/W when the flow is in equilibrium with the urban surface. The results suggest that using naphthalene sublimation with wind-tunnel models of urban surfaces can provide a direct measure of area-averaged scalar fluxes.

Observations of urban boundary layer structure during a strong urban heat island event

It has long been known that the urban surface energy balance is different to that of a rural surface, and that heating of the urban surface after sunset gives rise to the Urban Heat Island (UHI). Less well known is how flow and turbulence structure above the urban surface are changed during different phases of the urban boundary layer (UBL). This paper presents new observations above both an urban and rural surface and investigates how much UBL structure deviates from classical behaviour. A 5-day, low wind, cloudless, high pressure period over London, UK, was chosen for analysis, during which there was a strong UHI. Boundary layer evolution for both sites was determined by the diurnal cycle in sensible heat flux, with an extended decay period of approximately 4 h for the convective UBL. This is referred to as the “Urban Convective Island” as the surrounding rural area was already stable at this time. Mixing height magnitude depended on the combination of regional temperature profiles and surface temperature. Given the daytime UHI intensity of 1.5∘C, combined with multiple inversions in the temperature profile, urban and rural mixing heights underwent opposite trends over the period, resulting in a factor of three height difference by the fifth day. Nocturnal jets undergoing inertial oscillations were observed aloft in the urban wind profile as soon as the rural boundary layer became stable: clear jet maxima over the urban surface only emerged once the UBL had become stable. This was due to mixing during the Urban Convective Island reducing shear. Analysis of turbulent moments (variance, skewness and kurtosis) showed “upside-down” boundary layer characteristics on some mornings during initial rapid growth of the convective UBL. During the “Urban Convective Island” phase, turbulence structure still resembled a classical convective boundary layer but with some influence from shear aloft, depending on jet strength. These results demonstrate that appropriate choice of Doppler lidar scan patterns can give detailed profiles of UBL flow. Insights drawn from the observations have implications for accuracy of boundary conditions when simulating urban flow and dispersion, as the UBL is clearly the result of processes driven not only by local surface conditions but also regional atmospheric structure.

Secondary cyclogenesis along an occluded front leading to damaging wind gusts: windstorm Kyrill, January 2007

Windstorm Kyrill affected large parts of Europe in January 2007 and caused widespread havoc and loss of life. In this study the formation of a secondary cyclone, Kyill II, along the occluded front of the mature cyclone Kyrill and the occurrence of severe wind gusts as Kyrill II passed over Germany are investigated with the help of high-resolution regional climate model simulations. Kyrill underwent an explosive cyclogenesis south of Greenland as the storm crossed polewards of an intense upper-level jet stream. Later in its life cycle secondary cyclogenesis occurred just west of the British Isles. The formation of Kyrill II along the occluded front was associated (a) with frontolytic strain and (b) with strong diabatic heating in combination with a developing upper-level shortwave trough. Sensitivity studies with reduced latent heat release feature a similar development but a weaker secondary cyclone, revealing the importance of diabatic processes during the formation of Kyrill II. Kyrill II moved further towards Europe and its development was favored by a split jet structure aloft, which maintained the cyclone’s exceptionally deep core pressure (below 965 hPa) for at least 36 hours. The occurrence of hurricane force winds related to the strong cold front over North and Central Germany is analyzed using convection-permitting simulations. The lower troposphere exhibits conditional instability, a turbulent flow and evaporative cooling. Simulation at high spatio-temporal resolution suggests that the downward mixing of high momentum (the wind speed at 875 hPa widely exceeded 45 m s-1) accounts for widespread severe surface wind gusts, which is in agreement with observed widespread losses.

Dawn-dusk asymmetry in particles of solar wind origin within the magnetosphere

Solar wind/magnetosheath plasma in the magnetosphere can be identified using a component that has a higher charge state, lower density and, at least soon after their entry into the magnetosphere, lower energy than plasma from a terrestrial source. We survey here observations taken over 3 years of He2+ ions made by the Magnetospheric Ion Composition Sensor (MICS) of the Charge and Mass Mgnetospheric Ion Composition Experiment (CAMMICE) instrument aboard POLAR. The occurrence probability of these solar wind ions is then plotted as a function of Magnetic Local Time (MLT) and invariant latitude (3) for various energy ranges. For all energies observed by MICS (1.8–21.4 keV) and all solar wind conditions, the occurrence probabilities peaked around the cusp region and along the dawn flank. The solar wind conditions were filtered to see if this dawnward asymmetry is controlled by the Svalgaard-Mansurov effect (and so depends on the BY component of the interplanetary magnetic field, IMF) or by Fermi acceleration of He2+ at the bow shock (and so depends on the IMF ratio BX/BY ). It is shown that the asymmetry remained persistently on the dawn flank, suggesting it was not due to effects associated with direct entry into the magnetosphere. This asymmetry, with enhanced fluxes on the dawn flank, persisted for lower energy ions (below a “cross-over” energy of about 23 keV) but reversed sense to give higher fluxes on the dusk flank at higher energies. This can be explained by the competing effects of gradient/curvature drifts and the convection electric field on ions that are convecting sunward on re-closed field lines. The lower-energy He2+ ions E × B drift dawnwards as they move earthward, whereas the higher energy ions curvature/gradient drift towards dusk. The convection electric field in the tail is weaker for northward IMF. Ions then need less energy to drift to the dusk flank, so that the cross-over energy, at which the asymmetry changes sense, is reduced.