The minor planets on orbits that are dynamically stable in Neptune's 1:1 resonance on Gyr timescales were likely emplaced by Neptune's outward migration. We explore the intrinsic libration amplitude, eccentricity, and inclination distribution of Neptune's stable Trojans, using the detections and survey efficiency of the Outer Solar System Origins Survey (OSSOS) and Pan-STARRS1. We find that the libration amplitude of the stable Neptunian Trojan population can be well modeled as a Rayleigh distribution with a libration amplitude width $\sigma_{A_phi}$ of 15$^\circ$. When taken as a whole, the Neptune Trojan population can be acceptably modeled with a Rayleigh eccentricity distribution of width $\sigma_e$ of 0.045 and a typical sin(i) x gaussian inclination distribution with a width $\sigma_i$ of 14 +/- 2 degrees. However, these distributions are only marginally acceptable. This is likely because, even after accounting for survey detection biases, the known large Hr < 8 and small Hr >= 8 Neptune Trojans appear to have markedly different eccentricities and inclinations. We propose that like the classical Kuiper belt, the stable intrinsic Neptunian Trojan population have dynamically `hot' and dynamically `cold' components to its eccentricity/inclination distribution, with $\sigma_{e-cold}$ ~ 0.02 / $\sigma_{i-cold}$ ~ 6$^\circ$ and $\sigma_{e-hot}$~ 0.05 / $\sigma_{i-hot}$ ~ 18$^\circ$. In this scenario, the `cold' L4 Neptunian Trojan population lacks the Hr >= 8 member and has 13 +11/-6 `cold' Trojans with Hr < 8. On the other hand, the `hot' L4 Neptunian Trojan population has 136 +57/-48 Trojans with Hr < 10 -- a population 2.4 times greater than that of the L4 Jovian Trojans in the same luminosity range.

The high-inclination population in the 4:7 mean motion resonance (MMR) with Neptune has also substantial eccentricities ($e\gtrsim0.1$), with more inclined objects tending to occupy more eccentric orbits. For this high-order resonance, there are two different resonant modes. The principal one is the eccentricity-type mode, and we find that libration is permissible for orbits with $e\ge e_c^0$, where the critical eccentricity $e_c^0$ increases as a function of increasing inclination $i$. Correspondingly, we introduce a limiting curve $e_c^0(i)$, which puts constraints on the $(e, i)$ distribution of possible 4:7 resonators. We then perform numerical simulations on the sweep-up capture and long-term stability of the 4:7 MMR, and the results show that the simulated resonators are well-constrained by this theoretical limiting curve. The other 4:7 resonant mode is the mixed-$(e, i)$-type, and we show that stable resonators should exist at $i\gtrsim20^{\circ}$. We predict that the intrinsic number of these mixed-$(e, i)$-type resonators may provide a new clue into the Solar system's evolution, but, so far, only one real object has been observed resonating in this mode.

Currently, postdoctoral fellow (PDF) researchers in Canada face challenges due to the precarious nature of their employment and their overall low compensation and benefits coverage. This report presents three themes, written as statements of need, to support an inclusive and thriving PDF community. These themes are the need for better terms of employment and conditions, the need for access to grants by non-permanent research staff, and the need for a sustainable PDF hiring model that considers the outcomes for the PDFs. We make six recommendations: R1. PDFs should be hired and compensated as skilled experts in their areas, not as trainees. R2. Standard PDF hiring practices should be revised to be more inclusive of different life circumstances. - R2.1 Allow PDFs the option of part-time employment. - R2.2 Remove years-since-PhD time limits from PDF jobs. - R2.3 Financially support PDF hires for relocation and visa expenses. R3. CASCA should form a committee to advocate for and provide support to astronomy PDFs in Canada. R4. CASCA should encourage universities to create offices dedicated to their PDFs. R5. PDFs and other PhD-holding term researchers with a host institution should be able to compete for and win grants to self-fund their own research. R6. Astronomy in Canada should hire general-purpose continuing support scientist positions instead of term PDFs to fill project or mission-specific requirements. In short, we ask for prioritization of people over production of papers.

Resonant dynamics plays a significant role in the past evolution and current state of our outer Solar System. The population ratios and spatial distribution of Neptune's resonant populations are direct clues to understanding the history of our planetary system. The orbital structure of the objects in Neptune's 2:1 mean-motion resonance (\emph{twotinos}) has the potential to be a tracer of planetary migration processes. Different migration processes produce distinct architectures, recognizable by well-characterized surveys. However, previous characterized surveys only discovered a few twotinos, making it impossible to model the intrinsic twotino population. With a well-designed cadence and nearly 100\% tracking success, the Outer Solar System Origins Survey (OSSOS) discovered 838 trans-Neptunian objects, of which 34 are securely twotinos with well-constrained libration angles and amplitudes. We use the OSSOS twotinos and the survey characterization parameters via the OSSOS Survey Simulator to inspect the intrinsic population and orbital distributions of twotino. The estimated twotino population, 4400$^{+1500}_{-1100}$ with $H_r<8.66$ (diameter$\sim$100km) at 95\% confidence, is consistent with the previous low-precision estimate. We also constrain the width of the inclination distribution to a relatively narrow value of $\sigma_i$=6$^\circ$$^{+1}_{-1}$, and find the eccentricity distribution is consistent with a Gaussian centered on $e_\mathrm{c}=0.275$ with a width $e_\mathrm{w}=0.06$. We find a single-slope exponential luminosity function with $\alpha=0.6$ for the twotinos. Finally, we for the first time meaningfully constrain the fraction of symmetric twotinos, and the ratio of the leading asymmetric islands; both fractions are in a range of 0.2--0.6. These measurements rule out certain theoretical models of Neptune's migration history.

There is a vibrant and effective planetary science community in Canada. We do research in the areas of meteoritics, asteroid and trans-Neptunian object orbits and compositions, and space weather, and are involved in space probe missions to study planetary surfaces and interiors. For Canadian planetary scientists to deliver the highest scientific impact possible, we have several recommendations. Our top recommendation is to join LSST and gain access to the full data releases by hosting a data centre, which could be done by adding to the CADC, which is already highly involved in hosting planetary data and supporting computational modelling for orbital studies. We also support MSE, which can provide spectroscopy and thus compositional information for thousands of small bodies. We support a Canadian-led microsatellite, POEP, which will provide small body sizes by measuring occultations. We support the idea of piggybacking space weather instruments on other astronomical space probes to provide data for the space weather community. Many Canadian planetary scientists are involved in space probe missions, but through haphazard and temporary arrangements like co-appointments at US institutions, so we would like the community to support Canadian researchers to participate in these large, international missions.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

(Abridged) Canadian astronomy has, for decades, benefited from access to observatories and participating in international consortia on one of the best astronomical sites in the world: Maunakea. However, Maunakea is part of the unceded territory of the Native Hawaiian peoples and has always been of special significance to Hawaiian culture. The use of the summit and its science reserve has created tensions in the past decade, particularly with the development of the Thirty Meter Telescope. A meaningful and respectful response from the International astronomy community is still lacking. It is expected that the LRP 2020 will continue to support Canadian astronomy on Maunakea so a better official statement on the position and involvement of CASCA should be prepared. In this paper we present recommendations, based on the United Nation Declaration for the Rights of Indigenous Peoples, for the Canadian astronomical community to better support Indigenous rights on Maunakea and Hawaii while providing clear guidelines for the astronomical community to participate in activities conducted on Indigenous land. This framework is designed to motivate conversations with Indigenous communities regarding our place on Indigenous lands and our roles, and responsibilities toward the communities we are working with. Furthermore, we propose this framework as a basis for engaging with communities around the world regarding consent for astronomical facilities.

(Abridged) The Truth and Reconciliation Commission of Canada published its calls to action in 2015 with 94 recommendations. Many of these 94 recommendations are directly related to education, language, and culture, some of which the Canadian Astronomy community can address and contribute to as part of reconciliation. The Canadian Astronomy community has an additional obligation since it benefits from facilities on Indigenous territories across Canada and the world. Furthermore, Indigenous people are still underrepresented at all levels in Canadian astronomy. The purpose of this Community Paper is to develop recommendations for the Canadian astronomy community to support Indigenous inclusion in the science community, support Indigenous learning by developing Indigenous-based learning materials and facilitate access to professionals and science activities, and to recognize and acknowledge the great contributions of Indigenous communities to our science activities. As part of this work we propose the ten following recommendations for CASCA as an organization and throughout this Community Paper we will include additional recommendations for individuals: astronomers, students and academics.

Looking at the orbits of small bodies with large semimajor axes, we are compelled to see patterns. Some of these patterns are noted as strong indicators of new or hidden processes in the outer Solar System, others are substantially generated by observational biases, and still others may be completely overlooked. We can gain insight into the current and past structure of the outer Solar System through a careful examination of these orbit patterns. In this chapter, we discuss the implications of the observed orbital distribution of distant trans-Neptunian objects (TNOs). We start with some cautions on how observational biases must affect the known set of TNO orbits. Some of these biases are intrinsic to the process of discovering TNOs, while others can be reduced or eliminated through careful observational survey design. We discuss some orbital element correlations that have received considerable attention in the recent literature. We examine the known TNOs in the context of the gravitational processes that the known Solar System induces in orbital distributions. We discuss proposed new elements of the outer Solar System, posited ancient processes, and the types of TNO orbital element distributions that they predict to exist. We conclude with speculation.

Most known trans-Neptunian objects (TNOs) gravitationally scattering off the giant planets have orbital inclinations consistent with an origin from the classical Kuiper belt, but a small fraction of these "scattering TNOs" have inclinations that are far too large (i > 45 deg) for this origin. These scattering outliers have previously been proposed to be interlopers from the Oort cloud or evidence of an undiscovered planet. Here we test these hypotheses using N-body simulations and the 69 centaurs and scattering TNOs detected in the Outer Solar Systems Origins Survey and its predecessors. We confirm that observed scattering objects cannot solely originate from the classical Kuiper belt, and we show that both the Oort cloud and a distant planet generate observable highly inclined scatterers. Although the number of highly inclined scatterers from the Oort Cloud is ~3 times less than observed, Oort cloud enrichment from the Sun's galactic migration or birth cluster could resolve this. Meanwhile, a distant, low-eccentricity 5 Earth-mass planet replicates the observed fraction of highly inclined scatterers, but the overall inclination distribution is more excited than observed. Furthermore, the distant planet generates a longitudinal asymmetry among detached TNOs that is less extreme than often presumed, and its direction reverses across the perihelion range spanned by known TNOs. More complete models that explore the dynamical origins of the planet are necessary to further study these features. With observational biases well-characterized, our work shows that the orbital distribution of detected scattering bodies is a powerful constraint on the unobserved distant solar system.

The migration of Neptune's resonances through the proto-Kuiper belt has been imprinted in the distribution of small bodies in the outer Solar System. Here we analyze five published Neptune migration models in detail, focusing on the high pericenter distance (high-q) trans-Neptunian Objects (TNOs) near Neptune's 5:2 and 3:1 mean-motion resonances, because they have large resonant populations, are outside the main classical belt, and are relatively isolated from other strong resonances. We compare the observationally biased output from these dynamical models with the detected TNOs from the Outer Solar System Origins Survey, via its Survey Simulator. All of the four new OSSOS detections of high-q non-resonant TNOs are on the Sunward side of the 5:2 and 3:1 resonances. We show that even after accounting for observation biases, this asymmetric distribution cannot be drawn from a uniform distribution of TNOs at 2sigma confidence. As shown by previous work, our analysis here tentatively confirms that the dynamical model that uses grainy slow Neptune migration provides the best match to the real high-q TNO orbital data. However, due to extreme observational biases, we have very few high-q TNO discoveries with which to statistically constrain the models. Thus, this analysis provides a framework for future comparison between the output from detailed, dynamically classified Neptune migration simulations and the TNO discoveries from future well-characterized surveys. We show that a deeper survey (to a limiting r-magnitude of 26.0) with a similar survey area to OSSOS could statistically distinguish between these five Neptune migration models.

The Outer Solar System Origins Survey (OSSOS), a wide-field imaging program in 2013-2017 with the Canada-France-Hawaii Telescope, surveyed 155 deg$^{2}$ of sky to depths of $m_r = 24.1$-25.2. We present 838 outer Solar System discoveries that are entirely free of ephemeris bias. This increases the inventory of trans-Neptunian objects (TNOs) with accurately known orbits by nearly 50%. Each minor planet has 20-60 Gaia/Pan-STARRS-calibrated astrometric measurements made over 2-5 oppositions, which allows accurate classification of their orbits within the trans-Neptunian dynamical populations. The populations orbiting in mean-motion resonance with Neptune are key to understanding Neptune's early migration. Our 313 resonant TNOs, including 132 plutinos, triple the available characterized sample and include new occupancy of distant resonances out to semi-major axis $a \sim 130$ au. OSSOS doubles the known population of the non-resonant Kuiper belt, providing 436 TNOs in this region, all with exceptionally high-quality orbits of $a$ uncertainty $\sigma_{a} \leq 0.1\%$; they show the belt exists from $a \gtrsim 37$ au, with a lower perihelion bound of $35$ au. We confirm the presence of a concentrated low-inclination $a\simeq 44$ au "kernel" population and a dynamically cold population extending beyond the 2:1 resonance. We finely quantify the survey's observational biases. Our survey simulator provides a straightforward way to impose these biases on models of the trans-Neptunian orbit distributions, allowing statistical comparison to the discoveries. The OSSOS TNOs, unprecedented in their orbital precision for the size of the sample, are ideal for testing concepts of the history of giant planet migration in the Solar System.

We discuss the detection in the Outer Solar System Origins Survey (OSSOS) of two objects in Neptune's distant 9:1 mean motion resonance at semimajor axis $a\approx~130$~au. Both objects are securely resonant on 10~Myr timescales, with one securely in the 9:1 resonance's leading asymmetric libration island and the other in either the symmetric or trailing asymmetric island. These objects are the largest semimajor axis objects with secure resonant classifications, and their detection in a carefully characterized survey allows for the first robust resonance population estimate beyond 100~au. The detection of these objects implies a 9:1 resonance population of $1.1\times10^4$ objects with $H_r<8.66$ ($D~\gtrsim~100$~km) on similar orbits (95\% confidence range of $\sim0.4-3\times10^4$). Integrations over 4~Gyr of an ensemble of clones spanning these objects' orbit fit uncertainties reveal that they both have median resonance occupation timescales of $\sim1$~Gyr. These timescales are consistent with the hypothesis that these objects originate in the scattering population but became transiently stuck to Neptune's 9:1 resonance within the last $\sim1$~Gyr of solar system evolution. Based on simulations of a model of the current scattering population, we estimate the expected resonance sticking population in the 9:1 resonance to be 1000-4500 objects with $H_r<8.66$; this is marginally consistent with the OSSOS 9:1 population estimate. We conclude that resonance sticking is a plausible explanation for the observed 9:1 population, but we also discuss the possibility of a primordial 9:1 population, which would have interesting implications for the Kuiper belt's dynamical history.

All surveys include observational biases, which makes it impossible to directly compare properties of discovered trans-Neptunian Objects (TNOs) with dynamical models. However, by carefully keeping track of survey pointings on the sky, detection limits, tracking fractions, and rate cuts, the biases from a survey can be modelled in Survey Simulator software. A Survey Simulator takes an intrinsic orbital model (from, for example, the output of a dynamical Kuiper belt emplacement simulation) and applies the survey biases, so that the biased simulated objects can be directly compared with real discoveries. This methodology has been used with great success in the Outer Solar System Origins Survey (OSSOS) and its predecessor surveys. In this chapter, we give four examples of ways to use the OSSOS Survey Simulator to gain knowledge about the true structure of the Kuiper Belt. We demonstrate how to statistically compare different dynamical model outputs with real TNO discoveries, how to quantify detection biases within a TNO population, how to measure intrinsic population sizes, and how to use upper limits from non-detections. We hope this will provide a framework for dynamical modellers to statistically test the validity of their models.

The scattering trans-Neptunian Objects (TNOs) can be measured to smaller sizes than any other distant small-body population. We use the largest sample yet obtained, 68 discoveries, primarily by the Outer Solar System Origins Survey (OSSOS), to constrain the slope of its luminosity distribution, with sensitivity to much fainter absolute $H$ magnitudes than previous work. Using the analysis technique in Shankman et al. (2016), we confirm that a single slope for the $H$-distribution is not an accurate representation of the scattering TNOs and Centaurs, and that a break in the distribution is required, in support of previous conclusions. A bright-end slope of $\alpha_b=0.9$ transitioning to a faint-end slope $\alpha_f$ of 0.4-0.5 with a differential number contrast $c$ from 1 (a knee) to 10 (a divot) provides an acceptable match to our data. We find that break magnitudes $H_b$ of 7.7 and 8.3, values both previously suggested for dynamically hot Kuiper belt populations, are equally non-rejectable for a range of $\alpha_f$ and $c$ in our statistical analysis. Our preferred divot $H$-distribution transitions to $\alpha_f=0.5$ with a divot of contrast $c=3$ at $H_b=8.3$, while our preferred knee $H$-distribution transitions to $\alpha_f=0.4$ at $H_b=7.7$. The intrinsic population of scattering TNOs required to match the OSSOS detections is $3\times10^6$ for $H_r<12$, and $9\times10^4$ for $H_r<8.66$ ($D\gtrsim100$~km), with Centaurs having an intrinsic population two orders of magnitude smaller.

We analyze a detailed Nice model simulation of Kuiper Belt emplacement from Brasser & Morbidelli (2013), where Neptune undergoes a high eccentricity phase and migrates outward. In this work, which follows from Pike et al. (2017), we specifically focus on the details of structures within Neptune's mean motion resonances and in the high pericenter population of simulated trans-Neptunian Objects (TNOs). We find several characteristics of these populations which should be observable in the distant Solar System in future large-scale TNO surveys as a diagnostic of whether or not this mode of Neptune migration occurred in the early Solar System. We find that the leading asymmetric libration islands of the $n$:1 resonances are generally much more populated than the trailing islands. We also find the non-resonant high-$q$ population of TNOs should have higher inclinations than the low-$q$ population due to the importance of Kozai cycling during their emplacement histories. Finally, high-$q$ TNOs should be present in roughly equal numbers on either side of distant mean-motion resonances. These predictions contrast with predictions from other Kuiper Belt emplacement simulations, and will be testable by upcoming surveys.

The accumulating, but small, set of large semi-major axis trans-Neptunian objects (TNOs) shows an apparent clustering in the orientations of their orbits. This clustering must either be representative of the intrinsic distribution of these TNOs, or else arise as a result of observation biases and/or statistically expected variations for such a small set of detected objects. The clustered TNOs were detected across different and independent surveys, which has led to claims that the detections are therefore free of observational bias. This apparent clustering has led to the so-called "Planet 9" hypothesis that a super-Earth currently resides in the distant solar system and causes this clustering. The Outer Solar System Origins Survey (OSSOS) is a large program that ran on the Canada-France-Hawaii Telescope from 2013--2017, discovering more than 800 new TNOs. One of the primary design goals of OSSOS was the careful determination of observational biases that would manifest within the detected sample. We demonstrate the striking and non-intuitive biases that exist for the detection of TNOs with large semi-major axes. The eight large semi-major axis OSSOS detections are an independent dataset, of comparable size to the conglomerate samples used in previous studies. We conclude that the orbital distribution of the OSSOS sample is consistent with being detected from a uniform underlying angular distribution.

Debris discs are evidence of the ongoing destructive collisions between planetesimals, and their presence around stars also suggests that planets exist in these systems. In this paper, we present submillimetre images of the thermal emission from debris discs that formed the SCUBA-2 Observations of Nearby Stars (SONS) survey, one of seven legacy surveys undertaken on the James Clerk Maxwell telescope between 2012 and 2015. The overall results of the survey are presented in the form of 850 microns (and 450 microns, where possible) images and fluxes for the observed fields. Excess thermal emission, over that expected from the stellar photosphere, is detected around 49 stars out of the 100 observed fields. The discs are characterised in terms of their flux density, size (radial distribution of the dust) and derived dust properties from their spectral energy distributions. The results show discs over a range of sizes, typically 1-10 times the diameter of the Edgeworth-Kuiper Belt in our Solar System. The mass of a disc, for particles up to a few millimetres in size, is uniquely obtainable with submillimetre observations and this quantity is presented as a function of the host stars' age, showing a tentative decline in mass with age. Having doubled the number of imaged discs at submillimetre wavelengths from ground-based, single dish telescope observations, one of the key legacy products from the SONS survey is to provide a comprehensive target list to observe at high angular resolution using submillimetre/millimetre interferometers (e.g., ALMA, SMA).

This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed populate the outer resonances with orbital distributions consistent with trans-Neptunian object (TNO) detections in semi-major axis, inclination, and eccentricity, while capture into the closest resonances is too efficient. The relative populations of the simulated scattering objects and resonant objects in the 3:1 and 4:1 resonances are also consistent with observed populations based on debiased TNO surveys, but the 5:1 resonance is severely underpopulated compared to population estimates from survey results. Scattering emplacement results in the expected orbital distribution for the majority of the TNO populations, however the origin of the large observed population in the 5:1 resonance remains unexplained.

The orbital element distribution of trans-Neptunian objects (TNOs) with large pericenters has been suggested to be influenced by the presence of an undetected, large planet at >200 AU from the Sun. To find additional observables caused by this scenario, we here present the first detailed emplacement simulation in the presence of a massive ninth planet on the distant Kuiper Belt. We perform 4 Gyr N-body simulations with the currently known Solar System planetary architecture, plus a 10 Earth mass planet with similar orbital parameters to those suggested by Trujillo & Sheppard (2014) or Batygin & Brown (2016), and 10^5 test particles in an initial planetesimal disk. We find that including a distant superearth-mass planet produces a substantially different orbital distribution for the scattering and detached TNOs, raising the pericenters and inclinations of moderate semimajor axis (50<a<500 AU) objects. We test whether this signature is detectable via a simulator with the observational characteristics of four precisely characterized TNO surveys. We find that the qualitatively very distinct Solar System models that include a ninth planet are essentially observationally indistinguishable from an outer Solar System produced solely by the four giant planets. We also find that the mass of the Kuiper Belt's current scattering and detached populations is required to be 3-10 times larger in the presence of an additional planet. We do not find any evidence for clustering of orbital angles in our simulated TNO population. Wide-field, deep surveys targeting inclined high-pericenter objects will be required to distinguish between these different scenarios.

We explore the distant giant planet hypothesis by integrating the large semi-major axis, large pericenter Trans-Neptunian Objects (TNOs) in the presence of the giant planets and an external perturber whose orbit is consistent with the proposed distant, eccentric, and inclined giant planet, so called planet 9. We find that TNOs with semi-major axes greater than 250 au experience some longitude of perihelion shepherding, but that a generic outcome of such evolutions is that the TNOs evolve to larger pericenter orbits, and commonly get raised to retrograde inclinations. This pericenter and inclination evolution requires a massive disk of TNOs (tens of M$_\Earth$) in order to explain the detection of the known sample today. Some of the highly inclined orbits produced by the examined perturbers will be inside of the orbital parameter space probed by prior surveys, implying a missing signature of the 9th planet scenario. The distant giant planet scenarios explored in this work do not reproduce the observed signal of simultaneous clustering in argument of pericenter, longitude of the ascending node, and longitude of perihelion in the region of the known TNOs.

We report the discovery and orbit of a new dwarf planet candidate, 2015 RR$_{245}$, by the Outer Solar System Origins Survey (OSSOS). 2015 RR$_{245}$'s orbit is eccentric ($e=0.586$), with a semi-major axis near 82 au, yielding a perihelion distance of 34 au. 2015 RR$_{245}$ has $g-r = 0.59 \pm 0.11$ and absolute magnitude $H_{r} = 3.6 \pm 0.1$; for an assumed albedo of $p_V = 12$% the object has a diameter of $\sim670$ km. Based on astrometric measurements from OSSOS and Pan-STARRS1, we find that 2015 RR$_{245}$ is securely trapped on ten-Myr timescales in the 9:2 mean-motion resonance with Neptune. It is the first TNO identified in this resonance. On hundred-Myr timescales, particles in 2015 RR$_{245}$-like orbits depart and sometimes return to the resonance, indicating that 2015 RR$_{245}$ likely forms part of the long-lived metastable population of distant TNOs that drift between resonance sticking and actively scattering via gravitational encounters with Neptune. The discovery of a 9:2 TNO stresses the role of resonances in the long-term evolution of objects in the scattering disk, and reinforces the view that distant resonances are heavily populated in the current Solar System. This object further motivates detailed modelling of the transient sticking population.

We report the orbital distribution of the Trans-Neptunian objects (TNOs) discovered during the High Ecliptic Latitude (HiLat) extension of the Canada-France Ecliptic Plane Survey (CFEPS), conducted from June 2006 to July 2009. The HiLat component was designed to address one of the shortcomings of ecliptic surveys (like CFEPS), their lack of sensitivity to high-inclination objects. We searched 701~deg$^2$ of sky ranging from 12$^\circ$ to 85$^\circ$ ecliptic latitude and discovered \lKBO TNOs, with inclinations between 15$^\circ$ to 104$^\circ$. This survey places a very strong constraint on the inclination distribution of the hot component of the classical Kuiper Belt, ruling out any possibility of a large intrinsic fraction of highly inclined orbits. Using the parameterization of \citet{2001AJ....121.2804B}, the HiLat sample combined with CFEPS imposes a width $14^\circ \le \sigma \le 15.5^\circ$, with a best match for $\sigma = 14.5^\circ$. HiLat discovered the first retrograde TNO, 2008~KV$_{42}$, with an almost polar orbit with inclination 104$^\circ$, and (418993), a scattering object with perihelion in the region of Saturn's influence, with $a \sim 400$~AU and $i = 68^\circ$.

We present 1.3 mm observations of the Sun-like star $\tau$ Ceti with the Atacama Large Millimeter/submillimeter Array (ALMA) that probe angular scales of $\sim1$'' (4 AU). This first interferometric image of the $\tau$ Ceti system, which hosts both a debris disk and possible multiplanet system, shows emission from a nearly face-on belt of cold dust with a position angle of $90^\circ$ surrounding an unresolved central source at the stellar position. To characterize this emission structure, we fit parametric models to the millimeter visibilities. The resulting best-fit model yields an inner belt edge of $6.2^{+9.8}_{-4.6}$ AU, consistent with inferences from lower resolution, far-infrared Herschel observations. While the limited data at sufficiently short baselines preclude us from placing stronger constraints on the belt properties and its relation to the proposed five planet system, the observations do provide a strong lower limit on the fractional width of the belt, $\Delta R/R > 0.75$ with $99\%$ confidence. This fractional width is more similar to broad disks such as HD 107146 than narrow belts such as the Kuiper Belt and Fomalhaut. The unresolved central source has a higher flux density than the predicted flux of the stellar photosphere at 1.3 mm. Given previous measurements of an excess by a factor of $\sim2$ at 8.7 mm, this emission is likely due to a hot stellar chromosphere.

We report the discovery, tracking and detection circumstances for 85 trans-Neptunian objects (TNOs) from the first 42 deg$^{2}$ of the Outer Solar System Origins Survey (OSSOS). This ongoing $r$-band Solar System survey uses the 0.9 deg$^{2}$ field-of-view MegaPrime camera on the 3.6 m Canada-France-Hawaii Telescope. Our orbital elements for these TNOs are precise to a fractional semi-major axis uncertainty $<0.1\%$. We achieve this precision in just two oppositions, as compared to the normal 3-5 oppositions, via a dense observing cadence and innovative astrometric technique. These discoveries are free of ephemeris bias, a first for large trans-Neptunian surveys. We also provide the necessary information to enable models of TNO orbital distributions to be tested against our TNO sample. We confirm the existence of a cold "kernel" of objects within the main cold classical Kuiper belt, and infer the existence of an extension of the "stirred" cold classical Kuiper belt to at least several AU beyond the 2:1 mean motion resonance with Neptune. We find that the population model of Petit et al. (2011) remains a plausible representation of the Kuiper belt. The full survey, to be completed in 2017, will provide an exquisitely characterized sample of important resonant TNO populations, ideal for testing models of giant planet migration during the early history of the Solar System.

The first two observational sky "blocks" of the Outer Solar System Origins Survey (OSSOS) have significantly increased the number of well-characterized observed trans-Neptunian objects (TNOs) in Neptune's mean motion resonances. We describe the 31 securely resonant TNOs detected by OSSOS so far, and we use them to independently verify the resonant population models from the Canada-France Ecliptic Plane Survey (CFEPS; Gladman et al. 2012), with which we find broad agreement. We confirm that the 5:2 resonance is more populated than models of the outer Solar System's dynamical history predict; our minimum population estimate shows that the high eccentricity (e>0.35) portion of the resonance is at least as populous as the 2:1 and possibly as populated as the 3:2 resonance. One OSSOS block was well-suited to detecting objects trapped at low libration amplitudes in Neptune's 3:2 resonance, a population of interest in testing the origins of resonant TNOs. We detected three 3:2 objects with libration amplitudes below the cutoff modeled by CFEPS; OSSOS thus offers new constraints on this distribution. The OSSOS detections confirm that the 2:1 resonance has a dynamically colder inclination distribution than either the 3:2 or 5:2 resonances. Using the combined OSSOS and CFEPS 2:1 detections, we constrain the fraction of 2:1 objects in the symmetric mode of libration to be 0.2-0.85; we also constrain the fraction of leading vs. trailing asymmetric librators, which has been theoretically predicted to vary depending on Neptune's migration history, to be 0.05-0.8. Future OSSOS blocks will improve these constraints.

We present a new algorithm for space telescope high contrast imaging of close-to-face-on planetary disks called Optimized Spatially Filtered (OSFi) normalization. This algorithm is used on HR 8799 Hubble Space Telescope (HST) coronagraphic archival data, showing an over-luminosity after reference star point spread function (PSF) subtraction that may be from the inner disk and/or planetesimal belt components of this system. The PSF-subtracted radial profiles in two separate epochs from 2011 and 2012 are consistent with one another, and self-subtraction shows no residual in both epochs. We explore a number of possible false-positive scenarios that could explain this residual flux, including telescope breathing, spectral differences between HR 8799 and the reference star, imaging of the known warm inner disk component, OSFi algorithm throughput and consistency with the standard spider normalization HST PSF subtraction technique, and coronagraph misalignment from pointing accuracy. In comparison to another similar STIS dataset, we find that the over-luminosity is likely a result of telescope breathing and spectral difference between HR 8799 and the reference star. Thus, assuming a non-detection, we derive upper limits on the HR 8799 dust belt mass in small grains. In this scenario, we find that the flux of these micron-sized dust grains leaving the system due to radiation pressure is small enough to be consistent with measurements of other debris disk halos.

We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$\deg$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

The planet candidate Fomalhaut b is bright in optical light but undetected in longer wavelengths, requiring a large, reflective dust cloud. The most recent observations find an extremely eccentric orbit (e ~ 0.8), indicating that Fomalhaut b cannot be the planet that is constraining the system's eccentric debris ring. An irregular satellite swarm around a super-Earth has been proposed, however, explaining the well-constrained debris ring requires an additional planet on an orbit that crosses that of the putative super-Earth. This paper expands upon a second theory: Fomalhaut b is a transient dust cloud produced by a catastrophic collision between planetesimals in the disk. We perform collisional probability simulations of the Fomalhaut debris disk based on the structure of our Kuiper belt, finding that the catastrophic disruption rate of d ~ 100 km bodies in the high-eccentricity scattering component is several per decade. This model paints a picture of the Fomalhaut system as having recently (within ~10-100 Myr) experienced a dynamical instability within its planetary system, which scattered a massive number of planetesimals onto large, high eccentricity orbits similar to that of Fom b. If Fomalhaut b is indeed a dust cloud produced by such a collision, we should soon see another appear, while Fomalhaut b will expand until it is either resolved or becomes too faint to be seen.

$\tau$ Ceti is a nearby, mature G-type star very similar to our Sun, with a massive Kuiper Belt analogue (Greaves et al. 2004) and possible multiplanet system (Tuomi et al. 2013) that has been compared to our Solar System. We present Herschel Space Observatory images of the debris disk, finding the disk is resolved at 70 and 160 microns, and marginally resolved at 250 microns. The Herschel images and infrared photometry from the literature are best modelled using a wide dust annulus with an inner edge between 1-10 AU and an outer edge at ~55 AU, inclined from face-on by 35$\pm$10 degrees, and with no significant azimuthal structure. We model the proposed tightly-packed planetary system of five super-Earths and find that the innermost dynamically stable disk orbits are consistent with the inner edge found by the observations. The photometric modelling, however, cannot rule out a disk inner edge as close to the star as 1 AU, though larger distances produce a better fit to the data. Dynamical modelling shows that the 5 planet system is stable with the addition of a Neptune or smaller mass planet on an orbit outside 5 AU, where the Tuomi et al. analysis would not have detected a planet of this mass.

The dust measured in debris disks traces the position of planetesimal belts. In our Solar System, we are also able to measure the largest planetesimals directly and can extrapolate down to make an estimate of the dust. The zodiacal dust from the asteroid belt is better constrained than the only rudimentary measurements of Kuiper belt dust. Dust models will thus be based on the current orbital distribution of the larger bodies which provide the collisional source. The orbital distribution of many Kuiper belt objects is strongly affected by dynamical interactions with Neptune, and the structure cannot be understood without taking this into account. We present the debiased Kuiper belt as measured by the Canada-France Ecliptic Plane Survey (CFEPS). This model includes the absolute populations for objects with diameters >100 km, measured orbital distributions, and size distributions of the components of the Kuiper belt: the classical belt (hot, stirred, and kernel components), the scattering disk, the detached objects, and the resonant objects (1:1, 5:4, 4:3, 3:2 including Kozai subcomponent, 5:3, 7:4, 2:1, 7:3, 5:2, 3:1, and 5:1). Because a large fraction of known debris disks are consistent with dust at Kuiper belt distances from the host stars, the CFEPS Kuiper belt model provides an excellent starting point for a debris disk model, as the dynamical interactions with planets interior to the disk are well-understood and can be precisely modelled using orbital integrations.

Because of their proximity within the transneptunian region, the plutinos (objects in the 3:2 mean-motion resonance with Neptune) are numerous in flux-limited catalogs, and well-studied theoretically. We perform detailed modelling of the on-sky detection biases for plutinos, with special attention to those that are simultaneously in the Kozai resonance. In addition to the normal 3:2 resonant argument libration, Kozai plutinos also show periodic oscillations in eccentricity and inclination, coupled to the argument of perihelion (omega) oscillation. Due to the mean-motion resonance, plutinos avoid coming to pericenter near Neptune's current position in the ecliptic plane. Because Kozai plutinos are restricted to certain values of omega, perihelion always occurs out of the ecliptic plane, biasing ecliptic surveys against finding these objects. The observed Kozai plutino fraction (fkoz) has been measured by several surveys, finding values between 8% and 25%, while the true fkoz has been predicted to be between 10% and 30% by different giant planet migration simulations. We show that the observed fkoz varies widely depending on the ecliptic latitude and longitude of the survey, so debiasing to find the true ratio is complex. Even a survey that covers most or all of the sky will detect an apparent Kozai fraction that is different from the true fkoz. We present a map of the on-sky plutino Kozai fraction that would be detected by all-sky flux-limited surveys. This will be especially important for the Panoramic Survey Telescope & Rapid Response System (Pan-STARRS) and Large Synoptic Survey Telescope (LSST) projects, which may detect large numbers of plutinos as they sweep the sky. The Kozai fraction and the distribution of the orbital elements of Kozai plutinos may be a diagnostic of giant planet migration; future migration simulations should provide details on their resonant Kozai populations.

We present observations of the debris disk around gamma Doradus, an F1V star, from the Herschel Key Programme DEBRIS (Disc Emission via Bias-free Reconnaissance in the Infrared/Submillimetre). The disk is well-resolved at 70, 100 and 160 micron, resolved along its major axis at 250 micron, detected but not resolved at 350 micron, and confused with a background source at 500 micron. It is one of our best resolved targets and we find it to have a radially broad dust distribution. The modelling of the resolved images cannot distinguish between two configurations: an arrangement of a warm inner ring at several AU (best-fit 4 AU) and a cool outer belt extending from ~55 to 400 AU or an arrangement of two cool, narrow rings at ~70 AU and ~190 AU. This suggests that any configuration between these two is also possible. Both models have a total fractional luminosity of ~10^{-5} and are consistent with the disk being aligned with the stellar equator. The inner edge of either possible configuration suggests that the most likely region to find planets in this system would be within ~55 AU of the star. A transient event is not needed to explain the warm dust's fractional luminosity.

This paper describes Herschel observations of the nearby (8.5pc) G5V multi-exoplanet host star 61 Vir at 70-500micron carried out by the DEBRIS survey. These reveal emission that is extended out to >15arcsec with a morphology that can be fitted by a nearly edge-on (77deg inclination) radially broad (from 30AU to >100AU) debris disk of fractional luminosity 2.7x10^-5, with two unrelated sources nearby that are more prominent at longer wavelengths. Chance alignment with a background object seen at 1.4GHz provides potential for confusion, but the star's 1.4"/yr proper motion allows Spitzer 70micron images to confirm that what we are interpreting as disk emission really is circumstellar. Although the exact shape of the disk's inner edge is not well constrained, the region inside 30AU must be significantly depleted in planetesimals. This is readily explained if there are additional planets in the 0.5-30AU region, but is also consistent with collisional erosion. We also find tentative evidence that the presence of detectable debris around nearby stars correlates with the presence of the lowest mass planets that are detectable in current radial velocity surveys. Out of an unbiased sample of the nearest 60 G stars, 11 are known to have planets, of which 6 (including 61 Vir) have planets that are all less massive than Saturn, and 4 of these have evidence for debris. The debris toward one of these planet-hosts (HD20794) is reported here for the first time. This fraction (4/6) is higher than that expected for nearby field stars (15%), and implies that systems that form low-mass planets are also able to retain bright debris disks. We suggest that this correlation could arise because such planetary systems are dynamically stable and include regions that are populated with planetesimals in the formation process where the planetesimals can remain unperturbed over Gyr timescales.

The transneptunian objects (TNOs) trapped in mean-motion resonances with Neptune were likely emplaced there during planet migration late in the giant-planet formation process. We perform detailed modelling of the resonant objects detected in the Canada-France Ecliptic Plane Survey (CFEPS) in order to provide population estimates and, for some resonances, constrain the complex internal orbital element distribution. Detection biases play a critical role because phase relationships with Neptune make object discovery more likely at certain longitudes. This paper discusses the 3:2, 5:2, 2:1, 3:1, 5:1, 4:3, 5:3, 7:3, 5:4, and 7:4 mean-motion resonances, all of which had CFEPS detections, along with our upper limit on 1:1 Neptune Trojans (which is consistent with their small population estimated elsewhere). For the plutinos (TNOs in the 3:2 resonance) we refine the orbital element distribution given in Kavelaars et al. (2009) and show that steep H-magnitude distributions (N(H) proportional to 10aH, with a=0.8-0.9) are favoured in the range Hg=8-9, and confirm that this resonance does not share the inclination distribution of the classical Kuiper Belt. We give the first population estimate for the 5:2 resonance and find that, to within the uncertainties, the population is equal to that of the 3:2 (13,000 TNOs with Hg < 9.16), whereas the 2:1 population is smaller by a factor of 3-4 compared to the other two resonances. We also measure significant populations inhabiting the 4:3, 5:3, 7:3, 5:4, 7:4, 3:1, and 5:1 resonances, with Hg < 9.16 (D >100 km) populations in the thousands. We compare our intrinsic population and orbital-element distributions with several published models of resonant-TNO production; the most striking discrepancy is that resonances beyond the 2:1 are in reality more heavily populated than in published models.

The Kepler Mission recently identified 997 systems hosting candidate extrasolar planets, many of which are super-Earths. Realizing these planetary systems are candidates to host extrasolar asteroid belts, we use mid-infrared data from the Wide-field Infrared Survey Explorer (WISE) to search for emission from dust in these systems. We find excesses around eight stars, indicating the presence of warm to hot dust (~100-500 K), corresponding to orbital distances of 0.1-10 AU for these solar-type stars. The strongest detection, KOI 1099, demands ~500 K dust interior to the orbit of its exoplanet candidate. One star, KOI 904, may host very hot dust (~1200 K, corresponding to 0.02 AU). Although the fraction of these exoplanet-bearing stars with detectable warm excesses (~3%) is similar to that found by Spitzer surveys of solar-type field stars, the excesses detectable in the WISE data have much higher fractional luminosities (Ldust/L*) than most known debris disks, implying that the fraction with debris disks of comparable luminosity may actually be significantly higher. It is difficult to explain the presence of dust so close to the host stars, generally corresponding to dust rings at radii <0.3 AU; both the collisional and Poynting-Robertson drag timescales to remove dust from the system are hundreds of years or less at these distances. Assuming a steady-state for these systems implies large mass consumption rates with these short removal timescales, meaning that the dust production mechanism in these systems must almost certainly be episodic in nature.

In addition to the Sun, six other stars are known to harbor multiple planets and debris disks: HD 69830, HD 38529, HD 128311, HD 202206, HD 82943 and HR 8799. In this paper we set constraints on the location of the dust-producing planetesimals around the latter four systems. We use a radiative transfer model to analyze the spectral energy distributions of the dust disks (including two new Spitzer IRS spectra presented in this paper), and a dynamical model to assess the long-term stability of the planetesimals' orbits. As members of a small group of stars that show evidence of harboring a multiple planets and planetesimals, their study can help us learn about the diversity of planetary systems.

Na I D lines in the spectrum of the young binary KH 15D have been analyzed in detail. We find an excess absorption component that may be attributed to foreground interstellar absorption, and to gas possibly associated with the solids in the circumbinary disk. The derived column density is log N_NaI = 12.5 cm^-2, centered on a radial velocity that is consistent with the systemic velocity. Subtracting the likely contribution of the ISM leaves log N_NaI ~ 12.3 cm^-2. There is no detectable change in the gas column density across the "knife edge" formed by the opaque grain disk, indicating that the gas and solids have very different scale heights, with the solids being highly settled. Our data support a picture of this circumbinary disk as being composed of a very thin particulate grain layer composed of millimeter-sized or larger objects that are settled within whatever remaining gas may be present. This phase of disk evolution has been hypothesized to exist as a prelude to the formation of planetesimals through gravitational fragmentation, and is expected to be short-lived if much gas were still present in such a disk. Our analysis also reveals the presence of excess Na I emission relative to the comparison spectrum at the radial velocity of the currently visible star that plausibly arises within the magnetosphere of this still-accreting young star.

We have observed 152 nearby solar-type stars with the Infrared Spectrometer (IRS) on the Spitzer Space Telescope. Including stars that met our criteria but were observed in other surveys, we get an overall success rate for finding excesses in the long wavelength IRS band (30-34 micron) of 11.8% +/- 2.4%. The success rate for excesses in the short wavelength band (8.5-12 micron) is ~1% including sources from other surveys. For stars with no excess at 8.5-12 microns, the IRS data set 3 sigma limits of around 1,000 times the level of zodiacal emission present in our solar system, while at 30-34 microns set limits of around 100 times the level of our solar system. Two stars (HD 40136 and HD 10647) show weak evidence for spectral features; the excess emission in the other systems is featureless. If the emitting material consists of large (10 micron) grains as implied by the lack of spectral features, we find that these grains are typically located at or beyond the snow line, ~1-35 AU from the host stars, with an average distance of 14 +/- 6 AU; however smaller grains could be located at significantly greater distances from the host stars. These distances correspond to dust temperatures in the range ~50-450 K. Several of the disks are well modeled by a single dust temperature, possibly indicative of a ring-like structure. However, a single dust temperature does not match the data for other disks in the sample, implying a distribution of temperatures within these disks. For most stars with excesses, we detect an excess at both IRS and MIPS wavelengths. Only three stars in this sample show a MIPS 70 micron excess with no IRS excess, implying that very cold dust is rare around solar-type stars.

The Spitzer Space Telescope has advanced debris disk science tremendously with a wealth of information on debris disks around nearby A, F, G, K and M stars at 24 and 70 microns with the MIPS photometer and at 8-34 microns with IRS. Here we present 160 micron observations of a small sub-set of these stars. At this wavelength, the stellar photospheric emission is negligible and any detected emission corresponds to cold dust in extended Kuiper belt analogs. However, the Spitzer 160 micron observations are limited in sensitivity by the large beam size which results in significant ''noise'' due to cirrus and extragalactic confusion. In addition, the 160 micron measurements suffer from the added complication of a light leak next to the star's position whose flux is proportional to the near-infrared flux of the star. We are able to remove the contamination from the leak and report 160 micron measurements or upper limits for 24 stars. Three stars (HD 10647, HD 207129, and HD 115617) have excesses at 160 micron that we use to constrain the properties of the debris disks around them. A more detailed model of the spectral energy distribution of HD 10647 reveals that the 70 and 160 micron emission could be due to small water ice particles at a distance of 100 AU consistent with Hubble Space Telescope optical imaging of circumstellar material in the system.

We report the spectroscopic detection of mid-infrared emission from the transiting exoplanet HD 209458b. Using archive data taken with the Spitzer/IRS instrument, we have determined the spectrum of HD 209458b between 7.46 and 15.25 microns. We have used two independent methods to determine the planet spectrum, one differential in wavelength and one absolute, and find the results are in good agreement. Over much of this spectral range, the planet spectrum is consistent with featureless thermal emission. Between 7.5 and 8.5 microns, we find evidence for an unidentified spectral feature. If this spectral modulation is due to absorption, it implies that the dayside vertical temperature profile of the planetary atmosphere is not entirely isothermal. Using the IRS data, we have determined the broad-band eclipse depth to be 0.00315 +/- 0.000315, implying significant redistribution of heat from the dayside to the nightside. This work required development of improved methods for Spitzer/IRS data calibration that increase the achievable absolute calibration precision and dynamic range for observations of bright point sources.

We report the results of a spectroscopic search for debris disks surrounding 41 nearby solar type stars, including 8 planet-bearing stars, using the {\it Spitzer Space Telescope}. With accurate relative photometry using the Infrared Spectrometer (IRS) between 7-34 $\micron$ we are able to look for excesses as small as $\sim$2% of photospheric levels with particular sensitivity to weak spectral features. For stars with no excess, the $3\sigma$ upper limit in a band at 30-34 $\mu$m corresponds to $\sim$ 75 times the brightness of our zodiacal dust cloud. Comparable limits at 8.5-13 $\mu$m correspond to $\sim$ 1,400 times the brightness of our zodiacal dust cloud. These limits correspond to material located within the $<$1 to $\sim$5 AU region that, in our solar system, originates from debris associated with the asteroid belt. We find excess emission longward of $\sim$25 $\mu$m from five stars of which four also show excess emission at 70 $\mu$m. This emitting dust must be located around 5-10 AU. One star has 70 micron emission but no IRS excess. In this case, the emitting region must begin outside 10 AU; this star has a known radial velocity planet. Only two stars of the five show emission shortward of 25 $\micron$ where spectral features reveal the presence of a population of small, hot dust grains emitting in the 7-20 $\mu$m band. The data presented here strengthen the results of previous studies to show that excesses at 25 $\micron$ and shorter are rare: only 1 star out of 40 stars older than 1 Gyr or $\sim 2.5$% shows an excess. Asteroid belts 10-30 times more massive than our own appear are rare among mature, solar-type stars.

We report the results of a spectroscopic search for debris disks surrounding 41 nearby solar type stars, including 8 planet-bearing stars, using the {\it Spitzer Space Telescope}. With accurate relative photometry using the Infrared Spectrometer (IRS) between 7-34 $\micron$ we are able to look for excesses as small as $\sim$2% of photospheric levels with particular sensitivity to weak spectral features. For stars with no excess, the $3\sigma$ upper limit in a band at 30-34 $\mu$m corresponds to $\sim$ 75 times the brightness of our zodiacal dust cloud. Comparable limits at 8.5-13 $\mu$m correspond to $\sim$ 1,400 times the brightness of our zodiacal dust cloud. These limits correspond to material located within the $<$1 to $\sim$5 AU region that, in our solar system, originates from debris associated with the asteroid belt. We find excess emission longward of $\sim$25 $\mu$m from five stars of which four also show excess emission at 70 $\mu$m. This emitting dust must be located around 5-10 AU. One star has 70 micron emission but no IRS excess. In this case, the emitting region must begin outside 10 AU; this star has a known radial velocity planet. Only two stars of the five show emission shortward of 25 $\micron$ where spectral features reveal the presence of a population of small, hot dust grains emitting in the 7-20 $\mu$m band. The data presented here strengthen the results of previous studies to show that excesses at 25 $\micron$ and shorter are rare: only 1 star out of 40 stars older than 1 Gyr or $\sim 2.5$% shows an excess. Asteroid belts 10-30 times more massive than our own appear are rare among mature, solar-type stars.