Sensitivity and clay mineralogy of weathered tephra-derived soil materials in the Tauranga region

Soil sensitivity is defined as the ratio of peak to remoulded shear strength. Problem soil materials are those that show large strength losses on disruption, resulting in catastrophic failure, liquefaction and long run-out distances. This study focussed on sensitive, weathered, mainly tephra-fall derived soils of mid-Pleistocene age in the Tauranga region. The liquefiable character of these soils is well known, but little detailed study has been directed towards the reasons for sensitivity. The objective of this work was to examine soil sensitivity by investigating geomechanical properties, clay mineralogy, and microfabric, and to determine how these factors combine to develop sensitivity. To achieve these objectives a combination of both field and laboratory investigations was undertaken. Field investigations indicated that sensitive soils are common in the Tauranga region. Sampling was undertaken at sites in Tauriko and Otumoetai. Selected samples ranged across high (76) and low (≈8) field sensitivity. Stratigraphically, samples from Otumoetai lie below the Rangitawa Tephra (ca. 0.34 Ma), and those from Tauriko underlie the Te Ranga Ignimbrite (ca. 0.27 Ma). Geomechanical investigation revealed that the sensitive soils had high moisture contents (gt; 60 %), low dry bulk density (lt; 966kg m-3), and high porosity (gt; 60 %). Liquidity index values ranged between 0.27 and 2.41. Plasticity index values ranged from 13.2 % to 42.7 %, with all samples plotting below the A line. Strength tests indicated effective friction angles from 25.7 to 38.5 , effective cohesion from 4.7 kPa to 34.5 kPa, residual friction angles of 19.34 to 33.18 , and cohesions of 0 kPa to 4.87 kPa. Remoulded vane shear strengths ranged between 1 kPa and 36 kPa. Clay minerals were dominantly hydrated halloysite. Scanning electron microscopy indicated that clay morphology was in the form of hollow tubes, spheres, plates, and platy vermiforms ('books'). Tubes and spheres represent characteristic forms of halloysite in soils, plates are less common, and books have never previously been observed. Hence, these books represent a new morphology for halloysite. Individual plates in each of the books appear to show structural Fe enrichment (~5.2 %). This enrichment indicates that Fe had replaced Al in octahedral positions reducing the mismatch with the tetrahedral sheet, lessening layer curvature and thus generating flat plates. All microfabrics were continuous with larger sand and silt grains supported in a background of clay minerals. Microfabrics ranged from extremely open with components being loosely packed to those which were dense and tightly packed. A feature common to all structural types was an abundance of extremely small pores (lt; 20 μm) which are capable of tightly retaining water. The loosely packed microfabrics had void ratios that allowed moisture content to exceed liquid limits, producing a liquidity index gt; 1. These open microfabrics are probably a result of quick burial by subsequent pyroclastic beds; hence weathering to clays occurred as a process of subsurface diagenesis. Dense microfabrics with low void ratios and high liquid limits did not have liquidity indexes gt; 1. These dense microfabrics arose as a result of the deposits being at, or near, the land surface for a considerable time, thus allowing strong pedogenic processes to occur, which promoted clay formation and clay migration (illuviation) and reduced void ratios. Liquidity index was a major control on remoulded strength and sensitivity. Liquidity index is controlled by clay type and content, void ratio, and natural moisture content. When remoulded, structures with natural moisture contents exceeding the liquid limit release a large amount of water, which both dilutes the plasticity of binding clays and supports grains and broken aggregates of clay, allowing the material to flow. The development of sensitivity with low remoulded strength requires a number of factors. These include: a void ratio that is sufficiently high to allow natural moisture content to exceed the liquid limit; the presence of halloysite, which encourages samples to retain a coherent structure when saturated and to ensure the liquid limit remains sufficiently low so that it can be exceeded by natural moisture content; and a saturated environment, which ensures the liquid limit is exceeded

The Insignificance of P-R Drag in Detectable Extrasolar Planetesimal Belts

This paper considers a simple model in which dust produced in a planetesimal belt migrates in toward the star due to P-R drag suffering destructive collisions with other dust grains on the way. Assuming the dust is all of the same size, the resulting surface density distribution can be derived analytically and depends only on the parameter eta0=5000tau*sqrt(Mstar/r)/beta; this parameter can be determined observationally with the hypothesis that beta=0.5. For massive belts in which eta0>>1 dust is confined to the planetesimal belt, while the surface density of more tenuous belts, in which eta0<<1, is constant with distance from the star. The emission spectrum of dust from planetesimal belts at different distances from different mass stars shows that the dust belts which have been detected to date should have eta0>>1; dust belts with eta0<<1 are hard to detect as they are much fainter than the stellar photosphere. This is confirmed for a sample of 37 debris disk candidates for which eta0 was determined to be >10. This means that these disks are so massive that mutual collisions prevent dust from reaching the inner regions of these systems and P-R drag can be ignored when studying their dynamics. Models for the formation of structure in debris disks by the trapping of particles into planetary resonances by P-R drag should be reconsidered. However, since collisions do not halt 100% of the dust, this means that in the absence of planetary companions debris disk systems should be populated by small quantities of hot dust which may be detectable in the mid-IR. Even in disks with eta0<<1 the temperature of dust emission is shown to be a reliable tracer of the planetesimal distribution meaning that inner holes in the dust distribution imply a lack of colliding planetesimals in the inner regions.Comment: 6 pages. Accepted by A&

Landslides in sensitive soils, Tauranga, New Zealand.

In the Tauranga region sensitive soil failures commonly occur after heavy rainfall events, causing considerable infrastructure damage. Several notable landslides include a large failure at Bramley Drive, Omokoroa in 1979, the Ruahihi Canal collapse in 1981, and numerous landslides in May 2005; recently the Bramley Drive scarp was reactivated in 2011. These failures are associated with materials loosely classified as the Pahoia Tephras - a mixture of rhyolitic pyroclastic deposits of approximately 1 Ma. The common link with extreme rainfall events suggests a pore water pressure control on the initiation of these failures. Recent research on the structure of the soils shows a dominance of halloysite clay minerals packed loosely in arrangements with high porosity (51 – 77 %), but with almost entirely micropores. This leads us to conclude that the permeability is very low, and the materials remain continuously wet. The formation of halloysite is encouraged by a wet environment with no episodes of drying, supporting this assumption. A high-resolution CPT trace at Bramley Drive indicates induced pore water pressures rising steadily to a peak at approximately 25 m depth; this depth coincides with the base of the landslide scarp. We infer that elevated pore water pressures develop within this single, thick aquifer, triggering failure through reduced effective stresses. The inactive halloysite clay mineral results in low plasticity indices (13 – 44 %) and hence high liquidity indices (1.2 – 2.4) due to the saturated pore space; remoulding following failure is sudden and dramatic and results in large debris runout distances

Collisional Evolution of Irregular Satellite Swarms: Detectable Dust around Solar System and Extrasolar Planets

Since the 1980's it has been becoming increasingly clear that the Solar System's irregular satellites are collisionally evolved. We derive a general model for the collisional evolution of an irregular satellite swarm and apply it to the Solar System and extrasolar planets. Our model reproduces the Solar System's complement of observed irregulars well, and suggests that the competition between grain-grain collisions and Poynting-Robertson (PR) drag helps set the fate of the dust. Because swarm collision rates decrease over time the main dust sink can change with time, and may help unravel the accretion history of synchronously rotating regular satellites that show brightness asymmetries. Some level of dust must be present on AU scales around the Solar System's giant planets, which we predict may be at detectable levels. We also predict whether dust produced by extrasolar circumplanetary swarms can be detected. The coronagraphic instruments on JWST will have the ability to detect the dust generated by these swarms, which are most detectable around planets that orbit at tens of AU from the youngest stars. Because the collisional decay of swarms is relatively insensitive to planet mass, swarms can be much brighter than their host planets and allow discovery of Neptune-mass planets that would otherwise remain invisible. This dust may have already been detected. The observations of the planet Fomalhaut b can be explained as scattered light from dust produced by the collisional decay of an irregular satellite swarm around a 10 Earth-mass planet. Such a swarm comprises about 5 Lunar masses worth of irregular satellites. Finally, we consider what happens if Fomalhaut b passes through Fomalhaut's main debris ring, which allows the circumplanetary swarm to be replenished through collisions with ring planetesimals. (abridged)Comment: accepted to MNRA

Post-main sequence evolution of A star debris discs

While the population of main sequence debris discs is well constrained, little is known about debris discs around evolved stars. This paper provides a theoretical framework considering the effects of stellar evolution on debris discs, particularly the production and loss of dust within them. Here we repeat a steady state model fit to disc evolution statistics for main sequence A stars, this time using realistic grain optical properties, then evolve that population to consider its detectability at later epochs. Our model predicts that debris discs around giant stars are harder to detect than on the main sequence because radiation pressure is more effective at removing small dust around higher luminosity stars. Just 12% of first ascent giants within 100pc are predicted to have discs detectable with Herschel at 160um. However this is subject to the uncertain effect of sublimation on the disc, which we propose can thus be constrained with such observations. Our model also finds that the rapid decline in stellar luminosity results in only very young white dwarfs having luminous discs. As such systems are on average at larger distances they are hard to detect, but we predict that the stellar parameters most likely to yield a disc detection are a white dwarf at 200pc with cooling age of 0.1Myr, in line with observations of the Helix Nebula. Our model does not predict close-in (<0.01AU) dust, as observed for some white dwarfs, however we find that stellar wind drag leaves significant mass (~10^{-2}Msolar), in bodies up to ~10m in diameter, inside the disc at the end of the AGB phase which may replenish these discs

Are debris disks self-stirred?

This paper considers the evidence that debris disks are self-stirred by the formation of Plutos. A model for the dust produced during self-stirring is applied to statistics for A stars. As there is no significant difference between excesses of A-stars <50Myr old, we focus on reproducing the broad trends, the "rise and fall" of the fraction of stars with excesses. Using a population model, we find that the statistics and trends can be reproduced with a self-stirring model of planetesimal belts with radii distributed between 15-120AU. Disks must have this 15AU minimum radius to show a peak in disk fraction, rather than a monotonic decline. Populations of extended disks with fixed inner and/or outer radii fail to fit the statistics, due mainly to the slow 70um evolution as stirring moves further out in the disk. This conclusion, that debris disks are narrow belts, is independent of the significance of 24um trends for young A-stars. We show that the statistics can also be reproduced with a model in which disks are stirred by secular perturbations from a nearby eccentric planet. Detailed imaging is therefore the best way to characterise the stirring mechanism. From a more detailed look at beta Pictoris Moving Group and TW Hydrae Association A-stars we find that the disk around beta Pictoris is likely the result of secular stirring by the proposed planet at ~10AU; the structure of the HR 4796A disk also points to sculpting by a planet. The two other stars with disks, HR 7012 and eta Tel, possess transient hot dust, though the outer eta Tel disk is consistent with a self-stirred origin. Planet formation provides a natural explanation for the belt-like nature of debris disks, with inner regions cleared by planets that may also stir the disk, and the outer edges set by where planetesimals can form. [abridged]Comment: Accepted to MNRA

A general model of resonance capture in planetary systems: First and second order resonances

Mean motion resonances are a common feature of both our own Solar System and of extrasolar planetary systems. Bodies can be trapped in resonance when their orbital semi-major axes change, for instance when they migrate through a protoplanetary disc. We use a Hamiltonian model to thoroughly investigate the capture behaviour for first and second order resonances. Using this method, all resonances of the same order can be described by one equation, with applications to specific resonances by appropriate scaling. We focus on the limit where one body is a massless test particle and the other a massive planet. We quantify how the the probability of capture into a resonance depends on the relative migration rate of the planet and particle, and the particle's eccentricity. Resonant capture fails for high migration rates, and has decreasing probability for higher eccentricities. More massive planets can capture particles at higher eccentricities and migration rates. We also calculate libration amplitudes and the offset of the libration centres for captured particles, and the change in eccentricity if capture does not occur. Libration amplitudes are higher for larger initial eccentricity. The model allows for a complete description of a particle's behaviour as it successively encounters several resonances. We discuss implications for several scenarios: (i) Planet migration through gas discs trapping other planets or planetesimals in resonances. (ii) Planet migration through a debris disc. (iii) Dust migration through PR drag. The Hamiltonian model will allow quick interpretation of the resonant properties of extrasolar planets and Kuiper Belt Objects, and will allow synthetic images of debris disc structures to be quickly generated, which will be useful for predicting and interpreting disc images made with ALMA, Darwin/TPF or similar missions. [Abridged]Comment: 19 pages, 14 figures; accepted to MNRA

Mercurian impact ejecta: Meterorites and mantle

We have examined the fate of impact ejecta liberated from the surface of Mercury due to impacts by comets or asteroids, in order to study (1) meteorite transfer to Earth, and (2) re-accumulation of an expelled mantle in giant-impact scenarios seeking to explain Mercury's large core. In the context of meteorite transfer, we note that Mercury's impact ejecta leave the planet's surface much faster (on average) than other planet's in the Solar System because it is the only planet where impact speeds routinely range from 5-20 times the planet's escape speed. Thus, a large fraction of mercurian ejecta may reach heliocentric orbit with speeds sufficiently high for Earth-crossing orbits to exist immediately after impact, resulting in larger fractions of the ejecta reaching Earth as meteorites. We calculate the delivery rate to Earth on a time scale of 30 Myr and show that several percent of the high-speed ejecta reach Earth (a factor of -3 less than typical launches from Mars); this is one to two orders of magnitude more efficient than previous estimates. Similar quantities of material reach Venus. These calculations also yield measurements of the re-accretion time scale of material ejected from Mercury in a putative giant impact (assuming gravity is dominant). For mercurian ejecta escaping the gravitational reach of the planet with excess speeds equal to Mercury's escape speed, about one third of ejecta re-accretes in as little as 2 Myr. Thus collisional stripping of a silicate proto-mercurian mantle can only work effectively if the liberated mantle material remains in small enough particles that radiation forces can drag them into the Sun on time scale of a few million years, or Mercury would simply re-accrete the material.Comment: 14 pages. Submitted to Meteoritics and Planetary Scienc

Collisional dust avalanches in debris discs

We quantitatively investigate how collisional avalanches may developin debris discs as the result of the initial break-up of a planetesimal or comet-like object, triggering a collisional chain reaction due to outward escaping small dust grains. We use a specifically developed numerical code that follows both the spatial distribution of the dust grains and the evolution of their size-frequency distribution due to collisions. We investigate how strongly avalanche propagation depends on different parameters (e.g., amount of dust released in the initial break-up, collisional properties of dust grains and their distribution in the disc). Our simulations show that avalanches evolve on timescales of ~1000 years, propagating outwards following a spiral-like pattern, and that their amplitude exponentially depends on the number density of dust grains in the system. We estimate a probability for witnessing an avalanche event as a function of disc densities, for a gas-free case around an A-type star, and find that features created by avalanche propagation can lead to observable asymmetries for dusty systems with a beta Pictoris-like dust content or higher. Characteristic observable features include: (i) a brightness asymmetry of the two sides for a disc viewed edge-on, and (ii) a one-armed open spiral or a lumpy structure in the case of face-on orientation. A possible system in which avalanche-induced structures might have been observed is the edge-on seen debris disc around HD32297, which displays a strong luminosity difference between its two sides.Comment: 18 pages, 19 figures; has been accepted for publication in Astronomy and Astrophysics, section 6. Interstellar and circumstellar matter. The official date of acceptance is 29/08/200

Observing planet-disk interaction in debris disks

Context. Structures in debris disks induced by planetdisk interaction are promising to provide valuable constraints on the existence and properties of embedded planets. Aims. We investigate the observability of structures in debris disks induced by planet-disk interaction. Methods. The observability of debris disks with the Atacama Large Millimeter/submillimeter Array (ALMA) is studied on the basis of a simple analytical disk model. Furthermore, N-body simulations are used to model the spatial dust distribution in debris disks under the influence of planet-disk interaction. Images at optical scattered light to millimeter thermal re-emission are computed. Available information about the expected capabilities of ALMA and the James Webb Space Telescope (JWST) are used to investigate the observability of characteristic disk structures through spatially resolved imaging. Results. Planet-disk interaction can result in prominent structures. This provides the opportunity of detecting and characterizing extrasolar planets in a range of masses and radial distances from the star that is not accessible to other techniques. Facilities that will be available in the near future are shown to provide the capabilities to spatially resolve and characterize structures in debris disks. Limitations are revealed and suggestions for possible instrument setups and observing strategies are given. In particular, ALMA is limited by its sensitivity to surface brightness, which requires a trade-off between sensitivity and spatial resolution. Space-based midinfrared observations will be able to detect and spatially resolve regions in debris disks even at a distance of several tens of AU from the star, where the emission from debris disks in this wavelength range is expected to be low. [Abridged]Comment: 16 pages, 10 figures, accepted by A&