A general-purpose parallel raster processing programming library test application using a geographic cellular automata model

A general-purpose parallel raster processing programming library (pRPL) was developed and applied to speed up a commonly used cellular automaton model with known tractability limitations. The library is suitable for use by geographic information scientists with basic programming skills, but who lack knowledge and experience of parallel computing and programming. pRPL is a general-purpose programming library that provides generic support for raster processing, including local-scope, neighborhood-scope, regional-scope, and global-scope algorithms as long as they are parallelizable. The library also supports multilayer algorithms. Besides the standard data domain decomposition methods, pRPL provides a spatially adaptive quad-tree-based decomposition to produce more evenly distributed workloads among processors. Data parallelism and task parallelism are supported, with both static and dynamic load-balancing. By grouping processors, pRPL also supports data–task hybrid parallelism, i.e., data parallelism within a processor group and task parallelism among processor groups. pSLEUTH, a parallel version of a well-known cellular automata model for simulating urban land-use change (SLEUTH), was developed to demonstrate full utilization of the advanced features of pRPL. Experiments with real-world data sets were conducted and the performance of pSLEUTH measured. We conclude not only that pRPL greatly reduces the development complexity of implementing a parallel raster-processing algorithm, it also greatly reduces the computing time of computationally intensive raster-processing algorithms, as demonstrated with pSLEUTH.     

The inclusion of differentially assessed lands in urban growth model calibration: a comparison of two approaches using SLEUTH

This research sought to understand the role that differentially assessed lands (lands in the United States given tax breaks in return for their guarantee to remain in agriculture) play in influencing urban growth. Our method was to calibrate the SLEUTH urban growth model under two different conditions. The first used an excluded layer that ignored such lands, effectively rendering them available for development. The second treated those lands as totally excluded from development. Our hypothesis was that excluding those lands would yield better metrics of fit with past data. Our results validate our hypothesis since two different metrics that evaluate goodness of fit both yielded higher values when differentially assessed lands are treated as excluded. This suggests that, at least in our study area, differential assessment, which protects farm and ranch lands for tenuous periods of time, has indeed allowed farmland to resist urban development. Including differentially assessed lands also yielded very different calibrated coefficients of growth as the model tried to account for the same growth patterns over two very different excluded areas. Excluded layer design can greatly affect model behavior. Since differentially assessed lands are quite common through the United States and are often ignored in urban growth modeling, the findings of this research can assist other urban growth modelers in designing excluded layers that result in more accurate model calibration and thus forecasting.     

A bundle adjustment approach with inner constraints for the scaled orthographic projection

Bundle adjustment is a method for simultaneously calculating both the interior and exterior orientation parameters of a set of images, and the object-space coordinates of the observed points. In the case of long focal length lenses and narrow field-of-view (FOV) imaging situations, collinearity based (perspective projection) algorithms may result in linear dependencies between parameters that cause solution instability. The use of a scaled orthographic projection model based on linear algebraic formulations was therefore adopted to reduce this risk. Using quaternions, a new mathematical model is derived that includes the partial derivatives as well as the inner constraint equations for a scaled orthographic bundle adjustment. The model was then tested using two image sets of a single, small vessel (about 6 m length) with a cube target of known dimensions at two distinct ranges; perspective solutions were also calculated for comparison. RMS residual errors of 0.74-0.78 pixels associated with the new method compare favorably to a residual error range of 0.59-0.74 pixels using a perspective bundle adjustment of the same target points. Relative precisions (as a ratio of target size) of between 1:1650 and 1:750 have been achieved at ranges of 375 m and 850 m, respectively, given comparisons with the known cube dimensions. A third image dataset consisting of a network of 16 images was solved with a 1:2200 relative precision showing the new method can successfully handle high redundancy. For the experiments that were conducted, the new method was found to produce less precise results than the perspective bundle solution for a FOV of 0.50-0.65° where the object fills 5-8% of the image. However, it was found to match the precision of the perspective model (with an uncalibrated camera) for a FOV of 0.20-0.30° where the object of interest fills only 1-2% of the full image.     

Hurricane effects on subtropical pine rocklands of the Florida Keys

We investigate the effects of Hurricane Wilma??s storm surge (23?C24 October 2005) on the dominant tree Pinus elliottii var densa (South Florida slash pine) and rare plant species in subtropical pine rocklands of the Lower Florida Keys. We examine the role of elevation on species abundance in 1995 (Hurricane Betsy in 1965), 2005 (Hurricane Georges in 1998), and 2008 (Hurricane Wilma in 2005) to investigate if hurricanes influence abundance by eliminating plants at lower elevation on Big Pine Key, the largest island in the Lower Florida Keys. We compare densities before and after Hurricane Wilma over the 2005?C2008 sampling period and examine the role of elevation on changes in pine and rare species densities three years after Hurricane Wilma. We use elevation to assess the impact of hurricanes because elevation determined whether a location was influenced by storm surge (maximum surge of 2 m) in the Lower Florida Keys, where pine rocklands occur at a maximum elevation of 3 m. In 1995 (30 years after a major storm), elevation did not explain the abundance of South Florida slash pine or Chamaecrista lineata, but explained significant variation in abundance of Chamaesyce deltoidea. The latter two species are rare herbaceous plants restricted to pine rocklands. In 2008, 3 years after Hurricane Wilma, the positive relationship between elevation and abundance was strongest for South Florida slash pine, C. deltoidea, and C. lineata. Effects of Hurricane Wilma were not significant for rare species with wider distribution, occurring in plant communities adjacent to pine rocklands and in disturbed rocklands. Our results suggest that hurricanes drive population dynamics of South Florida slash pine and rare species that occur exclusively in pine rocklands at higher elevations. Rare species restricted to pine rocklands showed dramatic declines after Hurricane Wilma and were eliminated at elevations <0.5 m. Widely distributed rare species did not show significant changes in density after Hurricane Wilma. Abundance increased with elevation for South Florida slash pine and C. lineata after the hurricane. In an environment influenced by sea level rise, concrete plans to conserve pine ecosystems are warranted. Results from this study will help define conservation strategies by strengthening predictive understanding of plant responses to disturbance in the backdrop of sea level rise.     

Cytochrome P4501A biomarker indication of the timeline of chronic exposure of Barrow's goldeneyes to residual Exxon Valdez oil

We examined hepatic EROD activity, as an indicator of CYP1A induction, in Barrow’s goldeneyes captured in areas oiled during the 1989 Exxon Valdez spill and those from nearby unoiled areas. We found that average EROD activity differed between areas during 2005, although the magnitude of the difference was reduced relative to a previous study from 1996/1997, and we found that areas did not differ by 2009. Similarly, we found that the proportion of individuals captured from oiled areas with elevated EROD activity (?2 times unoiled average) declined from 41% in winter 1996/1997 to 10% in 2005 and 15% in 2009. This work adds to a body of literature describing the timelines over which vertebrates were exposed to residual Exxon Valdez oil and indicates that, for Barrow’s goldeneyes in Prince William Sound, exposure persisted for many years with evidence of substantially reduced exposure by 2 decades after the spill.     

On the use of IPCC-class models to assess the impact of climate on Living Marine Resources

The study of climate impacts on Living Marine Resources (LMRs) has increased rapidly in recent years with the availability of climate model simulations contributed to the assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Collaboration between climate and LMR scientists and shared understanding of critical challenges for such applications are essential for developing robust projections of climate impacts on LMRs. This paper assesses present approaches for generating projections of climate impacts on LMRs using IPCC-class climate models, recommends practices that should be followed for these applications, and identifies priority developments that could improve current projections. Understanding of the climate system and its representation within climate models has progressed to a point where many climate model outputs can now be used effectively to make LMR projections. However, uncertainty in climate model projections (particularly biases and inter-model spread at regional to local scales), coarse climate model resolution, and the uncertainty and potential complexity of the mechanisms underlying the response of LMRs to climate limit the robustness and precision of LMR projections. A variety of techniques including the analysis of multi-model ensembles, bias corrections, and statistical and dynamical downscaling can ameliorate some limitations, though the assumptions underlying these approaches and the sensitivity of results to their application must be assessed for each application. Developments in LMR science that could improve current projections of climate impacts on LMRs include improved understanding of the multi-scale mechanisms that link climate and LMRs and better representations of these mechanisms within more holistic LMR models. These developments require a strong baseline of field and laboratory observations including long time series and measurements over the broad range of spatial and temporal scales over which LMRs and climate interact. Priority developments for IPCC-class climate models include improved model accuracy (particularly at regional and local scales), inter-annual to decadal-scale predictions, and the continued development of earth system models capable of simulating the evolution of both the physical climate system and biosphere. Efforts to address these issues should occur in parallel and be informed by the continued application of existing climate and LMR models.     

Tectonic implications of some field relations of the Adelaidean Cooee Dolerite,Tasmania

Complex contact relations between thin sills of the Adelaidean Cooee Dolerite and the enclosing flysch sediments (Burnie Formation) indicate emplacement of the sills into substantially uncompacted quasi‐liquid sediments very shortly after deposition. The sills are consequently coeval with the sediments rather than synchronous with the first phase of folding of the sediments (the Penguin Orogeny).

The K/Ar date on the Cooee Dolerite, here revised as 725 ± 35 m.y. (J. Richards, pers. comm.), very probably provides an estimate of the age of the sediments, rather than the age of the orogeny. Age relationships of the Cooee Dolerite, its host sediments and their deformation proposed hitherto can no longer be regarded as settled. The Burnie Formation may be part of a passive continental margin to the Rocky Cape Block, which was deformed during the Cambrian as a result of subduction beneath the active margin of the Tyennan Block.     

Fore‐arc evolution in the Tasman Geosyncline: The origin of the southeast Australian continental crust

A new general model describing the extended evolution of fore‐arc terrains is used to analyse the evolution of the southern Tasman Geosyncline and the concomitant growth and kratonisation of the continental crust of southeast Australia during the Palaeozoic. The southern Tasman Geosyncline comprises ten arc terrains (here defined), most of which are east‐facing, and several features formed by crustal extension. Each arc terrain consists of several strato‐tectonic units: a volcanic arc, subduction complex and fore‐arc sequence formed during subduction; and an overlying post‐arc sequence which post‐dates subduction and is composed of flysch, neritic sediments or subaerial volcanics.

When these materials attained a thickness of c. 20 km their internal heat‐balance caused partial melting of the subduction complex and the hydrated oceanic lithosphere trapped beneath it, to yield S‐ and I‐type granitic magma. The magma rose, inducing pervasive deformation of each arc terrain and emplacement of granitoid plutons at high levels in the evolving crust. Transitional basins then developed in many terrains on top of their volcanic arcs or the thinner parts of the buried accretionary prisms. After deformation of the transitional sequences, platform cover accumulated, marking the completion of kratonisation.

Analysis of each arc terrain in terms of the above units leads to a predicted ‘stratigraphy’ for the continental crust of southeast Australia. The crust is complexly layered, with lateral discontinuities reflecting the boundaries of arc terrains which were successively accreted, principally back‐arc to fore‐arc, during crustal development.     

Permeameter data verify new turbulence process for MODFLOW

A sample of Key Largo Limestone from southern Florida exhibited turbulent flow behavior along three orthogonal axes as reported in recently published permeameter experiments. The limestone sample was a cube measuring 0.2 m on edge. The published nonlinear relation between hydraulic gradient and discharge was simulated using the turbulent flow approximation applied in the Conduit Flow Process (CFP) for MODFLOW-2005 mode 2, CFPM2. The good agreement between the experimental data and the simulated results verifies the utility of the approach used to simulate the effects of turbulent flow on head distributions and flux in the CFPM2 module of MODFLOW-2005.