1

## Numerical predictions of surface effects during the 2029 close

arXiv:1408.0168v1  [astro-ph.EP]  1 Aug 2014

## approach of asteroid 99942 Apophis

2

Yang Yu /star †

3

Derek C. Richardson †

4

Patrick Michel ‡

5

Stephen R. Schwartz ‡†

6

Ronald-Louis Ballouz †

7

/star School of Aerospace

8

Tsinghua University

9

Beijing 100084 China

10

† Department of Astronomy

11

University of Maryland

12

College Park MD 20740-2421 United States

13

‡ Lagrange Laboratory

14

University of Nice Sophia Antipolis, CNRS

15

Observatoire de la C o te dAzur, C.S. 34229, 06304 Nice Cedex 4, France

16

Printed August 4, 2014

17

Submitted to Icarus

18

61 manuscript pages

19

14 figures including 11 in color (online version only)

20

## Yang Yu et al.

- Proposed running page head: Surface effects on 99942 Apophis 21
- Please address all editorial correspondence and proofs to: Derek C. Richardson 22
- Derek C. Richardson 23
- Department of Astronomy 24
- Computer and Space Sciences Building 25
- Stadium Drive 26
- University of Maryland 27
- College Park MD 20742-2421 28
- Tel: 301-405-8786 29
- 30

Fax: 301-314-9067

- E-mail: dcr@astro.umd.edu 31

## Abstract 32

- Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small33 body encounter events in the near future and offers a good opportunity for in situ exploration 34 to determine the asteroid's surface properties and measure any tidal effects that might al35 ter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid 36 researchers due to its important implications for interpreting surface observations, including 37 space weathering effects. This paper provides a prediction for the tidal effects during the 38 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement 39 will occur. The potential shape modification of the object due to the tidal encounter is first 40 confirmed to be negligibly small with systematic simulations, thus only the external pertur41 bations are taken into account for this work (despite this, seismic shaking induced by shifting 42 blocks might still play a weak role and we will look into this mechanism in future work). 43 A two-stage approach is developed to model the responses of asteroid surface particles (the 44 regolith) based on the soft-sphere implementation of the parallel N -body gravity tree code 45 pkdgrav . A full-body model of Apophis is sent past the Earth on the predicted trajectory 46 to generate the data of all forces acting at a target point on the surface. A sandpile con47 structed in the local frame is then used to approximate the regolith materials; all the forces 48 the sandpile feels during the encounter are imposed as external perturbations to mimic the 49 regolith's behavior in the full scenario. The local mechanical environment on the asteroid 50 surface is represented in detail, leading to an estimation of the change in global surface 51 environment due to the encounter. Typical patterns of perturbation are presented that de52
- pend on the asteroid orientation and sense of rotation at perigee. We find that catastrophic 53

## Yang Yu et al.

- avalanches of regolith materials may not occur during the 2029 encounter due to the small 54 level of tidal perturbation, although slight landslides might still be triggered in positions 55 where a sandpile's structure is weak. Simulations are performed at different locations on 56 Apophis' surface and with different body- and spin-axis orientations; the results show that 57 the small-scale avalanches are widely distributed and manifest independently of the asteroid 58 orientation and the sandpile location. We also include simulation results of much closer 59 encounters of the Apophis with Earth than what is predicted to occur in 2029, showing that 60 much more drastic resurfacing takes place in these cases. 61
- Keywords : Asteroids, dynamics; asteroids, surfaces; near-Earth objects; regoliths; tides, 62 solid body 63

## 1 INTRODUCTION 64

In this paper, we investigate numerically the behavior of granular material at the surface of 65 an asteroid during close approach to the Earth. We focus on the specific case of Asteroid 66 (99942) Apophis, which will come as close as 5 6 Earth radii on April 13th, 2029. We provide . 67 predictions about possible reshaping and spin-alteration of-and surface effects on-Apophis 68 during this passage, as a function of plausible properties of the constituent granular mate69 rial. Studies of possible future space missions to Apophis are underway, including one by 70 the French space agency CNES calling for international partners (e.g., Michel et al. 2012), 71 with the aim of observing this asteroid during the 2029 close encounter and characterizing 72 whether reshaping, spin-alteration, and/or surface motion occur. The numerical investiga73 tions presented here allow for estimation of the surface properties that could lead to any 74 observed motion (or absence of motion) during the actual encounter. 75

Apophis made a passage to the Earth at ∼ 2300 Earth radii in early 2013. At that time, 76 the Herschel space telescope was used to refine the determination of the asteroid's albedo 77 and size (M¨ uller et al. 2013). According to these observations, the albedo is estimated 78 to be about 0.23 and the longest dimension about 325 ± 15 m, which is somewhat larger 79 than previous estimates (270 ± 60 m, according to Delbo' et al. 2007). Concurrent radar 80 observations improved the astrometry of the asteroid, ruling out the possibility of a collision 81 with the Earth in 2036 to better than 1 part in 10 . 6 However, Wlodarczyk et al. (2013) 82 presented a possible path of risk for 2068. This finding has put off any crisis by ∼ 32 years 83 and makes exploring Apophis in 2029 to be more for scientific interest. To date, nothing is 84

- known about the asteroid's surface mechanical properties, and this is why its close passage in 85

Yang Yu et al.

- 2029 offers a great opportunity to visit it with a spacecraft, determine its surface properties, 86
- and, for the first time, observe potential modifications of the surface due to tidal effects. 87

88

And as Apophis approaches, it is likely that international interest in a possible mission will

- increase, since such close approaches of a large object are relatively rare. 89

The case for tidally induced resurfacing was made by Binzel et al. (2010; also see DeMeo 90 et al. 2013) and discussed by Nesvorn´ et al. (2010) to explain the spectral properties of y 91 near-Earth asteroids (NEAs) belonging to the Q taxonomic type, which appear to have fresh 92 (unweathered) surface colors. Dynamical studies of these objects found that those bodies 93 had a greater tendency to come close to the Earth, within the Earth-Moon distance, than 94 bodies of other classes in the past 500 kyr. The authors speculated that tidal effects during 95 these passages could be at the origin of surface material disturbance leading to the renewed 96 exposure of unweathered material. We leave a more general and detailed investigation of 97 this issue for future work, but if this result is true for those asteroids, it may also be true 98 for Apophis, which will approach Earth on Friday, April 13, 2029 no closer than about 99 29,500 km from the surface (i.e., 4 6 Earth radii, or 5 6 Earth radii from the center of the . . 100 planet; Giorgini et al. 2008). It is predicted to go over the mid-Atlantic, appearing to the 101 naked eye as a moderately bright point of light moving rapidly across the sky. Our aim 102 is to determine whether, depending on assumed mechanical properties, it could experience 103 surface particulate motions, reshaping, or spin-state alteration due to tidal forces caused by 104 Earth's gravity field. The classical Roche limit for a cohesionless fluid body of bulk density 105 2 4 g . / cm 3 to not be disrupted by tidal forces is ∼ 3 22 Earth radii, so we do not expect any . 106 violent events to occur during the rocky asteroid's 2029 encounter at 5 6 Earth radii. . 107

- The presence of granular material (or regolith) and boulders at the surface of small bodies 108

## Yang Yu et al.

has been demonstrated by space missions that visited or flew by asteroids in the last few 109 decades (e.g., Veverka et al. 2000; Fujiwara et al. 2006). It appears that all encountered 110 asteroids to date, from the largest one, the main belt Asteroid (4) Vesta by the Dawn 111 mission, to the smallest one, the NEA (25143) Itokawa, sampled by the Hayabusa mission, 112 are covered with some sort of regolith. In fact, thermal infrared observations support the 113 idea that most asteroids are covered with regolith, given their preferentially low thermal 114 inertia (Delbo' et al. 2007). There even seems to be a trend as a function of the asteroid's 115 size based on thermal inertia measurements: larger objects are expected to have a surface 116 covered by a layer of fine regolith, while smaller ones are expected to have a surface covered 117 by a layer of coarse regolith (Clark et al. 2002). This trend is consistent with observations 118 by the NEAR-Shoemaker spacecraft of the larger ( ∼ 17 km mean diameter) Eros, whose 119 surface is covered by a deep layer of very fine grains, and by the Hayabusa spacecraft of the 120 much smaller (320 m mean diameter) Itokawa, whose surface is covered by a thin layer of 121 coarse grains. However, interpretation of thermal inertia measurements must be made with 122 caution, as we do not yet have enough comparisons with actual asteroid surfaces to verify 123 that the suggested trend is systematically correct. 124

Thus, we are left with a large parameter space to investigate possible surface motion 125 during an Earth close approach of an asteroid with unknown surface mechanical properties. 126 Our approach is to consider a range of simple and well-controlled cases that certainly do 127 not cover all possibilities regarding Apophis' surface mechanical properties, but rather aim 128 at demonstrating whether, even in a simple and possibly favorable case for surface motion, 129 some resurfacing event can be expected to occur during the passage. For instance, instead 130 of considering a flat granular surface, we consider a sandpile consisting of a size distribution 131

## Yang Yu et al.

of spherical grains (Section 2.1) and vary the grain properties in order to include more or 132 less favorable cases for motion (from a fluid-like case to a case involving rough particles). 133 Slight disturbances may manifest as very-small-scale avalanches in which grain connections 134 readjust slightly, for example. The forces acting on the sandpile are obtained by measuring 135 all 'external' perturbations during the encounter, including body spin magnitude and orien136 tation changes, for cases in which the global shape remains nearly fixed, and again assuming 137 simple and favorable configurations of the asteroid. Indirectly, the encounter may also lead 138 to internal reconfigurations of the asteroid, which in turn produce seismic vibrations that 139 could propagate to the surface and affect the regolith material. These secondary modifica140 tions are not modeled here, although it may be possible in future work to account for this 141 by shaking the surface in a prescribed manner. In any case, for this particular encounter, we 142 demonstrate (Section 3.1) that any global reconfiguration will likely be small to negligible 143 in magnitude. 144

In the following, we first present, in Section 2, the numerical method used to perform our 145 investigation, including the initial conditions of the sandpile adopted to investigate surface 146 motion, the representation of the encounter, and the mechanical environment. Results are 147 described in Section 3, including potential reshaping of the asteroid, tidal disturbances for 148 Apophis' encounter in 2029, which is a function of the sandpile properties, spin orientation 149 changes, and the dependency of the location of the sandpile on the asteroid to the outcome 150 of the encounter. We also show the responses of the sandpiles for artificially close approaches 151 (4 0 and 2 0 Earth radii) to demonstrate that our method does predict significant alteration . . 152 of the sandpiles when this is certainly expected to happen. The investigation is discussed in 153 Section 4 and conclusions are presented in Section 5. 154

Yang Yu et al.

## 2 NUMERICAL METHOD 155

We use pkdgrav , a parallel N -body gravity tree code (Stadel 2001) adapted for particle 156 collisions (Richardson et al. 2000; 2009; 2011). Originally collisions in pkdgrav were treated 157 as idealized single-point-of-contact impacts between rigid spheres. A soft-sphere option was 158 added recently (Schwartz et al. 2012); with this new functionality, particle contacts last 159 many timesteps, with reaction forces dependent on the degree of overlap (a proxy for surface 160 deformation) and contact history-this is appropriate for dense and/or near-static granular 161 systems with multiple persistent contacts per particle. The code uses a 2nd-order leapfrog in162 tegrator, with accelerations due to gravity and contact forces recomputed each step. Various 163 types of user-definable confining walls are available that can be combined to provide complex 164 boundary conditions for the simulations. The code also includes an optional variable gravity 165 field based on a user-specified set of rules. 166

The spring/dash-pot model used in pkdgrav 's soft-sphere implementation is described 167 fully in Schwartz et al. (2012). Briefly, a (spherical) particle overlapping with a neighbor 168 or confining wall feels a reaction force in the normal and tangential directions determined 169 by spring constants ( k n , k t ), with optional damping and effects that impose static, rolling, 170 and/or twisting friction. The damping parameters ( C n , C t ) are related to the conventional 171 normal and tangential coefficients of restitution used in hard-sphere implementations, ε n 172 and ε t . The static, rolling, and twisting friction components are parameterized by dimen173 sionless coefficients µ s , µ r , and µ t , respectively. Plausible values for these parameters are 174 obtained through comparison with laboratory experiments (also see Section 2.1). Careful 175 consideration of the soft-sphere parameters is needed to ensure internal consistency, partic176

## Yang Yu et al.

ularly with the choice of k n , k t , and timestep-a separate code is provided to assist the user 177 with configuring these parameters correctly. The numerical approach has been validated 178 through comparison with laboratory experiments; e.g., Schwartz et al. (2012) demonstrated 179 that pkdgrav correctly reproduces experiments of granular flow through cylindrical hoppers, 180 specifically the flow rate as a function of aperture size, and found the material properties 181 of the grains also affect the flow rate. Also simulated successfully with the soft-sphere code 182 in pkdgrav were laboratory impact experiments into sintered glass beads (Schwartz et al. 183 184

2013), and regolith, in support of asteroid sampling mechanism design (Schwartz et al. 2014). We use a two-stage approach to model the effect of a tidal encounter on asteroid surface 185 particles (regolith). First, a rigid (non-deformable) object approximating the size, shape, 186 and rotation state of the asteroid is sent past the Earth on a fixed trajectory (in the present 187 study, the trajectory is that expected of (99942) Apophis-see Section 2.2; note that the 188 actual shape of Apophis is poorly known beyond an estimate of axis ratios, so we assume an 189 idealized ellipsoid for this study). All forces acting at a target point designated on the object 190 surface are recorded (Section 2.2). Then, a second simulation is performed in the local frame 191 of the target point, allowing the recorded external forces to affect the motion of particles 192 arranged on the surface (in the present study we consider equilibrated sandpiles). This two193 stage approach is necessary due to the large difference in scale between the asteroid as a 194 whole and the tiny regolith particles whose reactive motion we are attempting to observe. 195 We approximate the regolith, which in reality likely consists of a mixture of powders and 196 gravel (Clark et al. 2002), by a size distribution of solid spheres. We mimic the properties 197 of different materials by adjusting the soft-sphere parameters (Section 2.1). The soft-sphere 198 approach permits simulation of the behavior of granular materials in the near-static regime, 199

Yang Yu et al.

appropriate for the present case in which the regolith particles generally remain stationary 200 for hours and suffer a rapid disturbance only during the moments of closest approach to the 201 planet. In particular, the model permits a detailed look at the responses of local individual 202 particles even when the tidal effects are too weak to cause any macroscopic surface or shape 203 changes; this gives insight into the limit of tidal resurfacing effects. 204

## 2.1 Sandpile Initial Conditions 205

In order to assess the effect of tidal encounters on a small surface feature, we carried out 206 numerical simulations in a local frame consisting of a flat horizontal surface (i.e., the local 207 plane tangential to the asteroid surface at the target point) with a 'sandpile' resting on 208 top. The sandpile consists of N = 1683 simulated spherical particles with radii drawn from 209 a power-law distribution of mean 1 46 cm and . ± 0 29 cm width truncated (minimum and . 210 maximum particle radii 1 17 . and 1 75 . cm, respectively). Between the sandpile and the 211 floor is a rigid monolayer of particles drawn from the same distribution and laid down in 212 a close-packed configuration in the shape of a flat disk. This rigid particle layer provides 213 a rough surface for the sandpile to reduce spreading (Fig. 1). Three different sandpiles 214 were constructed using the material parameters described below. Particles were dropped 215 though a funnel from a low height onto the rough surface to build up the sandpile. These 216 sandpiles were allowed to settle until the system kinetic energy dropped to near zero. This 217 approach eliminates any bias that might arise from simply arranging the spheres by hand in 218 a geometrical way; the result should better represent a natural sandpile. 219

## Yang Yu et al.

For this study, we compared three different sets of soft-sphere parameters for the sandpile 221 particles (Table 1). Our goal was to define a set of parameters that spans a plausible range of 222 material properties given that the actual mechanical properties of asteroid surface material 223 are poorly constrained. In the specific case of Apophis, very little is known beyond the 224 spectral type, Sq (Binzel et al. 2009). There are no measurements of thermal properties that 225 might give an indication of the presence or absence of regolith on Apophis. Consequently, 226 we chose three sets of material parameters that span a broad range of material properties. 227 The first set, denoted 'smooth' in the table, consists of idealized frictionless spheres with a 228 small amount of dissipation (5%, chosen to match the glass beads case). This is about as 229 close to the fluid case that a sandpile can achieve while still exhibiting shear strength arising 230 from the discrete nature of the particles (and the confining pressure of surface gravity). It 231 is assumed this set will respond most readily to tidal effects due to the absence of friction 232 between the particles and the non-uniform size distribution. The second set, 'glass beads,' 233 is modeled after actual glass beads being used in a set of laboratory experiments to calibrate 234 numerical simulations of granular avalanches (Richardson et al. 2012). In this case ε t was 235 measured directly, which informed our choice for C t (Schwartz et al. 2012), and µ s and µ r 236 were inferred from matching the simulations to the experiments. The glass beads provide 237 an intermediate case between the near-fluid smooth spheres and the third parameter set, 238 denoted 'gravel.' 239

240

## Table 1.

The gravel parameters were arrived at by carrying out simple avalanche experiments 241 using roughly equal-size rocks collected from a streambed. In these experiments, the rocks 242

Yang Yu et al.

(without sharp edges) were released near the top of a wooden board held at a 45 ◦ incline. 243 The dispersal pattern on the stone floor was measured, including the distance from the 244 end of the board to approximately half of the material, the furthest distance traveled by 245 a rock along the direction of the board, and the maximum angle of the dispersal pattern 246 relative to the end of the board (Fig. 2a-c). A series of numerical simulations was then 247 performed to reproduce the typical behavior by varying the soft-sphere parameters (Fig. 248 2d-f). The values used in Table 1 for 'gravel' were found in this way. The µ s and µ r values 249 are quite large, reflecting the fact that the actual particles being modeled were not spheres. 250 A correspondingly smaller timestep is needed to adequately sample the resulting forces on 251 the particles in the simulations. The value of ε n was measured by calculating the average 252 first rebound height of sample gravel pieces that were released from a certain height; ε t was 253 not measured but since the rocks were rough it was decided to simply set ε t = ε n . This 254 exercise was not meant to be an exhaustive or precise study; rather, we sought simply to 255 find representative soft-sphere parameters that can account plausibly for the irregularities 256 in the particle shapes. In any case, we will find that such rough particles are difficult to 257 displace using tidal forces in the parameter range explored here, so they provide a suitable 258 upper limit to the effects. Similarly, we do not consider cohesion in this numerical study, 259 which would further resist particle displacement due to tidal forces. 260

261

Yang Yu et al.

## 2.2 Representation of Mechanical Environment 262

A new adaptation of pkdgrav developed in this work allows for the simulation of sandpiles 263 located on an asteroid surface, based on the motion equations of granular material in the 264 noninertial frame fixed to the chosen spot. Detailed mechanics involved in the local motion of 265 the sandpile during a tidal encounter were considered and represented in the code, including 266 the contact forces (reaction, damping, and friction between particles and/or walls) and the 267 inertial forces derived from an analysis of the transport motion of a flyby simulation. Here 268 we provide a thorough derivation of the relevant external force expressions. 269

270

## Figure 3.

Figure 3 illustrates four frames used in the derivation, the space inertial frame (SPC), 271 mass center translating frame (MCT), body fixed frame (BDY) and local frame (LOC). Eq. 272 (1) gives the connection between the spatial position and the local position of an arbitrary 273 particle in the sandpile, for which R is the vector from SPC's origin to the particle, R C is 274 the vector from SPC to MCT/BDY, l is the vector from MCT/BDY to LOC, and r is the 275 vector from LOC to the particle: 276

$$R = R _ { C } + 1 + r.$$

Equations (2)-(3) are derived by calculating the first- and second-order time derivatives 277 of Eq. (1), which denote the connections of velocity and acceleration between the spatial and 278 local representations, respectively. ω indicates the angular velocity vector of the asteroid 279 in SPC or MCT. The operator d d t denotes the time derivative with respect to the inertial 280

Yang Yu et al.

frame SPC, while ˜ d d t denotes the time derivative with respect to the body-fixed frame BDY 281 or LOC: 282

$$\frac { d } { d t } \mathbf R = \frac { d } { d t } \mathbf R _ { C } + \omega \times ( 1 + \mathbf r ) + \frac { \tilde { d } } { d t } \mathbf r,$$

$$\frac { d ^ { 2 } } { d t ^ { 2 } } \mathbf R = \frac { d ^ { 2 } } { d t ^ { 2 } } \mathbf R _ { C } + \omega \times [ \omega \times ( 1 + \mathbf r ) ] + \frac { \tilde { d } } { d t } \omega \times ( 1 + \mathbf r ) + 2 \omega \times \frac { \tilde { d } } { d t } \mathbf r + \frac { \tilde { d } ^ { 2 } } { d t ^ { 2 } } \mathbf r.$$

The dynamical equation for an arbitrary particle in the inertial frame SPC can be written 283 as Eq. (4), which shows that the forces a particle in the sandpile feels can be grouped into 284 three categories: F A denotes the local gravity from the asteroid; F P denotes the attraction 285 from the planet; and F C represents all the contact forces coming from other particles and 286 walls: 287

$$\frac { d ^ { 2 } } { d t ^ { 2 } } R = \mathbf F _ { A } + \mathbf F _ { P } + \mathbf F _ { C }.$$

We get the corresponding dynamical equation in the local frame LOC (Eq. (5)) by sub288 stituting Eq. (3) into Eq. (4), 289

$$\frac { \tilde { d } ^ { 2 } } { d t ^ { 2 } } \mathbf r = \mathbf F _ { C } + \mathbf F _ { A } + \mathbf F _ { P } - \frac { d ^ { 2 } } { d t ^ { 2 } } \mathbf R _ { C } - \omega \times [ \omega \times ( 1 + \mathbf r ) ] - \frac { \tilde { d } } { d t } \omega \times ( 1 + \mathbf r ) - 2 \omega \times \frac { \tilde { d } } { d t } \mathbf r. \quad ( 5 )$$

The terms on the right-hand side of Eq. (5) represent the mechanical environment of a 290 particle in the local sandpile during a tidal encounter, and can be interpreted as follows. 291 F C is the reaction force from the walls and surrounding particles, which is directly exported 292 by pkdgrav . F A is approximated by a constant value at the LOC's origin that acts as 293

## Yang Yu et al.

the uniform gravity felt throughout the whole sandpile, since the latter's size is negligible 294 compared with the asteroid dimensions. The difference between the planetary attraction 295 and the translational transport inertia force, F P -d 2 d t 2 R C , is the tidal force that plays the 296 primary role of moving the surface materials. The term -ω × [ ω × ( l + )] is the centrifugal r 297 force due to the rotation of the asteroid, -˜ d d t ω × ( l + ) is the librational transport inertia r 298 force (LTIF) that results when the rotational state changes abruptly during the encounter, 299 and -2 ω × ˜ d d t r is the Coriolis effect that only acts when particle movement occurs. 300

301

## Figure 4.

In the numerical scheme, we derived the component dynamical equations by projecting 302 Eq. (5) to LOC and treating each term separately. The local gravity F A and LOC origin l are 303 constant vectors, which were initialized in the simulations as parameters. The complications 304 arising from the asteroid motion (tidal force, angular velocity, and acceleration) were solved 305 within a simulation in advance, using a rigid rubble pile as the asteroid model and exporting 306 the physical quantities required in each step. To specify Apophis's encounter with Earth in 307 2029, we computed the planet trajectory by tracing back from the perigee conditions, 5 6 . 308 R E distance and 8 4 km/s encounter speed (Giorgini et al. 2008), where R . E indicates the 309 average Earth radius (6371 km). As Fig. 4 illustrates, a tri-axial ellipsoid model with axis 310 ratio 1.4:1.0:0.8 was employed to approximate the overall shape of Apophis (Scheeres et al. 311 2005), corresponding to an equivalent radius of 162 5 m (M¨ller et al. 2013). . u A rubble pile 312 bounded by the tri-axial model was constructed for the flyby simulation. The rubble pile 313 model consisted of N = 2855 equal-size spherical particles in a close-packed configuration. 314

The total mass was set to 4

×

10

10

3

315

.

kg (bulk density 2 4 g

/

cm ) with an initial rotation around

Yang Yu et al.

the maximum principle axis of inertia of period 30 4 h (Tholen et al. 2013). . Three markers 316 (non-coplanar with the origin) were chosen from the surface particles to track the variation 317 of the rubble pile's attitude in SPC. A complete list of its motion states during the encounter 318 was exported to a data file, which was accessed by the subsequent simulations with local 319 sandpiles. 320

## 3 RESULTS 321

## 3.1 Evaluation of Reshaping Effects 322

Richardson et al. (1998) showed the complicated behavior of a rotating rubble pile due to a 323 tidal encounter. Generally, a variation in rotational state is usually induced because of the 324 coupling effect between librational motion and orbital motion. Especially for elongated bod325 ies like Apophis, the terrestrial torques during its flyby may force a strong alteration of the 326 rotational state in a short period (Scheeres et al. 2005). The modifications in structure of the 327 rubble pile are also notable. The tidal encounter outcomes show a strong dependence on the 328 progenitor elongation, the perigee distance, and the time spent within the Roche sphere. The 329 outcomes also depend on the spin orientation of the progenitor to the extent that retrograde 330 rotators always suffer less catastrophic consequences than prograde ones (Richardson et al. 331 1998). However, the parameter space considered in that study (Richardson et al. 1998) did 332 not include the Apophis scenario in 2029, since the focus was on the conditions for full-scale 333 distortion or disruption. For this study, we first redid the flyby simulations (Section 2.2) 334 using a true (non-rigid) rubble pile model with the Apophis encounter parameters of 5 6 . 335

Yang Yu et al.

R E and 8 4 km/s. . We quantified the tidal reshaping effect by measuring the maximum net 336 change of a normalized shape factor defined by Eq. (6) (Hu et al. 2004), in which I 1 , I 2 , I 3 337 are the principal inertia moments, with I 1 /lessorequalslant I 2 /lessorequalslant I 3 : 338

$$\sigma \equiv \frac { I _ { 2 } - I _ { 1 } } { I _ { 3 } - I _ { 1 } }.$$

The shape factor σ ∈ [0 , 1] roughly describes the mass distribution of the body (0=oblate, 339 1=prolate). The initial rubble-pile model (equilibrated, prior to the encounter) had σ ∼ 0 73, . 340 and its relative change, defined by ∆ σ/σ , was measured for several simulations parameterized 341 by bulk density and constituent material type (Table 2). The corresponding Roche limits 342 are provided in the table, with the values estimated for the case of a circular orbit. The 343 reshaping effects show consistent results for all three types of materials, varying between a 344 noticeable change in shape (magnitude 10 -2 ) and a negligible change in shape (magnitude 345 10 ), 5 with the sharpest transition occurring for a critical bulk density ∼ 0 4 . g / cm 3 with 346 corresponding Roche limit of 5 6 R . E . This agrees with the analysis of Apophis' disruption 347 limit using a continuum theory for a body made up of solid constituents (Holsapple et al. 348 2006). Note we find that even a bulk density as low as 0 1 g . / cm 3 did not dislodge any 349 particles enough for them to end up in a new geometrical arrangement, due to the very short 350 duration of the tidal encounter. 351

352

Figure 5.

353

Table 2.

Since idealized equal-size spheres can arrange into crystal-like structures that may arti354 ficially enhance the shear strength of the body (Tanga et al. 2009; Walsh et al. 2012), we 355

## Yang Yu et al.

carried out a second suite of simulations with rubble piles made up of a bimodal distribution 356 of spheres (Fig. 5). The model consisted of 181 big particles of radius 19 3 m and 1569 small . 357 particles of radius 9 7 m, bounded by the same tri-axial ellipsoid as shown in Fig. 4, and . 358 arranged randomly. As before, the bulk density and material type were varied, with notable 359 differences in outcome compared to the equal-size-sphere cases (Table 2). Particles using the 360 'smooth' parameter set are not able to hold the overall shape of Apophis (the 1.4:1.0:0.8 361 tri-axial ellipsoid); instead the rubble pile collapses and approaches a near-spherical shape. 362 For the other two parameter sets, however, the overall shape of the rubble pile is main363 tained. Figure 6 shows the relative change of shape factor σ as a function of bulk density 364 for both types of simulations. As shown, the results from the unimodal and bimodal rubble 365 piles present similar trends: before reaching the Roche limit density, the reshaping remains 366 small (the magnitudes are 10 -5 -10 -4 with a certain amount of stochasticity); after that, 367 the reshaping effect sharply increases to around 10 -1 . The bimodal rubble piles show less 368 resistance to the tidal disturbance, since they have less-well-organized crystalline structures 369 compared to the unimodal cases; they therefore exhibit larger shape modifications (Walsh 370 et al. 2012). 371

372

## Figure 6.

Regardless, catastrophic events were not detected for either type of rubble pile until the 373 bulk density was as small as 0 1 . g / cm , 3 for all material parameter sets. Therefore, the 374 reshaping effects on Apophis in 2029 should be negligibly small for bulk densities in the 375 likely range (2-3 g / cm ). 3 Even so, minor internal reconfigurations resulting from the tidal 376 encounter may produce seismic waves that could propagate and affect the configuration of 377

## Yang Yu et al.

surface regolith. Since we are only looking at localized areas on the surface, isolated from 378 the rest of the body, vibrations emanating from other regions, in or on the asteroid, are not 379 evaluated in the current work, although again we expect these to be small or non-existent for 380 the specific case of the Apophis 2029 encounter. Only the 'external' forces acting directly 381 on the surface particles in the considered localized region are taken into account. We will 382 look into the effects of seismic activity, which may stem from other regions of the asteroid, 383 on an actual high-resolution rubble pile in future work. In this study, we focus on external 384 forces, outside of the rubble pile itself; the rigid rubble pile model is employed as a reasonable 385 simplification for the purpose of measuring the external forces on a surface sandpile. 386

## 3.2 Global Change of Mechanical Environment 387

The right-hand side of Eq. (5) provides a description of the constituents of the mechanical 388 environment experienced by a sandpile particle, including local gravity, tidal force, centrifugal 389 force, LTIF, and the Coriolis effect. Note the Coriolis force does not play a role before particle 390 motion begins, so it can be ignored when examining the causes of avalanches/collapses. The 391 centrifugal force and LTIF show weak dependence on local position r , so these terms can be 392 simplified by substituting l + r ≈ l , since the dimension of the sandpile is much smaller than 393 that of the asteroid. This leads to an approximation of the resultant environmental force 394 (Eq. (7)), which provides a uniform expression of the field force felt throughout the sandpile: 395

$$\mathbf F _ { E } = \mathbf F _ { A } + \mathbf F _ { P } - \frac { \mathbf d ^ { 2 } } { \mathbf d t ^ { 2 } } \mathbf R _ { C } - \omega \times ( \omega \times 1 ) - \frac { \tilde { \mathbf d } } { \mathbf d t } \omega \times 1.$$

The variation of F E shows a common pattern for different encounter trajectories and 396

Yang Yu et al.

different locations on the asteroid: it stays nearly invariable for hours as the planet is ap397 proaching/departing and shows a rapid single perturbation around the rendezvous. This 398 feature enables us to evaluate the mechanical environment by measuring the short-term 399 change of F E , and moreover, to make a connection between these external stimuli with the 400 responses of sandpiles in the simulations. Spherical coordinates ( F, θ, α ) are used to rep401 resent the resultant force (Fig. 7), in which F is the magnitude and θ and α denote the 402 effective gravity slope angle and deflection angle, respectively. 403

404

## Figure 7.

The effects of the environment force must be considered in the context of the sandpile 405 modeling. A rough approximation of the catastrophic slope angle (Eq. (8)) can be derived 406 from the conventional theory of the angle of repose for conical piles (Brown et al. 1966), 407 namely that avalanches will occur when the effective gravity slope angle θ exceeds θ C : 408

$$\theta _ { C } = \tan ^ { - 1 } ( \mu _ { S } ) - \theta _ { S }.$$

Here θ S denotes the resting angle of the (conical) sandpile, assumed to be less than the 409 repose angle tan -1 ( µ S ), and θ C denotes the critical slope angle, which, if exceeded by the 410 effective gravity slope angle, results in structural failure of the sandpile, i.e., an avalanche. 411 Essentially, the avalanches of the sandpile depend primarily on the instantaneous change of 412 the slope angle θ of the resultant force and should be only weakly related to the magnitude 413 F and deflection angle α . The global distribution of changes of θ were examined for the 414 duration of the encounter (see below). It would be difficult to make an exhaustive search 415 over all asteroid encounter orientations due to the coupling effects between the asteroid's 416

## Yang Yu et al.

rotation and the orbital motion. Instead, we chose 12 representative trajectories along the 417 symmetry axes of the Apophis model (at perigee) for study, which serves as a framework for 418 understanding the influence of encounter orientation and the strength of tidal disturbance at 419 different locations. (The technique used to match the orientation of the asteroid at perigee is 420 to first run a simulation in reverse starting at perigee with the required orientation, thereby 421 providing the correct asteroid state for a starting position far away.) Figure 8 shows the 422 12 trajectories along three mutually perpendicular planes, including both the prograde and 423 retrograde cases. This choice is based on the symmetry of the dynamical system, and all 12 424 trajectories have a speed of 8 4 km/s at a perigee distance of 5 6 R . . E . 425

426

## Figure 8.

Figure 8 presents these representative trajectories in order, each corresponding to a spe427 cial possible orientation of the encounter. The effective gravity slope angle changes through428 out the surface of the tri-axial ellipsoid model were recorded for each trajectory, with partic429 ular attention to the location and time of the maximum change. We found that these maxi430 mum changes concentrate at several minutes around perigee. Figure 9 shows the maximum 431 values of slope angle change along these trajectories, indicating that the largest perturbation 432 on slope angle is less than 2 . ◦ We verified that the achievable range of the effective gravity 433 slope angle during the encounter (Section 2.2) is within the safe limit predicted by Eq. (8) 434 for all three sandpiles (Section 2.1); that is, a slope change below 2 ◦ is not enough to trigger 435 any massive avalanches throughout the sandpiles. 436

## Yang Yu et al.

Figure 9 also shows a general dependence of the tidal perturbation on the orientation of 438 the encounter trajectory, namely that the three step levels seen in the figure correspond to 439 three directions of the relative planet motion at perigee. Trajectories 1-4 correspond to the 440 direction of the body long axis, which leads to relatively strong tidal effects; trajectories 5-8 441 correspond to the direction of the intermediate axis, which leads to moderate effects; while 442 trajectories 9-12 correspond to the direction of the short axis, which leads to relatively weak 443 tidal effects. Figure 10 illustrates the distribution of effective gravity slope angle change at 444 perigee for the three trajectory sets, with A denoting trajectories 1-4, B denoting trajectories 445 5-8, and C denoting trajectories 9-12. We confirm that the four trajectories in each set 446 present visually the same distribution around perigee, therefore we chose the patterns due 447 to trajectories 1, 5, and 9 for demonstration of sets A, B, and C, and generalize these three 448 patterns as representative. Several points can be inferred from Fig. 10. First, the direction 449 to the planet at perigee largely determines the global distribution of tidal perturbation, that 450 is, the strongest effects tend to occur along the long axis and the weakest effects along the 451 short axis. Second, at perigee, the largest slope change occurs near areas surrounding the 452 pole for the most favorable orientation (set A), while the largest slope change at the pole 453 occurs several minutes before or after perigee (see Fig. 11). These maximum slope changes 454 are about equal in magnitude (the change is only slightly smaller at the pole compared to 455 the area immediately surrounding it, but not enough to make a difference to the avalanches), 456 so for simplicity we just use the poles themselves as our testing points. Third, the duration 457 of strong tidal effects depends on the eccentricity of the encounter trajectory. We checked 458 the trajectories shown in Fig. 8 and found the duration of strong perturbation, defined for 459 illustration as the period when the force magnitude stays above 90% of the peak value, is tens 460

## Yang Yu et al.

of minutes long; correspondingly, the responses of the surface material are also transitory 461 and weak. 462

463

## Figure 10.

The three poles p , p x y and p z (see Figure 8) were chosen as the test locations for the 464 sample sandpiles since the effective gravity slope angle is near zero at these locations, which 465 enables the sandpiles to hold their initial shapes. Figure 11 illustrates the time variation of 466 slope angle at the poles during the encounter. The results derived from all 12 trajectories in 467 Fig. 8 are shown and to the same scale. The most perturbed pole is p y on the medium axis, 468 which gains the largest slope angle change for most cases, except in trajectories { 3,4,9,10 } . 469 It is notable that trajectories in the same plane always share similar variation at the three 470 poles, such as { 1,2,5,6 } , { 3,4,9,10 } , and { 7,8,11,12 } . The curves at p x and p y show similar 471 doublet shapes since those poles are both located in the rotational plane. The variation at 472 p z is more complicated due to the fact that the spin pole is not perpendicular to the orbit 473 plane. 474

475

## Figure 11.

Figure 11 can be used to estimate the changes in the mechanical environment at these 476 pole locations. Since the slope angle change is the primary mechanism to drive avalanches 477 on the sandpiles, the results derived from different trajectories serve as a framework to 478 determine the magnitude of tidal perturbation at the candidate locations and locations in 479 between, which turns out to be quite small for the Apophis encounter (less than about 1 ◦ in 480 slope for these cases). 481

Yang Yu et al.

## 3.3 Tidal Disturbances in 2029 482

As analyzed above, catastrophic avalanches solely due to external forces acting directly on 483 the surface particles may never occur on Apophis during the 2029 encounter since the tidal 484 perturbation will be very weak, however it might still have the potential to trigger some local 485 tiny landslides of the surface materials. Generally, the regolith experiences more activity 486 than the constituents deep in the asteroid due to the dynamical environment on the surface, 487 including the microgravity, maximum tidal and centrifugal acceleration, and smaller damping 488 forces from the surroundings due to the smaller confining pressure, therefore we conclude 489 that resurfacing due to external perturbations should occur before wholesale reshaping during 490 the encounter (discussed in Section 3.1). In the following sections we detail our numerical 491 examination of the local sandpiles for the predicted Apophis encounter scenario to estimate 492 the limit and magnitude of the material responses in 2029. 493

- One problem that must be confronted in the numerical simulations is that the soft-sphere 494 sandpiles exhibit slow outward spreading due to the accumulation of velocity noise, and this 495 slow spreading may eventually lead to a collapse of the whole sandpile. Unfortunately, our 496 integration time required is long (hours) for a granular system, thus the numerical noise has 497 to be well limited in our method. Two techniques are used in this study. First, as described 498 in Section 2.1, a rough pallet made up of closely packed spheres in a disk fixed to the ground 499 is used to reduce spreading of the bottom particles. Second, we introduce a critical speed 500 in the code, below which all motions are considered to be noise and are forced to zero. We 501 found a speed threshold ∼ 10 -10 m/s by launching hours-long simulations of a static sandpile 502 for different material properties and finding the corresponding minimum value of the critical 503

## Yang Yu et al.

speed that limits the numerical spreading. 504

Simulations of sandpiles without a flyby were carried out first, serving as a reference for 505 subsequent flyby simulations. Sandpiles in this situation were confirmed to stay equilibrated 506 for a long period, which suggests the avalanches in our numerical experiments (if any) would 507 be attributed entirely to the effects of the tidal encounter. We performed local simulations 508 with the sandpiles generated in Section 2.1 to consider the response of different materials 509 to the Apophis encounter. The local frames at the three poles p , x p y and p z (Fig. 8) 510 were used to position the sample sandpiles, and the flyby simulation data derived from the 511 12 fiducial trajectories were employed as the source of external perturbations for the local 512 simulations. There are 108 possible combinations in total for different materials, locations, 513 and trajectory orientations, covering a wide range of these undetermined parameters. In all 514 the simulations, the start distance of Apophis from the center of the planet was set to 18 515 R E to ensure the sandpiles were equilibrated fully before it approached the perigee. For our 516 study, we concentrated on the connections between the sandpile's responses and significant 517 variables, e.g., spin orientation and sandpile locations, by which we will get a better sense of 518 the surface processes due to a tidal encounter. The sandpiles located at pole p y for trajectory 519 2 (Fig. 8) were examined in detail, which is the case suffering the strongest tidal perturbation 520 (Fig. 11). Figure 12 illustrates the disturbances on sandpiles of three different materials for 521 this configuration. Each panel includes a snapshot of the sandpile after the tidal encounter 522 with the displaced particles highlighted and a diagram showing the time variations of the 523 total kinetic energy with planetary distance. Intensive events can be identified by noting 524 the peaks of the kinetic energy curves, and the particles involved in these events are marked 525 with different colors to show the correspondence between the disturbed site and occurrence 526

527

## Yang Yu et al.

time.

528

## Figure 12.

The scale of the disturbances proves to be tiny even for the 'strongest' case in Fig. 529 12. As illustrated, the displaced particles are few in number and mostly distributed on the 530 surface of the sandpile. We found that the maximum displacement of a given particle is 531 less than 0 8 times its radius; that is, these disturbances only resulted from the collapse of . 532 some local weak structures and the small chain reaction among the surrounding particles. 533 Moreover, the 'gravel' sandpile proves to be unaffected by the tidal encounter, because 534 the static friction is quite large (Table 1) and all particles are locked in a stable configu535 ration, in which case the total kinetic energy of the sandpile remains very small ( ∼ 10 -14 536 J). The 'glass beads' sandpile suffered small but concentrated disturbances that involved 537 very few particles. Accordingly, the displacements of these particles are relatively large. The 538 'smooth' sandpile presents near-fuild properties with many surface particles experiencing 539 small-amplitude sloshing. The motion eventually damps out and the displacements of the 540 involved particles are tiny (smaller than that of 'glass beads'). 541

Based on the detailed analysis of this representative scenario, we performed simulations 542 at other poles and at other orientations, constructing a database to reveal any connections 543 between disturbances and these parameters. However, the results show little association 544 between these two: for the 'gravel' sandpile, no disturbances were detected for any location 545 and any orientation due to the large µ S ; for the 'glass beads' sandpile, some small-scale 546 avalanches can always be triggered, while the occurrence (site and time) of these avalanches 547 seems to be widely distributed and independent of location and orientation; and for the 548

## Yang Yu et al.

'smooth' sandpile, surface particles can feel the tidal perturbation and show slight sloshing, 549 but no visible avalanches eventually resulted, because the initial low slope angle in this 550 case imposes a relatively stable structure that is always able to recover from the external 551 perturbations. 552

## 3.4 Tide-induced Avalanches at Closer Approaches 553

To illustrate the effect of stronger perturbations on our model sandpiles, we carried out a 554 few simulations of closer approaches, specifically at 4 0 and 2 0 R . . E , for the same encounter 555 speed of 8 4 . km/s. These scenarios, though unlikely to occur in 2029 based on current 556 understanding of Apophis' orbit, are presented here to illustrate the magnitudes of tidal 557 resurfacing effects over a wider range of perturbation strengths. 558

Flyby simulations of rubble piles were first carried out to quantify the reshaping effects, 559 using the monodisperse and bidisperse particle rubble-pile models (see Section 3.1). The 560 bulk densities for both models were set to ∼ 2 4 . g/cc. Table 3 presents the relative net 561 changes of shape factor ∆ σ/σ for these rubble piles of different materials at different perigee 562 distances. 563

564

## Table 3.

Table 3 shows results consistent with Section 3.1, namely that the bimodal rubble pile 565 shows greater fluidity and larger shape changes during the encounter. To be more specific, 566 the magnitude of reshaping effects at perigee distance 4 0 R . E (still larger than the Roche 567 limit of 3 22 R . E for a fluid body) remains small ( ∼ 10 -4 ), and that at perigee distance 2 0 . 568 R E (smaller than the Roche limit) achieves a significant level of 10 -3 -10 -1 . It is notable 569

## Yang Yu et al.

that the rubble piles did not experience any irreversible distortion even for an approach 570 distance as close as 2 0 R . E , because the duration of strong tidal effects is relatively short 571 (see Section 3.1). Anyway, in this section we still adopt the assumption of rigidity to measure 572 the quantities required by the local simulations, which is acceptable since these scenarios are 573 fictitious and only designed to exhibit some massive resurfacing effects. 574

For the same reason, in this section we do not present systematic local simulations for 575 different sandpile locations and orientations as done for the real encounter scenario; instead, 576 we simply place the sample sandpiles (see Fig. 1) at the pole of the body long axis p x , and 577 choose trajectory 1 (see Fig. 8) as the encounter trajectory to give rise to the maximum tidal 578 effects. Figures 13 and 14 illustrate the responses of the three sandpiles for perigee distance 579 4 0 R . E and 2 0 R . E , respectively. Each panel includes snapshots of the sandpile before and 580 after the tidal encounter, and the overall shape of the sandpile is traced with a white border 581 for emphasis. Diagrams showing the time variation of the total kinetic energy with planetary 582 distance are also included, in which the intensive avalanches can be identified by noting the 583 peaks of the kinetic energy curves (note the different vertical scales). 584

585

Figure 13.

586

Figure 14.

As illustrated in Fig. 13, the shape changes of the three sandpiles are still small (but 587 visible) for the encounter at 4 0 . R E perigee, and involve many more particles than the 588 encounter at 5 6 R . E (Fig. 12). And accordingly, the magnitudes of the total kinetic energy 589 increase consistently for the three sandpiles of different materials, except for the gravel case 590 at 4 0 R . E , where it appears the frictional 'lock' established when the pile was first created is 591

Yang Yu et al.

disturbed enough to cause a stronger distortion than might otherwise be expected, causing 592 a sharp spike in the kinetic energy plot. Essentially, the 'gravel' sandpile experienced a 593 more significant collapse than the 'glass bead' sandpile. As stated in Section 2.1, the three 594 sandpiles were constructed in a manner that allowed for some inherent randomness, thus 595 their energy state and strength may differ to some degree, and appearently the 'gravel' 596 sandpile was closer to its failure limit than the other two. This adds an extra element of 597 stochasticity to the results, an aspect to be explored in future work. 598

Figure 14 shows the results from the encounter at 2 0 R . E perigee, for which the encounter 599 trajectory partly entered the Roche sphere (3 22 R . E ). In this case, the shapes of the three 600 sandpiles are highly distorted during the encounter, with the involved particles slumped 601 towards the direction where the planet is receding. The corresponding changes in total 602 kinetic energy become extremely large when the massive avalanches occur, especially for the 603 'smooth' particles. 604

- The results of this section suggest that a 4 0 R . E encounter may alter the regolith on 605 Apophis' surface slightly, and a 2 0 R . E encounter may produce a strong resurfacing effect (of 606 course, we would expect considerable global distortion as well, if the asteroid is a rubble pile). 607 The shape changes of the sandpiles depend on the orientation of the encounter trajectory; 608 609
- i.e., particles can be dragged away along the direction in which the planet recedes.

## 4 DISCUSSION 610

- The argument about whether a terrestrial encounter can reset the regolith of NEAs has been 611 discussed for a while (e.g., Binzel et al. 2010, Nesvorn´ et al. 2010). y Important evidence is 612

Yang Yu et al.

provided by measurements of the spectral properties, which suggest that asteroids with orbits 613 conducive to tidal encounters with planets show the freshest surfaces, while quantitative 614 evaluation of this mechanism is still rare and rough. The two-stage approach presented in 615 this paper enables the most detailed simulations to date of regolith migration due to tidal 616 effects, which we have applied to make a numerical prediction for the surface effects during 617 the 2029 approach of (99942) Apophis. 618

We confirm that the shape modification of Apophis due to this encounter will likely 619 be negligibly small based on the results of systematic simulations parameterized by bulk 620 density, internal structure, and material properties. The analysis of the global mechanical 621 environment over the surface of Apophis during the 2029 encounter reveals that the tidal 622 perturbation is even too small to result in any large-scale avalanches of the regolith materials, 623 based on a plausible range of material parameters and provided that external forces acting 624 directly on the surface dominate any surface effects due to seismic activity emanating from 625 other regions of the body. Future work will explore these second-order surface effects and 626 their relevance to this encounter, and to cases of tidal resurfacing of small bodies in general. 627 It is notable that our approach is capable of capturing slight changes in the regolith (modeled 628 as sandpiles); thus, through numerical simulation, we find that this weak tidal encounter does 629 trigger some local disturbances for appropriate material properties. 630

These possible disturbances, though local and tiny, are essentially related to the nature 631 of the sandpile, which is actually a bunch of rocks and powders in a jammed state. The 632 sandpile is formed by a competition between the constituents to fall under gravity and 633 eventually reach an equilibrium with mutual support. It can hold a stable shape under a 634 constant environment force at low sandpile densities (compared with close packing), which 635

## Yang Yu et al.

is primarily due to interior collective structures, called bridges, that are formed at the same 636 time as the sandpile and are distributed non-uniformly throughout the sandpile (Mehta 637 2007). The soft-sphere method provides a detailed approximation to a real sandpile in terms 638 of the granular shape and contact mechanics, which can well reproduce the bridges that 639 dominate the structure of the sandpile. Accordingly, we can describe the nuanced responses 640 due to the changes in ( F, θ, α )-see Section 3.2 and Fig. 7. Changing the slope angle θ 641 is the most efficient way to break the equilibrium of substructures in the sandpile, while 642 changing the magnitude F and deflection angle α may also play a role in the causes of 643 avalanches. Changes in F can readjust the interparticle overlaps and cause collapse of some 644 bridges, which is a principal mechanism of compaction by filling the voids during small 645 structural modifications (Mehta et al. 2004). This dependence was recently demonstrated in 646 experiments of measuring the angle of repose under reduced gravity (Kleinhans et al. 2011). 647 Similarly, the weak links in bridges of the sandpile may also be disturbed and broken during 648 the sweeping of α in the plane. In addition, Mehta et al. (1991) pointed out that even slight 649 vibration caused by the environment force may result in collapse of long-lived bridges, which 650 is also responsible for the sandpile landslides during the encounter. 651

Although the magnitudes of the avalanches are different for the three materials we tested, 652 they are nonetheless all very small. We note the occurrence of these local avalanches shows 653 little dependence on the orientation/location of the sandpile, because all the perturbations 654 are small in magnitude and the sandpile behavior, these small-scale avalanches, actually 655 depends on the presence of weak bridges inside the sandpile, which in turn could be sensitive 656 to the changes of the environmental force's magnitude, deflection angle, and slope angle in 657 a somewhat random way. 658

Yang Yu et al.

In any case, this numerical study shows that tidal resurfacing may not be particularly 659 effective at moderate encounter distances. We predict that overall resurfacing of Apophis, 660 regolith will not occur if the only source of disturbance is external perturbations. Mini661 landslides on the surface may still be observed by a visiting spacecraft with sufficiently 662 sensitive monitoring equipment. This provides a great motivation for in situ exploration of 663 Apophis in 2029 (Michel et al. 2012). 664

## 5 CONCLUSIONS 665

This paper provides a numerically derived prediction for the surface effects on (99942) 666 Apophis during its 2029 approach, which is likely to be one of the most significant as667 teroid encounter events in the near future. A two-stage scheme was developed based on the 668 soft-sphere code implementation in pkdgrav to mimic both a rubble pile's (rigid and flexi669 ble) responses to a planetary flyby and a sandpile's responses to all forms of perturbations 670 induced by the encounter. The flyby simulations with the rubble pile indicate that reshaping 671 effects due to the tidal force on Apophis in 2029 will be negligibly small for bulk densities 672 in the expected range (2-3 g / cm ). 3 The resultant environmental force felt by the sandpile 673 on the asteroid surface was approximated with a uniform analytical expression, which led 674 to an estimate of the changes in the global mechanical environment. Three typical patterns 675 of perturbation were presented based on the asteroid body and spin orientation at perigee, 676 showing a general dependence of the magnitude of tidal perturbation on the orientation of 677 the trajectory. Twelve fiducial trajectories were used to calculate the magnitude of the tidal 678 perturbation at three poles of the tri-axial ellipsoid model, indicating that the strongest 679

## Yang Yu et al.

tidal perturbation appears where the local slope is originally steep and that the duration of 680 the strong perturbation is short compared with the whole process. The tidal perturbation 681 on surface materials is confirmed to be quite weak for the 2029 encounter, therefore large682 scale avalanches may never occur. However, we showed this weak perturbation does trigger 683 some local tiny landslides on the sample sandpiles, though the involved particles are few 684 in number and are distributed on the surface of the sandpile. These small-scale avalanches 685 result from the breaking of weak substructures by slight external perturbations, therefore 686 the occurrence of these local avalanches is widely distributed and presents little dependence 687 on the encounter parameters. The simulations of closer approaches show that an encounter 688 at 2 0 R . E is capable of triggering some massive avalanches of the sandpiles, i.e., to alter 689 the regolith on Apophis' surface significantly (although the entire body would also undergo 690 significant shape change in this case). 691

Further research will be performed to generalize our work over a wide range of possible 692 asteroid-planet encounter conditions. And as stated above, we will also investigate whether 693 even slight internal perturbations in the asteroid during tidal encounters may contribute to 694 noticeable surface motions. 695

Yang Yu et al.

## Acknowledgments 696

This material is based on work supported by the US National Aeronautics and Space Ad697 ministration under Grant Nos. NNX08AM39G, NNX10AQ01G, and NNX12AG29G issued 698 through the Office of Space Science and by the National Science Foundation under Grant 699 No. AST1009579. Some simulations in this work were performed on the 'yorp' cluster 700 administered by the Center for Theory and Computation in the Department of Astron701 omy at the University of Maryland. Patrick Michel acknowledges support from the French 702 space agency CNES. Ray tracing for Figs.1, 2, 4, 5 and 12 was performed using POV-Ray 703 (http://www.povray.org/). 704

Yang Yu et al.

## References 705

Benz, W., Asphaug, E., 1999. Catastrophic disruptions revisited. Icarus 142, 5-20. 706

707

708

709

Binzel, R.P., Rivkin, A.S., Thomas, C.A., Vernazza, P., Burbine, T.H., DeMeo, F.E., Bus, S.J., Tokunaga, A.T., Birlan, M., 2009. Spectral properties and composition of potentially hazardous Asteroid (99942) Apophis. Icarus 200, 480-485.

710

711

712

713

714

715

716

717

718

719

720

721

722

723

Binzel, R.P., Morbidelli, A., Merouane, S., DeMeo, F.E., Birlan, M., Vernazza, P., Thomas, C.A., Rivkin, A.S., Bus, S.J., Tokunaga, A.T., 2010. Earth encounters as the origin of fresh surfaces on near-Earth asteroids. Nature 463, 331-334.

Brown, R.L., Richards, J.C., 1966. Principles of Powder Mechanics. Oxford: Pergamon.

Clark, B.E., Hapke, B., Pieters, C., Britt, D., 2002. Asteroid space weathering and regolith evolution. In: Bottke Jr., W.F., Cellino, A., Paolicchi, P., Binzel, R.P. (Eds.), Asteroids III. Univ. of Arizona Press, Tucson, pp. 485-500.

Delbo', M., Cellino, A., Tedesco, E.F., 2007. Albedo and size determination of potentially hazardous asteroids: (99942) Apophis. Icarus 188, 266-269.

DeMeo, F.E., Binzel, R.P., Lockhart, M., 2013. Mars encounters cause fresh surfaces on some near-Earth asteroids. Icarus 227, 112-122.

Fujiwara, A. Kawaguchi, J., Yeomans, D.K., Abe, M., Mukai, T., Okada, T., Saito, J., Yano, H., Yoshikawa, M., Scheeres, D.J., Barnouin-Jha, O., Cheng, A.F., Demura, H., Gaskell, R.W., Hirata, N., Ikeda, H., Kominato, T., Miyamoto, H., Nakamura,

724

725

726

727

728

729

730

731

732

733

734

735

736

737

738

739

740

- 741

742

743

Yang Yu et al.

- A.M., Nakamura, R., Sasaki, S., Uesugi, K., 2006. The rubble-pile asteroid Itokawa as observed by Hayabusa. Science 312, 1330-1334.

Giorgini, J.D., Benner, L.A.M., Ostro, S.J., Nolan, M.C., Busch, M.W., 2008. Predicting the Earth encounters of (99942) Apophis. Icarus 193, 1-19.

Holsapple, A., Michel, P., 2006. Tidal disruptions: A continuum theory for solid bodies. Icarus 183, 331-348.

Hu, W., Scheeres, D.J., 2004. Numerical determination of stability regions for orbital motion in uniformly rotating second degree and order gravity fields. Planetary and Space Science 52, 685-692.

Jutzi, M., Michel, P., Benz, W., Richardson, D.C., 2010. Fragment properties at the catastrophic disruption threshold: The effect of the parent body's internal structure. Icarus 207, 54-65.

Kleinhans, M.G., Markies, H., de Vet, S.J., in 't Veld, A.C., Postema, F.N., 2011. Static and dynamic angles of repose in loose granular materials under reduced gravity. Journal of Geophysical Research 116, E11004.

Nesvorn´, D., Bottke Jr., W.F., Vokrouhlick´, D., Chapman, C.R., Rafkin, S., 2010. y y Do planetary encounters reset surfaces of near Earth asteroids? Icarus 209, 510-519.

Mehta, A., 2007. Granular Physics. Cambridge University Press.

Mehta, A., Barker, G.C., 1991. Vibrated powders: A microscopic approach. Physics Review Letter 67, 394-397.

744

745

746

747

748

749

750

751

752

753

754

755

756

757

758

759

760

761

762

763

Yang Yu et al.

Mehta, A., Barker, G.C., Luck, J.M., 2004. Cooperativity in sandpiles: Statistics of bridge geometries. Journal of Statistical Mechanics: Theory and Experiment. P10014.

Michel, P., Benz, W., Tanga, P., Richardson, D.C., 2001. Collisions and gravitational reaccumulation: Forming asteroid families and satellites. Science 294, 1696-1700.

Michel, P., Prado, J.Y., Barucci, M.A., Bousquet, P., Chiavassa, F., Groussin, O., Herique, A., Hestroffer, D., Hinglais, E., Lopes, L., Martin, T., Mimoun, D., R´me, H., Thuillot, e W., 2012. Probing the interior of asteroid Apophis: A unique opportunity in 2029. Special Session 7 of the XXVIII IAU General Assembly, 29-31 August, 2012, Beijing (China), Abstract Book.

Richardson, D.C., Bottke, Jr., W.F., Love, S.G., 1998. Tidal distortion and disruption of Earth-crossing asteroids. Icarus 134, 47-76.

Richardson, D.C., Quinn, T., Stadel, J., Lake, G., 2000. Direct large-scale N -body simulations of planetesimal dynamics. Icarus 143, 45-59.

Richardson, D.C., Leinhardt, Z.M., Melosh, H.J., Bottke Jr., W.F., Asphaug, E., 2002. Gravitational aggregates: Evidence and evolution. In: Bottke Jr., W.F., Cellino, A., Paolicchi, P., Binzel, R.P. (Eds.), Asteroids III. Univ. of Arizona Press, Tucson, pp. 501-515.

Richardson, D.C., Michel, P., Walsh, K.J., Flynn, K.W., 2009. Numerical simulations of asteroids modeled as gravitational aggregates. Planetary and Space Science 57, 183192.

764

765

766

767

768

769

770

771

772

773

774

775

776

777

778

779

780

781

782

783

784

Yang Yu et al.

Richardson, D.C., Walsh, K.J., Murdoch, N., Michel, P., 2011. Numerical simulations of granular dynamics: I. Hard-sphere discrete element method and tests. Icarus 212, 427-437.

Richardson, D.C., Blum, J., Weinhart, T., Schwartz, S.R., Michel, P., Walsh, K.J., 2012. Numerical simulations of landslides calibrated against laboratory experiments for application to asteroid surface processes. American Astronomical Society, DPS meeting 44, 105.06.

Scheeres, D.J., Benner, L.A.M., Ostro, S.J., Rossi, A., Marzari, F., Washabaugh, P., 2005. Abrupt alteration of Asteroid 2004 MN4's spin state during its 2029 Earth flyby. Icarus 178, 281-283.

Schwartz, S.R., Richardson, D.C., Michel, P., 2012. An implementation of the soft-sphere discrete element method in a high-performance parallel gravity tree-code. Granular Matter 14, 363-380.

Schwartz, S.R., Michel, P., Richardson, D.C., 2013. Numerically simulating impact disruptions of cohesive glass bead agglomerates using the soft-sphere discrete element method. Icarus 226, 67-76.

Schwartz, S.R., Michel, P., Richardson, D.C., Yano, H., 2014. Low-speed impact simulations into regolith in support of asteroid sampling mechanism design I.: Comparison with 1-g Experiments. Planet. Space Sci., in press.

Stadel, J., 2001. Cosmological N -body simulations and their analysis. Thesis, University of Washington. 141 pp.

785

786

787

788

789

790

791

792

793

794

795

Yang Yu et al.

Tanga, P., Comito, C., Paolicchi, P., Hestroffer, D., Cellino, A., DellOro, A., Richardson, D. C., Walsh, K. J., Delbo', M., 2009. Rubble-pile reshaping reproduces overall asteroid shapes. The Astrophysical Journal 706, 197-202.

Tholen, D.J., Micheli, M., Elliott, G.T., 2013. Improved astrometry of (99942) Apophis. Acta Astronautica 90, 56-71.

Veverka, J., and 32 co-authors, 2000. NEAR at Eros: Imaging and spectral results. Science 289, 2088-2097.

Walsh, K. J., Richardson, D.C., Michel, P., 2012. Spin-up of rubble-pile asteroids: Disruption, satellite formation, and equilibrium shapes. Icarus 220, 514-529.

Wlodarczyk, I., 2013. The potentially dangerous asteroid (99942) Apophis. Monthly Notices of the Royal Astronomical Society 434(4), 3055-3060.

Table 1: Soft-sphere parameter sets used to represent three different material properties of sandpiles in the simulations. Here µ s is the coefficient of static friction (= tan -1 φ , where φ is the material friction angle), µ r is the coefficient of rolling friction, ε n is the normal coefficient of restitution (1 = elastic, 0 = plastic), ε t is the tangential coefficient of restitution (1 = no sliding friction), k n is the normal spring constant (kg / s 2 ), and k t is the tangential spring constant (also kg / s 2 ).

| Parameter   |   Gravel |   Glass Beads |   Smooth |
|-------------|----------|---------------|----------|
| µ s         |     1.31 |          0.43 |     0    |
| µ r         |     3    |          0.1  |     0    |
| ε n         |     0.55 |          0.95 |     0.95 |
| ε t         |     0.55 |          1    |     1    |
| k n         |    83.3  |         83.3  |    83.3  |
| k t         |    23.8  |         23.8  |    23.8  |

Table 2: Relative net change of shape factor ∆ σ/σ for the flexible rubble-pile model of monodisperse and bidispersed particles, parameterized by the bulk density and material property. The results of bidisperse smooth particles are omitted since the rubble pile in this case cannot hold the overall shape of a 1.4:1.0:0.8 tri-axial ellipsoid at Apophis' spin rate. The estimates of Apophis' Roche limit (second column) corresponding to the indicated bulk density are calculated assuming a fluid body in a circular orbit around Earth. For reference, the Apophis encounter distance in 2029 will be about 5.6 R E .

|   Density (g / cm 3 ) |   Roche Limit (R E ) | Dispersion       | Gravel                          | Glass Beads                     | Smooth            |
|-----------------------|----------------------|------------------|---------------------------------|---------------------------------|-------------------|
|                   2.4 |                 3.22 | Unimodal Bimodal | 6 . 62 × 10 - 5 8 . 36 × 10 - 4 | 7 . 76 × 10 - 5 6 . 00 × 10 - 4 | 3 . 26 × 10 - 5 - |
|                   2   |                 3.42 | Unimodal Bimodal | 5 . 19 × 10 - 5 9 . 67 × 10 - 4 | 8 . 84 × 10 - 5 6 . 06 × 10 - 4 | 8 . 67 × 10 - 5 - |
|                   1.5 |                 3.77 | Unimodal Bimodal | 2 . 75 × 10 - 5 4 . 73 × 10 - 4 | 9 . 54 × 10 - 5 3 . 51 × 10 - 4 | 3 . 60 × 10 - 5 - |
|                   1   |                 4.31 | Unimodal Bimodal | 1 . 88 × 10 - 5 9 . 07 × 10 - 4 | 4 . 95 × 10 - 5 5 . 57 × 10 - 4 | 4 . 68 × 10 - 5 - |
|                   0.5 |                 5.43 | Unimodal Bimodal | 1 . 99 × 10 - 5 2 . 82 × 10 - 3 | 1 . 45 × 10 - 5 2 . 31 × 10 - 4 | 2 . 16 × 10 - 5 - |
|                   0.1 |                 9.29 | Unimodal Bimodal | 2 . 37 × 10 - 2 6 . 54 × 10 - 2 | 3 . 12 × 10 - 3 1 . 31 × 10 - 1 | 2 . 47 × 10 - 2 - |

Table 3: Relative net change of shape factor ∆ σ/σ for the rubble piles of monodisperse and bidisperse particles, parameterized by the perigee distance and material property. As before, the results of bidisperse smooth particles are omitted since they cannot maintain Apophis' tri-axial ellipsoidal shape.

|   Perigee (R E ) | Dispersion       | Gravel                          | Glass Beads                     | Smooth            |
|------------------|------------------|---------------------------------|---------------------------------|-------------------|
|                2 | Unimodal Bimodal | 1 . 45 × 10 - 3 4 . 48 × 10 - 2 | 1 . 42 × 10 - 3 7 . 95 × 10 - 2 | 5 . 62 × 10 - 1 - |
|                4 | Unimodal Bimodal | 1 . 32 × 10 - 4 9 . 37 × 10 - 4 | 7 . 18 × 10 - 4 8 . 09 × 10 - 4 | 5 . 88 × 10 - 4 - |

## Figure Captions 796

797

798

799

800

801

Figure 1: Snapshots of sandpiles constructed using three different materials, which are (a) 'gravel,' (b) 'glass beads,' and (c) 'smooth,' with corresponding soft-sphere parameters listed in Table 1. The same constituent spheres are used in each sandpile both for the free particles (green) and the rigid particles (white). The values of average slope and pile height after equilibrium is achieved are indicated in the snapshots.

802

803

804

805

806

807

808

809

810

811

812

813

814

815

816

Figure 2: Snapshots of the avalanche experiment and corresponding simulations. In the experiments, ∼ 340 similar-size gravel pieces were piled up on the slope each time and released all at once by removing the supporting board. Snapshots from both the experiment and simulation include frames from the beginning (a, d), middle (b, e), and end (c, f) of the avalanche event, respectively.

Figure 3: Sketch of the frames used for deriving the motion equations of the local sandpile, indicated with different colors. SPC (black) is an inertial frame with the origin and axes set by pkdgrav . MCT (blue) is a noninertial translating frame with the origin at the mass center of the asteroid and x , y , z -axes parallel to SPC. BDY (green) is a noninertial frame that is fixed to the asteroid and initially coincides with MCT. LOC (red) is a frame fixed to the asteroid surface, taking the origin to be where the sandpile is located and choosing the z -axis to be normal to the surface at that point.

Figure 4: Diagram of Asteroid (99942) Apophis'. The tri-axial ellipsoid (light gray) denotes the overall shape model, while the arranged particles (khaki) denote the rubblepile model used for numerical simulations. The highlighted particles (red) are markers

817

818

819

820

821

822

823

824

825

826

827

828

829

830

831

832

833

834

835

836

837

## Yang Yu et al.

used for determining the attitude of the asteroid in SPC.

Figure 5: Snapshot of the rubble-pile model for Apophis with bimodal particles in irregular packing. The tri-axial ellipsoid (light gray) denotes the overall shape model, and the arranged particles (khaki) show the rubble-pile model used for numerical simulations.

Figure 6: The relative change of shape factor σ as a function of bulk density during the encounter. The solid line (with triangular points) shows the results for the original rubble pile with equal-sized particles. The dashed line (with square points) shows the results for the bimodal rubble pile with irregular packing.

Figure 7: Sketch of the environmental force that the sandpile feels in LOC. The plane (gray) indicates the local tangential plane. The red arrow indicates the environmental force given by Eq. (7). Label F indicates its magnitude, and θ and α indicate the slope angle (elevation) and deflection angle (azimuth) of the effective gravity, respectively.

Figure 8: Sketch of the fiducial trajectories in different orientations. These trajectories (blue numbered lines) are placed in the perpendicular planes (yellow) of the Apophis model (gray ellipsoid), matching the three axial directions at perigee. The blue arrows show the prograde and retrograde directions. The red arrow denotes the angular velocity. p , p , p x y z indicate the poles along the three axes x , y , z , respectively.

Figure 9: Maximum values of slope angle change for the fiducial trajectories. The stars denote the maximum values, and the dashed lines denote the levels of tidal perturbation (the average of the 4 values in each case), which divide the 12 trajectories into 3 basic catagories according to the spin orientation at perigee. See Fig. 8 for the trajectory

838

839

840

841

842

843

844

845

846

847

848

849

850

851

852

853

854

855

orientations.

Figure 10: Global distribution of slope angle change at perigee for the three trajectory sets A, B and C, each including 4 trajectories of the same spin orientation at perigee (see text for an explanation of each set). A uniform colormap ranging up to 1 8 . ◦ is used for the three plots.

Figure 11: Time variation of the effective gravity slope angle at the poles p x (dashed lines), p y (solid lines), and p z (dotted lines) during the encounter. Each subgraph corresponds to one fiducial trajectory marked with its number on the upper left (refer to Fig. 8).

Figure 12: Snapshots of sandpiles after the tide-induced avalanches in 2029, with corresponding time variations of their total kinetic energy. The three panels include the results of sandpiles constructed using the three materials (a) 'gravel,' (b) 'glass beads,' and (c) 'smooth,' respectively. Particles that moved more than 0 1 of a parti-. cle radius during the avalanches are highlighted. The peaks in the total kinetic energy curves indicate occurrence of avalanches; the peaks are labeled, and corresponding particles that moved appreciably are marked in different colors (red, blue, orange) in the snapshots. The dotted line in each graph shows the evolution of the planet-asteroid distance during the encounter.

Figure 13: Snapshots of the sample sandpiles before and after the avalanches induced 856 by the encounter at 4 0 R . E perigee, with corresponding time variations of their total 857 kinetic energy. The white lines in the snapshots indicate the original and final shapes 858

- of the sandpiles. The peaks in the total kinetic energy curves indicate occurrence of 859 avalanches, and the dotted line shows the variation of the planet-asteroid distance 860 during the encounter. 861
- Figure 14: Same as Fig. 13, but for the encounter at 2 0 R . E perigee. 862

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-b5c1d38f-e479-41c2-97f9-c28cf6bb7068.jpg)

(a) Sandpile with 'gravel' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-11435f76-fc4b-4b0f-89ff-4c9650065bad.jpg)

(b) Sandpile with 'glass beads' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-54349a1b-a5fb-4607-914f-03bdc1f2ec29.jpg)

- (c) Sandpile with 'smooth' parameter set.

Figure 1:

## Yang Yu et al.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-312627dd-423f-4bb0-b4b3-d6af8eb796e9.jpg)

- (a) Experiment: Beginning
- (d) Simulation: Beginning
- (b) Experiment: Middle
- (e) Simulation: Middle
- (c) Experiment: End
- (f) Simulation: End

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-0fb800b1-e110-415e-82a9-824e08ade312.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-fbb8d30b-9a6e-4db6-a0d5-fa44589f3aa2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-b3e6a772-d971-40d6-8b4f-c967746e7be0.jpg)

Figure 2:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-54f5e744-31ec-4357-ad6f-dd37ce050f86.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-64994703-6f2d-4e22-9baa-901a42c72aa8.jpg)

Figure 3:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-f5277091-ccfb-4b44-b25a-c38188c6557e.jpg)

Figure 4:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-6bdab53b-78c3-451d-ba07-c4662a35065d.jpg)

Figure 5:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-f3ad0b31-146a-42f6-8c0e-629241e1dea2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-98f77d5e-0106-4182-83c0-fc2634902aee.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-ace1c19f-d574-44ff-9929-5e15a018fd1e.jpg)

- (a) Rubble piles with 'gravel' parameter set.
- (b) Rubble piles with 'glass beads' parameter set.
- (c) Rubble piles with 'smooth' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-2b4f7c99-8620-41c6-9794-a769769bd241.jpg)

Figure 6:

Figure 7:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-a185ce5b-fcaa-4366-a62a-ea846815a4f3.jpg)

Figure 8:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-f25bcb44-e300-479c-90fc-01d9d7552670.jpg)

Figure 9:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-586ac1a0-98d8-4224-873e-e2cb28f9ed8d.jpg)

Figure 10:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-c51b53dd-178a-46f2-8e2b-13cdaf7c2e41.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-baa485d2-44e0-4227-89c3-3b75acc801aa.jpg)

Time (h)

Figure 11:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-27a7497a-7129-418a-9abb-53b9c400cabc.jpg)

Time (hr)

- (a) Avalanches with 'gravel' parameter set.
- (b) Avalanches with 'glass beads' parameter set.
- (c) Avalanches with 'smooth' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-a95a6616-af1e-4a07-b73a-ad008a559555.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-0cebf4d0-f3e6-41ed-ab4d-f9b0e0d9c422.jpg)

Figure 12:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-c28f57d5-ce0a-4ddd-8c27-68a958e97c9a.jpg)

- (a) Shape changes with 'gravel' parameter set.
- (b) Shape changes with 'glass beads' parameter set.
- (c) Shape changes with 'smooth' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-213ffaf1-679c-41bf-8fb9-ad5fa4dcd847.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-2792cea7-72f0-47ca-9d9c-13db138412c1.jpg)

Figure 13:

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-34c566c0-8e45-4492-989d-64bbd8aa308c.jpg)

- (a) Shape changes with 'gravel' parameter set.
- (b) Shape changes with 'glass beads' parameter set.
- (c) Shape changes with 'smooth' parameter set.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-d9d7bbfa-d4f0-4bb7-b1c7-059558a1dcc3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0168v1-figure-24c4e390-c6e8-470f-8410-ff31c209768d.jpg)

Figure 14:
Yang Yu Derek C. Richardson Patrick Michel Stephen R. Schwartz Ronald-Louis Ballouz
Numerical predictions of surface effects during the 2029 close approach of asteroid 99942 Apophis
Asteroid (99942) Apophis' close approach in 2029 will be one of the most
significant small-body encounter events in the near future and offers a good
opportunity for in situ exploration to determine the asteroid's surface
properties and measure any tidal effects that might alter its regolith
configuration. Resurfacing mechanics has become a new focus for asteroid
researchers due to its important implications for interpreting surface
observations, including space weathering effects. This paper provides a
prediction for the tidal effects during the 2029 encounter, with an emphasis on
whether surface refreshing due to regolith movement will occur. The potential
shape modification of the object due to the tidal encounter is first confirmed
to be negligibly small with systematic simulations, thus only the external
perturbations are taken into account for this work (despite this, seismic
shaking induced by shifting blocks might still play a weak role and we will
look into this mechanism in future work). A two-stage approach is developed to
model the responses of asteroid surface particles (the regolith) based on the
soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A
full-body model of Apophis is sent past the Earth on the predicted trajectory
to generate the data of all forces acting at a target point on the surface. A
sandpile constructed in the local frame is then used to approximate the
regolith materials; all the forces the sandpile feels during the encounter are
imposed as external perturbations to mimic the regolith's behavior in the full
scenario. The local mechanical environment on the asteroid surface is
represented in detail, leading to an estimation of the change in global surface
environment due to the encounter. Typical patterns of perturbation are
presented that depend on the asteroid orientation and sense of rotation at
perigee.