Binary near-Earth asteroid formation: Rubble pile model of tidal disruptions

We present numerical simulations of near-Earth asteroid (NEA) tidal disruption resulting in bound, mutually orbiting systems. Using a rubble pile model we have constrained the relative likelihoods for possible physical and dynamical properties of the binaries created. Overall 110,500 simulations were run, with each body consisting of ∼1000 particles. The encounter parameters of close approach distance and velocity were varied, as were the bodies' spin, elongation and spin axis direction. The binary production rate increases for closer encounters, at lower speeds, for more elongated bodies, and for bodies with greater spin. The semimajor axes for resultant binaries are peaked between 5 to 20 primary radii, and there is an overall trend for high eccentricity, with 97% of binaries having e > 0.1. The secondary-to-primary size ratios of the simulated binaries are peaked between 0.1 and 0.2, similar to trends among observed asteroid binaries. The spin rates of the primary bodies are narrowly distributed between 3.5- and 6-h periods, whereas the secondaries' periods are more evenly distributed and can exceed 15-h periods. The spin axes of the primary bodies are very closely aligned with the angular momenta of the binary orbits, whereas the secondary spin axes are nearly random. The shapes of the primaries show a large distribution of axis ratios, where those with low elongation (ratio of long and short axis) are both oblate and prolate, and nearly all with large elongation are prolate. This work presents results that suggest tidal disruption of gravitational aggregates can make binaries physically similar to those currently observed in the NEA population. As well, tidal disruption may create an equal number of binaries with qualities different from those observed, mostly binaries with large separation and with elongated primaries.     

A fast method for finding bound systems in numerical simulations: Results from the formation of asteroid binaries

We present a new code (companion) that identifies bound systems of particles in O(NlogN) time. Simple binaries consisting of pairs of mutually bound particles and complex hierarchies consisting of collections of mutually bound particles are identifiable with this code. In comparison, brute force binary search methods scale as O(N²) while full hierarchy searches can be as expensive as O(N³), making analysis highly inefficient for multiple data sets with N?¹⁰³. A simple test case is provided to illustrate the method. Timing tests demonstrating O(NlogN) scaling with the new code on real data are presented. We apply our method to data from asteroid satellite simulations [Durda et al., 2004. Icarus 167, 382-396; Erratum: Icarus 170, 242; reprinted article: Icarus 170, 243-257] and note interesting multi-particle configurations. The code is available at http://www.astro.umd.edu/zoe/companion/ and is distributed under the terms and conditions of the GNU Public License.     

Snow density variations: consequences for ground‐penetrating radar

Reliable hydrological forecasts of snowmelt runoff are of major importance for many areas. Ground‐penetrating radar (GPR) measurements are used to assess snowpack water equivalent for planning of hydropower production in northern Sweden. The travel time of the radar pulse through the snow cover is recorded and converted to snow water equivalent (SWE) using a constant snowpack mean density from the drainage basin studied. In this paper we improve the method to estimate SWE by introducing a depth‐dependent snowpack density. We used 6 years measurements of peak snow depth and snowpack mean density at 11 locations in the Swedish mountains. The original method systematically overestimates the SWE at shallow depths (+25% for 0·5 m) and underestimates the SWE at large depths (?35% for 2·0 m). A large improvement was obtained by introducing a depth–density relation based on average conditions for several years, whereas refining this by using separate relations for individual years yielded a smaller improvement. The SWE estimates were substantially improved for thick snow covers, reducing the average error from 162 ± 23 mm to 53 ± 10 mm for depth range 1·2–2·0 m. Consequently, the introduction of a depth‐dependent snow density yields substantial improvements of the accuracy in SWE values calculated from GPR data. Copyright © 2005 John Wiley & Sons, Ltd.     

Multibeam Observations of Mine Burial Near Clearwater, FL, Including Comparisons to Predictions of Wave-Induced Burial

A Kongsberg Simrad EM 3000 multibeam sonar (Kongsberg Simrad, Kongsberg, Norway) was used to conduct a set of six repeat high-resolution bathymetric surveys west of Indian Rocks Beach (IRB), just to the south of Clearwater, FL, between January and March 2003, to observe in situ scour and burial of instrumented inert mines and mine-like cylinders. Three closely located study sites were chosen: two fine-sand sites, a shallow one located in 13 m of water depth and a deep site located in 14 m of water depth; and a coarse-sand site in 13 m. Results from these surveys indicate that mines deployed in fine sand are nearly buried within two months of deployment (i.e., they sunk 74.5% or more below the ambient seafloor depth). Mines deployed in coarse sand showed a lesser amount of scour, burying until they present roughly the same hydrodynamic roughness as the surrounding rippled bedforms. These data were also used to test the validity of the Virginia Institute of Marine Science (VIMS, Gloucester Point, VA) 2-D burial model. The model worked well in areas of fine sand, sufficiently predicting burial over the course of the experiment. In the area of coarse sand, the model greatly overpredicted the amount of burial. This is believed to be due to the presence of rippled bedforms around the mines, which affect local bottom morphodynamics and are not accounted for in the model, an issue currently being addressed by the modelers. This paper focuses specifically on two instrumented mines: an acoustic mine located in fine sand and an optical instrumented mine located in coarse sand.     

Predicting Seabed Burial of Cylinders by Wave-Induced Scour: Application to the Sandy Inner Shelf Off Florida and Massachusetts

A simple parameterized model for wave-induced burial of mine-like cylinders as a function of grain-size, time-varying, wave orbital velocity and mine diameter was implemented and assessed against results from inert instrumented mines placed off the Indian Rocks Beach (IRB, FL), and off the Martha's vineyard coastal observatory (MVCO, Edgartown, MA). The steady flow scour parameters provided by Whitehouse (1998) for self-settling cylinders worked well for predicting burial by depth below the ambient seabed for (0.5 m) diameter mines in fine sand at both sites. By including or excluding scour pit infilling, a range of percent burial by surface area was predicted that was also consistent with observations. Rapid scour pit infilling was often seen at MVCO but never at IRB, suggesting that the environmental presence of fine sediment plays a key role in promoting infilling. Overprediction of mine scour in coarse sand was corrected by assuming a mine within a field of large ripples buries only until it generates no more turbulence than that produced by surrounding bedforms. The feasibility of using a regional wave model to predict mine burial in both hindcast and real-time forecast mode was tested using the National Oceanic and Atmospheric Administration (NOAA, Washington, DC) WaveWatch 3 (WW3) model. Hindcast waves were adequate for useful operational forcing of mine burial predictions, but five-day wave forecasts introduced large errors. This investigation was part of a larger effort to develop simple yet reliable predictions of mine burial suitable for addressing the operational needs of the U.S. Navy.     

Partitioning of rainfall in a eucalypt forest and pine plantation in southeastern australia: II stemflow and factors affecting stemflow in a dry sclerophyll eucalypt forest and a pinus radiata plantation

Stemflow of a dry sclerophyll eucalypt forest and a nearby Pinus radiata plantation was studied on a rainfall event basis. The stemflow yields of the forests are quantified, compared, and presented on an annual basis for four years. Yields of the individual eucalypt species are compared and the tree characteristics responsible for the yield differences are discussed. The influence of event size, type, and season on stemflow are also shown. Rainfall angle is shown to have a significant effect on stemflow yield.     

Dust motions in quasi-statically charged binary asteroid systems

In this paper, we discuss dust motion and investigate possible mass transfer of charged particles in a binary asteroid system, in which the asteroids are electrically charged due to solar radiation. The surface potential of the asteroids is assumed to be a piecewise function, with positive potential on the sunlit half and negative potential on the shadow half. We derive the nonautonomous equations of motion for charged particles and an analytic representation for their lofting conditions. Particle trajectories and temporary relative equilibria are examined in relation to their moving forbidden regions, a concept we define and discuss. Finally, we use a Monte Carlo simulation for a case study on mass transfer and loss rates between the asteroids.