Delay coordinates: a sensitive indicator of nonlinear dynamics in single charged particle motion in magnetic reversals

The delay coordinate technique is examined as an indicator of the regime of particle dynamics for the system of single charged particle motion in magnetic reversals. Examples of numerically integrated trajectories in both static (zero electric field) and time dependent (corresponding nonzero induction electric field) simple models for magnetic reversals are considered. In the static case, the dynamics can in principle be directly classified by constructing Poincaré surfaces of section; here we demonstrate that whilst the Poincaré surface contains the relevant information to classify the dynamics, the corresponding delay coordinate plot can provide a far more sensitive indication of the onset of nonregular behaviour. In the case of nonperiodic time dependence considered here Poincaré plots cannot in general be constructed directly. Nevertheless, delay coordinate plots can still reveal details of the phase space portrait of the system, and here are shown to indicate whether segments of stochastic motion exist in a given trajectory. It is anticipated that the delay coordinate plot technique as realized here will be a valuable tool in characterizing the behaviour in large numbers of trajectories that are evolved in time-dependent systems, thereby giving us insight into the evolution of the distribution function as a whole, either in prescribed fields or in self-consistent numerical simulations.     

A chemical survey of the Mississippi estuary

A “snap shot” survey of the Mississippi estuary was made during a period of low river discharge, when the estuarine mixing zone was within the deltaic channels. Concentrations of H⁺, Ca²⁺, inorganic phosphorus and inorganic carbon suggest that the waters of the river and the low salinity (<5‰) portion of the estuary are near saturation with respect to calcite and sedimentary calcium phosphate. An input of oxidized nitrogen species and N₂O was observed in the estuary between 0 and 4‰ salinity. The concentrations of dissolved NH₄ ⁺ and O₂, over most of the estuary, appeared to be influenced by decomposition of terrestrial organic matter in bottom sediments. The estuarine bottom also appears to be a source of CH₄ which has been suggested to originate from petroleum shipping and refining operations. Estuarine mixing with offshore Gulf waters was the dominant influence on distributions of dissolved species over most of the estuary (i.e., from salinities >5‰). The phytoplankton abundance (measured as chlorophylla) increased as the depth of the mixed layer decreased in a manner consistent with that expected for a light-limited ecosystem. Fluxes of NO₃ ^?+NO₂ ^? and soluble inorganic phosphorus to the Gulf of Mexico were estimated to be 3.4±0.2×10³ g N s^?1 and 1.9±0.2 g P s^?1 respectively, at the time of this study.     

Sorption of noble gases by solids,with reference to meteorites. I. Magnetite and carbon

To simulate trapping of meteoritic noble gases by solids, 18 samples of Fe₃O₄ were synthesized in a noble gas atmosphere at 350–720 K by the reactions: 3Fe + 4H₂O → Fe₃O₄ + 4H₂ (Ne, Ar, Kr, Xe) 3Fe + 4CO → Fe₃O₃ + 4C + carbides (Xe only) Phases were separated by selective solvents (HgCl₂, HCl). Noble gas contents were analyzed by mass spectrometry, or, in runs where 36 d Xe¹²⁷ tracer was used, by γ-counting. Surface areas, as measured by the BET method, ranged from 1 to 400 m²/g. Isotopic fractionations were below the detection limit of 0.5%/m.u.Sorption of Xe on Fe₃O₄ and C obeys Henry's Law between $\text{1 × 10}{}^{\mathrm{?8}}$ and $\text{4 × 10}{}^{\mathrm{?5}}$ atm, but shows only a slight temperature dependence between 650 and 720 K ($\text{ΔH}{}_{\mathrm{sol}}\text{= ?4 ± 2}\text{kcal/mole}$). The mean distribution coefficient K_Xe is $\text{0.28 ± 0.09}$ cc STP/g atm for Fe₃O₄ and only a factor of $\text{1.2 ± 0.4}$ greater for C; such similarity for two cogenetic phases was predicted by Lewis et al. (1977). Stepped heating and etching experiments show that 20–50% of the total Xe is physically adsorbed and about 20% is trapped in the solid. The rest is chemisorbed with $\text{ΔH}{}_{\mathrm{s}}\text{? ?13}\text{kcal/mole}$. The desorption or exchange half-time for the last two components is >10² yr at room temperature.Etching experiments showed a possible analogy to “Phase $\text{Q}$” in meteorites. A typical carbon + carbide sample, when etched with HNO₃, lost 47% of its Xe but only 0.9% of its mass, corresponding to a ~0.6 Å layer. Though this etchable, surficial gas component was more thermolabile than $\text{Q}$ (release $\text{T}$ below 1000°C, compared to 1200–1600°C), another experiment shows that the proportion of chemisorbed Xe increases upon moderate heating (1 hr at 450°C). Apparently adsorbed gases can become “fixed” to the crystal, by processes not involving volume diffusion (recrystallization, chemical reaction, migration to traps, etc.). Such mechanisms may have acted in the solar nebula, to strengthen the binding of adsorbed gases.Adsorbed atmospheric noble gases are present in all samples, and dominate whenever the noble gas partial pressure in the atmosphere is greater than that in the synthesis. Many of the results of Lancet and Anders (1973) seem to have been dominated by such an atmospheric component; others are suspect for other reasons, whereas still others seem reliable. When the doubtful samples of Lancet and Anders are eliminated or corrected, the fractionation pattern—as in our samples—no longer peaks at Ar, but rises monotonically from Ne to Xe. No clear evidence remains for the strong temperature dependence claimed by these authors.     

Water movement within lac du bonnet batholith as revealed by detailed thermal studies of three closely-spaced boreholes

Successive temperature logs have been obtained over a period of two years in three closely-spaced boreholes in the Lac du Bonnet batholith of the Superior Province of the Canadian Shield. Two of the boreholes, of depth 450 m and 830 m, intersect a dipping fracture zone at 435–450 m. In both holes water is flowing from near the surface to the fracture zone at approximately 1.5–1.9·10⁻⁵ m³ s⁻¹, the flow being inferred from analysis of the temperature logs. Below 25 m, temperatures in these two holes are 0.22–0.28 K lower than those in the third, 145 m, hole.The temperature data have been combined with over 200 thermal conductivity measurements on core samples to produce heat flow values. In the deepest hole heat flow above the fracture zone is 16% higher than that below the zone. This indicates that water is flowing up the fracture zone. The flow rate is approximately 0.3 g s⁻¹ m⁻¹, and the flow has existed for thousands of years.Observation of thermal effects of water flow in massive, relatively unfractured plutons in a region having little topographic relief causes one to be concerned about the reliability of heat flow values measured in similar environments.The regional heat flow is taken to be 50 mW m⁻² after correction for glaciation effects. The average value of 24 determinations of radioactive heat generation in granitic core samples is 5.23 ± 1.11 μW m⁻³, which is more than three times higher than expected for such a heat flow in the Superior Province. This implies that the layer of high radioactive heat generation is thin, being not more than 4 km and probably about 1.3 km thick.     

Impact of prescribed burning on Acacia nilotica seed banks and seedlings in the Astrebla grasslands of northern Australia

Acacia nilotica is a tree of international significance both as a beneficial plant and as a species prone to thicket formation and negative impacts on savannas throughout much of its range. While fire has been identified as a useful tool for controlling negative impacts of some Acacia species in Africa, A. nilotica adult trees are apparently fire tolerant. The effects of fire on seed bank mortality, viability and germination, and on juvenile plants of this species are unknown. We undertook a study to investigate the impacts of fire on seeds and first and second year seedlings in Australian Astrebla savannas. The primary objective was to determine the effects of control burns applied during the early stages of invasion and thicket formation. Because of high rates of seedling dessication post fire in 1997, supplementary watering was undertaken until fire mortality could be measured the following year. High seed mortality occurred for seeds located either in cattle dung (100%) or on the soil surface (c. 80%). Germination of those seeds that survived was enhanced following exposure to fire. Very high fire mortality (99%) was observed among first year seedlings and moderate mortality (60%) for second year seedlings even after a low intensity burn. Negative impacts on major pasture species were minimal. The implications of these findings are discussed in relation to A. nilotica fire ecology and to the efficacy of fire as a tool for managing this species.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Distorted Froude‐scaled flume analysis of large woody debris

This paper presents the results of a movable‐boundary, distorted, Froude‐scaled hydraulic model based on Abiaca Creek, a sand‐bedded channel in northern Mississippi. The model was used to examine the geomorphic and hydraulic impact of simplified large woody debris (LWD) elements. The theory of physical scale models is discussed and the method used to construct the LWD test channel is developed. The channel model had bed and banks moulded from 0·8 mm sand, and flow conditions were just below the threshold of motion so that any sediment transport and channel adjustment were the result of the debris element. Dimensions and positions of LWD elements were determined using a debris jam classification model. Elements were attached to a dynamometer to measure element drag forces, and channel adjustment was determined through detailed topographic surveys. The fluid drag force on the elements decreased asymptotically over time as the channel boundary eroded around the elements due to locally increased boundary shear stress. Total time for geomorphic adjustment computed for the prototype channel at the Q₂ discharge (discharge occurring once every two years on average) was as short as 45 hours. The size, depth and position of scour holes, bank erosion and bars created by flow acceleration past the elements were found to be related to element length and position within the channel cross‐section. Morphologies created by each debris element in the model channel were comparable with similar jams observed in the prototype channel. Published in 2001 John Wiley & Sons, Ltd.     

Interactive Web-Based Map: Applications to Large Data Sets in the Geosciences

Current advances in computer hardware, information technology and data collection techniques have produced very large data sets in a wide variety of scientific and engineering disciplines. We must harness this opportunity to visualize and extract useful information from geophysical and geological data. We have taken the task of data mining by implementing a map-like approach over a web server for interrogating the data, using a client-server paradigm. The spatial-data is mapped onto a two-dimensional grid from which the user (client) can quiz the data with the map-interface as a user extension. The data is stored on the server, while the computational gateway separating the client and the server can be the front-end of an electronic publication, electronic classroom, a survey, or an e-business. We have used a combination of Java, Java3D, and Perl for processing the data and communicating between the client and the server. The user can interrogate the geospatial data over any particular region with arbitrary dimensions and then receive back relevant statistical analysis, such as the histogram plots and local statistics. We have applied this method for the following data sets: (1.) distribution of prime numbers (2.) two-dimensional mantle convection (3.) three-dimensional mantle convection (4.) high-resolution satellite reflectance data over multiple wavelengths (5.) molecular dynamics describing the flow of blood in narrow vessels. Using the map-interface, one can actually interrogate this data over the Internet.     

Interactions of growing folds and coeval depositional systems

Responses of both modern and ancient fluvial depositional systems to growing folds can be interpreted in terms of interactions among competing controlling variables which can be incorporated into simple conceptual models. The ratio of the rate of sediment accumulation to the rate of structural uplift determines whether a fold develops a topographic expression above local base level. The balance between (a) stream power and rates of upstream deposition vs. (b) bedrock resistance and rates of crestal uplift and of fold widening determines whether an antecedent stream maintains its course or is defeated by a growing structure. Modern drainage configurations in actively folding landscapes can often be interpreted in terms of these competing variables, and through analysis of digital topography, detailed topographic characteristics of these folds can be quantified. Modern examples of growing folds display both defeated and persistent antecedent rivers, deflected drainages and laterally propagating structures. The topography associated with a defeated antecedent river at Wheeler Ridge, California, is consistent with a model in which defeat results from forced aggradation in the piggyback basin, without the need to vary discharge or uplift rate. Reconstruction of the long-term interplay between a depositional system and evolving folds requires a stratigraphic perspective, such as that provided by syntectonic strata which are directly juxtaposed with ancient folds and faults. Analysis of Palaeogene growth strata bounding the Catalan Coastal Ranges of NE Spain demonstrates the synchronous growth and the kinematic history of multiple folds and faults in the proximal foreland basin. Although dominated by transverse rivers which crossed fold crests, palaeovalleys, interfan lows, structural re-entrants and saddles, and rising anticlines diverted flow and influenced local deposition. In the ancient record, drainage-network events, such as avulsion or defeat of a transverse stream, usually cannot be unambiguously attributed to a single cause. Examination of ancient syntectonic strata from a geomorphological perspective, however, permits successive reconstructions of synorogenic topography, landscapes and depositional systems.     

The role of choke points in the ocean context

Choke points are narrow international waterways where three characteristics are present. The waterway should be narrow and capable of being closed off to both commercial and military snipping. There should be no readily available maritime route to utilize in the event of closure. Finally, the choke point should be of considerable significance to at least several States. This article identifies seven “primary” choke points which seem to satisfy these criteria. These are Gibraltar, Bab el Mandeb, Hormuz, the Danish and Turkish Straits, and the Suez and Panama Canals. It also identifies eleven “secondary” choke points, where at least one of the reouired characteristics is missing. The eleven include Dover, Bering and Magellan Straits, as well as Malacca-Singapore and a number of others in the Western Pacific. Turning to a consideration of the status of choke points in a Post-Cold War era, the article notes four basic assumptions: (1) the intense military rivalry between the United States and the Soviet Union will be reduced in coming years; (2) for many countries marine-related environmental concerns may be increasing considerably; (3) as the Cold War recedes, regional contests and confrontations will grow more intense; and (4) the overall densities of water-borne traffic will undergo change as the relative economic growth of regional centers changes. The article then considers the potential impact of these trends on the future role of choke points.