The surface properties of small asteroids: Peculiar Betulia—A case study

We present the results of extensive thermal-infrared observations of the C-type near-Earth Asteroid (1580) Betulia obtained in June 2002 with the 3-m NASA Infrared Telescope Facility (IRTF) on Mauna Kea, Hawaii. Betulia is a highly unusual object for which earlier radiometric observations, interpreted on the basis of simple thermal models, indicated a surface of high thermal inertia. A high thermal inertia implies a lack of thermally insulating regolith. Radiometric observations of other asteroids of comparable size indicate that regolith is present in nearly all cases. Knowledge of the surface thermal properties of small near-Earth asteroids is crucial for meaningful calculations of the Yarkovsky effect, which is invoked to explain the delivery of collisional fragments from the main belt into near-Earth orbits, and apparently has a significant influence on the orbital evolution of potentially hazardous near-Earth asteroids. Furthermore, apart from being an indicator of the presence of thermally insulating regolith on the surface of an asteroid, the thermal inertia determines the magnitude of the diurnal temperature variation and is therefore of great importance in the design of instrumentation for lander missions to small asteroids. In the case of Betulia our database is sufficiently broad to allow the use of more sophisticated thermal models than were available for earlier radiometric observations. The measured fluxes have been fitted with thermal-model emission continua to determine the asteroid's size and geometric albedo, pv. Fits obtained with a new thermophysical model imply an effective diameter of 4.57±0.46 km and an albedo of 0.077±0.015 and indicate a moderate surface thermal inertia of around 180 J m⁻² s^−0.5 K⁻¹. It is difficult to reconcile our results with earlier work, which indicate a larger diameter for Betulia and a high-thermal-inertia surface of bare rock.     

Archean high-Mg monzodiorite–syenite,epidote skarn,and biotite–sericite gold lodes in the Granny Smith–Wallaby district,Australia: U–Pb and Re–Os chronometry of two intrusion-related hydrothermal systems

The Granny Smith (37 t Au production) and Wallaby deposits (38 t out of a 180 t Au resource) are located northeast of Kalgoorlie, in 2.7 Ga greenstones of the Eastern Goldfields Province, the youngest orogenic belt of the Yilgarn craton, Western Australia. At Granny Smith, a zoned monzodiorite–granodiorite stock, dated by a concordant titanite–zircon U–Pb age of 2,665 ± 3 Ma, cuts across east-dipping thrust faults. The stock is fractured but not displaced and sets a minimum age for large-scale (1 km) thrust faulting (D2), regional folding (D1), and dynamothermal metamorphism in the mining district. The local gold–pyrite mineralization, controlled by fractured fault zones, is younger than 2,665 ± 3 Ma. In augite–hornblende monzodiorite, alteration progressed from a hematite-stained alkali feldspar–quartz–calcite assemblage and quartz–molybdenite–pyrite veins to a late reduced sericite–dolomite–albite assemblage. Gold-related monazite and xenotime define a U–Pb age of 2,660 ± 5 Ma, and molybdenite from veins a Re–Os isochron age of 2,661 ± 6 Ma, indicating that mineralization took place shortly after the emplacement of the main stock, perhaps coincident with the intrusion of late alkali granite dikes. At Wallaby, a NE-trending swarm of porphyry dikes comprising augite monzonite, monzodiorite, and minor kersantite intrudes folded and thrust-faulted molasse. The conglomerate and the dikes are overprinted by barren (<0.01 g/t Au) anhydrite-bearing epidote–actinolite–calcite skarn, forming a 600-m-wide and >1,600-m-long replacement pipe, which is intruded by a younger ring dike of syenite porphyry pervasively altered to muscovite + calcite + pyrite. Skarn and syenite are cut by pink biotite–calcite veins, containing magnetite + pyrite and subeconomic gold–silver mineralization (Au/Ag = 0.2). The veins are associated with red biotite–sericite–calcite–albite alteration in adjacent monzonite dikes. Structural relations and the concordant titanite U–Pb age of the skarn constrain intrusion-related mineralization to 2,662 ± 3 Ma. The main-stage gold–pyrite ore (Au/Ag >10) forms hematite-stained sericite–dolomite–albite lodes in stacked D2 reverse faults, which offset skarn, syenite, and the biotite–calcite veins by up to 25 m. The molybdenite Re–Os age (2,661 ± 10 Ma) of the ore suggests a genetic link to intrusive activity but is in apparent conflict with a monazite–xenotime U–Pb age (2,651 ± 6 Ma), which differs from that of the skarn at the 95% confidence level. The time relationships at both gold deposits are inconsistent with orogenic models invoking a principal role for metamorphic fluids released during the main phase of compression in the fold belt. Instead, mineralization is related in space and time to late-orogenic, magnetite-series, high-Mg monzodiorite–syenite intrusions of mantle origin, characterized by Mg/(Mg + Fe_TOTAL) = 0.31–0.57, high Cr (34–96 ppm), Ni (22–63 ppm), Ba (1,056–2,321 ppm), Sr (1,268–2,457 ppm), Th (15–36 ppm), and rare earth elements (total REE: 343–523 ppm). At Wallaby, shared Ca–K–CO₂ metasomatism and Th-REE enrichment (in allanite) link Au–Ag mineralization in biotite–calcite veins to the formation of the giant epidote skarn, implicating a Th + REE-rich syenite pluton at depth as the source of the oxidized hydrothermal fluid. At Granny Smith, lead isotope data and the Rb–Th–U signature of early hematite-bearing wall-rock alteration point to fluid released by the source pluton of the differentiated alkali granite dikes.     

Exurbia from the bottom-up: Confronting empirical challenges to characterizing a complex system

We describe empirical results from a multi-disciplinary project that support modeling complex processes of land-use and land-cover change in exurban parts of Southeastern Michigan. Based on two different conceptual models, one describing the evolution of urban form as a consequence of residential preferences and the other describing land-cover changes in an exurban township as a consequence of residential preferences, local policies, and a diversity of development types, we describe a variety of empirical data collected to support the mechanisms that we encoded in computational agent-based models. We used multiple methods, including social surveys, remote sensing, and statistical analysis of spatial data, to collect data that could be used to validate the structure of our models, calibrate their specific parameters, and evaluate their output. The data were used to investigate this system in the context of several themes from complexity science, including have (a) macro-level patterns; (b) autonomous decision making entities (i.e., agents); (c) heterogeneity among those entities; (d) social and spatial interactions that operate across multiple scales and (e) nonlinear feedback mechanisms. The results point to the importance of collecting data on agents and their interactions when producing agent-based models, the general validity of our conceptual models, and some changes that we needed to make to these models following data analysis. The calibrated models have been and are being used to evaluate landscape dynamics and the effects of various policy interventions on urban land-cover patterns.     

Estimating Deep Seafloor Interface and Volume Roughness Parameters Using the Multibeam-Hydrosweep System

Utilizing a hull-mounted, multinarrow beam echosounder onboard RV Polarstern, we measured variation of acoustic backscatter with incidence angles at two different sites in the Southern Oceans (Agulhas Plateau and the Riiser Larsen Sea). We modeled the data, using a composite roughness model, including water-sediment interface roughness and sediment volume roughness parameters. The model effectively uses the near normal incidence angle backscatter to determine the seafloor interface roughness parameters employing Helmholtz-Kirchhoff theory. Beyond 20° incidence angles, an application of Rayleigh-Rice theory is made by using a necessary splicing technique (combining both of the theories at 20° incidence angle). The estimated interface and volume roughness parameters are found to be in accordance with the known area geology.     

Bayesian estimation of seismic hazard for two sites in Switzerland

Seismic hazard analysis is based on data and models, which both are imprecise and uncertain. Especially the interpretation of historical information into earthquake parameters, e.g. earthquake size and location, yields ambiguous and imprecise data. Models based on probability distributions have been developed in order to quantify and represent these uncertainties. Nevertheless, the majority of the procedures applied in seismic hazard assessment do not take into account these uncertainties, nor do they show the variance of the results. Therefore, a procedure based on Bayesian statistics was developed to estimate return periods for different ground motion intensities (MSK scale).Bayesian techniques provide a mathematical model to estimate the distribution of random variables in presence of uncertainties. The developed method estimates the probability distribution of the number of occurrences in a Poisson process described by the parameter . The input data are the historical occurrences of intensities for a particular site, represented by a discrete probability distribution for each earthquake. The calculation of these historical occurrences requires a careful preparation of all input parameters, i.e. a modelling of their uncertainties. The obtained results show that the variance of the recurrence rate is smaller in regions with higher seismic activity than in less active regions. It can also be demonstrated that long return periods cannot be estimated with confidence, because the time period of observation is too short. This indicates that the long return periods obtained by seismic source methods only reflects the delineated seismic sources and the chosen earthquake size distribution law.     

Potential of High-Resolution Satellite Data in the Context of Vulnerability of Buildings

High-resolution space-borne remote sensing data are investigated for their potential to extract relevant parameters for a vulnerability analysis of buildings in European countries. For an evaluation of large earthquake scenarios, the number of parameters in models for vulnerability is reduced to a minimum of relevant information such as the type of building (age, material, number of storeys) and the geological and spatial context. Building-related parameters can be derived from remote sensing data either directly (e.g. height) or indirectly based on the recognition of the urban structure type in which the buildings are located. With the potential of a fully- or semi-automatic inventory of the buildings and their parameters, high-resolution satellite data and techniques for their processing are a useful supporting tool for the assessment of vulnerability.     

Spatial process and data models: Toward integration of agent-based models and GIS

The use of object-orientation for both spatial data and spatial process models facilitates their integration, which can allow exploration and explanation of spatial-temporal phenomena. In order to better understand how tight coupling might proceed and to evaluate the possible functional and efficiency gains from such a tight coupling, we identify four key relationships affecting how geographic data (fields and objects) and agent-based process models can interact: identity, causal, temporal and topological. We discuss approaches to implementing tight integration, focusing on a middleware approach that links existing GIS and ABM development platforms, and illustrate the need and approaches with example agent-based models.