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.

Ionospheric ion and electron heating at the poleward boundary of a poleward expanding substorm-disturbed region

We present observations of a poleward propagating substorm-disturbed region which was observed by the European Incoherent SCATter (EISCAT) radar and the Svalbard International Monitor for Auroral Geomagnetic Effects (IMAGE) magnetometers in the postmidnight sector. The expansion of the disturbance was launched by a substorm intensification which started similar to 25 min after the initial onset, and similar to 10 min before the disturbance arrived over Svalbard. In association with the magnetic disturbance, a poleward expanding enduring enhancement in the F region electron temperature was observed, indicative of soft electron precipitation, with a narrow band of enhanced ion temperature straddling its poleward edge, indicative of fast ion flows and ion-neutral collisional heating. This electron temperature boundary was coincident with the poleward propagating electrojet current system detected by the high-latitude IMAGE magnetometer stations and is taken to be a proxy for the observation of a substorm auroral bulge. The electron temperature boundary is inferred to have a width comparable or less than one radar range gate (similar to 60 km transverse to the magnetic field), while the region of high ion temperature was found to be approximately three gates wide, extending approximately two gates (similar to 120 km) poleward of the electron temperature boundary, and approximately one gate (similar to 60 km) equatorward. The two-beam radar line-of-sight velocity data are found to be consistent with the existence of a layer of high-speed flow in the boundary, peaking at values similar to1.5-3 km s(-1), roughly consistent with the ion temperature data. The flow is directed either east or west along the boundary depending on the direction of the flow in the poleward region. We infer that the flow is deflected along and around the boundary of the substorm-disturbed region due to the high conductivity of the latter. Variations in the flow poleward of the boundary produced no discernible magnetic effects on the ground, confirming the low conductivity of the preboundary ionosphere.

Meteorological forcing data for urban outdoor thermal comfort models from a coupled convective boundary layer and surface energy balance scheme

Site-specific meteorological forcing appropriate for applications such as urban outdoor thermal comfort simulations can be obtained using a newly coupled scheme that combines a simple slab convective boundary layer (CBL) model and urban land surface model (ULSM) (here two ULSMs are considered). The former simulates daytime CBL height, air temperature and humidity, and the latter estimates urban surface energy and water balance fluxes accounting for changes in land surface cover. The coupled models are tested at a suburban site and two rural sites, one irrigated and one unirrigated grass, in Sacramento, U.S.A. All the variables modelled compare well to measurements (e.g. coefficient of determination = 0.97 and root mean square error = 1.5 °C for air temperature). The current version is applicable to daytime conditions and needs initial state conditions for the CBL model in the appropriate range to obtain the required performance. The coupled model allows routine observations from distant sites (e.g. rural, airport) to be used to predict air temperature and relative humidity in an urban area of interest. This simple model, which can be rapidly applied, could provide urban data for applications such as air quality forecasting and building energy modelling, in addition to outdoor thermal comfort.

A Summary of the NATO ASI on Polar Cap Boundary Phenomena

The polar cap boundary is a subject of central importance to current magnetosphere-ionosphere research and its applications in “space weather” activities. The problems are that it has a number of definitions, and that the most physically meaningful definition (namely the open-closed field line boundary) is very difficult to identify in observations. New understanding of the importance of the structure and dynamics of the boundary region made the time right for a meeting reviewing our knowledge in this area. The Advanced Study Institute (ASI) on Svalbard in June 1997 discussed the boundary on both the dayside and the nightside, mapping magnetically to the dayside magnetopause and to tail plasma sheet/lobe interface, respectively. We held a “brainstorming” session, in which different ideas which arose from the presented papers were discussed and developed, and a summary session, in which session convenors gave a personal view of progress that has been made and problems which still need solving. Both were designed as ways of promoting further discussion. This paper attempts to distil some of the themes that emerged from these discussions.

Auroral and plasma flow transients at magnetic noon

We present observations of a transient event in the dayside auroral ionosphere at magnetic noon. F-region plasma convection measurements were made by the EISCAT radar, operating in the beamswinging “Polar” experiment mode, and simultaneous observations of the dayside auroral emissions were made by optical meridian-scanning photometers and all-sky TV cameras at Ny Ålesund, Spitzbergen. The data were recorded on 9 January 1989, and a sequence of bursts of flow, with associated transient aurora, were observed between 08:45 and 11:00 U.T. In this paper we concentrate on an event around 09:05 U.T. because that is very close to local magnetic noon. The optical data show a transient intensification and widening (in latitude) of the cusp/cleft region, as seen in red line auroral emissions. Over an interval of about 10 min, the band of 630 nm aurora widened from about 1.5° of invariant latitude to over 5° and returned to its original width. Embedded within the widening band of 630 nm emissions were two intense, active 557.7 nm arc fragments with rays which persisted for about 2 min each. The flow data before and after the optical transient show eastward flows, with speeds increasing markedly with latitude across the band of 630 nm aurora. Strong, apparently westward, flows appeared inside the band while it was widening, but these rotated round to eastward, through northward, as the band shrunk to its original width. The observed ion temperatures verify that the flow speeds during the transient were, to a large extent, as derived using the beamswinging technique; but they also show that the flow increase initially occurred in the western azimuth only. This spatial gradient in the flow introduces ambiguity in the direction of these initial flows and they could have been north-eastward rather than westward. However, the westward direction derived by the beamswinging is consistent with the motion of the colocated and coincident active 557.7 nm arc fragment, A more stable transient 557.7 nm aurora was found close to the shear between the inferred westward flows and the persisting eastward flows to the North. Throughout the transient, northward flow was observed across the equatorward boundary of the 630 nm aurora. Interpretation of the data is made difficult by lack of IMF data, problems in distinguishing the cusp and cleft aurora and uncertainty over which field lines are open and which are closed. However, at magnetic noon there is a 50% probability that we were observing the cusp, in which case from its southerly location we infer that the IMF was southward and many features are suggestive of time-varying reconnection at a single X-line on the dayside magnetopause. This IMF orientation is also consistent with the polar rain precipitation observed simultaneously by the DMSP-F9 satellite in the southern polar cap. There is also a 25% chance that we were observing the cleft (or the mantle poleward of the cleft). In this case we infer that the IMF was northward and the transient is well explained by reconnection which is not only transient in time but occurs at various sites located randomly on the dayside magnetopause (i.e. patchy in space). Lastly, there is a 25% chance that we were observing the cusp poleward of the cleft, in which case we infer that IMF Bz was near zero and the transient is explained by a mixture of the previous two interpretations.

The neurological underpinnings of cluttering: some initial findings

Background Cluttering is a fluency disorder characterised by overly rapid or jerky speech patterns that compromise intelligibility. The neural correlates of cluttering are unknown but theoretical accounts implicate the basal ganglia and medial prefrontal cortex. Dysfunction in these brain areas would be consistent with difficulties in selection and control of speech motor programs that are characteristic of speech disfluencies in cluttering. There is a surprising lack of investigation into this disorder using modern imaging techniques. Here, we used functional MRI to investigate the neural correlates of cluttering. Method We scanned 17 adults who clutter and 17 normally fluent control speakers matched for age and sex. Brain activity was recorded using sparse-sampling functional MRI while participants viewed scenes and either (i) produced overt speech describing the scene or (ii) read out loud a sentence provided that described the scene. Speech was recorded and analysed off line. Differences in brain activity for each condition compared to a silent resting baseline and between conditions were analysed for each group separately (cluster-forming threshold Z > 3.1, extent p < 0.05, corrected) and then these differences were further compared between the two groups (voxel threshold p < 0.01, extent > 30 voxels, uncorrected). Results In both conditions, the patterns of activation in adults who clutter and control speakers were strikingly similar, particularly at the cortical level. Direct group comparisons revealed greater activity in adults who clutter compared to control speakers in the lateral premotor cortex bilaterally and, as predicted, on the medial surface (pre-supplementary motor area). Subcortically, adults who clutter showed greater activity than control speakers in the basal ganglia. Specifically, the caudate nucleus and putamen were overactive in adults who clutter for the comparison of picture description with sentence reading. In addition, adults who clutter had reduced activity relative to control speakers in the lateral anterior cerebellum bilaterally. Eleven of the 17 adults who clutter also stuttered. This comorbid diagnosis of stuttering was found to contribute to the abnormal overactivity seen in the group of adults who clutter in the right ventral premotor cortex and right anterior cingulate cortex. In the remaining areas of abnormal activity seen in adults who clutter compared to controls, the subgroup who clutter and stutter did not differ from the subgroup who clutter but do not stutter. Conclusions Our findings were in good agreement with theoretical predictions regarding the neural correlates of cluttering. We found evidence for abnormal function in the basal ganglia and their cortical output target, the medial prefrontal cortex. The findings are discussed in relation to models of cluttering that point to problems with motor control of speech.

A multi-model evaluation of aerosols over South Asia: common problems and possible causes

Atmospheric pollution over South Asia attracts special attention due to its effects on regional climate, water cycle and human health. These effects are potentially growing owing to rising trends of anthropogenic aerosol emissions. In this study, the spatio-temporal aerosol distributions over South Asia from seven global aerosol models are evaluated against aerosol retrievals from NASA satellite sensors and ground-based measurements for the period of 2000–2007. Overall, substantial underestimations of aerosol loading over South Asia are found systematically in most model simulations. Averaged over the entire South Asia, the annual mean aerosol optical depth (AOD) is underestimated by a range 15 to 44% across models compared to MISR (Multi-angle Imaging SpectroRadiometer), which is the lowest bound among various satellite AOD retrievals (from MISR, SeaWiFS (Sea-Viewing Wide Field-of-View Sensor), MODIS (Moderate Resolution Imaging Spectroradiometer) Aqua and Terra). In particular during the post-monsoon and wintertime periods (i.e., October–January), when agricultural waste burning and anthropogenic emissions dominate, models fail to capture AOD and aerosol absorption optical depth (AAOD) over the Indo–Gangetic Plain (IGP) compared to ground-based Aerosol Robotic Network (AERONET) sunphotometer measurements. The underestimations of aerosol loading in models generally occur in the lower troposphere (below 2 km) based on the comparisons of aerosol extinction profiles calculated by the models with those from Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) data. Furthermore, surface concentrations of all aerosol components (sulfate, nitrate, organic aerosol (OA) and black carbon (BC)) from the models are found much lower than in situ measurements in winter. Several possible causes for these common problems of underestimating aerosols in models during the post-monsoon and wintertime periods are identified: the aerosol hygroscopic growth and formation of secondary inorganic aerosol are suppressed in the models because relative humidity (RH) is biased far too low in the boundary layer and thus foggy conditions are poorly represented in current models, the nitrate aerosol is either missing or inadequately accounted for, and emissions from agricultural waste burning and biofuel usage are too low in the emission inventories. These common problems and possible causes found in multiple models point out directions for future model improvements in this important region.

On the structure of infinity-harmonic maps

Let H ∈ C 2(ℝ N×n ), H ≥ 0. The PDE system arises as the Euler-Lagrange PDE of vectorial variational problems for the functional E ∞(u, Ω) = ‖H(Du)‖ L ∞(Ω) defined on maps u: Ω ⊆ ℝ n → ℝ N . (1) first appeared in the author's recent work. The scalar case though has a long history initiated by Aronsson. Herein we study the solutions of (1) with emphasis on the case of n = 2 ≤ N with H the Euclidean norm on ℝ N×n , which we call the “∞-Laplacian”. By establishing a rigidity theorem for rank-one maps of independent interest, we analyse a phenomenon of separation of the solutions to phases with qualitatively different behaviour. As a corollary, we extend to N ≥ 2 the Aronsson-Evans-Yu theorem regarding non existence of zeros of |Du| and prove a maximum principle. We further characterise all H for which (1) is elliptic and also study the initial value problem for the ODE system arising for n = 1 but with H(·, u, u′) depending on all the arguments.

Sparse space-time boundary element methods for the heat equation

The goal of this work is the efficient solution of the heat equation with Dirichlet or Neumann boundary conditions using the Boundary Elements Method (BEM). Efficiently solving the heat equation is useful, as it is a simple model problem for other types of parabolic problems. In complicated spatial domains as often found in engineering, BEM can be beneficial since only the boundary of the domain has to be discretised. This makes BEM easier than domain methods such as finite elements and finite differences, conventionally combined with time-stepping schemes to solve this problem. The contribution of this work is to further decrease the complexity of solving the heat equation, leading both to speed gains (in CPU time) as well as requiring smaller amounts of memory to solve the same problem. To do this we will combine the complexity gains of boundary reduction by integral equation formulations with a discretisation using wavelet bases. This reduces the total work to O(h

A sodar for profiling in a spatially inhomogeneous urban environment

The urban boundary layer, above the canopy, is still poorly understood. One of the challenges is obtaining data by sampling more than a few meters above the rooftops, given the spatial and temporal inhomogeneities in both horizontal and vertical. Sodars are generally useful tools for ground-based remote sensing of winds and turbulence, but rely on horizontal homogeneity (as do lidars) in building up 3-component wind vectors from sampling three or more spatially separated volumes. The time taken for sound to travel to a typical range of 200 m and back is also a limitation. A sodar of radically different design is investigated, aimed at addressing these problems. It has a single vertical transmitted sound pulse. Doppler shifted signals are received from a number of volumes around the periphery of the transmitted beam with microphones which each having tight angular sensitivity at zenith angles slightly off-vertical. The spatial spread of sampled volumes is therefore smaller. By having more receiver microphones than a conventional sodar, the effect of smaller zenith angle is offset. More rapid profiling is also possible with a single vertical transmitted beam, instead of the usual multiple beams.A prototype design is described, together with initial field measurements. It is found that the beam forming using a single dish antenna and the drift of the sound pulse downwind both give rise to reduced performance compared with expectations. It is concluded that, while the new sodar works in principle, the compromises arising in the design mean that the expected advantages have not been realized

Advancements in decadal climate predictability: the role of nonoceanic drivers

We review recent progress in understanding the role of sea ice, land surface, stratosphere, and aerosols in decadal-scale predictability and discuss the perspectives for improving the predictive capabilities of current Earth system models (ESMs). These constituents have received relatively little attention because their contribution to the slow climatic manifold is controversial in comparison to that of the large heat capacity of the oceans. Furthermore, their initialization as well as their representation in state-of-the-art climate models remains a challenge. Numerous extraoceanic processes that could be active over the decadal range are proposed. Potential predictability associated with the aforementioned, poorly represented, and scarcely observed constituents of the climate system has been primarily inspected through numerical simulations performed under idealized experimental settings. The impact, however, on practical decadal predictions, conducted with realistically initialized full-fledged climate models, is still largely unexploited. Enhancing initial-value predictability through an improved model initialization appears to be a viable option for land surface, sea ice, and, marginally, the stratosphere. Similarly, capturing future aerosol emission storylines might lead to an improved representation of both global and regional short-term climatic changes. In addition to these factors, a key role on the overall predictive ability of ESMs is expected to be played by an accurate representation of processes associated with specific components of the climate system. These act as “signal carriers,” transferring across the climatic phase space the information associated with the initial state and boundary forcings, and dynamically bridging different (otherwise unconnected) subsystems. Through this mechanism, Earth system components trigger low-frequency variability modes, thus extending the predictability beyond the seasonal scale.

Langmuir turbulence and surface heating in the ocean surface boundary layer

This study uses large-eddy simulation to investigate the structure of the ocean surface boundary layer (OSBL) in the presence of Langmuir turbulence and stabilizing surface heat fluxes. The OSBL consists of a weakly stratified layer, despite a surface heat flux, above a stratified thermocline. The weakly stratified (mixed) layer is maintained by a combination of a turbulent heat flux produced by the wave-driven Stokes drift and downgradient turbulent diffusion. The scaling of turbulence statistics, such as dissipation and vertical velocity variance, is only affected by the surface heat flux through changes in the mixed layer depth. Diagnostic models are proposed for the equilibrium boundary layer and mixed layer depths in the presence of surface heating. The models are a function of the initial mixed layer depth before heating is imposed and the Langmuir stability length. In the presence of radiative heating, the models are extended to account for the depth profile of the heating.

Resonant Wave Interactions in the Equatorial Waveguide

Weakly nonlinear interactions among equatorial waves have been explored in this paper using the adiabatic version of the equatorial beta-plane primitive equations in isobaric coordinates. Assuming rigid lid vertical boundary conditions, the conditions imposed at the surface and at the top of the troposphere were expanded in a Taylor series around two isobaric surfaces in an approach similar to that used in the theory of surface-gravity waves in deep water and capillary-gravity waves. By adopting the asymptotic method of multiple time scales, the equatorial Rossby, mixed Rossby-gravity, inertio-gravity, and Kelvin waves, as well as their vertical structures, were obtained as leading-order solutions. These waves were shown to interact resonantly in a triad configuration at the O(epsilon) approximation. The resonant triads whose wave components satisfy a resonance condition for their vertical structures were found to have the most significant interactions, although this condition is not excluding, unlike the resonant conditions for the zonal wavenumbers and meridional modes. Thus, the analysis has focused on such resonant triads. In general, it was found that for these resonant triads satisfying the resonance condition in the vertical direction, the wave with the highest absolute frequency always acts as an energy source (or sink) for the remaining triad components, as usually occurs in several other physical problems in fluid dynamics. In addition, the zonally symmetric geostrophic modes act as catalyst modes for the energy exchanges between two dispersive waves in a resonant triad. The integration of the reduced asymptotic equations for a single resonant triad shows that, for the initial mode amplitudes characterizing realistic magnitudes of atmospheric flow perturbations, the modes in general exchange energy on low-frequency (intraseasonal and/or even longer) time scales, with the interaction period being dependent upon the initial mode amplitudes. Potential future applications of the present theory to the real atmosphere with the inclusion of diabatic forcing, dissipation, and a more realistic background state are also discussed.

Continuity of the dynamics in a localized large diffusion problem with nonlinear boundary conditions

This paper is concerned with singular perturbations in parabolic problems subjected to nonlinear Neumann boundary conditions. We consider the case for which the diffusion coefficient blows up in a subregion Omega(0) which is interior to the physical domain Omega subset of R(n). We prove, under natural assumptions, that the associated attractors behave continuously as the diffusion coefficient blows up locally uniformly in Omega(0) and converges uniformly to a continuous and positive function in Omega(1) = (Omega) over bar\Omega(0). (C) 2009 Elsevier Inc. All rights reserved.

Coincidence theorems and its applications to equilibrium problems

Given two maps h : X x K -> R and g : X -> K such that, for all x is an element of X, h(x, g(x)) = 0, we consider the equilibrium problem of finding (x) over tilde is an element of X such that h((x) over tilde, g(x)) >= 0 for every x is an element of X. This question is related to a coincidence problem.