We conduct a search for strong gravitational lenses in the Dark Energy Survey (DES) Year 6 imaging data. We implement a pre-trained Vision Transformer (ViT) for our machine learning (ML) architecture and adopt interactive machine learning to construct a training sample with multiple classes to address common types of false positives. Our ML model reduces similar to 236 million DES cutout images to 22,564 targets of interest, including similar to 85% of previously reported galaxy-galaxy lens candidates discovered in DES. These targets were visually inspected by citizen scientists, who ruled out similar to 90% as false positives. Of the remaining 2618 candidates, 149 were expert-classified as "definite" lenses and 516 as "probable" lenses, for a total of 665 systems, with 147 of these candidates being newly identified. Additionally, we trained a second ViT to find double-source plane lens systems, finding at least one double-source system. Our main ViT excels at identifying galaxy-galaxy lenses, consistently assigning high scores to candidates with high expert assessments. The top 800 ViT-scored images include similar to 100 of our "definite" lens candidates. This selection is an order of magnitude higher in purity than previous convolutional neural-network-based lens searches and demonstrates the feasibility of applying our methodology for discovering large samples of lenses in future surveys.

Characterization of the redshift distribution of ensembles of galaxies is pivotal for large scale structure cosmological studies. In this work, we focus on improving the Self-Organizing Map (SOM) methodology for photometric redshift estimation (SOMPZ), specifically in anticipation of the Dark Energy Survey Year 6 (DES Y6) data. This data set, featuring deeper and fainter galaxies than DES Year 3 (DES Y3), demands adapted techniques to ensure accurate recovery of the underlying redshift distribution. We investigate three strategies for enhancing the existing SOM-based approach used in DES Y3: 1) Replacing the Y3 SOM algorithm with one tailored for redshift estimation challenges; 2) Incorporating $\textit{g}$-band flux information to refine redshift estimates (i.e. using $\textit{griz}$ fluxes as opposed to only $\textit{riz}$); 3) Augmenting redshift data for galaxies where available. These methods are applied to DES Y3 data, and results are compared to the Y3 fiducial ones. Our analysis indicates significant improvements with the first two strategies, notably reducing the overlap between redshift bins. By combining strategies 1 and 2, we have successfully managed to reduce redshift bin overlap in DES Y3 by up to 66$\%$. Conversely, the third strategy, involving the addition of redshift data for selected galaxies as an additional feature in the method, yields inferior results and is abandoned. Our findings contribute to the advancement of weak lensing redshift characterization and lay the groundwork for better redshift characterization in DES Year 6 and future stage IV surveys, like the Rubin Observatory.

We describe the photometric dataset assembled from the full 6 yr of observations by the Dark Energy Survey (DES) in support of static-sky cosmology analyses. DES Y6 Gold is a curated dataset derived from DES Data Release 2 (DR2) that incorporates improved measurement, photometric calibration, object classification and value-added information. Y6 Gold comprises nearly 5000 deg ^2 of grizY imaging in the south Galactic cap and includes 669 million objects with a depth of i _AB ∼ 23.4 mag at a signal-to-noise ratio ∼ 10 for extended objects and a top-of-the-atmosphere photometric uniformity <2 mmag. Y6 Gold augments DES DR2 with simultaneous fits to multiepoch photometry for more robust galaxy shapes, colors, and photometric redshift estimates. Y6 Gold features improved morphological star–galaxy classification with an efficiency of 98.6% and a contamination of 0.8% for galaxies with 17.5 < i _AB < 22.5. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used for cosmology analyses. After quality selections, benchmark samples contain 448 million galaxies and 120 million stars. This publication is complemented by data access and documentation.

Modeling the intrinsic alignment (IA) of galaxies poses a challenge to weak lensing analyses. Using the Dark Energy Survey Year 3 shape catalog, we expect less impact from IA when we limit the sample to blue, star-forming galaxies. The cosmological parameter constraints from this BLUE cosmic shear sample are stable to IA model choice, unlike passive galaxies in the full DES Y3 sample, the goodness-of-fit is improved and the Omega(m) and S-8 better agree with the observations from Planck on the cosmic microwave background. Mitigating IA with sample selection in DES, rather than flexible model choices, can reduce uncertainty in S-8 by a factor of 1.5.

Using 1533 type Ia supernovae (SNe Ia) from the five-year sample of the Dark Energy Survey (DES), we investigate the relationship between the projected galactocentric separation of the SNe and their host galaxies and their light curves and standardization. We show, for the first time, that the difference in SN Ia post-standardization brightnesses between high and low-mass hosts reduces from 0.078±0.011 mag in the full sample to 0.036 ± 0.018 mag for SNe Ia located in the outer regions of their host galaxies, while increasing to 0.100 ± 0.014 mag for SNe in the inner regions. The difference in the size of the mass step between inner and outer regions is 0.064±0.023 mag. In these inner regions, the step can be reduced (but not removed) using a model where the R_V of dust along the line-of-sight to the SN changes as a function of galaxy properties. We investigate the remaining difference using the distributions of the SN Ia stretch parameter to test the inferred age of SN progenitors. Comparing red (older) environments only, outer regions have a higher proportion of high-stretch SNe and a more homogeneous stretch distribution. However, this effect cannot explain the reduction in significance of any Hubble residual step in outer regions. We conclude that the standardized distances of SNe Ia located in the outer regions of galaxies are less affected by their global host galaxy properties than those in the inner regions.

The present state of cosmology is facing a crisis where there is a fundamental disagreement in measurements of the Hubble constant (H-0), with significant tension between the early and late Universe methods. Type Ia supernovae (SNe Ia) are important to measuring H-0 through the astronomical distance ladder. However, there remains potential to better standardize SN Ia light curves by using known dependencies on host galaxy properties after the standard light curve width and colour corrections have been applied to the peak SN Ia luminosities. To explore this, we use the 5-yr photometrically identified SNe Ia sample obtained by the Dark Energy Survey, along with host galaxy spectra obtained by the Australian Dark Energy Survey. Using host galaxy spectroscopy, we find a significant trend with the equivalent width (EW) of the [O II] lambda lambda 3727, 29 doublet, a proxy for specific star formation rate, and Hubble residuals. We find that the correlation with [O II] EW is a powerful alternative to the commonly used mass step after initial light-curve corrections. Applying this [O II] EW correction to 20 SNe Ia in calibrator galaxies observed with WiFeS, we examined the impact on SN Ia absolute magnitudes and H-0. Our [O II] EW corrections result in H-0 values ranging between 73.04 and 73.51 kms(-1)Mpc(-1), with a combined statistical and systematic uncertainty of similar to 1.31kms(-1)Mpc(-1). However, even with this additional correction, the impact of host galaxy properties in standardizing SNe Ia appears limited in reducing the current tension (similar to 5 sigma) with the cosmic microwave background result for H-0.

We measure the current expansion rate of the Universe, Hubble's constant H-0, by calibrating the absolute magnitudes of supernovae to distances measured by baryon acoustic oscillations (BAO). This 'inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning 0.01 < z < 1.13 anchored to the recent baryon acoustic oscillation measurements from Dark Energy Spectroscopic Instrument (DESI) spanning 0.30 < z(eff) < 2.33. To trace cosmology to z = 0, we use the third-, fourth-, and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find H-0 = 67. 19(-0.64)(+0.66) km s(-1) Mpc(-1), in agreement with the value found by Planck without the need to assume Flat-Lambda CDM. However, the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

We present a study of the weak lensing inferred matter profiles Delta Sigma(R) of 698 South Pole Telescope (SPT) thermal Sunyaev-Zel'dovich effect (tSZE) selected and MCMF optically confirmed galaxy clusters in the redshift range 0.25 < z < 0.94 that have associated weak gravitational lensing shear profiles from the Dark Energy Survey (DES). Rescaling these profiles to account for the mass dependent size and the redshift dependent density produces average rescaled matter profiles Delta Sigma(R/R-200c)/(rho R-crit(200c)) with a lower dispersion than the unscaled Delta Sigma(R) versions, indicating a significant degree of self-similarity. Galaxy clusters from hydrodynamical simulations also exhibit matter profiles that suggest a high degree of self-similarity, with RMS variation among the average rescaled matter profiles with redshift and mass falling by a factor of approximately six and 23, respectively, compared to the unscaled average matter profiles. We employed this regularity in a new Bayesian method for weak lensing mass calibration that employs the so-called cluster mass posterior P(M-200|zeta, lambda, z), which describes the individual cluster masses given their tSZE (zeta) and optical (lambda, z) observables. This method enables simultaneous constraints on richness lambda-mass and tSZE detection significance zeta-mass relations using average rescaled cluster matter profiles. We validated the method using realistic mock datasets and present observable-mass relation constraints for the SPTxDES sample, where we constrained the amplitude, mass trend, redshift trend, and intrinsic scatter. Our observable-mass relation results are in agreement with the mass calibration derived from the recent cosmological analysis of the SPTxDES data based on a cluster-by-cluster lensing calibration. Our new mass calibration technique offers a higher efficiency when compared to the single cluster calibration technique. We present new validation tests of the observable-mass relation that indicate the underlying power-law form and scatter are adequate to describe the real cluster sample but that also suggest a redshift variation in the intrinsic scatter of the lambda-mass relation may offer a better description. In addition, the average rescaled matter profiles offer high signal-to-noise ratio (S/N) constraints on the shape of real cluster matter profiles, which are in good agreement with available hydrodynamical Lambda CDM simulations. This high S/N profile contains information about baryon feedback, the collisional nature of dark matter, and potential deviations from general relativity.

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the Universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3 x 2pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining Lambda cold dark matter (Lambda CDM) parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure Omega(m) = 0.300 +/- 0.017 and sigma(8) = 0.797 +/- 0.026. Compared to constraints from Planck primary cosmic microwave background (CMB) anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 (1.2 sigma) for the two-parameter difference. We further obtain S-8 equivalent to sigma(8)(Omega(m)/0.3)(0.5) = 0.796 +/- 0.013 which is lower than the Planck measurement at the 1.6 sigma level. The combined SPT cluster, DES 3 x 2pt, and Planck datasets mildly prefer a nonzero positive neutrino mass, with a 95% upper limit Sigma m(nu) < 0.25 eV on the sum of neutrino masses. Assuming a wCDM model, we constrain the dark energy equation of state parameter w = -1.15(-0.17)(+0.23) and when combining with Planck primary CMB anisotropies, we recover w = -1.20(-0.09)(+0.15), a 1.7 sigma difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology and our analysis paves the way for upcoming joint analyses of next-generation datasets.

Using a sample of 2800 galaxy clusters identified in the Dark Energy Survey across the redshift range $0.20 < z < 0.60$, we characterize the hierarchical assembly of Bright Central Galaxies (BCGs) and the surrounding intracluster light (ICL). To quantify hierarchical formation we use the stellar mass - halo mass (SMHM) relation, comparing the halo mass, estimated via the mass-richness relation, to the stellar mass within the BCG+ICL system. Moreover, we incorporate the magnitude gap (M14), the difference in brightness between the BCG (measured within 30 kpc) and 4th brightest cluster member galaxy within 0.5 $R_{200,c}$, as a third parameter in this linear relation. The inclusion of M14, which traces BCG hierarchical growth, increases the slope and decreases the intrinsic scatter, highlighting that it is a latent variable within the BCG+ICL SMHM relation. Moreover, the correlation with M14 decreases at large radii. However, the stellar light within the BCG+ICL transition region (30 kpc - 80 kpc) most strongly correlates with halo mass and has a statistically significant correlation with M14. Since the transition region and M14 are independent measurements, the transition region may grow due to the BCG's hierarchical formation. Additionally, as M14 and ICL result from hierarchical growth, we use a stacked sample and find that clusters with large M14 values are characterized by larger ICL and BCG+ICL fractions, which illustrates that the merger processes that build the BCG stellar mass also grow the ICL. Furthermore, this may suggest that M14 combined with the ICL fraction can identify dynamically relaxed clusters.

We present cosmological constraints from the Dark Energy Camera All Data Everywhere (DECADE) cosmic shear analysis. This work uses shape measurements for 107 million galaxies measured through Dark Energy Camera (DECam) imaging of $5,\!412$ deg$^2$ of sky that is outside the Dark Energy Survey (DES) footprint. We derive constraints on the cosmological parameters $S_8 = 0.791^{+0.027}_{-0.032}$ and $\Omega_{\rm m} =0.269^{+0.034}_{-0.050}$ for the $\Lambda$CDM model, which are consistent with those from other weak lensing surveys and from the cosmic microwave background. We combine our results with cosmic shear results from DES Y3 at the likelihood level, since the two datasets span independent areas on the sky. The combined measurements, which cover $\approx\! 10,\!000$ deg$^2$, prefer $S_8 = 0.791 \pm 0.023$ and $\Omega_{\rm m} = 0.277^{+0.034}_{-0.046}$ under the $\Lambda$CDM model. These results are the culmination of a series of rigorous studies that characterize and validate the DECADE dataset and the associated analysis methodologies (Anbajagane et. al 2025a,b,c). Overall, the DECADE project demonstrates that the cosmic shear analysis methods employed in Stage-III weak lensing surveys can provide robust cosmological constraints for fairly inhomogeneous datasets. This opens the possibility of using data that have been previously categorized as ``unusable'' for cosmic shear analyses, thereby increasing the statistical power of upcoming weak lensing surveys.

Synthetic source injection (SSI), the insertion of sources into pixel-level on-sky images, is a powerful method for characterizing object detection and measurement in wide-field, astronomical imaging surveys. Within the Dark Energy Survey (DES), SSI plays a critical role in characterizing all necessary algorithms used in converting images to catalogs, and in deriving quantities needed for the cosmology analysis, such as object detection rates, galaxy redshift estimation, galaxy magnification, star-galaxy classification, and photometric performance. We present here a source injection catalog of 146 million injections spanning the entire 5000 deg^2 DES footprint, generated using the Balrog SSI pipeline. Through this sample, we demonstrate that the DES Year 6 (Y6) image processing pipeline provides accurate estimates of the object properties, for both galaxies and stars, at the percent-level, and we highlight specific regimes where the accuracy is reduced. We then show the consistency between SSI and data catalogs, for all galaxy samples developed within the weak lensing and galaxy clustering analyses of DES Y6. The consistency between the two catalogs also extends to their correlations with survey observing properties (seeing, airmass, depth, extinction, etc.). Finally, we highlight a number of applications of this catalog to the DES Y6 cosmology analysis. This dataset is the largest SSI catalog produced at this fidelity and will serve as a key testing ground for exploring the utility of SSI catalogs in upcoming surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time.

The number density , redshift evolution of optically selected galaxy clusters offer an independent measurement of the amplitude of matter fluctuations, S 8 . However, recent results have shown that clusters chosen by the redMaPPer algorithm show richness-dependent biases that affect the weak lensing signals and number densities of clusters, increasing uncertainty in the cluster mass calibration and reducing their constraining power. In this work, we evaluate an alternative cluster proxy, outskirt stellar mass, M out , defined as the total stellar mass within a [50, 100] kpc envelope centered on a massive galaxy. This proxy exhibits scatter comparable to redMaPPer richness, A, but is less likely to be subject to projection effects. We compare the Dark Energy Survey Year 3 redMaPPer cluster catalog with a M out selected cluster sample from the Hyper-Suprime Camera survey. We use weak lensing measurements to quantify and compare the scatter of M out and A with halo mass. Our results show M out has a scatter consistent with A, with a similar halo mass dependence , that both proxies contain unique information about the underlying halo mass. We find A-selected samples introduce features into the measured Delta Sigma signal that are not well fit by a log- normal scatter only model, absent in M out selected samples. Our findings suggest that M out offers an alternative for cluster selection with more easily calibrated selection biases, at least at the generally lower richnesses probed here. Combining both proxies may yield a mass proxy with a lower scatter and more tractable selection biases, enabling the use of lower mass clusters in cosmology. Finally, we find the scatter and slope in the A - M out scaling relation to be 0.49 + 0.02 and 0.38 + 0.09.

We present the pipeline for the cosmic shear analysis of the Dark Energy Camera All Data Everywhere (DECADE) weak lensing dataset: a catalog consisting of 107 million galaxies observed by the Dark Energy Camera (DECam) in the northern Galactic cap. The catalog derives from a large number of disparate observing programs and is therefore more inhomogeneous across the sky compared to existing lensing surveys. First, we use simulated data-vectors to show the sensitivity of our constraints to different analysis choices in our inference pipeline, including sensitivity to residual systematics. Next we use simulations to validate our covariance modeling for inhomogeneous datasets. Finally, we show that our choices in the end-to-end cosmic shear pipeline are robust against inhomogeneities in the survey, by extracting relative shifts in the cosmology constraints across different subsets of the footprint/catalog and showing they are all consistent within 1σ to 2σ. This is done for forty-six subsets of the data and is carried out in a fully consistent manner: for each subset of the data, we re-derive the photometric redshift estimates, shear calibrations, survey transfer functions, the data vector, measurement covariance, and finally, the cosmological constraints. Our results show that existing analysis methods for weak lensing cosmology can be fairly resilient towards inhomogeneous datasets. This also motivates exploring a wider range of image data for pursuing such cosmological constraints.

Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev-Zel'dovich (tSZ) maps to measure the anisotropic distribution of hot gas $2.5-40$ Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of $\sim15$ Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a $\sim3\times$ more significant measurement of the extended gas compared to the technique's proof-of-concept. Decomposing stacks into cosine multipoles of order $m$, we detect a dipole ($m=1$) and quadrupole ($m=2$) at $8-10\sigma$, as well as evidence for $m=4$ signal at up to $6\sigma$, indicating sensitivity to late-time non-Gaussianity. We compare to the Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

Our understanding of the γ -ray sky has improved dramatically in the past decade, however, the unresolved γ -ray background (UGRB) still has a potential wealth of information about the faintest γ -ray sources pervading the Universe. Statistical cross-correlations with tracers of cosmic structure can indirectly identify the populations that most characterize the γ -ray background. In this study, we analyze the angular correlation between the γ -ray background and the matter distribution in the Universe as traced by gravitational lensing, leveraging more than a decade of observations from the Fermi -Large Area Telescope (LAT) and 3 years of data from the Dark Energy Survey (DES). We detect a correlation at signal-to-noise ratio of 8.9. Most of the statistical significance comes from large scales, demonstrating, for the first time, that a substantial portion of the UGRB aligns with the mass clustering of the Universe as traced by weak lensing. Blazars provide a plausible explanation for this signal, especially if those contributing to the correlation reside in halos of large mass (∼ 10 14 M ⊙ ) and account for approximately 30–40% of the UGRB above 10 GeV. Additionally, we observe a preference for a curved γ -ray energy spectrum, with a log-parabolic shape being favored over a power-law. We also discuss the possibility of modifications to the blazar model and the inclusion of additional γ -ray sources, such as star-forming galaxies, misalinged active galactic nuclei, or particle dark matter.

We present the Metadetection weak lensing galaxy shape catalogue from the six-year Dark Energy Survey (DES Y6) imaging data. This dataset is the final release from DES, spanning 4422 deg$^2$ of the southern sky. We describe how the catalogue was constructed, including the two new major processing steps, cell-based image coaddition and shear measurements with Metadetection. The DES Y6 Metadetection weak lensing shape catalogue consists of 151,922,791 galaxies detected over riz bands, with an effective number density of $n_{\rm eff}$ =8.22 galaxies per arcmin$^2$ and shape noise of $\sigma_e$ = 0.29. We carry out a suite of validation tests on the catalogue, including testing for PSF leakage, testing for the impact of PSF modeling errors, and testing the correlation of the shear measurements with galaxy, PSF, and survey properties. In addition to demonstrating that our catalogue is robust for weak lensing science, we use the DES Y6 image simulation suite (Mau, Becker et al. 2025) to estimate the overall multiplicative shear bias of our shear measurement pipeline. We find no detectable multiplicative bias at the roughly half-percent level, with m = (3.4 $\pm$ 6.1) x $10^{-3}$, at 3$\sigma$ uncertainty. This is the first time both cell-based coaddition and Metadetection algorithms are applied to observational data, paving the way to the Stage-IV weak lensing surveys.

We present the point-spread function (PSF) modeling for weak lensing shear measurement using the full six years of the Dark Energy Survey (DES Y6) data. We review the PSF estimation procedure using the PIFF (PSFs In the Full FOV) software package and describe the key improvements made to PIFF and modeling diagnostics since the DES year three (Y3) analysis: (i) use of external Gaia and infrared photometry catalogs to ensure higher purity of the stellar sample used for model fitting, (ii) addition of color-dependent PSF modeling, the first for any weak lensing analysis, and (iii) inclusion of model diagnostics inspecting fourth-order moments, which can bias weak lensing measurements to a similar degree as second-order modeling errors. Through a comprehensive set of diagnostic tests, we demonstrate the improved accuracy of the Y6 models evident in significantly smaller systematic errors than those of the Y3 analysis, in which all g band data were excluded due to insufficiently accurate PSF models. For the Y6 weak lensing analysis, we include g band photometry data in addition to the riz bands, providing a fourth band for photometric redshift estimation. Looking forward to the next generation of wide-field surveys, we describe several ongoing improvements to PIFF, which will be the default PSF modeling software for weak lensing analyses for the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

Cosmology from weak gravitational lensing has been limited by astrophysical uncertainties in baryonic feedback and intrinsic alignments. By calibrating these effects using external data, we recover non-linear information, achieving a 2

For the 696 trans-Neptunian objects (TNOs) with absolute magnitudes 5.5 < H _r < 8.2 detected in the Dark Energy Survey, we characterize the relationships between their dynamical state and physical properties—namely H _r , indicating size; colors, indicating surface composition; and flux variation semiamplitude A , indicating asphericity and surface inhomogeneity. We seek “birth” physical distributions that can recreate these parameters in every dynamical class. We show that the observed colors of these TNOs are consistent with two Gaussian distributions in griz space, “near-infrared bright” (NIRB) and “near-infrared faint” (NIRF), presumably an inner and outer birth population, respectively. We find a model in which both the NIRB and NIRF H _r and A distributions are independent of current dynamical states, supporting their assignment as birth populations. All objects are consistent with a common rolling p ( H _r ), but NIRF objects are significantly more variable. Cold classicals (CCs) are purely NIRF, while hot classical (HC), scattered, and detached TNOs are consistent with ≈ 70% NIRB and the resonance NIRB fractions show significant variation. The NIRB components of the HCs and of some resonances have broader inclination distributions than the NIRFs, i.e. their current dynamics retains information about birth location. We find evidence for radial stratification within the birth NIRB population, in that HC NIRBs are on average redder than detached or scattered NIRBs; a similar effect distinguishes CCs from other NIRFs. We estimate total object counts and masses of each class within our H _r range. These results will strongly constrain models of the outer solar system.