Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.     

Shape and distribution analysis of Merensky Reef potholing, Northam Platinum Mine, western Bushveld Complex: implications for pothole formation and growth

Syn-magmatic removal of the cumulate pile during the formation of the Bushveld Complex resulted in “potholes”. Erosion progressed downward in the cumulate pile, resulting in a series of steep, transgressive contacts between locally conformable potholed reefs in the regional pothole sub-facies of the Swartklip Facies in the western limb of the Bushveld Complex. The deepest of these potholes, “third-order” or “FWP2” potholing, occurs where the base of the Merensky Cyclic Unit transgresses the Upper Pseudo-Reef Chromitite marker horizon. The base of a FWP2 pothole on Northam Platinum Mine consists of an unconformable stringer Merensky Chromitite overlain by a medium-grained, poikilitic orthopyroxenite and underlain by either a pegmatitic harzburgite or the medium-grained Lower Pseudo-Reef Anorthosite. Detailed shape and distribution analysis of FWP2 potholes reveals underlying patterns in their shape and distribution which, in turn, suggest a structural control. The ratio between pothole short vs long axes is 0.624 (N=1,385), although the ratio increases from 0.48 to 0.61 in the long axis range 10 to 60 m, then decreases from 0.61 to 0.57 from 61 to 100 m, increasing again from 0.57 to 0.61 from 101 to 400 m, suggesting that there is not a simple relationship between pothole shape and size. Shape (circularity, eccentricity, and dendricity) analysis of a subset of 638 potholes indicates that potholes with long axes <100 m have an elliptical, average normalized shape, elongate on a 120–150° orientation. Potholes with long axis lengths >100 m have an average normalized shape that is bilobate and elongate on a 120° orientation. The average aspect ratio (short axis length divided by long axis length) of potholes is highest for potholes with long axis lengths >100 m and lowest for potholes with long axis lengths between 35 and 60 m. The most common long axis orientation for potholes with long axis lengths <100 m is 150° but 120° for long axis lengths >100 m. Fractal analysis indicates that the distribution of pothole centers is controlled neither by a single nor several interacting fractal dimensions. Autocorrelation (Fry) analysis of the distribution of pothole centers shows recurring pothole distribution trends at 038, 070, and 110° for potholes over the full range of long axis lengths, while the trends of 008 and 152° occur in potholes with long axes lengths between 60 and 100 m. Chi-squared (X ²) analysis of the locations of pothole centers suggests that the distribution of small potholes is highly non-uniform but becomes exponentially more uniform with increasing pothole size. The model which best fits the observed shape and distribution analysis is a combination of protracted independent growth and “nearest neighbor” merging along specific orientations. For instance, the clustered distribution of original pothole centers resulted in merged potholes with long axes lengths of up to 60 m, exhibiting short vs long axes ratios of 0.61, preferred orientations of 150°, and alignment along 010 and 150° trends. Further independent growth allowed for merging of similar-sized (and smaller) neighboring potholes, generating potholes with long axes of up to 100 m in length, a preferred long axis orientation of 150°, and alignment along 010, 040, 075, and 150°. Subsequent preferential merging occurred along a 120° trend, thereby preserving a bilobate form. This implies that while pothole initiation and enlargement may be driven by a “top-down” (i.e., possibly thermomechanical) process, an underlying linear or structural catalyst/control is revealed in changes in pothole shape during enlargement and, furthermore, in the preferred trends along which potholes merged over a considerable period, possibly concomitant with adjustment of major structures in the footwall to the Bushveld Complex and pulses into the magma chamber.     

Cross-Calibration and Long-Term Monitoring of the Microwave Radiometers of ERS, TOPEX, GFO, Jason, and Envisat

The radiometers on board the satellites ERS-1, TOPEX/Poseidon, ERS-2, GFO, Jason-1, and Envisat measure brightness temperatures at two or three different frequencies to determine the total columnal water vapor content and wet tropospheric path delay, a major correction to the altimeter range measurements. In order to asses the long-term stability of the path delay, the radiometers are calibrated against vicarious cold and hot references, against each other, and against several atmospheric models. Four of these radiometers exhibit significant drifts in at least one of the channels, resulting in yet unmodeled errors in path delay of up to 1 mm/year, thus limiting the accuracy at which global sea level rise can be inferred from the altimeter range measurements.     

Chemical variations and genetic relationships between the Hess and Foy offset dikes at the Sudbury impact structure

Offset dikes are found concentrically around—and extending radially outward from—the Sudbury Igneous Complex (SIC), which represents an ~3 km thick differentiated impact melt sheet. The dikes are typically composed of an inclusion‐rich, so‐called quartz diorite (IQD) in the center of the dike, and an inclusion‐poor quartz diorite (QD) along the margins of the dike. New exposures of the intersection between the concentric Hess and radial Foy offset dikes provide an excellent opportunity to understand the relationship between the radial and concentric offset dikes and their internal phases. The goal was to constrain the timing of the dike emplacements relative to the impact and formation of the SIC. Results herein suggest that (1) the timing between the emplacement of the QD and IQD melts was geologically short, (2) the Hess and Foy dikes coexisted as melts at the same time and the intersection between them represents a mixture of the two, (3) the Foy dike has a slightly more evolved chemical composition than the Hess dike, and (4) the IQD melt from the Foy dike underwent some degree of chemical fractionation after its initial emplacement.     

Impulsively generated MHD waves and their detectability in solar coronal loops

Impulsively generated magnetohydrodynamic waves in solar coronal loops, with arbitrary plasma , are studied numerically by a flux-corrected transport algorithm. Numerical results show that the total reflection which occurs in the region of low Alfvén speed leads to trapped fast kink magnetosonic waves. These waves propagate along the slab and exhibit periodic, quasi-periodic, and decay phases. As a consequence of the difference in wave propagation speeds, the time signatures of the slow magnetosonic waves are delayed in time in comparison to the time signatures of the fast magnetosonic and Alfvén waves. An interaction between the waves can generate a longer lasting and complex quasi-periodic phase of the fast wave. We discuss also the observational detectability of such MHD waves in optical, radio, and soft X-ray wavelenghts.     

Transport of ionizing radiation in terrestrial-like exoplanet atmospheres

The propagation of ionizing radiation through model atmospheres of terrestrial-like exoplanets is studied for a large range of column densities and incident photon energies using a Monte Carlo code we have developed to treat Compton scattering and photoabsorption. Incident spectra from parent star flares, supernovae, and gamma-ray bursts are modeled and compared to energetic particles in importance. Large irradiation events with fluences of 10⁶-10⁹ erg cm⁻² at the conventional habitable zone can occur at a rate from many per day (flares from young low-mass parent stars) to ∼100 per Gyr (supernovae and gamma-ray bursts). We find that terrestrial-like exoplanets with atmospheres thinner than about 100 g cm⁻² block nearly all X-rays, but transmit and reprocess a significant fraction of incident γ-rays, producing a characteristic, flat surficial spectrum. Thick atmospheres (?100 g cm⁻²) efficiently block even γ-rays, but nearly all the incident energy is redistributed into diffuse UV and visible aurora-like emission, increasing the effective atmospheric transmission by many orders of magnitude. Depending on the presence of molecular UV absorbers and atmospheric thickness, up to 10% of the incident energy can reach the surface as UV reemission. For the Earth, between 2×10⁻³ and 4×10⁻² of the incident flux reaches the ground in the biologically effective 200-320 nm range, depending on O₂/O₃ shielding. For atmospheres thicker than ∼50 g cm⁻² in the case of pure Rayleigh scattering and ∼100 g cm⁻² in the case of O₂/O₃ absorption, the UV reemission exceeds the surficial transmitted ionizing radiation. We also discuss the effects of angle of incidence and derive a modified two-stream approximation solution for the UV transfer. Finally, we suggest that transient atmospheric ionization layers can be frequently created at altitudes lower than the equilibrium layers that result from steady irradiation and winds from the parent star. We suggest that these events can produce frequent fluctuations in atmospheric ionization levels and surficial UV fluxes on terrestrial-like planets.     

Detectors For Space Science

The design of detector systems for flight applications requires the consideration of a number of issues unique to space instrumentation. Flight detectors must endure hostile radiation environments and thermal extremes. Paramount importance is given to reliability since inflight replacement is at best difficult and usually impossible. Flight detectors are also significant cost and design drivers since they often determine key requirements for flight instruments such as volume, mass, power consumption, heat dissipation and communications budgets. In this paper we describe the primary concerns in developing flight detector systems, and review the challenges posed by future NASA and ESA space science missions for detector development.     

Hydrated sulphuric acid in dense molecular clouds

We consider sulphur depletion in dense molecular clouds, and suggest hydrated sulphuric acid, H₂SO₄ ·  n H₂O, as a component of interstellar dust in icy mantles. We discuss the formation of hydrated sulphuric acid in collapsing clouds and its instability in heated regions in terms of the existing hot core models and observations. We also show that some features of the infrared spectrum of hydrated sulphuric acid have correspondence in the observed spectra of young stellar objects.     

NIGHTGLOW: an instrument to measure the Earth’s nighttime ultraviolet glow—results from the first engineering flight

We have designed and built an instrument to measure and monitor the “nightglow” of the Earth’s atmosphere in the near ultraviolet (NUV). In this paper we describe the design of this instrument, called NIGHTGLOW. NIGHTGLOW is designed to be flown from a high altitude research balloon, and circumnavigate the globe. NIGHTGLOW is a NASA, University of Utah, and New Mexico State University project. A test flight took place from Palestine, Texas on July 5, 2000, lasting about 8 h. The instrument performed well and landed safely in Stiles, Texas with little damage. The resulting measurements of the NUV nightglow are compared with previous measurements from sounding rockets and balloons.