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.

Plasma structure within poleward-moving cusp/cleft auroral transients: EISCAT Svalbard radar observations and an explanation in terms of large local time extent of events

We report high-resolution observations of the southward-IMF cusp/cleft ionosphere made on December 16th 1998 by the EISCAT (European incoherent scatter) Svalbard radar (ESR), and compare them with observations of dayside auroral luminosity, as seen at a wavelength of 630 nm by a meridian scanning photometer at Ny Alesund, and of plasma flows, as seen by the CUTLASS (co-operative UK twin location auroral sounding system) Finland HF radar. The optical data reveal a series of poleward-moving transient red-line (630 nm) enhancements, events that have been associated with bursts in the rate of magnetopause reconnection generating new open flux. The combined observations at this time have strong similarities to predictions of the effects of soft electron precipitation modulated by pulsed reconnection, as made by Davis and Lockwood (1996); however, the effects of rapid zonal flow in the ionosphere, caused by the magnetic curvature force on the newly opened field lines, are found to be a significant additional factor. In particular, it is shown how enhanced plasma loss rates induced by the rapid convection can explain two outstanding anomalies of the 630 nm transients, namely how minima in luminosity form between the poleward-moving events and how events can re-brighten as they move poleward. The observations show how cusp/cleft aurora and transient poleward-moving auroral forms appear in the ESR data and the conditions which cause enhanced 630 nm emission in the transients: they are an important first step in enabling the ESR to identify these features away from the winter solstice when supporting auroral observations are not available.

Relationship of dayside auroral precipitations to the open-closed separatrix and the pattern of convective flow

The implications are discussed of acceleration of magnetospheric ions by reflection off two magnetopause Alfvén waves, launched by the reconnection site into the inflow regions on both sides of the boundary. The effects of these waves on the ion populations, predicted using the model described by Lockwood et al. [1996], offer a physical interpretation of all the various widely used classifications of precipitation into the dayside ionosphere, namely, central plasma sheet, dayside boundary plasma sheet (BPS), void, low-latitude boundary layer (LLBL), cusp, mantle, and polar cap. The location of the open-closed boundary and the form of the convection flow pattern are discussed in relation to the regions in which these various precipitations are typically found. Specifically, the model predicts that both the LLBL and the dayside BPS precipitations are on newly opened field lines and places the convection reversal within the LLBL, as is often observed. It is shown that this offers solutions to a number of paradoxes and problems that arise if the LLBL and BPS precipitations are thought of as being on closed field lines. This model is also used to make quantitive predictions of the longitudinal extent and latitudinal width of the cusp, as a function of solar wind density.

Time-dependent flows in the coupled solar wind-magnetosphere-ionosphere system

Concepts of time-dependent flow in the coupled solar wind-magnetosphere-ionosphere system are discussed and compared with the frequently-adopted steady-state paradigm. Flows are viewed as resulting from departures of the system from equilibrium excited by dayside and nightside reconnection processes, with the flows then taking the system back towards a new equilibrium configuration. The response of the system to reconnection impulses, continuous but unbalanced reconnection and balanced steady-state reconnection are discussed in these terms. It is emphasized that in the time-dependent case the ionospheric and interplanetary electric fields are generally inductively decoupled from each other; a simple mapping of the interplanetary electric field along equipotential field lines into the ionosphere occurs only in the electrostatic steady-state case.

Dynamic auroral structure in the vicinity of the polar cusp - multipoint observations during southward and northward IMF

The low- and high-latitude boundary layers of the earth's magnetosphere [low-latitude boundary layer (LLBL) and mantle] play important roles in transferring momentum and energy from the solar wind to the magnetosphere-ionosphere system. Particle precipitation, field-aligned current, auroral emission, ionospheric ion drift and ground magnetic perturbations are among the low-altitude parameters that show signatures of various plasma processes in the LLBL and the magnetopause current layer. Magnetic merging events, Kelvin-Helmholtz waves, and pressure pulses excited by the variable solar wind/magnetosheath plasma are examples of boundary phenomena that may be coupled to the ionosphere via field-aligned currents. Optical auroral observation, by photometry and all-sky TV cameras, is a unique technique for investigating the spatial and temporal structure of the electron precipitation associated with such phenomena. However, the distinction between the different boundary layer plasma populations cannot in general be unambiguously determined by optics alone. Additional information, such as satellite observations of particle boundaries and field-aligned currents, is needed in order to identify the plasma source(s) and the magnetosphere-ionosphere coupling mode(s). Two categories of auroral activity/structure in the vicinity of the polar cusp are discussed in this paper, based on combined ground and satellite data. In one case, the quasi-periodic sequence of auroral events at the polar cap boundary involves accelerated electrons (< 1 keV) moving poleward (< 1 km s-1) and azimuthally along the persistent cusp/cleft arc poleward boundary with velocities (< 4 km s-1), comparable to the local ionospheric ion drift during periods of southward IMF. A critical question is whether or not the optical events signify a corresponding plasma flow across the open/closed field line boundary in such cases. Near-simultaneous observations of magnetopause flux transfer events (FTEs) and such optical/ion drift events are reported. The reverse pattern of motion of discrete auroral forms is observed during positive interplanetary magnetic field (IMF) B(Z), i.e. equatorward motion into the cusp/cleft background arc from the poleward edge. Combined satellite and ground-based information for the latter cases indicate a source mechanism, poleward of the cusp at the high-latitude magnetopause or plasma mantle, giving rise to strong momentum transfer and electron precipitation structures within a approximately 200 km-wide latitudinal zone at the cusp/cleft poleward boundary. The striking similarities of auroral electrodynamics in the cleft/mantle region during northward and southward IMF indicate that a qualitatively similar solar wind-magnetosphere coupling mode is operating. It is suggested that, in both cases, the discrete auroral forms represent temporal/spatial structure of larger-scale convection over the polar magnetosphere.