Characterising the variability and extremes of the stratospheric polar vortices using 2D moment analysis

The mean state, variability and extreme variability of the stratospheric polar vortices, with an emphasis on the Northern Hemisphere vortex, are examined using 2-dimensional moment analysis and Extreme Value Theory (EVT). The use of moments as an analysis to ol gives rise to information about the vortex area, centroid latitude, aspect ratio and kurtosis. The application of EVT to these moment derived quantaties allows the extreme variability of the vortex to be assessed. The data used for this study is ECMWF ERA-40 potential vorticity fields on interpolated isentropic surfaces that range from 450K-1450K. Analyses show that the most extreme vortex variability occurs most commonly in late January and early February, consistent with when most planetary wave driving from the troposphere is observed. Composites around sudden stratospheric warming (SSW) events reveal that the moment diagnostics evolve in statistically different ways between vortex splitting events and vortex displacement events, in contrast to the traditional diagnostics. Histograms of the vortex diagnostics on the 850K (∼10hPa) surface over the 1958-2001 period are fitted with parametric distributions, and show that SSW events comprise the majority of data in the tails of the distributions. The distribution of each diagnostic is computed on various surfaces throughout the depth of the stratosphere, and shows that in general the vortex becomes more circular with higher filamentation at the upper levels. The Northern Hemisphere (NH) and Southern Hemisphere (SH) vortices are also compared through the analysis of their respective vortex diagnostics, and confirm that the SH vortex is less variable and lacks extreme events compared to the NH vortex. Finally extreme value theory is used to statistically mo del the vortex diagnostics and make inferences about the underlying dynamics of the polar vortices.

Long-term acclimation of leaf production, development, longevity and quality following 3 yr exposure to free-air CO2 enrichment during canopy closure in Populus

The effects of elevated CO2 on leaf development in three genotypes of Populus were investigated during canopy closure, following exposure to elevated CO2 over 3 yr using free-air enrichment.• Leaf quality was altered such that nitrogen concentration per unit d. wt (Nmass) declined on average by 22 and 13% for sun and shade leaves, respectively, in elevated CO2. There was little evidence that this was the result of ‘dilution’ following accumulation of nonstructural carbohydrates. Most likely, this was the result of increased leaf thickness. Specific leaf area declined in elevated CO2 on average by 29 and 5% for sun and shade leaves, respectively.• Autumnal senescence was delayed in elevated CO2 with a 10% increase in the number of days at which 50% leaf loss occurred in elevated as compared with ambient CO2.• These data suggest that changes in leaf quality may be predicted following long-term acclimation of fast-growing forest trees to elevated CO2, and that canopy longevity may increase, with important implications for forest productivity.

Bistability through triadic closure

We propose and analyse a class of evolving network models suitable for describing a dynamic topological structure. Applications include telecommunication, on-line social behaviour and information processing in neuroscience. We model the evolving network as a discrete time Markov chain, and study a very general framework where, conditioned on the current state, edges appear or disappear independently at the next timestep. We show how to exploit symmetries in the microscopic, localized rules in order to obtain conjugate classes of random graphs that simplify analysis and calibration of a model. Further, we develop a mean field theory for describing network evolution. For a simple but realistic scenario incorporating the triadic closure effect that has been empirically observed by social scientists (friends of friends tend to become friends), the mean field theory predicts bistable dynamics, and computational results confirm this prediction. We also discuss the calibration issue for a set of real cell phone data, and find support for a stratified model, where individuals are assigned to one of two distinct groups having different within-group and across-group dynamics.

Harringtonian virtue: Harrington, Machiavelli, and the method of the 'Moment'

This article presents a reinterpretation of James Harrington's writings. It takes issue with J. G. A. Pocock's reading, which treats him as importing into England a Machiavellian ‘language of political thought’. This reading is the basis of Pocock's stress on the republicanism of eighteenth-century opposition values. Harrington's writings were in fact a most implausible channel for such ideas. His outlook owed much to Stoicism. Unlike the Florentine, he admired the contemplative life; was sympathetic to commerce; and was relaxed about the threat of ‘corruption’ (a concept that he did not understand). These views can be associated with his apparent aims: the preservation of a national church with a salaried but politically impotent clergy; and the restoration of the royalist gentry to a leading role in English politics. Pocock's hypothesis is shown to be conditioned by his method; its weaknesses reflect some difficulties inherent in the notion of ‘languages of thought’.

First order k-th moment finite element analysis of nonlinear operator equations with stochastic data

We develop and analyze a class of efficient Galerkin approximation methods for uncertainty quantification of nonlinear operator equations. The algorithms are based on sparse Galerkin discretizations of tensorized linearizations at nominal parameters. Specifically, we consider abstract, nonlinear, parametric operator equations J(\alpha ,u)=0 for random input \alpha (\omega ) with almost sure realizations in a neighborhood of a nominal input parameter \alpha _0. Under some structural assumptions on the parameter dependence, we prove existence and uniqueness of a random solution, u(\omega ) = S(\alpha (\omega )). We derive a multilinear, tensorized operator equation for the deterministic computation of k-th order statistical moments of the random solution's fluctuations u(\omega ) - S(\alpha _0). We introduce and analyse sparse tensor Galerkin discretization schemes for the efficient, deterministic computation of the k-th statistical moment equation. We prove a shift theorem for the k-point correlation equation in anisotropic smoothness scales and deduce that sparse tensor Galerkin discretizations of this equation converge in accuracy vs. complexity which equals, up to logarithmic terms, that of the Galerkin discretization of a single instance of the mean field problem. We illustrate the abstract theory for nonstationary diffusion problems in random domains.

3D cloud reconstructions: evaluation of scanning radar scan strategy with a view to surface shortwave radiation closure

The ability of six scanning cloud radar scan strategies to reconstruct cumulus cloud fields for radiation study is assessed. Utilizing snapshots of clean and polluted cloud fields from large eddy simulations, an analysis is undertaken of error in both the liquid water path and monochromatic downwelling surface irradiance at 870 nm of the reconstructed cloud fields. Error introduced by radar sensitivity, choice of radar scan strategy, retrieval of liquid water content (LWC), and reconstruction scheme is explored. Given an in␣nitely sensitive radar and perfect LWC retrieval, domain average surface irradiance biases are typically less than 3 W m␣2 ␣m␣1, corresponding to 5–10% of the cloud radiative effect (CRE). However, when using a realistic radar sensitivity of ␣37.5 dBZ at 1 km, optically thin areas and edges of clouds are dif␣cult to detect due to their low radar re-ectivity; in clean conditions, overestimates are of order 10 W m␣2 ␣m␣1 (~20% of the CRE), but in polluted conditions, where the droplets are smaller, this increases to 10–26 W m␣2 ␣m␣1 (~40–100% of the CRE). Drizzle drops are also problematic; if treated as cloud droplets, reconstructions are poor, leading to large underestimates of 20–46 W m␣2 ␣m␣1 in domain average surface irradiance (~40–80% of the CRE). Nevertheless, a synergistic retrieval approach combining the detailed cloud structure obtained from scanning radar with the droplet-size information and location of cloud base gained from other instruments would potentially make accurate solar radiative transfer calculations in broken cloud possible for the first time.

The full infinite dimensional moment problem on semi-algebraic sets of generalized functions

We consider a generic basic semi-algebraic subset S of the space of generalized functions, that is a set given by (not necessarily countably many) polynomial constraints. We derive necessary and sufficient conditions for an infinite sequence of generalized functions to be realizable on S, namely to be the moment sequence of a finite measure concentrated on S. Our approach combines the classical results about the moment problem on nuclear spaces with the techniques recently developed to treat the moment problem on basic semi-algebraic sets of Rd. In this way, we determine realizability conditions that can be more easily verified than the well-known Haviland type conditions. Our result completely characterizes the support of the realizing measure in terms of its moments. As concrete examples of semi-algebraic sets of generalized functions, we consider the set of all Radon measures and the set of all the measures having bounded Radon–Nikodym density w.r.t. the Lebesgue measure.

Reconfiguration and closure of lobe flux by reconnection during northward IMF: possible evidence for signatures in cusp/cleft auroral emissions

Observations are presented of the response of the dayside cusp/cleft aurora to changes in both the clock and elevation angles of the interplanetary magnetic field (IMF) vector, as monitored by the WIND spacecraft. The auroral observations are made in 630 nm light at the winter solstice near magnetic noon, using an all-sky camera and a meridian-scanning photometer on the island of Spitsbergen. The dominant change was the response to a northward turning of the IMF which caused a poleward retreat of the dayside aurora. A second, higher-latitude band of aurora was seen to form following the northward turning, which is interpreted as the effect of lobe reconnection which reconfigures open flux. We suggest that this was made possible in the winter hemisphere, despite the effect of the Earth's dipole tilt, by a relatively large negative X component of the IMF. A series of five events then formed in the poleward band and these propagated in a southwestward direction and faded at the equatorward edge of the lower-latitude band as it migrated poleward. It is shown that the auroral observations are consistent with overdraped lobe flux being generated by lobe reconnection in the winter hemisphere and subsequently being re-closed by lobe reconnection in the summer hemisphere. We propose that the balance between the reconnection rates at these two sites is modulated by the IMF elevation angle, such that when the IMF points more directly northward, the summer lobe reconnection site dominates, re-closing all overdraped lobe flux and eventually becoming disconnected from the Northern Hemisphere.

A moment problem for random discrete measures

Let X be a locally compact Polish space. A random measure on X is a probability measure on the space of all (nonnegative) Radon measures on X. Denote by K(X) the cone of all Radon measures η on X which are of the form η =