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

Qingfeng Guan  Keith C. Clarke 《International journal of geographical information science》2013,27(5):695-722

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     

Jeff Onsted  Keith C. Clarke 《International journal of geographical information science》2013,27(5):881-898

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.  相似文献   

Hubble observations of Jupiter’s north–south conjugate ultraviolet aurora     

J.-C. Gérard  D. Grodent  A. Radioti  B. Bonfond  J.T. Clarke 《Icarus》2013

Comparisons of the northern and southern far ultraviolet (UV) auroral emissions of Jupiter from the Hubble Space Telescope (HST) or any other ultraviolet imager have mostly been made so far on a statistical basis or were not obtained with high sensitivity and resolution. Such observations are important to discriminate between different mechanisms responsible for the electron acceleration of the different components of the aurora such as the satellite footprints, the «main oval» or the polar emissions. The field of view of the ACS and STIS cameras on board HST is not wide enough to provide images of the full jovian disk. We thus compare the morphology of the north and south aurora observed 55 min apart and we point out similarities and differences. On one occasion HST pointed successively the two polar regions and auroral images were seen separated by only 3 min. This makes it possible to compare the emission structure and the emitted FUV power of corresponding regions. We find that most morphological features identified in one hemisphere have a conjugate counterpart in the other hemisphere. However, the power associated with conjugate regions of the main oval, diffuse or discrete equatoward emission observed quasi-simultaneously may be different in the two hemispheres. It is not directly nor inversely proportional to the strength of the B-field as one might expect for diffuse precipitation or field-aligned acceleration with equal ionospheric electron density in both hemispheres. Finally, the lack of symmetry of some polar emissions suggests that some of them could be located on open magnetic field lines.  相似文献   

A numerical modelling and experimental study of flow width dynamics on alluvial fans     

A. P. Nicholas  L. Clarke  T. A. Quine 《地球表面变化过程与地形》2009,34(15):1985-1993

Alluvial fans are dynamic landforms, the evolution of which is controlled by both external environmental forcing (climate, tectonics and base level change) and internal process‐form feedbacks. The latter include changes in flow configuration (between sheetflow and channelized flow states), driven by aggradation and degradation, which may in turn promote changes in sediment transport capacity. Recent numerical modelling indicates that such feedbacks may lead to dramatic and persistent fan entrenchment in the absence of external forcing. However, the parameterization of flow width within such models is untested to date and is subject to considerable uncertainty. This paper presents results from an experimental study of flow width dynamics on an aggrading fan in which spatial and temporal patterns of fan inundation are monitored continuously using analysis of digital vertical photography. Observed flow widths are compared with results from a simple theoretical model developed for non‐equilibrium (aggradational) conditions. Results demonstrate that the theoretical model is capable of capturing the first‐order characteristics of width adjustment over the course of the experiment, and indicate that flow width is a function of fan aggradation rate. This illustrates that models of alluvial flow width derived for equilibrium conditions may have limited utility in non‐equilibrium situations, despite their widespread use to date. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Estimating suspended sediment concentrations from turbidity measurements and the calibration problem   总被引：1，自引：0，他引：1  

Jean P. G. Minella  Gustavo H. Merten  José M. Reichert  Robin T. Clarke 《水文研究》2008,22(12):1819-1830

In situ turbidity meters are being increasingly used to generate continuous records of suspended sediment concentration in rivers. However, the usefulness of the information obtained depends heavily on the existence of a close relationship between fluctuations in suspended sediment concentration and turbidity and the calibration procedure that relates suspended sediment concentration to the turbidity meter's signal. This study assesses the relationship between suspended sediment concentration and turbidity for a small (1·19 km²) rural catchment in southern Brazil and evaluates two calibration methods by comparing the estimates of suspended sediment concentration obtained from the calibrated turbidity readings with direct measurements obtained using a USDH 48 suspended sediment sampler. With the first calibration method, the calibration relationship is derived by relating the turbidity readings to simultaneous measurements of concentration obtained from suspended sediment samples collected from the vicinity of the turbidity probe during flood events. With the second method, the calibration is based on the readings obtained from the turbidity meter when the probe immersed in samples of known concentration prepared using soils collected from the catchment. Overall, there was a close link between fluctuations in suspended sediment concentration and turbidity in the stream at the outlet of the catchment, and the estimates of sediment concentration obtained using the first calibration method corresponded closely with the conventionally measured sediment concentrations. However, use of the second calibration method introduced appreciable errors. When the estimated sediment concentrations were compared with the measured values, the mean errors were ± 122 mg l^?1 and + 601 mg l^?1 for the first and second calibration procedures respectively. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

The architecture, eruptive history, and evolution of the Table Rock Complex, Oregon: From a Surtseyan to an energetic maar eruption     

Brittany D. Brand  Amanda B. Clarke   《Journal of Volcanology and Geothermal Research》2009,180(2-4):203

The Table Rock Complex (TRC; Pliocene–Pleistocene), first documented and described by Heiken [Heiken, G.H., 1971. Tuff rings; examples from the Fort Rock-Christmas Lake valley basin, south-central Oregon. J. Geophy. Res. 76, 5615-5626.], is a large and well-exposed mafic phreatomagmatic complex in the Fort Rock–Christmas Lake Valley Basin, south-central Oregon. It spans an area of approximately 40 km², and consists of a large tuff cone in the south (TRC1), and a large tuff ring in the northeast (TRC2). At least seven additional, smaller explosion craters were formed along the flanks of the complex in the time between the two main eruptions. The first period of activity, TRC1, initiated with a Surtseyan-style eruption through a 60–70 m deep lake. The TRC1 deposits are dominated by multiple, 1-2 m thick, fining upward sequences of massive to diffusely-stratified lapilli tuff with intermittent zones of reverse grading, followed by a finely-laminated cap of fine-grained sediment. The massive deposits are interpreted as the result of eruption-fed, subaqueous turbidity current deposits; whereas, the finely laminated cap likely resulted from fallout of suspended fine-grained material through a water column. Other common features are erosive channel scour-and-fill deposits, massive tuff breccias, and abundant soft sediment deformation due to rapid sediment loading. Subaerial TRC1 deposits are exposed only proximal to the edifice, and consist of cross-stratified base-surge deposits. The eruption built a large tuff cone above the lake surface ending with an effusive stage, which produced a lava lake in the crater (365 m above the lake floor). A significant repose period occurred between the TRC1 and TRC2 eruptions, evidenced by up to 50 cm of diatomitic lake sediments at the contact between the two tuff sequences. The TRC2 eruption was the last and most energetic in the complex. General edifice morphology and a high percentage of accidental material suggest eruption through saturated TRC1 deposits and/or playa lake sediments. TRC2 deposits are dominated by three-dimensional dune features with wavelengths 200–500 m perpendicular to the flow, and 20–200 m parallel to the direction of flow depending on distance from source. Large U-shaped channels (10–32 m deep), run-up features over obstacles tens of meters high, and a large (13 m) chute-and-pool feature are also identified. The TRC2 deposits are interpreted as the products of multiple, erosive, highly-inflated pyroclastic surges resulting from collapse of an unusually high eruption column relative to previously documented mafic phreatomagmatic eruptions.  相似文献   

The lower mass function of young open clusters: clues to (sub)stellar formation     

E. Moraux  J. Bouvier  C. Clarke 《Astronomische Nachrichten》2005,326(10):985-990

We report the present day mass functions (PDMFs) of 4 young open clusters over a mass range from 30 Jupiter masses to 3M_ . Three of these clusters have been chosen to have a similar age of ∼100 Myr. Their PDMFs are remarkably similar and are comparable to the field mass function. This suggests little impact of initial conditions (stellar density, metallicity) on the mass distribution and raises some issues concerning the currently debated star and brown dwarf formation theories. The fourth cluster is older (600 Myr) which allow us to investigate the effect of the cluster dynamical evolution on the shape of the mass function. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

Guest Editorial: Selected Papers from the Fourth International Conference on Integrating Geographic Information Systems and Environmental Modeling (GIS/EM4)     

Keith Clarke  Brad Parks  & Michael Crane 《Transactions in GIS》2000,4(3):177-180

  相似文献   

A Study of DMS Oxidation in the Tropics: Comparison of Christmas Island Field Observations of DMS, SO2, and DMSO with Model Simulations   总被引：2，自引：0，他引：2  

G. Chen  D. D. Davis  P. Kasibhatla  A. R. Bandy  D. C. Thornton  B. J. Huebert  A. D. Clarke  B. W. Blomquist 《Journal of Atmospheric Chemistry》2000,37(2):137-160

This study reports comparisonsbetween model simulations, based on current sulfurmechanisms, with the DMS, SO₂ and DMSOobservational data reported by Bandy et al.(1996) in their 1994 Christmas Island field study. For both DMS and SO₂, the model results werefound to be in excellent agreement with theobservations when the observations were filtered so asto establish a common meteorological environment. Thisfiltered DMS and SO₂ data encompassedapproximately half of the total sampled days. Basedon these composite profiles, it was shown thatoxidation of DMS via OH was the dominant pathway withno more than 5 to 15% proceeding through Cl atoms andless than 3% through NO₃. This analysis wasbased on an estimated DMS sea-to-air flux of 3.4 ×10⁹ molecs cm^-2 s^-1. The dominant sourceof BL SO₂ was oxidation of DMS, the overallconversion efficiency being evaluated at 0.65 ± 0.15. The major loss of SO₂ was deposition to theocean's surface and scavenging by aerosol. Theresulting combined first order k value was estimated at 1.6 × 10^-5 s^-1. In contrast to the DMSand SO₂ simulations, the model under-predictedthe observed DMSO levels by nearly a factor of 50. Although DMSO instrument measurement problems can notbe totally ruled out, the possibility of DMSO sourcesother than gas phase oxidation of DMS must beseriously considered and should be explored in futurestudies.  相似文献   

Watson: A new link in the HE iron chain     

Edward Olsen  Andrew Davis  Roy S. Clarke  Ludolf Schultz  Hartwig W. Weber  Robert Clayton  Toshiko Mayeda  Eugene Jarosewich  Paul Sylvester  Lawrence Grossman  Ming-Sheng Wang  Michael E. Lipschutz  Ian M. Steele  James Schwade 《Meteoritics & planetary science》1994,29(2):200-213

Abstract— Watson, which was found in 1972 in South Australia, contains the largest single silicate rock mass seen in any known iron meteorite. A comprehensive study has been completed on this unusual meteorite: petrography, metallography, analyses of the silicate inclusion (whole rock chemical analysis, INAA, RNAA, noble gases, and oxygen isotope analysis) and mineral compositions (by electron microprobe and ion microprobe). The whole rock has a composition of an H-chondrite minus the normal H-group metal and troilite content. The oxygen isotope composition is that of the silicates in the HE iron meteorites and lies along an oxygen isotope fractionation line with the H-group chondrites. Trace elements in the metal confirm Watson is a new HE iron. Whole rock Watson silicate shows an enrichment in K and P (each ～2X H-chondrites). The silicate inclusion has a highly equilibrated igneous (peridotite-like) texture with olivine largely poikilitic within low-Ca pyroxene: olivine (Fa₂₀), opx (Fs₁₇Wo₃), capx (Fs₉Wo₄₁) (with very fine exsolution lamellae), antiperthite feldspar (An_1–Or₅) with <1 μm exsolution lamellae (An_1–3Or_>40), shocked feldspar with altered stoichiometry, minor whitlockite (also a poorly characterized interstitial phosphate-rich phase) and chromite, and only traces of metal and troilite. The individual silicate minerals have normal chondritic REE patterns, but whitlockite has a remarkable REE pattern. It is very enriched in light REE (La is 720X C1, and Lu is 90X C1, as opposed to usual chonditic values of ～300X and 100–150X, respectively) with a negative Eu anomaly. The enrichment of whole rock K is expressed both in an unusually high mean modal Or content of the feldspar, Or₁₃, and in the presence of antiperthite. Whole rock trace element data for the silicate mass support the petrography. Watson silicate was an H-chondrite engulfed by metal and melted at > 1550 °C. A flat refractory lithophile and flat REE pattern (at ～1x average H-chondrites) indicate that melting took place in a relatively closed system. Immiscible metal and sulfide were occluded into the surrounding metal host. Below 1100 °C, the average cooling rate is estimated to have been ～1000 °C/Ma; Widmanstätten structure formed, any igneous zoning in the silicates was equilibrated, and feldspar and pyroxene exsolution took place. Cooling to below 300 °C was completed by 3.5 Ga B. P. At 8 Ma, a shock event took place causing some severe metal deformation and forming local melt pockets of schreibersite/metal. This event likely caused the release of Watson into interplanetary space. The time of this event, 8Ma, corresponds to the peak frequency of exposure ages of the H-chondrites. This further confirms the link between HE irons and the H-chondrites, a relationship already indicated by their common oxygen isotope source. Watson metal structures are very similar to those in Kodaikanal. Watson, Kodaikanal and Netschaëvo form the young group of HE meteorites (ages 3.7 ± 0.2 Ga). They appear to represent steps in a chain of events that must have taken place repeatedly on the HE parent body or bodies from which they came: chondrite engulfed in metal (Netschaëvo); chondrite melted within metal (Watson); and finally melted silicate undergoing strong fractionation with the fractionated material emplaced as globules within metal (Kodaikanal). Watson fills an important gap in understanding the sequence of events that took place in the evolution of the IIE-H parent body(ies). This association of H-chondrite with HE metal suggests a surface, or near surface process-a suggestion made by several other researchers.  相似文献