Least squares best-fit circles (with applications to Mohr''s diagram)

Procedures are outlined for the selection of a least squares best-fit circle to data points defined by rectangular Cartesian coordinates. Equations are derived to allow fitting of circles centered on the x-axis as well as off-axis Mohr circles. These procedures are applicable to the estimation of second-order tensors such as stress and strain by means of Mohr's diagram.     

Stable boron isotope fractionation between dissolved B(OH)3 and B(OH)4

The stable boron isotope ratio (¹¹B/¹⁰B) in marine carbonates is used as a paleo-pH recorder and is one of the most promising paleo-carbonate chemistry proxies. Understanding the thermodynamic basis of the proxy is of fundamental importance, including knowledge on the equilibrium fractionation factor between dissolved boric acid, B(OH)₃, and borate ion, B(OH)₄⁻ (, hereafter α_(B3-B4)). However, this factor has hitherto not been determined experimentally and a theoretically calculated value (Kakihana and Kotaka, 1977, hereafter KK77) has therefore been widely used. I examine the calculations underlying this value. Using the same spectroscopic data and methods as KK77, I calculate the same α_(B3−B4) = 1.0193 at 300 K. Unfortunately, it turns out that in general the result is sensitive to the experimentally determined vibrational frequencies and the theoretical methods used to calculate the molecular forces. Using analytical techniques and ab initio molecular orbital theory, the outcome for α_(B3-B4) varies between ∼1.020 and ∼1.050 at 300 K. However, several arguments suggest that α_(B3-B4) ? 1.030. Measured isotopic shifts in various ¹⁰B-, ²D-, and ¹⁸O-labeled isotopomers do not provide a constraint on stable boron isotope fractionation. I conclude that in order to anchor the fundamentals of the boron pH proxy, experimental work is required. The critics of the boron pH proxy should note, however, that uncertainties in α_(B3-B4) do not bias pH reconstructions provided that organism-specific calibrations are used.     

Combining in situ isotopic, trace element and textural analyses of quartz from four magmatic-hydrothermal ore deposits

This study couples in situ ¹⁶O, ¹⁷O and ¹⁸O isotope and in situ trace element analyses to investigate and characterize the geochemical and textural complexity of magmatic-hydrothermal quartz crystals. Euhedral quartz crystals contemporaneous with mineralization were obtained from four magmatic-hydrothermal ore deposits: El Indio Au–Ag–Cu deposit; Summitville Au–Ag–Cu deposit; North Parkes Cu–Au deposit and Kingsgate quartz-Mo–Bi–W deposit. The internal features of the crystals were imaged using cathodoluminescence and qualitative electron microprobe maps. Quantitative isotopic data were collected in situ using 157 nm laser ablation inductively coupled plasma mass spectrometry (for 40 trace elements in quartz) and sensitive high-resolution ion microprobe (for 3 isotopes in quartz). Imaging revealed fine oscillatory zoning, sector zoning, complex “macromosaic” textures and hidden xenocrystic cores. In situ oxygen isotope analyses revealed a δ¹⁸O range of up to 12.4 ± 0.3 ‰ in a single crystal—the largest isotopic range ever ascribed to oscillatory zonation in quartz. Some of these crystals contain a heavier δ¹⁸O signature than expected by existing models. While sector-zoned crystals exhibited strong trace element variations between faces, no evidence for anisotropic isotope fractionation was found. We found: (1) isotopic heterogeneity in hydrothermal quartz crystals is common and precludes provenance analysis (e.g., δD–δ¹⁸O) using bulk analytical techniques, (2) the trace element signature of quartz is not an effective pathfinder toward noble metal mineralization and (3) in three of the four samples, both textural and isotopic data indicate non-equilibrium deposition of quartz.     

Winter Arctic Amplification at the synoptic timescale, 1979–2018, its regional variation and response to tropical and extratropical variability

We investigate winter Arctic Amplification (AA) on synoptic timescales and at regional scales using a daily version of the Arctic Amplification Index (AAI) and examine causes on a synoptic scale. The persistence, frequency and intensity of high AAI events show significant increases over the Arctic. Similarly, low AAI events are decreasing in frequency, persistence and intensity. In both cases, there are regional variations in these trends, in terms of significance and timing. Significant trends in increasing persistence, frequency and intensity of high AAI events in winter are concentrated in the period 2000–2009, with few significant trends before and after this. There are some decreases in sea-ice concentration in response to synoptic-scale AA events and these AA events can contribute to the decadal trends in AA found in other studies. A sectoral analysis of the Arctic indicates that in the Beaufort–Chukchi and East Siberian–Laptev Seas, synoptic scale high AAI events can be driven by tropical teleconnections while in other Arctic sectors, it is the intrusion of moisture-transporting synoptic cyclones into the Arctic that is most important in synoptic-scale AA. The presence of Rossby wave breaking during high AAI events is indicative of forcing from lower latitudes, modulated by variations in the jet stream. An important conclusion is that the increased persistence, frequency and intensity of synoptic-scale high AAI events make significant contributions to the interannual trend in AA.

     

Oxygen isotopic constraints on the origin and parent bodies of eucrites, diogenites, and howardites

A few eucrites have anomalous oxygen isotopic compositions. To help understand their origin and identify additional samples, we have analyzed the oxygen isotopic compositions of 18 eucrites and four diogenites. Except for five eucrites, these meteorites have Δ¹⁷O values that lie within 2σ of their mean value viz., −0.242 ± 0.016‰, consistent with igneous isotopic homogenization of Vesta. The five exceptional eucrites—NWA 1240, Pasamonte (both clast and matrix samples), PCA 91007, A-881394, and Ibitira—have Δ¹⁷O values that lie, respectively, 4σ, 5σ, 5σ, 15σ, and 21σ away from this mean value. NWA 1240 has a δ¹⁸O value that is 5σ below the mean eucrite value. Four of the five outliers are unbrecciated and unshocked basaltic eucrites, like NWA 011, the first eucrite found to have an anomalous oxygen isotopic composition. The fifth outlier, Pasamonte, is composed almost entirely of unequilibrated basaltic clasts. Published chemical data for the six eucrites with anomalous oxygen isotopic compositions (including NWA 011) exclude contamination by chondritic projectiles as a source of the oxygen anomalies. Only NWA 011 has an anomalous Fe/Mn ratio, but several anomalous eucrites have exceptional Na, Ti, or Cr concentrations. We infer that the six anomalous eucrites are probably derived from five distinct Vesta-like parent bodies (Pasamonte and PCA 91007 could come from one body). These anomalous eucrites, like the isotopically normal, unbrecciated eucrites with 4.48 Gyr Ar-Ar ages, are probably deficient in brecciation and shock effects because they were sequestered in small asteroids (10 km diameter) during the Late Heavy Bombardment following ejection from Vesta-like bodies. The preservation of Vesta’s crust and the lack of deeply buried samples from the hypothesized Vesta-like bodies are consistent with the removal of these bodies from the asteroid belt by gravitational perturbations from planets and protoplanets, rather than by collisional grinding.     

The nature of cosmic ray sources

The limitations on the nature of cosmic ray acceleration regions and processes, as deduced from cosmic ray measurements and propagation studies, are reviewed. The power requirements for these acceleration regions are estimated from measurements of the local cosmic ray energy density, anisotropy and spallation-deduced pathlength. Possible constraints on the acceleration spectrum of the cosmic rays and on a charge dependence of the acceleration process, implied by the measured cosmic ray spectrum and composition, are considered. Various suggested sources and processes of cosmic ray acceleration are discussed in the light of these limitations.Astrophysics and Space Science Review Paper.     

The formation of Mg i 4571 Å in the solar atmosphere

An analysis of the 4571 Å line of neutral magnesium is presented in which one-dimensional macroscopic velocity fields are included. It is shown that gradients over restricted heights in the vertical and horizontal components of the velocity field of order -0.005 s^–1 and -0.004 s^–1 (such that velocity towards the observer decreases as height increases), respectively, result in asymmetries in the computed line profile similar to those observed. The heights in the solar atmosphere at which these velocity gradients exist are shown to be very critical in reproducing the observations. It was found that the best results were obtained when the gradients existed in the height range from 200 km to 300 km below the temperature minimum. The results indicate that for the Mg i 4571 Å line model calculations that do not include one-dimensional flow velocities may safely be compared with frequency-averaged observations.     

The opposition and tilt effects of Saturn’s rings from HST observations

The two major factors contributing to the opposition brightening of Saturn’s rings are (i) the intrinsic brightening of particles due to coherent backscattering and/or shadow hiding on their surfaces, and (ii) the reduced interparticle shadowing when the solar phase angle α → 0°. We utilize the extensive set of Hubble Space Telescope observations (Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 199–223) for different elevation angles B and wavelengths λ to disentangle these contributions. We assume that the intrinsic contribution is independent of B, so that any B dependence of the phase curves is due to interparticle shadowing, which must also act similarly for all λ’s. Our study complements that of Poulet et al. (Poulet, F., Cuzzi, J.N., French, R.G., Dones, L. [2002]. Icarus 158, 224), who used a subset of data for a single B ∼ 10°, and the French et al. (French, R.G., Verbiscer, A., Salo, H., McGhee, C.A., Dones, L. [2007b] PASP 119, 623–642) study for the B ∼ 23° data set that included exact opposition. We construct a grid of dynamical/photometric simulation models, with the method of Salo and Karjalainen (Salo and Karjalainen [2003]. Icarus 164, 428–460), and use these simulations to fit the elevation-dependent part of opposition brightening. Eliminating the modeled interparticle component yields the intrinsic contribution to the opposition effect: for the B and A rings it is almost entirely due to coherent backscattering; for the C ring, an intraparticle shadow hiding contribution may also be present.Based on our simulations, the width of the interparticle shadowing effect is roughly proportional to B. This follows from the observation that as B decreases, the scattering is primarily from the rarefied low filling factor upper ring layers, whereas at larger B’s the dense inner parts are visible. Vertical segregation of particle sizes further enhances this effect. The elevation angle dependence of interparticle shadowing also explains most of the B ring tilt effect (the increase of brightness with elevation). From comparison of the magnitude of the tilt effect at different filters, we show that multiple scattering can account for at most a 10% brightness increase as B → 26°, whereas the remaining 20% brightening is due to a variable degree of interparticle shadowing. The negative tilt effect of the middle A ring is well explained by the the same self-gravity wake models that account for the observed A ring azimuthal brightness asymmetry (Salo, H., Karjalainen, R., French, R.G. [2004]. Icarus 170, 70–90; French, R.G., Salo, H., McGhee, C.A., Dones, L. [2007]. Icarus 189, 493–522).     

The Fossil Starburst in M82

We present high-resolution optical and near-infrared HST observations of two adjacent regions in the fossil starburst region inM82, M82 B1 and B2. The presence of both the active and the fossil starburst in M82 provides a unique physical environment to study the stellar and dynamical evolution of star cluster systems. The cluster population in B2 is more heavily affected by internal extinction than that in B1, amounting to an excess extinction in B2 of A_V,excess≃1.1±0.3 mag. Preliminary age estimates date the cluster population in the fossil starburst between ∼2× 10⁸ and ∼ 10⁹ years. The radial luminosity profiles of the brightest clusters are more closely approximated by power laws than by a Gaussian model, in particular in their wings, which favours a slow star formation scenario. This revised version was published online in July 2006 with corrections to the Cover Date.