On Hamiltonian balanced dynamics and the slowest invariant manifold

The concept of a slowest invariant manifold is investigated for the five-component model of Lorenz under conservative dynamics. It is shown that Lorenz's model is a two-degree-of-freedom canonical Hamiltonian system, consisting of a nonlinear vorticity-triad oscillator coupled to a linear gravity wave oscillator, whose solutions consist of regular and chaotic orbits. When either the Rossby number or the rotational Froude number is small, there is a formal separation of timescales, and one can speak of fast and slow motion. In the same regime, the coupling is weak, and the Kolmogorov–Arnold-Moser theorem is shown to apply. The chaotic orbits are inherently unbalanced and are confined to regions sandwiched between invariant tori consisting of quasi-periodic regular orbits. The regular orbits generally contain free fast motion, but a slowest invariant manifold may be geometrically defined as the set of all slow cores of invariant tori (defined by zero fast action) that are smoothly related to such cores in the uncoupled system. This slowest invariant manifold is not global; in fact, its structure is fractal; but it is of nearly full measure in the limit of weak coupling. It is also nonlinearly stable. As the coupling increases, the slowest invariant manifold shrinks until it disappears altogether. The results clarify previous definitions of a slowest invariant manifold and highlight the ambiguity in the definition of “slowness.” An asymptotic procedure, analogous to standard initialization techniques, is found to yield nonzero free fast motion even when the core solutions contain none. A hierarchy of Hamiltonian balanced models preserving the symmetries in the original low-order model is formulated; these models are compared with classic balanced models, asymptotically initialized solutions of the full system and the slowest invariant manifold defined by the core solutions. The analysis suggests that for sufficiently small Rossby or rotational Froude numbers, a stable slowest invariant manifold can be defined for this system, which has zero free gravity wave activity, but it cannot be defined everywhere. The implications of the results for more complex systems are discussed.

Antarctic ocean and sea ice response to ozone depletion: a two timescale problem

The response of the Southern Ocean to a repeating seasonal cycle of ozone loss is studied in two coupled climate models and found to comprise both fast and slow processes. The fast response is similar to the inter-annual signature of the Southern Annular Mode (SAM) on Sea Surface Temperature (SST), on to which the ozone-hole forcing projects in the summer. It comprises enhanced northward Ekman drift inducing negative summertime SST anomalies around Antarctica, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year-round. The enhanced northward Ekman drift, however, results in upwelling of warm waters from below the mixed layer in the region of seasonal sea ice. With sustained bursts of westerly winds induced by ozone-hole depletion, this warming from below eventually dominates over the cooling from anomalous Ekman drift. The resulting slow-timescale response (years to decades) leads to warming of SSTs around Antarctica and ultimately a reduction in sea-ice cover year-round. This two-timescale behavior - rapid cooling followed by slow but persistent warming - is found in the two coupled models analysed, one with an idealized geometry, the other a complex global climate model with realistic geometry. Processes that control the timescale of the transition from cooling to warming, and their uncertainties are described. Finally we discuss the implications of our results for rationalizing previous studies of the effect of the ozone-hole on SST and sea-ice extent. %Interannual variability in the Southern Annular Mode (SAM) and sea ice covary such that an increase and southward shift in the surface westerlies (a positive phase of the SAM) coincides with a cooling of Sea Surface Temperature (SST) around 70-50$^\circ$S and an expansion of the sea ice cover, as seen in observations and models alike. Yet, in modeling studies, the Southern Ocean warms and sea ice extent decreases in response to sustained, multi-decadal positive SAM-like wind anomalies driven by 20th century ozone depletion. Why does the Southern Ocean appear to have disparate responses to SAM-like variability on interannual and multidecadal timescales? Here it is demonstrated that the response of the Southern Ocean to ozone depletion has a fast and a slow response. The fast response is similar to the interannual variability signature of the SAM. It is dominated by an enhanced northward Ekman drift, which transports heat northward and causes negative SST anomalies in summertime, earlier sea ice freeze-up the following winter, and northward expansion of the sea ice edge year round. The enhanced northward Ekman drift causes a region of Ekman divergence around 70-50$^\circ$S, which results in upwelling of warmer waters from below the mixed layer. With sustained westerly wind enhancement in that latitudinal band, the warming due to the anomalous upwelling of warm waters eventually dominates over the cooling from the anomalous Ekman drift. Hence, the slow response ultimately results in a positive SST anomaly and a reduction in the sea ice cover year round. We demonstrate this behavior in two models: one with an idealized geometry and another, more detailed, global climate model. However, the models disagree on the timescale of transition from the fast (cooling) to the slow (warming) response. Processes that controls this transition and their uncertainties are discussed.

The quasi-nondispersive regimes of long extratropical baroclinic rossby waves over (slowly varying) topography

Actual energy paths of long, extratropical baroclinic Rossby waves in the ocean are difficult to describe simply because they depend on the meridional-wavenumber-to-zonal-wavenumber ratio tau, a quantity that is difficult to estimate both observationally and theoretically. This paper shows, however, that this dependence is actually weak over any interval in which the zonal phase speed varies approximately linearly with tau, in which case the propagation becomes quasi-nondispersive (QND) and describable at leading order in terms of environmental conditions (i.e., topography and stratification) alone. As an example, the purely topographic case is shown to possess three main kinds of QND ray paths. The first is a topographic regime in which the rays follow approximately the contours f/h(alpha c) = a constant (alpha(c) is a near constant fixed by the strength of the stratification, f is the Coriolis parameter, and h is the ocean depth). The second and third are, respectively, "fast" and "slow" westward regimes little affected by topography and associated with the first and second bottom-pressure-compensated normal modes studied in previous work by Tailleux and McWilliams. Idealized examples show that actual rays can often be reproduced with reasonable accuracy by replacing the actual dispersion relation by its QND approximation. The topographic regime provides an upper bound ( in general a large overestimate) of the maximum latitudinal excursions of actual rays. The method presented in this paper is interesting for enabling an optimal classification of purely azimuthally dispersive wave systems into simpler idealized QND wave regimes, which helps to rationalize previous empirical findings that the ray paths of long Rossby waves in the presence of mean flow and topography often seem to be independent of the wavenumber orientation. Two important side results are to establish that the baroclinic string function regime of Tyler and K se is only valid over a tiny range of the topographic parameter and that long baroclinic Rossby waves propagating over topography do not obey any two-dimensional potential vorticity conservation principle. Given the importance of the latter principle in geophysical fluid dynamics, the lack of it in this case makes the concept of the QND regimes all the more important, for they are probably the only alternative to provide a simple and economical description of general purely azimuthally dispersive wave systems.

Probing the large-scale topology of the heliospheric magnetic field using Jovian electrons

Jupiter’s magnetosphere acts as a point source of near-relativistic electrons within the heliosphere. In this study, three solar cycles of Jovian electron data in near-Earth space are examined. Jovian electron intensity is found to peak for an ideal Parker spiral connection, but with considerable spread about this point. Assuming the peak in Jovian electron counts indicates the best magnetic connection to Jupiter, we find a clear trend for fast and slow solar wind to be over- and under-wound with respect to the ideal Parker spiral, respectively. This is shown to be well explained in terms of solar wind stream interactions. Thus, modulation of Jovian electrons by corotating interaction regions (CIRs) may primarily be the result of changing magnetic connection, rather than CIRs acting as barriers to cross-field diffusion. By using Jovian electrons to remote sensing magnetic connectivity with Jupiter’s magnetosphere, we suggest that they provide a means to validate solar wind models between 1 and 5 AU, even when suitable in situ solar wind observations are not available. Furthermore, using Jovian electron observations as probes of heliospheric magnetic topology could provide insight into heliospheric magnetic field braiding and turbulence, as well as any systematic under-winding of the heliospheric magnetic field relative to the Parker spiral from footpoint motion of the magnetic field.

Understanding electron heat flux signatures in the solar wind

Suprathermal electrons (E > 80 eV) carry heat flux away from the Sun. Processes controlling the heat flux are not well understood. To gain insight into these processes, we model heat flux as a linear dependence on two independent parameters: electron number flux and electron pitch angle anisotropy. Pitch angle anisotropy is further modeled as a linear dependence on two solar wind components: magnetic field strength and plasma density. These components show no correlation with number flux, reinforcing its independence from pitch angle anisotropy. Multiple linear regression applied to 2 years of Wind data shows good correspondence between modeled and observed heat flux and anisotropy. The results suggest that the interplay of solar wind parameters and electron number flux results in distinctive heat flux dropouts at heliospheric features like plasma sheets but that these parameters continuously modify heat flux. This is inconsistent with magnetic disconnection as the primary cause of heat flux dropouts. Analysis of fast and slow solar wind regimes separately shows that electron number flux and pitch angle anisotropy are equally correlated with heat flux in slow wind but that number flux is the dominant correlative in fast wind. Also, magnetic field strength correlates better with pitch angle anisotropy in slow wind than in fast wind. The energy dependence of the model fits suggests different scattering processes in fast and slow wind.

Shear banding in molecular dynamics of polymer melts

In order to establish constitutive equations for a viscoelastic fluid uniform shear flow is usually required. However, in the last 10 years S. Q. Wang and co-workers have demonstrated that some entangled polymers do not flow with the uniform shear rate as usually assumed, but instead choose to separate into fast and slow flowing regions. This phenomenon, known as shear banding, causes flow instabilities and in principle invalidates all rheological measurements when it occurs. In this Letter we report the first observation of shear banding in molecular dynamics simulations of entangled polymer melts. We show that our observations are in a very good agreement with the phenomenology developed by Fielding and Olmsted. Our findings provide a simple way of validating the empirical macroscopic phenomenology of shear banding. © 2012 American Physical Society

A low-order model investigation of the analysis of gravity waves in the ensemble Kalman filter.

The behavior of the ensemble Kalman filter (EnKF) is examined in the context of a model that exhibits a nonlinear chaotic (slow) vortical mode coupled to a linear (fast) gravity wave of a given amplitude and frequency. It is shown that accurate recovery of both modes is enhanced when covariances between fast and slow normal-mode variables (which reflect the slaving relations inherent in balanced dynamics) are modeled correctly. More ensemble members are needed to recover the fast, linear gravity wave than the slow, vortical motion. Although the EnKF tends to diverge in the analysis of the gravity wave, the filter divergence is stable and does not lead to a great loss of accuracy. Consequently, provided the ensemble is large enough and observations are made that reflect both time scales, the EnKF is able to recover both time scales more accurately than optimal interpolation (OI), which uses a static error covariance matrix. For OI it is also found to be problematic to observe the state at a frequency that is a subharmonic of the gravity wave frequency, a problem that is in part overcome by the EnKF.However, error in themodeled gravity wave parameters can be detrimental to the performance of the EnKF and remove its implied advantages, suggesting that a modified algorithm or a method for accounting for model error is needed.

Balance model for equatorial long waves

Geophysical fluid models often support both fast and slow motions. As the dynamics are often dominated by the slow motions, it is desirable to filter out the fast motions by constructing balance models. An example is the quasi geostrophic (QG) model, which is used widely in meteorology and oceanography for theoretical studies, in addition to practical applications such as model initialization and data assimilation. Although the QG model works quite well in the mid-latitudes, its usefulness diminishes as one approaches the equator. Thus far, attempts to derive similar balance models for the tropics have not been entirely successful as the models generally filter out Kelvin waves, which contribute significantly to tropical low-frequency variability. There is much theoretical interest in the dynamics of planetary-scale Kelvin waves, especially for atmospheric and oceanic data assimilation where observations are generally only of the mass field and thus do not constrain the wind field without some kind of diagnostic balance relation. As a result, estimates of Kelvin wave amplitudes can be poor. Our goal is to find a balance model that includes Kelvin waves for planetary-scale motions. Using asymptotic methods, we derive a balance model for the weakly nonlinear equatorial shallow-water equations. Specifically we adopt the ‘slaving’ method proposed by Warn et al. (Q. J. R. Meteorol. Soc., vol. 121, 1995, pp. 723–739), which avoids secular terms in the expansion and thus can in principle be carried out to any order. Different from previous approaches, our expansion is based on a long-wave scaling and the slow dynamics is described using the height field instead of potential vorticity. The leading-order model is equivalent to the truncated long-wave model considered previously (e.g. Heckley & Gill, Q. J. R. Meteorol. Soc., vol. 110, 1984, pp. 203–217), which retains Kelvin waves in addition to equatorial Rossby waves. Our method allows for the derivation of higher-order models which significantly improve the representation of Rossby waves in the isotropic limit. In addition, the ‘slaving’ method is applicable even when the weakly nonlinear assumption is relaxed, and the resulting nonlinear model encompasses the weakly nonlinear model. We also demonstrate that the method can be applied to more realistic stratified models, such as the Boussinesq model.

The heliospheric magnetic field

The heliospheric magnetic field (HMF) is the extension of the coronal magnetic field carried out into the solar system by the solar wind. It is the means by which the Sun interacts with planetary magnetospheres and channels charged particles propagating through the heliosphere. As the HMF remains rooted at the solar photosphere as the Sun rotates, the large-scale HMF traces out an Archimedean spiral. This pattern is distorted by the interaction of fast and slow solar wind streams, as well as the interplanetary manifestations of transient solar eruptions called coronal mass ejections. On the smaller scale, the HMF exhibits an array of waves, discontinuities, and turbulence, which give hints to the solar wind formation process. This review aims to summarise observations and theory of the small- and large-scale structure of the HMF. Solar-cycle and cycle-to-cycle evolution of the HMF is discussed in terms of recent spacecraft observations and pre-spaceage proxies for the HMF in geomagnetic and galactic cosmic ray records.

Impact of mass lumping on gravity and Rossby waves in 2D finite-element shallow-water models

The goal of this study is to evaluate the effect of mass lumping on the dispersion properties of four finite-element velocity/surface-elevation pairs that are used to approximate the linear shallow-water equations. For each pair, the dispersion relation, obtained using the mass lumping technique, is computed and analysed for both gravity and Rossby waves. The dispersion relations are compared with those obtained for the consistent schemes (without lumping) and the continuous case. The P0-P1, RT0 and P-P1 pairs are shown to preserve good dispersive properties when the mass matrix is lumped. Test problems to simulate fast gravity and slow Rossby waves are in good agreement with the analytical results.

Magnetic cloud distortion resulting from propagation through a structured solar wind: Models and observations

Numerical simulations of magnetic clouds (MCs) propagating through a structured solar wind suggest that MC-associated magnetic flux ropes are highly distorted by inhomogeneities in the ambient medium. In particular, a solar wind configuration of fast wind from high latitudes and slow wind at low latitudes, common at periods close to solar minimum, should distort the cross section of magnetic clouds into concave-outward structures. This phenomenon has been reported in observations of shock front orientations, but not in the body of magnetic clouds. In this study an analytical magnetic cloud model based upon a kinematically distorted flux rope is modified to simulate propagation through a structured medium. This new model is then used to identify specific time series signatures of the resulting concave-outward flux ropes. In situ observations of three well studied magnetic clouds are examined with comparison to the model, but the expected concave-outward signatures are not present. Indeed, the observations are better described by the convex-outward flux rope model. This may be due to a sharp latitudinal transition from fast to slow wind, resulting in a globally concave-outward flux rope, but with convex-outward signatures on a local scale.

Classification-based probabilistic modeling of texture transition for fast line search tracking and delineation

We introduce a classification-based approach to finding occluding texture boundaries. The classifier is composed of a set of weak learners, which operate on image intensity discriminative features that are defined on small patches and are fast to compute. A database that is designed to simulate digitized occluding contours of textured objects in natural images is used to train the weak learners. The trained classifier score is then used to obtain a probabilistic model for the presence of texture transitions, which can readily be used for line search texture boundary detection in the direction normal to an initial boundary estimate. This method is fast and therefore suitable for real-time and interactive applications. It works as a robust estimator, which requires a ribbon-like search region and can handle complex texture structures without requiring a large number of observations. We demonstrate results both in the context of interactive 2D delineation and of fast 3D tracking and compare its performance with other existing methods for line search boundary detection.

Pesticides and other micro-organic contaminants in freshwater sedimentary environments - a review

A wide range of issues relating to the presence and fate of pesticides and other micro-organic contaminants (MOCs) in surface freshwater sedimentary environments is reviewed. These issues include the sources, transport and occurrence of MOCs in freshwater environments; their ecological effects; their interaction with sedimentary material; and a range of processes related to their fate, including degradation, diffusion in bed sediments, bioturbation and slow contaminant release. An emphasis is placed on those processes-chemical, physical or biological-in which sediments play a role in determining the fate of micro-organics in freshwater environments. The issues of occurrence, source and transport, and the ecological effects of micro-organics are introduced more briefly, the focus where these aspects are concerned being largely on pesticides. In the concluding section, key points and issues relating to the study of micro-organics in freshwater environments are summarised and areas where initial or further research would be welcome are highlighted. It is hoped that this paper will both form a useful reference for workers in the field of micro-organic contaminants, and also stimulate new work in the freshwater environment and beyond. (C) 2003 Elsevier Science Ltd. All rights reserved.

Pesticides and other micro-organic contaminants in freshwater sedimentary environments - a review

A wide range of issues relating to the presence and fate of pesticides and other micro-organic contaminants (MOCs) in surface freshwater sedimentary environments is reviewed. These issues include the sources, transport and occurrence of MOCs in freshwater environments; their ecological effects; their interaction with sedimentary material; and a range of processes related to their fate, including degradation, diffusion in bed sediments, bioturbation and slow contaminant release. An emphasis is placed on those processes-chemical, physical or biological-in which sediments play a role in determining the fate of micro-organics in freshwater environments. The issues of occurrence, source and transport, and the ecological effects of micro-organics are introduced more briefly, the focus where these aspects are concerned being largely on pesticides. In the concluding section, key points and issues relating to the study of micro-organics in freshwater environments are summarised and areas where initial or further research would be welcome are highlighted. It is hoped that this paper will both form a useful reference for workers in the field of micro-organic contaminants, and also stimulate new work in the freshwater environment and beyond. (C) 2003 Elsevier Science Ltd. All rights reserved.

A surface energy perspective on climate change

A surface forcing response framework is developed that enables an understanding of time-dependent climate change from a surface energy perspective. The framework allows the separation of fast responses that are unassociated with global-mean surface air temperature change (ΔT), which is included in the forcing, and slow feedbacks that scale with ΔT. The framework is illustrated primarily using 2 × CO2 climate model experiments and is robust across the models. For CO2 increases, the positive downward radiative component of forcing is smaller at the surface than at the tropopause, and so a rapid reduction in the upward surface latent heat (LH) flux is induced to conserve the tropospheric heat budget; this reduces the precipitation rate. Analysis of the time-dependent surface energy balance over sea and land separately reveals that land areas rapidly regain energy balance, and significant land surface warming occurs before global sea temperatures respond. The 2 × CO2 results are compared to a solar increase experiment and show that some fast responses are forcing dependent. In particular, a significant forcing from the fast hydrological response found in the CO2 experiments is much smaller in the solar experiment. The different fast response explains why previous equilibrium studies found differences in the hydrological sensitivity between these two forcings. On longer time scales, as ΔT increases, the net surface longwave and LH fluxes provide positive and negative surface feedbacks, respectively, while the net surface shortwave and sensible heat fluxes change little. It is found that in contrast to their fast responses, the longer-term response of both surface energy fluxes and the global hydrological cycle are similar for the different forcing agents.