Four-dimensional data assimilation and balanced dynamics

The problem of spurious excitation of gravity waves in the context of four-dimensional data assimilation is investigated using a simple model of balanced dynamics. The model admits a chaotic vortical mode coupled to a comparatively fast gravity wave mode, and can be initialized such that the model evolves on a so-called slow manifold, where the fast motion is suppressed. Identical twin assimilation experiments are performed, comparing the extended and ensemble Kalman filters (EKF and EnKF, respectively). The EKF uses a tangent linear model (TLM) to estimate the evolution of forecast error statistics in time, whereas the EnKF uses the statistics of an ensemble of nonlinear model integrations. Specifically, the case is examined where the true state is balanced, but observation errors project onto all degrees of freedom, including the fast modes. It is shown that the EKF and EnKF will assimilate observations in a balanced way only if certain assumptions hold, and that, outside of ideal cases (i.e., with very frequent observations), dynamical balance can easily be lost in the assimilation. For the EKF, the repeated adjustment of the covariances by the assimilation of observations can easily unbalance the TLM, and destroy the assumptions on which balanced assimilation rests. It is shown that an important factor is the choice of initial forecast error covariance matrix. A balance-constrained EKF is described and compared to the standard EKF, and shown to offer significant improvement for observation frequencies where balance in the standard EKF is lost. The EnKF is advantageous in that balance in the error covariances relies only on a balanced forecast ensemble, and that the analysis step is an ensemble-mean operation. Numerical experiments show that the EnKF may be preferable to the EKF in terms of balance, though its validity is limited by ensemble size. It is also found that overobserving can lead to a more unbalanced forecast ensemble and thus to an unbalanced analysis.

Parallel pipelined histogram architecture via c-slow retiming

A parallel pipelined array of cells suitable for realtime computation of histograms is proposed. The cell architecture builds on previous work to now allow operating on a stream of data at 1 pixel per clock cycle. This new cell is more suitable for interfacing to camera sensors or to microprocessors of 8-bit data buses which are common in consumer digital cameras. Arrays using the new proposed cells are obtained via C-slow retiming techniques and can be clocked at a 65% faster frequency than previous arrays. This achieves over 80% of the performance of two-pixel per clock cycle parallel pipelined arrays.

Very high resolution bioclimatic zoning of Portuguese wine regions: present and future scenarios

Wine production is strongly affected by weather and climate and thus highly vulnerable to climate change. In Portugal, viticulture and wine production are an important economic activity. In the present study, current bioclimatic zoning in Portugal (1950–2000) and its projected changes under future climate conditions (2041–2070) are assessed through the analysis of an aggregated, categorized bioclimatic index (CatI) at a very high spatial resolution (near 1 km). CatI incorporates the most relevant bioclimatic characteristics of a given region, thus allowing the direct comparison between different regions. Future viticultural zoning is achieved using data from 13 climate model transient experiments following the A1B emission scenario. These data are downscaled using a two-step method of spatial pattern downscaling. This downscaling approach allows characterizing mesoclimatic influences on viticulture throughout Portugal. Results for the recent past depict the current spatial variability of Portuguese viticultural regions. Under future climate conditions, the current viticultural zoning is projected to undergo significant changes, which may represent important challenges for the Portuguese winemaking sector. The changes are quite robust across the different climate models. A lower bioclimatic diversity is also projected, resulting from a more homogeneous warm and dry climate in most of the wine regions. This will lead to changes in varietal suitability and wine characteristics of each region.

Comments on 'Balance and the slow quasimanifold: some explicit results'

The concept of slow vortical dynamics and its role in theoretical understanding is central to geophysical fluid dynamics. It leads, for example, to “potential vorticity thinking” (Hoskins et al. 1985). Mathematically, one imagines an invariant manifold within the phase space of solutions, called the slow manifold (Leith 1980; Lorenz 1980), to which the dynamics are constrained. Whether this slow manifold truly exists has been a major subject of inquiry over the past 20 years. It has become clear that an exact slow manifold is an exceptional case, restricted to steady or perhaps temporally periodic flows (Warn 1997). Thus the concept of a “fuzzy slow manifold” (Warn and Ménard 1986) has been suggested. The idea is that nearly slow dynamics will occur in a stochastic layer about the putative slow manifold. The natural question then is, how thick is this layer? In a recent paper, Ford et al. (2000) argue that Lighthill emission—the spontaneous emission of freely propagating acoustic waves by unsteady vortical flows—is applicable to the problem of balance, with the Mach number Ma replaced by the Froude number F, and that it is a fundamental mechanism for this fuzziness. They consider the rotating shallow-water equations and find emission of inertia–gravity waves at O(F2). This is rather surprising at first sight, because several studies of balanced dynamics with the rotating shallow-water equations have gone beyond second order in F, and found only an exponentially small unbalanced component (Warn and Ménard 1986; Lorenz and Krishnamurthy 1987; Bokhove and Shepherd 1996; Wirosoetisno and Shepherd 2000). We have no technical objection to the analysis of Ford et al. (2000), but wish to point out that it depends crucially on R 1, where R is the Rossby number. This condition requires the ratio of the characteristic length scale of the flow L to the Rossby deformation radius LR to go to zero in the limit F → 0. This is the low Froude number scaling of Charney (1963), which, while originally designed for the Tropics, has been argued to be also relevant to mesoscale dynamics (Riley et al. 1981). If L/LR is fixed, however, then F → 0 implies R → 0, which is the standard quasigeostrophic scaling of Charney (1948; see, e.g., Pedlosky 1987). In this limit there is reason to expect the fuzziness of the slow manifold to be “exponentially thin,” and balance to be much more accurate than is consistent with (algebraic) Lighthill emission.

C-slow retimed parallel histogram architectures for consumer imaging devices

A parallel pipelined array of cells suitable for real-time computation of histograms is proposed. The cell architecture builds on previous work obtained via C-slow retiming techniques and can be clocked at 65 percent faster frequency than previous arrays. The new arrays can be exploited for higher throughput particularly when dual data rate sampling techniques are used to operate on single streams of data from image sensors. In this way, the new cell operates on a p-bit data bus which is more convenient for interfacing to camera sensors or to microprocessors in consumer digital cameras.

Retrieval characteristics of non-linear sea surface temperature from the Advanced Very High Resolution Radiometer

Criteria are proposed for evaluating sea surface temperature (SST) retrieved from satellite infra-red imagery: bias should be small on regional scales; sensitivity to atmospheric humidity should be small; and sensitivity of retrieved SST to surface temperature should be close to 1 K K−1. Their application is illustrated for non-linear sea surface temperature (NLSST) estimates. 233929 observations from the Advanced Very High Resolution Radiometer (AVHRR) on Metop-A are matched with in situ data and numerical weather prediction (NWP) fields. NLSST coefficients derived from these matches have regional biases from −0.5 to +0.3 K. Using radiative transfer modelling we find that a 10% increase in humidity alone can change the retrieved NLSST by between −0.5 K and +0.1 K. A 1 K increase in SST changes NLSST by <0.5 K in extreme cases. The validity of estimates of sensitivity by radiative transfer modelling is confirmed empirically.

Relationship between classroom authority and epistemological beliefs as espoused by primary school mathematics teachers from the very high and very low socio-economic regions in Thailand

This article presents findings of a larger single-country comparative study which set out to better understand primary school teachers’ mathematics education-related beliefs in Thailand. By combining the interview and observation data collected in the initial stage of this study with data gathered from the relevant literature, the 8-belief / 22-item ‘Thai Teachers’ Mathematics Education-related Beliefs’ (TTMEB) Scale was developed. The results of the Mann-Whitney U Test showed that Thai teachers in the two examined socio-economic regions espouse statistically different beliefs concerning the source and stability of mathematical knowledge, as well as classroom authority. Further, these three beliefs are found to be significantly and positively correlated.

Potential of the slow component of quartz OSL for age determination of sedimentary samples

The slow component of quartz OSL exhibits a high thermal stability, and, in many of the samples studied, a high dose saturation level (several hundreds or, even thousands, of Grays). These properties suggest that the slow component has potential as a long-range dating tool. Initial attempts have been made to obtain equivalent doses for a number of sedimentary samples. Single- and multiple-aliquot techniques were modified for use with the slow component. The results indicate that there is a good potential for sediment dating, particularly for samples of significant age. Experiments concerning the optical resetting of the slow component suggest that, given its slow optical depletion rate, dating may be restricted to aeolian sediments.

Global patterns of vegetation response to millennial-scale variability and rapid climate change during the last glacial period

Ninety-four sites worldwide have sufficient resolution and dating to document the impact of millennial-scale climate variability on vegetation and fire regimes during the last glacial period. Although Dansgaard–Oeschger (D–O) cycles all show a basically similar gross structure, they vary in the magnitude and the length of the warm and cool intervals. We illustrate the geographic patterns in the climate-induced changes in vegetation by comparing D–O 6, D–O 8 and D–O 19. There is a strong response to both D–O warming events and subsequent cooling, most marked in the northern extratropics. Pollen records from marine cores from the northern extratropics confirm that there is no lag between the change in climate and the vegetation response, within the limits of the dating resolution (50–100 years). However, the magnitude of the change in vegetation is regionally specific and is not a simple function of either the magnitude or the duration of the change in climate as registered in Greenland ice cores. Fire regimes also show an initial immediate response to climate changes, but during cooling intervals there is a slow recovery of biomass burning after the initial reduction, suggesting a secondary control through the recovery of vegetation productivity. In the extratropics, vegetation changes are largely determined by winter temperatures while in the tropics they are largely determined by changes in plant-available water. Tropical vegetation records show changes corresponding to Heinrich Stadials but the response to D–O warming events is less marked than in the northern extratropics. There are very few high-resolution records from the Southern Hemisphere extratropics, but these records also show both a vegetation and fire response to millennial-scale climate variability. It is not yet possible to determine unequivocally whether terrestrial records reflect the asynchroneity apparent in the ice-core records.

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.

Strong armies, slow adaptation: civil-military relations and diffusion of military power

Why are some states more willing to adopt military innovations than others? Why, for example, were the great powers of Europe able to successfully reform their military practices to better adapt to and participate in the so-called military revolution of the sixteenth and seventeenth centuries while their most important extra-European competitor, the Ottoman Empire, failed to do so? This puzzle is best explained by two factors: civil-military relations and historical timing. In the Ottoman Empire, the emergence of an institutionally strong and internally cohesive army during the early stages of state formation—in the late fourteenth century—equipped the military with substantial bargaining powers. In contrast, the great powers of Europe drew heavily on private providers of military power during the military revolution and developed similar armies only by the second half of the seventeenth century, limiting the bargaining leverage of European militaries over their rulers. In essence, the Ottoman standing army was able to block reform efforts that it believed challenged its parochial interests. Absent a similar institutional challenge, European rulers initiated military reforms and motivated officers and military entrepreneurs to participate in the ongoing military revolution.

“There’s something very familiar about all this”: generic play and performance in the Back to the Future trilogy (Robert Zemeckis, 1985, 1989, 1990)

The Back to the Future Trilogy incorporates several different generic elements, including aspects of the fifties teen movie, science fiction, comedy and the western. These different modes playfully intertwine with each other creating a complex world of repetitions, echoes and modulations. This essay seeks to interrogate the construction of generic elements and the play between them through a close analysis of a repeated performance. Genre is signalled through various strategies employed within the construction of mise-en-scène, a significant portion of this, as I would like to argue, is transmitted through performance. The material detail of a performance – incorporating gesture, movement, voice, and even surrounding elements such as costume – as well as the way it its presented within a film is key to the establishment, invocation and coherence of genre. Furthermore, attention to the complexity of performance details, particularly in the manner in which they reverberate across texts, demonstrates the intricacy of genre and its inherent mutability. The Back to the Future trilogy represents a specific interest in the flexibility of genre. Within each film, and especially across all three, aspects of various genres are interlaced through both visual and narrative detail, thus constructing a dense layer of references both within and without the texts. To explore this patterning in more detail I will interrogate the contribution of performance to generic play through close analysis of Thomas F. Wilson’s performance of Biff/Griff/Burford Tannen and his central encounter with Marty McFly (Michael J. Fox) in each film. These moments take place in a fifties diner, a 1980s retro diner and a saloon respectively, each space contributing the similarities and differences in each repetition. Close attention to Wilson’s performance of each related character, which contains both modulations and repetitions used specifically to place each film’s central generic theme, demonstrates how embedded the play between genres and their flexibility is within the trilogy.

Instantaneous gelation in Smoluchowski’s coagulation equation revisited

We study the solutions of the Smoluchowski coagulation equation with a regularization term which removes clusters from the system when their mass exceeds a specified cutoff size, M. We focus primarily on collision kernels which would exhibit an instantaneous gelation transition in the absence of any regularization. Numerical simulations demonstrate that for such kernels with monodisperse initial data, the regularized gelation time decreasesas M increases, consistent with the expectation that the gelation time is zero in the unregularized system. This decrease appears to be a logarithmically slow function of M, indicating that instantaneously gelling kernels may still be justifiable as physical models despite the fact that they are highly singular in the absence of a cutoff. We also study the case when a source of monomers is introduced in the regularized system. In this case a stationary state is reached. We present a complete analytic description of this regularized stationary state for the model kernel, K(m1,m2)=max{m1,m2}ν, which gels instantaneously when M→∞ if ν>1. The stationary cluster size distribution decays as a stretched exponential for small cluster sizes and crosses over to a power law decay with exponent ν for large cluster sizes. The total particle density in the stationary state slowly vanishes as [(ν−1)logM]−1/2 when M→∞. The approach to the stationary state is nontrivial: Oscillations about the stationary state emerge from the interplay between the monomer injection and the cutoff, M, which decay very slowly when M is large. A quantitative analysis of these oscillations is provided for the addition model which describes the situation in which clusters can only grow by absorbing monomers.

Variational data assimilation for very large environmental problems

Variational data assimilation is commonly used in environmental forecasting to estimate the current state of the system from a model forecast and observational data. The assimilation problem can be written simply in the form of a nonlinear least squares optimization problem. However the practical solution of the problem in large systems requires many careful choices to be made in the implementation. In this article we present the theory of variational data assimilation and then discuss in detail how it is implemented in practice. Current solutions and open questions are discussed.