Stochastic diffusion in inverse square fields: General formulation for interplanetary gas

The Fokker-Planck equation for small stochastic changes to particles in Kepler orbits has to be formulated in terms of the integrals of motion. We generalize the modelling of proton and electron collisional perturbations to gas particles on trajectories through the solar system in order to include both spatial and velocity diffusion. The general solution is obtained in terms of a 4-dimensional normal distribution. Treatment of the singularity in the Fokker-Planck operator reduces the dimensionality by one. In addition to extending earlier results for anisotropic collisional heating in the thermal approximation, the present formulation gives the changes in density due to the mean repulsive force and to perturbations of trajectories (spatial diffusion). The net diffusion is almost everywhere towards the sun and the density increase is significant in the downstream hydrogen wake, particularly where destructive depletion is strong and gravitational focussing weak.     

Methods for rapidly processing angular masks of next-generation galaxy surveys

As galaxy surveys become larger and more complex, keeping track of the completeness, magnitude limit and other survey parameters as a function of direction on the sky becomes an increasingly challenging computational task. For example, typical angular masks of the Sloan Digital Sky Survey contain about   N = 300 000  distinct spherical polygons. Managing masks with such large numbers of polygons becomes intractably slow, particularly for tasks that run in time     with a naive algorithm, such as finding which polygons overlap each other. Here we present a 'divide-and-conquer' solution to this challenge: we first split the angular mask into pre-defined regions called 'pixels', such that each polygon is in only one pixel, and then perform further computations, such as checking for overlap, on the polygons within each pixel separately. This reduces     tasks to     , and also reduces the important task of determining in which polygon(s) a point on the sky lies from     to     , resulting in significant computational speedup. Additionally, we present a method to efficiently convert any angular mask to and from the popular healpix format. This method can be generically applied to convert to and from any desired spherical pixelization. We have implemented these techniques in a new version of the mangle software package, which is freely available at http://space.mit.edu/home/tegmark/mangle/ , along with complete documentation and example applications. These new methods should prove quite useful to the astronomical community, and since mangle is a generic tool for managing angular masks on a sphere, it has the potential to benefit terrestrial mapmaking applications as well.     

The geology of Howth and Ireland's Eye,Co. Dublin

The Cambrian sediments in Howth and Ireland's Eye to the northeast of Dublin have been mapped and the submarine deltaic Censure Group and the overlying Nose of Howth Group, which consists of interbedded turbidites and olistostromes, have been distinguished. Large scale sedimentary disruption has accompanied detachment of quartzitic sandstone beds in the Censure Group, where the bases of detached blocks have slump features while their tops generally do not. This contrasts with the chaotic sedimentary breccias in the olistostromes of the Nose of Howth Group. Blocks in the Censure and Nose of Howth Groups are up to 450 m and 700 m across respectively. The turbidies have come from the southeast and east whereas the slumping appears to have been generated from the southwest. Slumping during sedimentation was followed closely by large-scale folding and sliding of sedimentary or tectono-sedimentary origin, and the bulk appearance of the sediments is due to these sedimentary and tectono-sedimentary mechanisms. The main structure is an easterly facing, steep plunging syncline which locally is associated with a poorly preserved cleavage. This structure predated the main penetrative cleavage which is not axial planar to it.     

On Rates and Acceleration Trends of Global Glacier Mass Changes

Worldwide glacier mass changes are considered to represent natural key variables within global climate-related monitoring programmes, especially with respect to strategies concerning early detection of enhanced greenhouse effects on climate. This is due to the fact that glacier mass changes provide important quantitative information on rates of change, acceleration tendencies and pre-industrial variability relating to energy exchange at the earth/athmosphere interface. During the coming decades, excess radiation income and sensible heat (a few watts per square metre) as calculated with numerical climate models are both estimated to increase by a factor of about two to four as compared to the mean of the 20th century. The rate of average annual mass loss (a few decimetres per year) measured today on mountain glaciers in various parts of the world now appears to accelerate accordingly, even though detailed interpretation of the complex processes involved remains difficult. Within the framework of secular glacier retreat and Holocene glacier fluctuations, similar rates of change and acceleration must have taken place before, i.e. during times of weak anthropogenic forcing. However, the anthropogenic influences on the atmosphere could now and for the first time represent a major contributing factor to the observed glacier shrinkage at a global scale. Problems with such assessments mainly concern aspects of statistical averaging, regional climate variability, strong differences in glacier sensitivity and relations between mass balance and cumulative glacier length change over decadal to secular time scales. Considerable progress has recently been achieved in these fields of research.     

Australia's Box–Ironbark Forests and Woodlands: saving the fragments of a threatened ecosystem

Australia's box–ironbark forests and woodlands once covered about 14 per cent of the State of Victoria on the riverine plains and foothills of the Great Dividing Range. But approximately 83 per cent of the total original habitat has been destroyed and what remains of this significant ecosystem is now highly fragmented and vulnerable to further degradation. Moreover, only 14 per cent of the area remaining is on public land. A 10 year campaign on the part of the environmental movement eventually succeeded in forcing the State government to conduct an independent inquiry into this ecosystem and make recommendations on future management. This paper outlines the innovative public participation process adopted by the Victorian State government and the outcomes of the inquiry. A subsequent compensation package for commercial operations disadvantaged by the proclamation of a series of new national parks is also discussed, as are the shortcomings of a process that can have little or no impact on what happens on private land.     

A Temporal Variant‐Invariant Validation Approach for Agent‐based Models of Landscape Dynamics

Agent‐based modeling provides a means for addressing the way human and natural systems interact to change landscapes over time. Until recently, evaluation of simulation models has focused on map comparison techniques that evaluate the degree to which predictions match real‐world observations. However, methods that change the focus of evaluation from patterns to processes have begun to surface; that is, rather than asking if a model simulates a correct pattern, models are evaluated on their ability to simulate a process of interest. We build on an existing agent‐based modeling validation method in order to present a temporal variant‐invariant analysis (TVIA). The enhanced method, which focuses on analyzing the uncertainty in simulation results, examines the degree to which outcomes from multiple model runs match some reference to how land use parcels make the transition from one land use class to another over time. We apply TVIA to results from an agent‐based model that simulates the relationships between landowner decisions and wildfire risk in the wildland‐urban interface of the southern Willamette Valley, Oregon, USA. The TVIA approach demonstrates a novel ability to examine uncertainty across time to provide an understanding of how the model emulates the system of interest.     

Ra-partitioning between phlogopite and silicate melt and 226Ra/Ba–230Th/Ba isochrons

In this study we experimentally determine phlogopite/melt partition coefficients of Ra and other trace elements in a lamproitic system. This work was achieved using an analytical technique (LA-ICP-MS) with low detection limits (~ 0.01 fg) permitting the measurement of the very low Ra concentrations feasible in experiments (~ 1 ppb). D_Ra^{phlogopite/melt} was determined to 2.28 ± 0.44 and 2.84 ± 0.47 in two experiments, the ratio D_Ra/D_Ba is around 1.6. The compatibility of Ra in phlogopite results from an ionic radius being close to the apex of the lattice strain parabola for earth alkalis in the large XII-coordinated interlayer site of phlogopite. A re-evaluation of D_Ra and D_Ra/D_Ba for magmatic minerals containing appreciable Ra, yields D_Ra^mineral/melt ranging from ~ 2.6 for phlogopite down to 2–3 ? 10^? 5 for pyroxenes, and D_Ra/D_Ba^mineral/melt from ~ 4 for leucite to 2 ? 10^? 2 for orthopyroxene. The influence of melt composition on D_Ra/D_Ba is less than 10%. All investigated minerals have different D_Ra/D_Ba, strongly fractionating Ra from Ba. Thus, for magmatic systems, (²²⁶Ra)/Ba in the various minerals is not constant, these minerals do not form a straight line in the (²²⁶Ra)/Ba–(²³⁰Th)/Ba system at the time of crystallization and thus, there is no (²²⁶Ra)/Ba–(²³⁰Th)/Ba isochron at t₀. ²²⁶Ra–²³⁰Th–Ba mineral dating is thus applicable only to model ages calculated from mineral–glass pairs with known D_Ra.     

Magnetic crustal character in Central Ireland

The simplest magnetic model for the upper 20 km of crust in central Ireland is of two layers. A thin, weakly magnetic upper layer carrying a number of magnetic bodies that can be related to surface geology overlies a magnetic crustal zone whose gross magnetization increases towards the north. The deeper trace of the lapetus Suture in central Ireland forms the boundary between a more magnetic northern crustal area and a less magnetic southern area. A pre-Carboniferous major dextral offset of this suture boundary occurs in the Limerick area. The suture is traced beneath the Dingle peninsula in the far south-west of Ireland. In east-central Ireland, between the Iapetus Suture and Southern Upland Line, the crust contains a large, south-dipping magnetic body with no surface expression, which produces the Virginia magnetic anomaly. The development of magnetic bodies — in particular, volcanic rocks of Carboniferous age — has been controlled by regional stress patterns. The major volcanic centre at Limerick lies at the intersection of two major Caledonian block fault systems just to the north of the surface trace of the Hercynian Front.