Researcher profile

Nicolas F. Martin

Nicolas F. Martin contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
18works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2023arXiv

The Pristine survey -- XX: GTC follow-up observations of extremely metal-poor stars identified from Pristine and LAMOST

Ultra metal-poor stars ([Fe/H] < -4.0) are very rare, and finding them is a challenging task. Both narrow-band photometry and low-resolution spectroscopy have been useful tools for identifying candidates, and in this work we combine both approaches. We cross-matched metallicity-sensitive photometry from the Pristine survey with the low-resolution spectroscopic LAMOST database, and re-analysed all LAMOST spectra with [Fe/H]_Pristine < -2.5. We find that ~1/3rd of this sample (selected without [Fe/H]_Pristine quality cuts) also have spectroscopic [Fe/H] < -2.5. From this sample, containing many low signal-to-noise (S/N) spectra, we selected eleven stars potentially having [Fe/H] < -4.0 or [Fe/H] < -3.0 with very high carbon abundances, and we performed higher S/N medium-resolution spectroscopic follow-up with OSIRIS on the 10.4m Gran Telescopio Canarias (GTC). We confirm their extremely low metallicities, with a mean of [Fe/H] = -3.4 and the most metal-poor star having [Fe/H]= -3.8. Three of these are clearly carbon-enhanced metal-poor (CEMP) stars with +1.65 < [C/Fe] < +2.45. The two most carbon-rich stars are either among the most metal-poor CEMP-s stars or the most carbon-rich CEMP-no stars known, the third is likely a CEMP-no star. We derived orbital properties for the OSIRIS sample and find that only one of our targets can be confidently associated with known substructures/accretion events, and that three out of four inner halo stars have prograde orbits. Large spectroscopic surveys may contain many hidden extremely and ultra metal-poor stars, and adding additional information from e.g. photometry as in this work can uncover them more efficiently and confidently.

preprint2022arXiv

A stellar stream remnant of a globular cluster below the metallicity floor

Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo.

preprint2022arXiv

On the inconsistency of [C/Fe] abundances and the fractions of carbon-enhanced metal-poor stars among various stellar surveys

Carbon-enhanced metal-poor (CEMP) stars are a unique resource for Galactic archaeology because they probe the properties of the First Stars, early chemical evolution and binary interactions at very low metallicity. Comparing the fractions and properties of CEMP stars in different Galactic environments can provide us with unique insights into the formation and evolution of the Milky Way halo and its building blocks. In this work, we investigate whether directly comparing fractions of CEMP stars from different literature samples of very metal-poor ([Fe/H] < -2.0) stars is valid. We compiled published CEMP fractions and samples of Galactic halo stars from the past 25 years, and find that they are not all consistent with each other. Focusing on giant stars, we find significant differences between various surveys when comparing their trends of [Fe/H] versus [C/Fe] and their distributions of CEMP stars. To test the role of the analysis pipelines for low-resolution spectroscopic samples, we re-analysed giant stars from various surveys with the SSPP and FERRE pipelines. We found systematic differences in [C/Fe] of ~0.1-0.4 dex, partly independent of degeneracies with the stellar atmospheric parameters. These systematics are likely due to the different pipeline approaches, different assumptions in the employed synthetic grids, and/or the comparison of different evolutionary phases. We conclude that current biases in (the analysis of) very metal-poor samples limit the conclusions one can draw from comparing different surveys. We provide some recommendations and suggestions that will hopefully aid the community to unlock the full potential of CEMP stars for Galactic archaeology.

preprint2022arXiv

Optimizing Nitrogen Management with Deep Reinforcement Learning and Crop Simulations

Nitrogen (N) management is critical to sustain soil fertility and crop production while minimizing the negative environmental impact, but is challenging to optimize. This paper proposes an intelligent N management system using deep reinforcement learning (RL) and crop simulations with Decision Support System for Agrotechnology Transfer (DSSAT). We first formulate the N management problem as an RL problem. We then train management policies with deep Q-network and soft actor-critic algorithms, and the Gym-DSSAT interface that allows for daily interactions between the simulated crop environment and RL agents. According to the experiments on the maize crop in both Iowa and Florida in the US, our RL-trained policies outperform previous empirical methods by achieving higher or similar yield while using less fertilizers

preprint2022arXiv

Stellar kinematics of dwarf galaxies from multi-epoch spectroscopy: application to Triangulum II

We present new MMT/Hectochelle spectroscopic measurements for 257 stars observed along the line of sight to the ultra-faint dwarf galaxy Triangulum II. Combining with results from previous Keck/DEIMOS spectroscopy, we obtain a sample that includes 16 likely members of Triangulum II, with up to 10 independent redshift measurements per star. To this multi-epoch kinematic data set we apply methodology that we develop in order to infer binary orbital parameters from sparsely sampled radial velocity curves with as few as two epochs. For a previously-identified (spatially unresolved) binary system in Tri~II, we infer an orbital solution with period $296.0_{-3.3}^{+3.8} \rm~ days$ , semi-major axis $1.12^{+0.41}_{-0.24}\rm~AU$, and a systemic velocity $ -380.0 \pm 1.7 \rm~km ~s^{-1}$ that we then use in the analysis of Tri~II&#39;s internal kinematics. Despite this improvement in the modeling of binary star systems, the current data remain insufficient to resolve the velocity dispersion of Triangulum II. We instead find a 95% confidence upper limit of $σ_{v} \lesssim 3.4 \rm ~km~s^{-1}$.

preprint2022arXiv

The Complexity of the Cetus Stream Unveiled from the Fusion of STREAMFINDER and StarGO

We combine the power of two stream-searching tools, STREAMFINDER and StarGO applied to the Gaia EDR3 data, to detect stellar debris belonging to the Cetus stream system that forms a complex, nearly polar structure around the Milky Way. In this work, we find the southern extensions of the northern Cetus stream as the Palca stream and a new southern stream, which overlap on the sky but have different distances. These two stream wraps extend over more than $\sim100°$ on the sky ($-60°<δ<+40°$). The current N-body model of the system reproduces both wraps in the trailing arm. We also show that the Cetus system is confidently associated with the Triangulum/Pisces, Willka Yaku, and the recently discovered C-20 streams. The association with the ATLAS-Aliqa Uma stream is much weaker. All of these stellar debris are very metal-poor, comparable to the average metallicity of the southern Cetus stream with [Fe/H] = $-2.17\pm0.2$. The estimated stellar mass of the Cetus progenitor is at least $10^{5.6}\rm\,M_\odot$, compatible with Ursa Minor or Draco dwarf galaxies. The associated globular cluster with similar stellar mass, NGC 5824 very possibly was accreted in the same group infall. The multi-wrap Cetus stream is a perfect example of a dwarf galaxy that has undergone several periods of stripping, leaving behind debris at multiple locations in the halo. The full characterization of such systems is crucial to unravel the assembly history of the Milky Way and, as importantly, to provide nearby fossils to study ancient low-mass dwarf galaxies.

preprint2022arXiv

The Global Dynamical Atlas of the Milky Way mergers: Constraints from Gaia EDR3 based orbits of globular clusters, stellar streams and satellite galaxies

The Milky Way halo was predominantly formed by the merging of numerous progenitor galaxies. However, our knowledge of this process is still incomplete, especially in regard to the total number of mergers, their global dynamical properties and their contribution to the stellar population of the Galactic halo. Here, we uncover the Milky Way mergers by detecting groupings of globular clusters, stellar streams and satellite galaxies in action ($\mathbf{J}$) space. While actions fully characterize the orbits, we additionally use the redundant information on their energy ($\textit{E}$) to enhance the contrast between groupings. For this endeavour, we use $\textit{Gaia}$ EDR3 based measurements of $170$ globular clusters, $41$ streams and $46$ satellites to derive their $\mathbf{J}$ and $\textit{E}$. To detect groups, we use the $\texttt{ENLINK}$ software, coupled with a statistical procedure that accounts for the observed phase-space uncertainties of these objects. We detect a total of $N=6$ groups, including the previously known mergers $\textit{Sagittarius}$, $\textit{Cetus}$, $\textit{Gaia-Sausage/Enceladus}$, $\textit{LMS-1/Wukong}$, $\textit{Arjuna/Sequoia/I&#39;itoi}$ and one new merger that we call $\textit{Pontus}$. All of these mergers, together, comprise $62$ objects ($\approx 25\%$ of our sample). We discuss their members, orbital properties and metallicity distributions. We find that the three most metal-poor streams of our Galaxy -- &#34;C-19&#34; ([Fe/H]$=-3.4$ dex), &#34;Sylgr&#34; ([Fe/H]$=-2.9$ dex) and &#34;Phoenix&#34; ([Fe/H]$=-2.7$ dex) -- are associated with $\textit{LMS-1/Wukong}$; showing it to be the most metal-poor merger. The global dynamical atlas of Milky Way mergers that we present here provides a present-day reference for galaxy formation models.

preprint2022arXiv

The PAndAS View of the Andromeda Satellite System. III. Dwarf galaxy detection limits

We determine the detection limits of the search for dwarf galaxies in the Pan-Andromeda Archaeological Survey (PAndAS) using the algorithm developed by the PAndAS team. The recovery fractions of artificial dwarf galaxies are, as expected, a strong function of physical size and luminosity and, to a lesser extent, distance. We show that these recovery fractions vary strongly with location in the surveyed area because of varying levels of contamination from both the Milky Way foreground stars and the stellar halo of Andromeda. We therefore provide recovery fractions that are a function of size, luminosity, and location within the survey on a scale of 1 square degree. Overall, the effective surface brightness for a 50-percent detection rate range between 28 and 30 mag per square arcsecond. This is in line with expectations for a search that relies on photometric data that are as deep as the PAndAS survey. The derived detection limits are an essential ingredient on the path to constraining the global properties of Andromeda&#39;s system of satellite dwarf galaxies and, more broadly, to provide constraints on dwarf galaxy formation and evolution in a cosmological context.

preprint2022arXiv

The Pristine survey -- XVII. The C-19 stream is dynamically hot and more extended than previously thought

The C-19 stream is the most metal poor stellar system ever discovered, with a mean metallicity $[Fe/H] = -3.38\pm0.06$. Its low metallicity dispersion ($σ_{\rm [Fe/H]}$ $<$ 0.18 at the 95\% confidence level) as well as variations in sodium abundances strongly suggest a globular cluster origin. In this work, we use VLT/UVES spectra of seven C-19 stars to derive more precise velocity measurements for member stars, and to identify two new members with radial velocities and metallicities consistent with the stream&#39;s properties. One of these new member stars is located 30 $°$ away from the previously identified body of C-19, implying that the stream is significantly more extended than previously known and that more members likely await discovery. In the main part of C-19, we measure a radial velocity dispersion $σ_v$ = 6.2$^{+2.0}_{-1.4}$ km s$^{-1}$ from nine members, and a stream width of 0.56$°\pm0.08°$, equivalent to $\sim$158 pc at a heliocentric distance of 18 kpc. These confirm that C-19 is comparatively hotter, dynamically, than other known globular cluster streams and shares the properties of faint dwarf galaxy streams. On the other hand, the variations in the Na abundances of the three newly observed bright member stars, the variations in Mg and Al for two of them, and the normal Ba abundance of the one star where it can be measured provide further evidence for a globular cluster origin. The tension between the dynamical and chemical properties of C-19 suggests that its progenitor experienced a complex birth environment or disruption history.

preprint2020arXiv

SMASHing the low surface brightness SMC

The periphery of the Small Magellanic Cloud (SMC) can unlock important information regarding galaxy formation and evolution in interacting systems. Here, we present a detailed study of the extended stellar structure of the SMC using deep colour-magnitude diagrams (CMDs), obtained as part of the Survey of the MAgellanic Stellar History (SMASH). Special care was taken in the decontamination of our data from MW foreground stars, including from foreground globular clusters NGC 362 and 47 Tuc. We derived the SMC surface brightness using a ``conservative&#39;&#39; approach from which we calculated the general parameters of the SMC, finding a staggered surface brightness profile. We also traced the fainter outskirts by constructing a stellar density profile. This approach, based on stellar counts of the oldest main sequence turn-off (MSTO) stars, uncovered a tidally disrupted stellar feature that reaches as far out as 12 degrees from the SMC centre. We also serendipitously found a faint feature of unknown origin located at $\sim 14$ degrees from the centre of the SMC and that we tentatively associated to a more distant structure. We compared our results to in-house simulations of a $1\times10^{9} M_\odot$ SMC, finding that its elliptical shape can be explained by its tidal disruption under the combined presence of the MW and the LMC. Finally, we found that the older stellar populations show a smooth profile while the younger component presents a jump in the density followed by a flat profile, confirming the heavily disturbed nature of the SMC.

preprint2020arXiv

The Pristine Inner Galaxy Survey (PIGS) II: Uncovering the most metal-poor populations in the inner Milky Way

Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy.

preprint2020arXiv

The Second Data Release of the Survey of the MAgellanic Stellar History (SMASH)

The Large and Small Magellanic Clouds (LMC and SMC) are the largest satellite galaxies of the Milky Way and close enough to allow for a detailed exploration of their structure and formation history. The Survey of the MAgellanic Stellar History (SMASH) is a community Dark Energy Camera (DECam) survey of the Magellanic Clouds using $\sim$50 nights to sample over $\sim$2400 deg$^2$ centered on the Clouds at $\sim$20% filling factor (but with contiguous coverage in the central regions) and to depths of $\sim$24th mag in $ugriz$. The primary goals of SMASH are to map out the extended stellar peripheries of the Clouds and uncover their complicated interaction and accretion history as well as to derive spatially-resolved star formation histories of the central regions and create a &#34;movie&#34; of their past star formation. Here we announce the second SMASH public data release (DR2), which contains all 197 fully-calibrated DECam fields including the main body fields in the central regions. The DR2 data are available through the Astro Data Lab hosted by the NSF&#39;s National Optical-Infrared Astronomy Research Laboratory. We highlight three science cases that make use of the SMASH DR2 data and will be published in the future: (1) preliminary star formation histories of the LMC; (2) the search for Magellanic star clusters using citizen scientists; and, (3) photometric metallicities of Magellanic Cloud stars using the DECam $u$-band.

preprint2019arXiv

A detailed study of Andromeda XIX, an extreme local analogue of ultra diffuse galaxies

With a central surface brightness of $μ_0=29.3$ mag. per sq. arcsec, and half-light radius of $r_{\rm half}=3.1^{+0.9}_{-1.1}$~kpc, Andromeda XIX (And XIX) is an extremely diffuse satellite of Andromeda. We present spectra for $\sim100$ red giant branch stars in this galaxy, plus 16 stars in a nearby stellar stream. With this exquisite dataset, we re-derive the properties of And XIX, measuring a systemic velocity of $v_r=-109.0\pm1.6$ km/s and a velocity dispersion of $σ_v = 7.8^{+1.7}_{-1.5}$ km/s (higher than derived in our previous work). We marginally detect a velocity gradient along the major axis of ${\rm d}v/{\rm d}r = -2.1\pm1.8$ km/s kpc$^{-1}$. We find its mass-to-light ratio is higher than galaxies of comparable stellar mass ($[M/L]_{\rm half} = 278^{+146}_{-198}M_\odot/L_\odot$), but its dynamics place it in a halo with a similar total mass to these galaxies. This could suggest that And XIX is a &#34;puffed up&#34; dwarf galaxy, whose properties have been altered by tidal processes, similar to its Milky Way counterpart, Antlia II. For the nearby stream, we measure $v_r=-279.2\pm3.7$ km/s, and $σ_v=13.8^{+3.5}_{-2.6}$ km/s. We measure its metallicity, and find it to be more metal rich than And XIX, implying that the two features are unrelated. Finally, And XIX&#39;s dynamical and structural properties imply it is a local analogue to ultra diffuse galaxies (UDGs). Its complex dynamics suggest that the masses of distant UDGs measured from velocity dispersions alone should be carefully interpreted

preprint2019arXiv

Comparing the Quenching Times of Faint M31 and Milky Way Satellite Galaxies

We present the star formation histories (SFHs) of 20 faint M31 satellites ($-12 \lesssim M_V \lesssim -6$) that were measured by modeling sub-horizontal branch (HB) depth color-magnitude diagrams constructed from Hubble Space Telescope (HST) imaging. Reinforcing previous results, we find that virtually all galaxies quenched between 3 and 9 Gyr ago, independent of luminosity, with a notable concentration $3-6$ Gyr ago. This is in contrast to the Milky Way (MW) satellites, which are generally either faint with ancient quenching times or luminous with recent ($<3$ Gyr) quenching times. We suggest that systematic differences in the quenching times of M31 and MW satellites may be a reflection of the varying accretion histories of M31 and the MW. This result implies that the formation histories of low-mass satellites may not be broadly representative of low-mass galaxies in general. Among the M31 satellite population we identify two distinct groups based on their SFHs: one with exponentially declining SFHs ($τ\sim 2$ Gyr) and one with rising SFHs with abrupt quenching. We speculate how these two groups could be related to scenarios for a recent major merger involving M31. The Cycle 27 HST Treasury survey of M31 satellites will provide well-constrained ancient SFHs to go along with the quenching times we measure here. The discovery and characterization of M31 satellites with $M_V \gtrsim -6$ would help quantify the relative contributions of reionization and environment to quenching of the lowest-mass satellites.

preprint2019arXiv

Phase-space correlation in stellar streams of the Milky Way halo: The clash of Kshir and GD-1

We report the discovery of a 70-degree long stellar stream in the Milky Way halo, which criss-crosses the well known &#34;GD-1&#34; stream. We show that this new stellar structure (&#34;Kshir&#34;) and GD-1 lie at similar distance, and are remarkably correlated in kinematics. We propose several explanations for the nature of this new structure and its possible association with GD-1. However, a scenario in which these two streams were accreted onto the Milky Way within the same dark matter sub-halo seems to provide a natural explanation for their phase-space entanglement, and other complexities of this coupled-system.

preprint2019arXiv

The Canada-France Imaging Survey: Reconstructing the Milky Way Star Formation History from its White Dwarf Population

As the remnants of stars with initial masses $\lesssim$ 8 M$_{\odot}$, white dwarfs contain valuable information on the formation histories of stellar populations. In this paper, we use deep, high-quality, u-band photometry from the Canada France Imaging Survey (CFIS), griz photometry from Pan-STARRS 1 (PS1), as well as proper motions from Gaia DR2, to select 25,156 white dwarf candidates over $\sim$4500 deg$^2$ using a reduced proper motion diagram. We develop a new white dwarf population synthesis code that returns mock observations of the Galactic field white dwarf population for a given star formation history, while simultaneously taking into account the geometry of the Milky Way, survey parameters, and selection effects. We use this model to derive the star formation histories of the thin disk, thick disk, and stellar halo. Our results show that the Milky Way disk began forming stars (11.3 $\pm$ 0.5) Gyr ago, with a peak rate of (8.8 $\pm$ 1.4) M$_{\odot}$yr$^{-1}$ at (9.8 $\pm$ 0.4) Gyr, before a slow decline to a constant rate until the present day --- consistent with recent results suggesting a merging event with a satellite galaxy. Studying the residuals between the data and best-fit model shows evidence for a slight increase in star formation over the past 3 Gyr. We fit the local fraction of helium-atmosphere white dwarfs to be (21 $\pm$ 3) %. Incorporating this methodology with data from future wide-field surveys such as LSST, Euclid, CASTOR, and WFIRST should provide an unprecedented view into the formation of the Milky Way at its earliest epoch through its white dwarfs.

preprint2019arXiv

The Pristine survey X: a large population of low-metallicity stars permeates the Galactic disk

The orbits of the least chemically enriched stars open a window on the formation of our Galaxy when it was still in its infancy. The common picture is that these low-metallicity stars are distributed as an isotropic, pressure-supported component since these stars were either accreted from the early building blocks of the assembling Milky Way, or were later brought by the accretion of faint dwarf galaxies. Combining the metallicities and radial velocities from the Pristine and LAMOST surveys and Gaia DR2 parallaxes and proper motions for an unprecedented large and unbiased sample of very metal-poor stars at $[Fe/H]\leq-2.5$ we show that this picture is incomplete. This sample shows strong statistical evidence (at the $5.0σ$ level) of asymmetry in their kinematics, favouring prograde motion. Moreover, we find that $31\%$ of the stars that currently reside in the disk do not venture outside of the disk plane throughout their orbit. The discovery of this population implies that a significant fraction of stars with iron abundances $[Fe/H]\leq-2.5$ formed within or concurrently with the Milky Way disk and that the history of the disk was quiet enough to allow them to retain their disk-like orbital properties.

preprint2010arXiv

Mapping the stellar structure of the Milky Way thick disk and halo using SEGUE photometry

We map the stellar structure of the Galactic thick disk and halo by applying color-magnitude diagram (CMD) fitting to photometric data from the SEGUE survey, allowing, for the first time, a comprehensive analysis of their structure at both high and low latitudes using uniform SDSS photometry. Incorporating photometry of all relevant stars simultaneously, CMD fitting bypasses the need to choose single tracer populations. Using old stellar populations of differing metallicities as templates we obtain a sparse 3D map of the stellar mass distribution at |Z|>1 kpc. Fitting a smooth Milky Way model comprising exponential thin and thick disks and an axisymmetric power-law halo allows us to constrain the structural parameters of the thick disk and halo. The thick-disk scale height and length are well constrained at 0.75+-0.07 kpc and 4.1+-0.4 kpc, respectively. We find a stellar halo flattening within ~25 kpc of c/a=0.88+-0.03 and a power-law index of 2.75+-0.07 (for 7<R_{GC}<~30 kpc). The model fits yield thick-disk and stellar halo densities at the solar location of rho_{thick,sun}=10^{-2.3+-0.1} M_sun pc^{-3} and rho_{halo,sun}=10^{-4.20+-0.05} M_sun pc^{-3}, averaging over any substructures. Our analysis provides the first clear in situ evidence for a radial metallicity gradient in the Milky Way&#39;s stellar halo: within R<~15 kpc the stellar halo has a mean metallicity of [Fe/H]=-1.6, which shifts to [Fe/H]=-2.2 at larger radii. Subtraction of the best-fit smooth and symmetric model from the overall density maps reveals a wealth of substructures at all latitudes, some attributable to known streams and overdensities, and some new. A simple warp cannot account for the low latitude substructure, as overdensities occur simultaneously above and below the Galactic plane. (abridged)