Asteroid families are formed as the result of collisions. Large fragments are ejected with speeds of the order of the escape velocity from the parent body. After the family formation, the fragments' orbits evolve in the space of proper elements because of gravitational and non-gravitational perturbations, such as the Yarkovsky effect. Disentangling the contribution to the current orbital position of family members caused by the initial ejection velocity field and the subsequent orbital evolution is usually a difficult task. Among the more than 100 asteroid families currently known, some interact with secular resonances. Linear secular resonances occur when there is a commensurability between the precession frequency of the longitude of the pericenter (g) or of the longitude of node (s) of an asteroid and a planet, or a massive asteroid. The linear secular resonance most effective in increasing an asteroid eccentricity is the ${\nu}_6$, that corresponds to a commensurability between the precession frequency g of an asteroid and Saturn's $g_6$. Non-linear secular resonances involve commensurabilities of higher order, and can often be expressed as combinations of linear secular resonances. This is the case, for instance, of the $z_k=k(g-g_6)+(s-s_6)$ resonances. Asteroid families that are crossed by secular resonances are of particular interest in dynamical astronomy. First, they often provide a clear evidence of asteroid orbit evolution due to the Yarkovsky effect. Second, conserved quantities of secular dynamics can be used to set valuable constraints on the magnitude of the original ejection velocity field. Finally, by changing the value of inclination of family members nodal secular resonances with massive asteroids or dwarf planets can cause the i distribution to become more and more leptokurtic (i.e., more peaked and with larger tails than that of a Gaussian distribution).

The Nele group, formerly known as the Iannini family, is one of the youngest asteroid families in the main belt. Previously, it has been noted that the pericenter longitudes $\varpi$ and nodal longitudes $\Omega$ of its largest member asteroids are clustered at the present time, therefore suggesting that the collisional breakup of parent body must have happened recently. Here we verify this conclusion by detailed orbit-propagation of a synthetic Nele family and show that the current level of clustering of secular angles of the largest Nele family members requires an approximate age limit of $4.5$ Myr. Additionally, we make use of an updated and largely extended Nele membership to obtain, for the first time, an age estimate of this family using the Backward Integration Method (BIM). Convergence of the secular angles in a purely gravitational model and in a model including the non-gravitational forces caused by the Yarkovsky effect are both compatible with an age younger than $7$ Myr. More accurate determination of the Nele family age would require additional data about the spin state of its members.

Asteroid families are valuable source of information to many asteroid-related re- searches, assuming a reliable list of their members could be obtained. However, as the number of known asteroids increases fast it becomes more and more difficult to obtain robust list of members of an asteroid family. Here we are proposing a new approach to deal with the problem, based on the well known Hierarchical Clustering Method (HCM). An additional step in the whole procedure is introduced in order to reduce a so-called chaining effect. The main idea is to prevent chaining through an al- ready identified interloper. We show that in this way a number of potential interlopers among family members is significantly reduced. Moreover, we developed an automatic on-line based portal to apply this procedure, i.e to generate a list of family members as well as a list of potential interlopers. The Asteroid Families Portal (AFP) is freely available to all interested researchers.

The age of a young asteroid family can be determined by tracking the orbits of family members backward in time and showing that they converge at some time in the past. Here we consider the Veritas family. We find that the membership of the Veritas family increased enormously since the last detailed analysis of the family. Using backward integration, we confirm the convergence of nodal longitudes $\Omega$, and, for the first time, also obtain a simultaneous convergence of pericenter longitudes $\varpi$. The Veritas family is found to be $8.23^{+0.37}_{-0.31}$~Myr old. To obtain a tight convergence of $\Omega$ and $\varpi$, as expected from low ejection speeds of fragments, the Yarkovsky effect needs to be included in the modeling of the past orbital histories of Veritas family members. Using this method, we compute the Yarkovsky semi-major axis drift rates, ${\rm d}a/{\rm d}t$, for 274 member asteroids. The distribution of ${\rm d}a/{\rm d}t$ values is consistent with a population of C-type objects with low densities and low thermal conductivities. The accuracy of individual ${\rm d}a/{\rm d}t$ measurements is limited by the effect of close encounters of member asteroids to (1) Ceres and other massive asteroids, which cannot be evaluated with confidence.

V-type asteroids are a taxonomic class whose surface is associated to a basaltic composition. The only known source of V-type asteroids in the Main Asteroid Belt is (4) Vesta, that is located in the inner part of the belt. However, many V-type asteroids cannot be dynamically linked to Vesta., in particular, those asteroids located in the middle and outer parts of the Belt. Previous works have failed to find mechanisms to transport V-type asteroids from the inner to the middle and outer belt. In this work we propose a dynamical mechanism that could have acted on primordial asteroid families. We consider a model of the giant planets migration known as the jumping Jupiter model with five planets. Our study is focused on the period of 10 Myr that encompasses the instability phase of the giant planets. We show that, for different hypothetical Vesta-like paleo-families in the inner belt, the perturbations caused by the ice giant that is scattered into the asteroid belt before being ejected from the solar system, are able to scatter V-type asteroids to the middle and outer belt. Based on the orbital distribution of V-type candidates identified from the Sloan Digital Sky Survey and the VISTA Survey colours, we show that this mechanism is efficient enough provided that the hypothetical paleo-family originated from a 100 to 500 km crater excavated on the surface of (4) Vesta. This mechanism is able to explain the currently observed V-type asteroids in the middle and outer belt, with the exception of (1459) Magnya.

The Hoffmeister family is a C-type group located in the central main belt. Dynamically, it is important because of its interaction with the ${\nu}_{1C}$ nodal secular resonance with Ceres, that significantly increases the dispersion in inclination of family members at lower semi-major axis. As an effect, the distribution of inclination values of the Hoffmeister family at semi-major axis lower than its center is significantly leptokurtic, and this can be used to set constraints on the terminal ejection velocity field of the family at the time it was produced. By performing an analysis of the time behaviour of the kurtosis of the $v_W$ component of the ejection velocity field (${\gamma}_2(v_W)$), as obtained from Gauss' equations, for different fictitious Hoffmeister families with different values of the ejection velocity field, we were able to exclude that the Hoffmeister family should be older than 335 Myr. Constraints from the currently observed inclination distribution of the Hoffmeister family suggest that its terminal ejection velocity parameter $V_{EJ}$ should be lower than 25~m/s. Results of a Yarko-YORP Monte Carlo method to family dating, combined with other constraints from inclinations and ${\gamma}_2(v_W)$, indicate that the Hoffmeister family should be $220^{+60}_{-40}$~Myr old, with an ejection parameter $V_{EJ}= 20\pm5$~m/s.

One of the most accurate models currently used to represent the gravity field of irregular bodies is the polyhedral approach. In this model, the mass of the body is assumed to be homogeneous, which may not be true for a real object. The main goal of the present paper is to study the dynamical effects induced by three different internal structures (uniform, three- and four-layers) of asteroid (21) Lutetia, an object that recent results from space probe suggest being at least partially differentiated. The Mascon gravity approach used in the present work, consists of dividing each tetrahedron into eight parts to calculate the gravitational field around the asteroid. The zero-velocity curves show that the greatest displacement of the equilibrium points occurs in the position of the E4 point for the four-layers structure and the smallest one occurs in the position of the E3 point for the three-layers structure. Moreover, stability against impact shows that the planar limit gets slightly closer to the body with the four-layered structure. We then investigated the stability of orbital motion in the equatorial plane of (21) Lutetia and propose numerical stability criteria to map the region of stable motions. Layered structures could stabilize orbits that were unstable in the homogeneous model.

$v_W$ leptokurtic asteroid families are families for which the distribution of the normal component of the terminal ejection velocity field $v_W$ is characterized by a positive value of the ${\gamma}_2$ Pearson kurtosis, i.e., they have a distribution with a more concentrated peak and larger tails than the Gaussian one. Currently, eight families are known to have ${\gamma}_2(v_W) > 0.25$. Among these, three are highly inclined asteroid families, the Hansa, Barcelona, and Gallia families. As observed for the case of the Astrid family, the leptokurtic inclination distribution seems to be caused by the interaction of these families with node secular resonances. In particular, the Hansa and Gallia family are crossed by the $s-s_V$ resonance with Vesta, that significantly alters the inclination of some of their members. In this work we use the time evolution of ${\gamma}_2(v_W)$ for simulated families under the gravitational influence of all planets and the three most massive bodies in the main belt to assess the dynamical importance (or lack of) node secular resonances with Ceres, Vesta, and Pallas for the considered families, and to obtain independent constraints on the family ages. While secular resonances with massive bodies in the main belt do not significantly affect the dynamical evolution of the Barcelona family, they significantly increase the ${\gamma}_2(v_W)$ values of the simulated Hansa and Gallia families. Current values of the ${\gamma}_2(v_W)$ for the Gallia family are reached over the estimated family age only if secular resonances with Vesta are accounted for.

Among asteroid families, the Astrid family is peculiar because of its unusual inclination distribution. Objects at $a\simeq$~2.764 au are quite dispersed in this orbital element, giving the family a "crab-like" appearance. Recent works showed that this feature is caused by the interaction of the family with the $s-s_C$ nodal secular resonance with Ceres, that spreads the inclination of asteroids near its separatrix. As a consequence, the currently observed distribution of the $v_W$ component of terminal ejection velocities obtained from inverting Gauss equation is quite leptokurtic, since this parameter mostly depends on the asteroids inclination. The peculiar orbital configuration of the Astrid family can be used to set constraints on key parameters describing the strength of the Yarkovsky force, such as the bulk and surface density and the thermal conductivity of surface material. By simulating various fictitious families with different values of these parameters, and by demanding that the current value of the kurtosis of the distribution in $v_W$ be reached over the estimated lifetime of the family, we obtained that the thermal conductivity of Astrid family members should be $\simeq$ 0.001 W/m/K, and that the surface and bulk density should be higher than 1000 kg/m$^{3}$. Monte Carlo methods simulating Yarkovsky and stochastic YORP evolution of the Astrid family show its age to be $T$ = 140$\pm$30 Myr old, in good agreement with estimates from other groups. Its terminal ejection velocity parameter is in the range $V_{EJ}= 5^{+17}_{-5}$~m/s.] Values of $V_{EJ}$ larger than 25 m/s are excluded from constraints from the current inclination distribution.

The Karin cluster is a young asteroid family thought to have formed only $\simeq 5.75$~My ago. The young age can be demonstrated by numerically integrating the orbits of Karin cluster members backward in time and showing the convergence of the perihelion and nodal longitudes (as well as other orbital elements). Previous work has pointed out that the convergence is not ideal if the backward integration only accounts for the gravitational perturbations from the Solar System planets. It improves when the thermal radiation force known as the Yarkovsky effect it is accounted for. This argument can be used to estimate the spin obliquities of the Karin cluster members. Here we take advantage of the fast growing membership of the Karin cluster and show that the obliquity distribution of diameter $D\simeq 1-2$ km Karin asteroids is bimodal, as expected if the YORP effect acted to move obliquities toward the extreme values ($0^\circ$ or $180^\circ$). The measured magnitude of the effect is consistent with the standard YORP model. The surface thermal conductivity is inferred to be $0.07$-0.2 W m$^{-1}$ K$^{-1}$ (thermal inertia $\simeq 300-500$ J m$^{-2}$ K$^{-1}$s$^{-1/2}$). We find that the strength of the YORP effect is roughly $\simeq 0.7$ of the nominal strength obtained for a collection of random Gaussian spheroids. These results are consistent with a surface composed of rough, rocky regolith. The obliquity values predicted here for 480 members of the Karin cluster can be validated by the light-curve inversion method.

Asteroid families are groups of minor bodies produced by high-velocity collisions. After the initial dispersions of the parent bodies fragments, their orbits evolve because of several gravitational and non-gravitational effects,such as diffusion in mean-motion resonances, Yarkovsky and YORP effects, close encounters of collisions, etc. The subsequent dynamical evolution of asteroid family members may cause some of the original fragments to travel beyond the conventional limits of the asteroid family. Eventually, the whole family will dynamically disperse and no longer be recognizable. A natural question that may arise concerns the timescales for dispersion of large families. In particular, what is the oldest still recognizable family in the main belt? Are there any families that may date from the late stages of the Late Heavy Bombardment and that could provide clues on our understanding of the primitive Solar System? In this work, we investigate the dynamical stability of seven of the allegedly oldest families in the asteroid main belt. Our results show that none of the seven studied families has a nominally mean estimated age older than 2.7 Gyr, assuming standard values for the parameters describing the strength of the Yarkovsky force. Most "paleo-families" that formed between 2.7 and 3.8 Gyr would be characterized by a very shallow size-frequency distribution, and could be recognizable only if located in a dynamically less active region (such as that of the Koronis family). V-type asteroids in the central main belt could be compatible with a formation from a paleo-Eunomia family.

Ceres is the largest and most massive body in the asteroid main belt. Observational data from the Dawn spacecraft reveal the presence of at least two impact craters about 280~km in diameter on the Ceres surface, that could have expelled a significant number of fragments. Yet, standard techniques for identifying dynamical asteroid families have not detected any Ceres family. In this work, we argue that linear secular resonances with Ceres deplete the population of objects near Ceres. Also, because of the high escape velocity from Ceres, family members are expected to be very dispersed, with a considerable fraction of km-sized fragments that should be able to reach the pristine region of the main belt, the area between the 5J:-2A and 7J:-3A mean-motion resonances, where the observed number of asteroids is low. Rather than looking for possible Ceres family members near Ceres, here we propose to search in the pristine region. We identified 156 asteroids whose taxonomy, colors, albedo could be compatible with being fragments from Ceres. Remarkably, most of these objects have inclinations near that of Ceres itself.

An asteroid family forms as a result of a collision between an impactor and a parent body. The fragments with ejection speeds higher than the escape velocity from the parent body can escape its gravitational pull. The cloud of escaping debris can be identified by the proximity of orbits in proper element, or frequency, domains. Obtaining estimates of the original ejection speed can provide valuable constraints on the physical processes occurring during collision, and used to calibrate impact simulations. Unfortunately, proper elements of asteroids families are modified by gravitational and non-gravitational effects, such as resonant dynamics, encounters with massive bodies, and the Yarkovsky effect, such that information on the original ejection speeds is often lost, especially for older, more evolved families. It has been recently suggested that the distribution in proper inclination of the Koronis family may have not been significantly perturbed by local dynamics, and that information on the component of the ejection velocity that is perpendicular to the orbital plane ($v_W$), may still be available, at least in part. In this work we estimate the magnitude of the original ejection velocity speeds of Koronis members using the observed distribution in proper eccentricity and inclination, and accounting for the spread caused by dynamical effects. Our results show that i) the spread in the original ejection speeds is, to within a 15% error, inversely proportional to the fragment size, and ii) the minimum ejection velocity is of the order of 50 m/s, with larger values possible depending on the orbital configuration at the break-up.

Asteroid families form as a result of large-scale collisions among main belt asteroids. The orbital distribution of fragments after a family-forming impact could inform us about their ejection velocities. Unfortunately, however, orbits dynamically evolve by a number of effects, including the Yarkovsky drift, chaotic diffusion, and gravitational encounters with massive asteroids, such that it is difficult to infer the ejection velocities eons after each family's formation. Here we analyze the inclination distribution of asteroid families, because proper inclination can remain constant over long time intervals, and could help us to understand the distribution of the component of the ejection velocity that is perpendicular to the orbital plane ($v_{W}$). From modeling the initial breakup, we find that the distribution of $v_{W}$ of the fragments, which manage to escape the parent body's gravity, should be more peaked than a Gaussian distribution (i.e., be leptokurtic) even if the initial distribution was Gaussian. We surveyed known asteroid families for signs of a peaked distribution of $v_{W}$ using a statistical measure of the distribution peakedness or flatness known as kurtosis. We identified eight families whose $v_{W}$ distribution is significantly leptokurtic. These cases (e.g. the Koronis family) are located in dynamically quiet regions of the main belt, where, presumably, the initial distribution of $v_{W}$ was not modified by subsequent orbital evolution. We suggest that, in these cases, the inclination distribution can be used to obtain interesting information about the original ejection velocity field.

The Erigone family is a C-type group in the inner main belt. Its age has been estimated by several researchers to be less then 300 My, so it is a relatively young cluster. Yarko-YORP Monte Carlo methods to study the chronology of the Erigone family confirm results obtained by other groups. The Erigone family, however, is also characterized by its interaction with the $z_2$ secular resonance. While less than 15% of its members are currently in librating states of this resonance, the number of objects, members of the dynamical group, in resonant states is high enough to allow to use the study of dynamics inside the $z_2$ resonance to set constraints on the family age. Like the ${\nu}_{6}$ and $z_1$ secular resonances, the $z_2$ resonance is characterized by one stable equilibrium point at $\sigma = 180^{\circ}$ in the $z_2$ resonance plane $(\sigma, \frac{d\sigma}{dt})$, where $\sigma$ is the resonant angle of the $z_2$ resonance. Diffusion in this plane occurs on timescales of $\simeq 12$ My, which sets a lower limit on the Erigone family age. Finally, the minimum time needed to reach a steady-state population of $z_2$ librators is about 90 My, which allows to impose another, independent constraint on the group age.

The Euphrosyne asteroid family is uniquely situated at high inclination in the outer Main Belt, bisected by the nu_6 secular resonance. This large, low albedo family may thus be an important contributor to specific subpopulations of the near-Earth objects. We present simulations of the orbital evolution of Euphrosyne family members from the time of breakup to the present day, focusing on those members that move into near-Earth orbits. We find that family members typically evolve into a specific region of orbital element-space, with semimajor axes near ~3 AU, high inclinations, very large eccentricities, and Tisserand parameters similar to Jupiter family comets. Filtering all known NEOs with our derived orbital element limits, we find that the population of candidate objects is significantly lower in albedo than the overall NEO population, although many of our candidates are also darker than the Euphrosyne family, and may have properties more similar to comet nuclei. Followup characterization of these candidates will enable us to compare them to known family properties, and confirm which ones originated with the breakup of (31) Euphrosyne.

The Cybele region, located between the 2J:-1A and 5J:-3A mean-motion resonances, is adjacent and exterior to the asteroid main belt. An increasing density of three-body resonances makes the region between the Cybele and Hilda populations dynamically unstable, so that the Cybele zone could be considered the last outpost of an extended main belt. The presence of binary asteroids with large primaries and small secondaries suggested that asteroid families should be found in this region, but only relatively recently the first dynamical groups were identified in this area. Among these, the Sylvia group has been proposed to be one of the oldest families in the extended main belt. In this work we identify families in the Cybele region in the context of the local dynamics and non-gravitational forces such as the Yarkovsky and stochastic YORP effects. We confirm the detection of the new Helga group at $\simeq$3.65~AU, that could extend the outer boundary of the Cybele region up to the 5J:-3A mean-motion resonance. We obtain age estimates for the four families, Sylvia, Huberta, Ulla and Helga, currently detectable in the Cybele region, using Monte Carlo methods that include the effects of stochastic YORP and variability of the Solar luminosity. The Sylvia family should be $T = 1220 \pm 40$ Myr old, with a possible older secondary solution. Any collisional Cybele group formed prior to the late heavy bombardment would have been most likely completely dispersed in the jumping Jupiter scenario of planetary migration.

Asteroids formed in a dynamically quiescent disk but their orbits became gravitationally stirred enough by Jupiter to lead to high-speed collisions. As a result, many dozen large asteroids have been disrupted by impacts over the age of the Solar System, producing groups of fragments known as asteroid families. Here we explain how the asteroid families are identified, review their current inventory, and discuss how they can be used to get insights into long-term dynamics of main belt asteroids. Electronic tables of the membership for 122 notable families are reported on the Planetary Data System node.

V-type asteroids, characterized by two absorption bands at 1.0 and 2.0 $\mu m$, are usually thought to be portions of the crust of differentiated or partially differentiated bodies. Most V-type asteroids are found in the inner main belt and are thought to be current or past members of the Vesta dynamical family. Recently, several V-type photometric candidates have been identified in the central and outer main belt. While the dynamical evolution of V-type photometric candidates in the central main belt has been recently investigated, less attention has been given to the orbital evolution of basaltic material in the outer main belt as a whole. Here we identify known and new V-type photometric candidates in this region, and study their orbital evolution under the effect of gravitational and non-gravitational forces. A scenario in which a minimum of three local sources, possibly associated with the parent bodies of (349) Dembowska, (221) Eos, and (1459) Magnya, could in principle explain the current orbital distribution of V-type photometric candidates in the region.

V-type asteroids are associated with basaltic composition, and are supposed to be fragments of crust of differentiated objects. Most V-type asteroids in the main belt are found in the inner main belt, and are either current members of the Vesta dynamical family (Vestoids), or past members that drifted away. However, several V-type photometric candidates have been recently identified in the central and outer main belt. The origin of this large population of V-type objects is not well understood. Since it seems unlikely that Vestoids crossing the 3J:-1A mean-motion resonance with Jupiter could account for the whole population of V-type asteroids in the central and outer main belt, origin from local sources, such as the parent bodies of the Eunomia, and of the Merxia and Agnia asteroid families, has been proposed as an alternative mechanism. In this work we investigated the dynamical evolution of the V-type photometric candidates in the central main belt, under the effect of gravitational and non-gravitational forces. Our results show that dynamical evolution from the parent bodies of the Eunomia and Merxia/Agnia families on timescales of 2 Byr or more could be responsible for the current orbital location of most of the low-inclined V-type asteroids.

The asteroid (10) Hygiea is the fourth largest asteroid of the Main Belt, by volume and mass, and it is the largest member of its own family. Previous works investigated the long-term effects of close encounters with (10) Hygiea of asteroids in the orbital region of the family, and analyzed the taxonomical and dynamical properties of members of this family. In this paper we apply the high-quality SDSS-MOC4 taxonomic scheme of DeMeo and Carry (2013) to members of the Hygiea family core and halo, we obtain an estimate of the minimum time and number of encounter necessary to obtain a $3\sigma$ (or 99.7%) compatible frequency distribution function of changes in proper $a$ caused by close encounters with (10) Hygiea, we study the behavior of asteroids near secular resonance configurations, in the presence and absence of the Yarkovsky force, and obtain a first estimate of the age of the family based on orbital diffusion by the Yarkovsky and YORP effects with two methods. The Hygiea family is at least 2 Byr old, with an estimated age of $T = 3200^{+380}_{-120}$ Myr and a relatively large initial ejection velocity field, according to the approach of Vokrouhlick\'{y} et al. (2006a, b). Surprisingly, we found that the family age can be shortened by $\simeq$ 25% if the dynamical mobility caused by close encounters with (10) Hygiea is also accounted for, which opens interesting new research lines for the dynamical evolution of families associated with massive bodies. In our taxonomical analysis of the Hygiea asteroid family, we also identified a new V-type candidate: the asteroid (177904) (2005 SV5). If confirmed, this could be the fourth V-type object ever to be identified in the outer main belt.

Previous works have identified families halos by an analysis in proper elements domains, or by using Sloan Digital Sky Survey-Moving Object Catalog data, fourth release (SDSS-MOC4) multi-band photometry to infer the asteroid taxonomy, or by a combination of the two methods. The limited number of asteroids for which geometric albedo was known until recently discouraged in the past the extensive use of this additional parameter, which is however of great importance in identifying an asteroid taxonomy. The new availability of geometric albedo data from the Wide-field Infrared Survey Explorer (WISE) mission for about 100,000 asteroids significantly increased the sample of objects for which such information, with some errors, is now known. In this work we proposed a new method to identify families halos in a multi-domain space composed by proper elements, SDSS-MOC4 (a*,i-z) colors, and WISE geometric albedo for the whole main belt (and the Hungaria and Cybele orbital regions). Assuming that most families were created by the breakup of an undifferentiated parent body, they are expected to be homogeneous in colors and albedo. The new method is quite effective in determining objects belonging to a family halo, with low percentages of likely interlopers, and results that are quite consistent in term of taxonomy and geometric albedo of the halo members.

Among the largest objects in the main belt, asteroid 4 Vesta is unique in showing a basaltic crust. It is also the biggest member of the Vesta family, which is supposed to originate from a large cratering event about 1 Gyr ago (Marzari et al. 1996). Most of the members of the Vesta family for which a spectral classification is available show a V-type spectra. Before the discovery of 1459 Magnya (Lazzaro et al. 2000) and of several V-type NEA (Xu 1995), all the known V-type asteroids were members of the Vesta family. Recently two V-type asteroids, 809 Lundia and 956 Elisa, (Florczak et al. 2002) have been discovered well outside the limits of the family, near the Flora family. We currently know 22 V-type asteroids outside the family, in the inner asteroid belt. In this work we investigate the possibility that these objects are former family members that migrated to their current positions via the interplay of Yarkovsky effect and nonlinear secular resonances. The main dynamical feature of 956 Elisa and 809 Lundia is that they are currently inside the 2(g-g6)+s-s6 (z2 by Milani and Knezevic, 1993) secular resonance. Our investigations show that members of the Vesta dynamical family may drift in three-body and weak secular resonances until they are captured in the strong z2 secular resonance. Only asteroids with diameters larger than 16 km can remain in one of the three-body or secular resonances long enough to reach the region of the z2 resonance. This two-step mechanism of capture into the z2 resonance could explain: i) the current resonant orbits of 956 Elisa and 809 Lundia, ii) why their size is significantly larger than that of the typical member of the Vesta family, and iii) provide a lower limit on the Vesta family age.

Among the many new irregular satellites that have been discovered in the last five years, at least six are in the so-called Kozai resonance. Due to solar perturbations, the argument of pericenter of a satellite usually precesses from 0 to 360 degrees. However, at inclinations higher than 39.3 degrees and lower than 140.7 degrees a new kind of behavior occurs for which the argument of pericenter oscillates around +/-90 degrees. In this work we will concentrate on the orbital history of the saturnian satellite S/2000 S5 Kiviuq, one of the satellites currently known to be in such resonance Kiviuq's orbit is very close to the separatrix of the Kozai resonance. Due to perturbations from the other jovian planets, it is expected that orbits near the Kozai separatrix may show significant chaotic behavior. This is important because chaotic diffusion may transfer orbits from libration to circulation, and vice versa. To identify chaotic orbits we used two well-known methods: the Frequency Analysis Method (Laskar 1990) and Maximum Lyapunov Exponents (Benettin et al. 1980). Our results show that the Kozai resonance is crossed by a web of secondary resonances, whose arguments involve combinations of the argument of pericenter, the argument of the Great Inequality, longitude of the node, and other terms related to the secular frequencies g5, g6, and s6. Many test orbits whose precession period is close to the period of the Great Inequality (883 yrs), or some of its harmonics, are trapped by these secondary resonances, and show significant chaotic behavior. Planetary migration, by moving the locations of these secondary resonances, may have depleted an original population of Kozai resonators.