Optimal power management strategy for energy storage with stochastic loads

In this paper, a power management strategy (PMS) has been developed for the control of energy storage in a system subjected to loads of random duration. The PMS minimises the costs associated with the energy consumption of specific systems powered by a primary energy source and equipped with energy storage, under the assumption that the statistical distribution of load durations is known. By including the variability of the load in the cost function, it was possible to define the optimality criteria for the power flow of the storage. Numerical calculations have been performed obtaining the control strategies associated with the global minimum in energy costs, for a wide range of initial conditions of the system. The results of the calculations have been tested on a MATLAB/Simulink model of a rubber tyre gantry (RTG) crane equipped with a flywheel energy storage system (FESS) and subjected to a test cycle, which corresponds to the real operation of a crane in the Port of Felixstowe. The results of the model show increased energy savings and reduced peak power demand with respect to existing control strategies, indicating considerable potential savings for port operators in terms of energy and maintenance costs.

Stability of an ice sheet on an elastic bed

The stability of stationary flow of a two-dimensional ice sheet is studied when the ice obeys a power flow law (Glen's flow law). The mass accumulation rate at the top is assumed to depend on elevation and span and the bed supporting the ice sheet consists of an elastic layer lying on a rigid surface. The normal perturbation of the free surface of the ice sheet is a singular eigenvalue problem. The singularity of the perturbation at the front of the ice sheet is considered using matched asymptotic expansions, and the eigenvalue problem is seen to reduce to that with fixed ice front. Numerical solution of the perturbation eigenvalue problem shows that the dependence of accumulation rate on elevation permits the existence of unstable solutions when the equilibrium line is higher than the bed at the ice divide. Alternatively, when the equilibrium line is lower than the bed, there are only stable solutions. Softening of the bed, expressed through a decrease of its elastic modulus, has a stabilising effect on the ice sheet.

Scattered power from non-thermal, F-region plasma observed by EISCAT—evidence for coherent echoes?

Three rapid, poleward bursts of plasma flow, observed by the U.K.-POLAR EISCAT experiment, are studied in detail. In all three cases the large ion velocities (> 1 kms−1) are shown to drive the ion velocity distribution into a non-Maxwellian form, identified by the characteristic shape of the observed spectra and the fact that analysis of the spectra with the assumption of a Maxwellian distribution leads to excessive rises in apparent ion temperature, and an anticorrelation of apparent electron and ion temperatures. For all three periods the total scattered power is shown to rise with apparent ion temperature by up to 6 dB more than is expected for an isotropic Maxwellian plasma of constant density and by an even larger factor than that expected for non-thermal plasma. The anomalous increases in power are only observed at the lower altitudes (< 300 km). At greater altitudes the rise in power is roughly consistent with that simulated numerically for homogeneous, anisotropic, non-Maxwellian plasma of constant density, viewed using the U.K.-POLAR aspect angle. The spectra at times of anomalously high power are found to be asymmetric, showing an enhancement near the downward Doppler-shifted ion-acoustic frequency. Although it is not possible to eliminate completely rapid plasma density fluctuations as a cause of these power increases, such effects cannot explain the observed spectra and the correlation of power and apparent ion temperature without an unlikely set of coincidences. The observations are made along a beam direction which is as much as 16.5° from orthogonality with the geomagnetic field. Nevertheless, some form of coherent-like echo contamination of the incoherent scatter spectrum is the most satisfactory explanation of these data.

Inter-hemispheric mesospheric coupling in a comprehensive middle atmosphere model

Observations of noctilucent clouds have revealed a surprising coupling between the winter stratosphere and the summer polar mesopause region. In spite of the great distance involved, this inter-hemispheric link has been suggested to be the principal reason for both the year-to-year variability and the hemispheric differences in the frequency of occurrence of these high-altitude clouds. In this study, we investigate the dynamical influence of the winter stratosphere on the summer mesosphere using simulations from the vertically extended version of the Canadian Middle Atmosphere Model (CMAM). We find that for both Northern and Southern Hemispheres, variability in the summer polar mesopause region from one year to another can be traced back to the planetary-wave flux entering the winter stratosphere. The teleconnection pattern is the same for both positive and negative wave-flux anomalies. Using a composite analysis to isolate the events, it is argued that the mechanism for interhemispheric coupling is a feedback between summer mesosphere gravity-wave drag (GWD) and zonal wind, which is induced by an anomaly in mesospheric cross-equatorial flow, the latter arising from the anomaly in winter hemisphere GWD induced by the anomaly in stratospheric conditions.

Detection of solar dimming and brightening effects on Northern Hemisphere river flow

Anthropogenic aerosols in the atmosphere have the potential to affect regional-scale land hydrology through solar dimming. Increased aerosol loading may have reduced historical surface evaporation over some locations, but the magnitude and extent of this effect is uncertain. Any reduction in evaporation due to historical solar dimming may have resulted in an increase in river flow. Here we formally detect and quantify the historical effect of changing aerosol concentrations, via solar radiation, on observed river flow over the heavily industrialized, northern extra-tropics. We use a state-of-the-art estimate of twentieth century surface meteorology as input data for a detailed land surface model, and show that the simulations capture the observed strong inter-annual variability in runoff in response to climatic fluctuations. Using statistical techniques, we identify a detectable aerosol signal in the observed river flow both over the combined region, and over individual river basins in Europe and North America. We estimate that solar dimming due to rising aerosol concentrations in the atmosphere around 1980 led to an increase in river runoff by up to 25% in the most heavily polluted regions in Europe. We propose that, conversely, these regions may experience reduced freshwater availability in the future, as air quality improvements are set to lower aerosol loading and solar dimming.

Becoming a region, becoming global, becoming imperceptable: territorialising salmon in Chilean Patagonia

Our article focuses on the region of Chilean Patagonia and considers how it has developed as a leading producer of salmon for global food markets. It addresses the problem of how to decentre conventional views of the forces driving regional development that give primacy to the role of capital and technology, instead giving due recognition to the knowledge and practices of situated actors and to the relationships that form between human and non-human entities in food producing regions. As an alternative, we ask whether an assemblage approach can improve our understanding of regional transformation. To explore this question, we present original ethnographic data on constitutive practices that have transformed the Patagonian region, from the territorialization of Salmonidae species to experimentation in ocean ranching and sea water fish farming, and finally the development of a global industry. The evidence leads us to argue that in a complex globalised world, assemblage theory offers a valuable approach for understanding how regional potential is realised. In the case of Chilean Patagonia, it is apparent that forms of bio-power generate new relations between life, agency and nature, stimulating contemporary regional transformations in ways overlooked by the lineal logic of capital objectification discourses. Applying an assemblage approach enables the significance of new contemporary human – non-human relationships and inter-subjectivities to come to the fore, keeping the social in view as potential for regional transformation and new power asymmetries continuously emerge.

Tropical cyclones in a T159 Resolution Global Climate Model: comparison with observations and re-analyses

Tropical cyclones have been investigated in a T159 version of the MPI ECHAM5 climate model using a novel technique to diagnose the evolution of the 3-dimensional vorticity structure of tropical cyclones, including their full life cycle from weak initial vortex to their possible extra-tropical transition. Results have been compared with reanalyses (ERA40 and JRA25) and observed tropical storms during the period 1978-1999 for the Northern Hemisphere. There is no indication of any trend in the number or intensity of tropical storms during this period in ECHAM5 or in re-analyses but there are distinct inter-annual variations. The storms simulated by ECHAM5 are realistic both in space and time, but the model and even more so the re-analyses, underestimate the intensities of the most intense storms (in terms of their maximum wind speeds). There is an indication of a response to ENSO with a smaller number of Atlantic storms during El Niño in agreement with previous studies. The global divergence circulation responds to El Niño by setting up a large-scale convergence flow, with the center over the central Pacific with enhanced subsidence over the tropical Atlantic. At the same time there is an increase in the vertical wind shear in the region of the tropical Atlantic where tropical storms normally develop. There is a good correspondence between the model and ERA40 except that the divergence circulation is somewhat stronger in the model. The model underestimates storms in the Atlantic but tends to overestimate them in the Western Pacific and in the North Indian Ocean. It is suggested that the overestimation of storms in the Pacific by the model is related to an overly strong response to the tropical Pacific SST anomalies. The overestimation in 2 the North Indian Ocean is likely to be due to an over prediction in the intensity of monsoon depressions, which are then classified as intense tropical storms. Nevertheless, overall results are encouraging and will further contribute to increased confidence in simulating intense tropical storms with high-resolution climate models.

Impact of Tibetan Orography and Heating on the Summer Flow over Asia

The influence on the summer flow over Asia of both the orographic and thermal forcing of the Tibetan Plateau is investigated using a sequence of idealised experiments with a global primitive equation model. The zonally averaged flow is prescribed and both realistic and idealised orography and heating are used. There is some similarity between the responses to the two forcings when applied separately. The upper tropospheric Tibetan anticyclone is predominantly forced by the heating but also weakly by the orography. Below this, both forcings tend to give air descending in an equatorward anticyclonic circulation down the isentropes to the west and rising in a similar poleward circulation to the east. However the heating-only response has a strong ascending southwesterly flow that is guided around the south and south-east of the orography when it is included. On the northern side, the westerly flow over the orography gives ascent on the upslope and descent on the downslope. It is found that heating over the Plateau leads to a potential vorticity (PV) minimum and that if it is sufficiently strong the flow is unstable, producing a quasi-biweekly oscillation. During this oscillation the Tibetan anticyclone changes between a single centre over the southwestern side of the Plateau and a split/double structure with centres over China and the Middle East. These characteristics are similar to observed variability in the region. Associated with this quasi-biweekly oscillation are significant variations in the strength of the ascent over the Plateau and the Rossby wave pattern over the North Pacific. The origin of the variability is instability associated with the zonally extended potential vorticity PV minimum on a θ-surface, as proposed by Hsu and Plumb (2000). This minimum is due to the tendency to reduce the PV above the heating over the Plateau and to advection by the consequent anticyclone of high PV around from the east and low PV to the west. The deep convection to the south and southeast of the Plateau tends to suppress the quasi-biweekly oscillation because the low PV produced above it acts to reduce the meridional PV gradient reversal. The occurrence of the oscillation depends on the relative magnitude of the heating in the two regions.

Spatial variability of flow statistics within regular building arrays

Turbulence statistics obtained by direct numerical simulations are analysed to investigate spatial heterogeneity within regular arrays of building-like cubical obstacles. Two different array layouts are studied, staggered and square, both at a packing density of λp=0.25 . The flow statistics analysed are mean streamwise velocity ( u− ), shear stress ( u′w′−−−− ), turbulent kinetic energy (k) and dispersive stress fraction ( u˜w˜ ). The spatial flow patterns and spatial distribution of these statistics in the two arrays are found to be very different. Local regions of high spatial variability are identified. The overall spatial variances of the statistics are shown to be generally very significant in comparison with their spatial averages within the arrays. Above the arrays the spatial variances as well as dispersive stresses decay rapidly to zero. The heterogeneity is explored further by separately considering six different flow regimes identified within the arrays, described here as: channelling region, constricted region, intersection region, building wake region, canyon region and front-recirculation region. It is found that the flow in the first three regions is relatively homogeneous, but that spatial variances in the latter three regions are large, especially in the building wake and canyon regions. The implication is that, in general, the flow immediately behind (and, to a lesser extent, in front of) a building is much more heterogeneous than elsewhere, even in the relatively dense arrays considered here. Most of the dispersive stress is concentrated in these regions. Considering the experimental difficulties of obtaining enough point measurements to form a representative spatial average, the error incurred by degrading the sampling resolution is investigated. It is found that a good estimate for both area and line averages can be obtained using a relatively small number of strategically located sampling points.

Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

The impact of targeted observations made during the Greenland Flow Distortion Experiment

The impact of targeted sonde observations on the 1-3 day forecasts for northern Europe is evaluated using the Met Office four-dimensional variational data assimilation scheme and a 24 km gridlength limited-area version of the Unified Model (MetUM). The targeted observations were carried out during February and March 2007 as part of the Greenland Flow Distortion Experiment, using a research aircraft based in Iceland. Sensitive area predictions using either total energy singular vectors or an ensemble transform Kalman filter were used to predict where additional observations should be made to reduce errors in the initial conditions of forecasts for northern Europe. Targeted sonde data was assimilated operationally into the MetUM. Hindcasts show that the impact of the sondes was mixed. Only two out of the five cases showed clear forecast improvement; the maximum forecast improvement seen over the verifying region was approximately 5% of the forecast error 24 hours into the forecast. These two cases are presented in more detail: in the first the improvement propagates into the verification region with a developing polar low; and in the second the improvement is associated with an upper-level trough. The impact of cycling targeted data in the background of the forecast (including the memory of previous targeted observations) is investigated. This is shown to cause a greater forecast impact, but does not necessarily lead to a greater forecast improvement. Finally, the robustness of the results is assessed using a small ensemble of forecasts.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

A Glacial Geomorphological Map of Victoria Island, Canadian Arctic

Victoria Island lies at the north-western extremity of the region covered by the vast North American Laurentide Ice Sheet (LIS) in the Canadian Arctic Archipelago. This area is significant because it linked the interior of the LIS to the Arctic Ocean, probably via a number of ice streams. Victoria Island, however, exhibits a remarkably complex glacial landscape, with several successive generations of ice flow indicators superimposed on top of each other and often at abrupt (90 degrees) angles. This complexity represents a major challenge to those attempting to produce a detailed reconstruction of the glacial history of the region. This paper presents a map of the glacial geomorphology of Victoria Island. The map is based on analysis of Landsat Enhanced Thematic Plus (ETM+) satellite imagery and contains over 58,000 individual glacial features which include: glacial lineations, moraines (terminal, lateral, subglacial shear margin), hummocky moraine, ribbed moraine, eskers, glaciofluvial deposits, large meltwater channels, and raised shorelines. The glacial features reveal marked changes in ice flow direction and vigour over time. Moreover, the glacial geomorphology indicates a non-steady withdrawal of ice during deglaciation, with rapidly flowing ice streams focussed into the inter-island troughs and several successively younger flow patterns superimposed on older ones. It is hoped that detailed analysis of this map will lead to an improved reconstruction of the glacial history of this area which will provide other important insights, for example, with respect to the interactions between ice streaming, deglaciation and Arctic Ocean meltwater events.

The unsteady flow of a weakly compressible fluid in a thin porous layer. I: Two-dimensional theory

We consider the problem of determining the pressure and velocity fields for a weakly compressible fluid flowing in a two-dimensional reservoir in an inhomogeneous, anisotropic porous medium, with vertical side walls and variable upper and lower boundaries, in the presence of vertical wells injecting or extracting fluid. Numerical solution of this problem may be expensive, particularly in the case that the depth scale of the layer h is small compared to the horizontal length scale l. This is a situation which occurs frequently in the application to oil reservoir recovery. Under the assumption that epsilon=h/l<<1, we show that the pressure field varies only in the horizontal direction away from the wells (the outer region). We construct two-term asymptotic expansions in epsilon in both the inner (near the wells) and outer regions and use the asymptotic matching principle to derive analytical expressions for all significant process quantities. This approach, via the method of matched asymptotic expansions, takes advantage of the small aspect ratio of the reservoir, epsilon, at precisely the stage where full numerical computations become stiff, and also reveals the detailed structure of the dynamics of the flow, both in the neighborhood of wells and away from wells.