This is an open-access article distributed under the terms of the Creative Commons Attribution CC-BY 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The Mastcam CCD cameras on the Mars Science Laboratory Curiosity Rover each use an 8-position filter wheel in acquiring up to 1600 × 1200 pixel images. The filter set includes a broadband near-infrared cutoff filter for RGB Bayer imaging on each camera and 12 narrow-band geology filters distributed between the two cameras, spanning the wavelength range 445–1013 nm. This wavelength region includes the relatively broad charge-transfer and crystal-field absorption bands that are most commonly due to the presence of iron-bearing minerals. To identify such spectral features, sequences of images taken with identical pointings through different filters have been calibrated to relative reflectance using pre-flight calibration coefficients and in-flight measurements of an onboard calibration target. Within the first 1000 sols of the mission, Mastcam observed a spectrally diverse set of materials displaying absorption features consistent with the presence of iron-bearing silicate, iron oxide, and iron sulfate minerals. Dust-coated surfaces as well as soils possess a strong positive reflectance slope in the visible, consistent with the presence of nanophase iron oxides, which have long been considered the dominant visible-wavelength pigmenting agent in weathered martian surface materials. Fresh surfaces, such as tailings produced by the drill tool and the interiors of rocks broken by the rover wheels, are grayer in visible wavelengths than their reddish, dust-coated surfaces but possess reflectance spectra that vary considerably between sites. To understand the mineralogical basis of observed Mastcam reflectance spectra, we focus on a subset of the multispectral data set for which additional constraints on the composition of surface materials are available from other rover instruments, with an emphasis on sample sites for which detailed mineralogy is provided by the results of CheMin X-ray diffraction analyses. We also discuss the results of coordinated observations with the ChemCam instrument, whose passive mode of operation is capable of acquiring reflectance spectra over wavelengths that considerably overlap the range spanned by the Mastcam filter set (Johnson et al. 2016). Materials that show a distinct 430 nm band in ChemCam data also are observed to have a strong near-infrared absorption band in Mastcam spectral data, consistent with the presence of a ferric sulfate mineral. Long-distance Mastcam observations targeted toward the flanks of the Gale crater central mound are in agreement with both ChemCam spectra and orbital results, and in particular exhibit the spectral features of a crystalline hematite layer identified in MRO/CRISM data. Variations observed in Mastcam multi-filter images acquired to date have shown that multispectral observations can discriminate between compositionally different materials within Gale Crater and are in qualitative agreement with mineralogies from measured samples and orbital data.

Introduction

The Mars Science Laboratory (MSL) Curiosity Rover has been exploring the environment of Gale Crater since its successful landing in August 2012. The scientific payload was selected to investigate the potential for past habitable environments through a detailed investigation of the layered sedimentary units of the central mound, informally known as Mt. Sharp, and the surrounding plains (Grotzinger et al. 2012). A history of aqueous alteration within the crater is indicated by spectral and geomorphologic evidence identified from orbit (e.g., Anderson and Bell 2010; Milliken et al. 2010; Thomson et al. 2011) and is hypothesized to represent a global transition from clay- to sulfate-dominated alteration mineralogy that took place early in the planet's history (e.g., Bibring et al. 2006; Milliken et al. 2010). Multispectral imaging capability on Curiosity is provided by the Mast Camera (Mastcam) instrument suite (Malin et al. 2010, 2016 in review1; Bell et al. 2012, 2016 in review1), which comprises two 1600 × 1200 pixel Bayer-patterned CCD cameras located ~2 m above the surface on the rover's remote sensing mast, along with an accompanying calibration target mounted on the rover deck. Each camera is equipped with an 8-position filter wheel designed to characterize the visible to near-infrared reflectance spectra of surface materials at up to 12 unique wavelengths from 445–1013 nm, including broadband imaging over Bayer filter red, green, and blue (RGB) wavelengths. This wavelength range includes the positions of numerous absorption features of both primary iron-bearing basaltic minerals as well as certain iron-bearing alteration products (e.g., Hunt et al. 1974; Hunt and Ashley 1979; Burns 1993; Clark et al. 2007).

Multispectral observations supplement information on morphology and stratigraphic relationships provided by broadband RGB stereo Mastcam or single-band engineering stereo camera images, which together provide geologic context for other instruments. The importance of this contextual information cannot be overstated: many of the scientific instruments aboard the rover perform measurements with relatively small spot sizes (APXS, MAHLI, ChemCam) or require material to be transferred internally via the Sample Acquisition, Processing, and Handling (SA/SPaH) subsystem (CheMin, SAM) (Grotzinger et al. 2012). Mastcam multispectral observations can document compositional variability within and across geologic units on a broad scale, including distant targets that are out of range to most other rover instruments, to the extent that such variability is observable by the filter set. In the near field they can serve as reconnaissance to identify spectrally distinct materials for follow-up analyses by other instruments. Operational and data volume constraints limit the number of multispectral survey sequences that can be acquired, however, and therefore the development of imaging strategies that make use of a subset of the full pixel array, or of the full filter set, is an ongoing consideration.

As was the case for similar limited-filter visible to near-infrared multispectral imaging on the prior Mars Pathfinder, Mars Exploration Rover, and Phoenix Lander missions (Bell et al. 2000, 2004a, 2004b; Farrand et al. 2007, 2008, 2016; Blaney et al. 2009), even complete Mastcam 12-point spectra are generally not sufficiently diagnostic to provide unique mineralogical interpretations, at least not without separate supporting information. Fortunately, analyses by other on-board instruments can help to constrain the interpretation of Mastcam spectral features. For this reason, we focus in this paper on a subset of multispectral observations acquired in conjunction with multiple other instruments to better understand the mineralogy underlying the observed spectral characteristics. Of most relevance to the interpretation of Mastcam multispectral data are elemental and mineralogic analyses by ChemCam (Maurice et al. 2012; Wiens et al. 2012), elemental analyses by the Alpha Particle X-ray Spectrometer (APXS; Gellert et al. 2009), and mineralogic analyses by the CheMin X-ray diffraction (XRD) instrument (Blake et al. 2012). In the case of the ChemCam instrument, Johnson et al. (2015) have shown that ChemCam passive observations (acquired when the LIBS laser is not active) can be used to generate relative reflectance spectra in the 400–840 nm wavelength range. This range overlaps significantly with that of the Mastcam filter set, and hence inter-comparisons between the two data sets can provide an important check on the identification of spectral features within the range of overlap (as well as important cross-calibration information for both instrument investigations). Many of the observations presented herein are also the subject of ChemCam passive observations, and detailed interpretations of those data are presented by Johnson et al. (2016).

Because the most diagnostic information on mineralogy is provided by the CheMin XRD instrument, the set of Mastcam spectra presented here primarily focuses on multispectral observations of soil scoops and drill fines that have also been examined by CheMin. While drilling has been the predominant method of sample acquisition, several scoops of soil were processed early in the mission at the Rocknest location (Blake et al. 2013; Anderson et al. 2015b), uncovering fresh material for multispectral analysis in an aeolian sand ripple. The drill, part of the Powder Acquisition Drill System (PADS) portion of the SA/SPaH subsystem, produces a 1.6 cm diameter hole in the surface of up to approximately 5 cm depth (Anderson et al. 2012). The upper ~1.5 cm of material is not collected by the sampling subsystem but is instead distributed on the surface as a tailings pile. Multispectral image sequences of these drill fines allow Mastcam to observe surfaces largely uncontaminated by reddish, Fe3+-bearing dust, which influences the visible to near-infrared spectra of practically all surfaces to a varying degree [including, to a lesser extent, surfaces brushed by the Dust Removal Tool (DRT)]. In addition, drill fines are typically subject to analyses by ChemCam and APXS as well, providing a comparatively comprehensive multiple-instrument data set for these targets.

Background

The Mastcam left (M-34, 34 mm focal length) and right (M-100, 100 mm focal length, thus ~3× better spatial resolution than the M-34) cameras each possess a filter wheel holding eight different optical filters for multispectral imaging. One position on each camera is occupied by a broadband infrared-cutoff filter (“filter zero”) for RGB color imaging, making use of the 2 × 2 unit-cell Bayer pattern bonded directly to the detector to acquire three broadband visible wavelength channels. A second filter slot is taken by a narrow-band, 10-5 neutral-density-coated filter designed for direct solar imaging. The remaining 12 filter positions are occupied by narrow-band filters selected to characterize the visible to near-infrared reflectance spectra of rock, soil, and other targets in the wavelength range of 445–1013 nm (see Table 1 and Fig. 1 for filter bandpass characteristics). The combination of filters between the two cameras provides 12 unique (differing by more than a few nanometers; Bell et al. 2016 in review1) center wavelengths for multispectral analysis, including the three RGB Bayer bands. Owing to the presence of the Bayer pattern, filters below 850 nm have different throughput between nonequivalent Bayer pixels (Malin et al. 2010; Bell et al. 2012, 2016 in review1). When lossy JPEG compression is used for downlink of observations made with these narrow-band filters (which it is for most cases), the flight software produces the full-size image by bilinear interpolation from the Bayer color pixel with a wavelength response closest to that of the filter itself, discarding the other Bayer pixels [see also the Software Interface Specification (SIS) document for the instrument, Malin et al. (2013)]. For filters L2 and R2 the Bayer blue pixels are used, for L1 and R1 the greens, and for L3, L4, and R3 the red pixels, effectively decreasing the spatial resolution at these wavelengths by a factor of about 1.4 (for L1/R1) or 2 (for the other short-wavelength filters). The broadband filters L0 and R0 are demosaicked by the algorithm of Malvar et al. (2004) to produce the individual red, green, and blue color images.

Spectral features in the wavelength range of the camera are predominantly due to the crystal-field and charge-transfer absorptions of iron-bearing minerals (e.g., Burns 1993), while most vibrational features lie beyond the sensitivity range of the cameras' silicon CCD detectors. One exception is a narrow H2O vibrational overtone/combination band that, in certain hydrated minerals, coincides approximately with the longest wavelength Mastcam filter. Detection of this feature has been reported in Mars Exploration Rover (MER) Pancam observations (e.g., Wang et al. 2008; Rice et al. 2010) and its detection in Mast-cam spectra is being explored by Rice et al. (2013a, 2013b). Broad electronic absorption bands that have been identified on Mars in this wavelength region are attributed to the presence of iron-bearing silicate, oxide, and sulfate minerals (e.g., Bell 1996, 2008). Ferrous silicates such as iron-bearing pyroxenes and olivines possess an absorption band near 1000 nm, a result of a spin-allowed transition of Fe2+ in octahedral coordination (e.g., Hunt 1977). Reflectance spectra of pyroxenes vary systematically with composition (e.g., Adams 1974), with the 900–1000 nm absorption band tending to shorter wavelengths with lower calcium content. Ferric minerals generally have several crystal-field bands in this range, as well as a strong charge-transfer band extending from the ultraviolet into the visible (e.g., Hunt and Ashley 1979). As an example, hematite has a distinct band near 860 nm, a shoulder near 630 nm, and an intense absorption wing extending from the ultraviolet into the visible to about 550 nm, formed from overlapping crystal-field and charge-transfer bands (Morris et al. 1985). Nanophase hematite, an X-ray amorphous material in which the particle size is less than approximately 10 nm, lacks distinct crystal-field bands but still possesses a strong iron-oxygen charge transfer absorption edge through the visible wavelengths (Morris et al. 1989). Such nanophase oxides are believed to be primarily responsible for the reddish color of martian dust and soil (Morris et al. 1993, 1997). Representative library spectra of several of these and other iron-bearing minerals that will be mentioned below are shown in Figure 2, along with their values as convolved to the Mastcam filter bandpasses.

Variability in Mastcam spectral data acquired within Gale Crater was anticipated on the basis of orbital observations. Specifically, previous analyses of orbital spectral data have shown evidence for the existence of nontronite, magnesium sulfates, and crystalline hematite in the lowermost layers of Mt. Sharp (Milliken et al. 2010; Thomson et al. 2011). A hematite-rich layer near the base of the mound has been mapped to the uppermost stratum of an erosion-resistant ridge (Fraeman et al. 2013) and should be easily identifiable by Mastcam, whose 867 nm filter is located near the center of a crystalline hematite absorption band. Nontronite possesses features similar to other Fe3+-rich minerals, but is distinct from hematite by the longer wavelength position of its broad near-infrared band, centered around 950 nm (e.g., Singer 1982; Bishop et al. 2008). Magnesium sulfates, on the other hand, lack broad absorption features at Mastcam wavelengths (although hydrous varieties possess an H2O vibrational band near 1000 nm (e.g., Drake 1995) that may be detectable by the cameras). Iron-bearing varieties of sulfate were not detected from orbit, although this does not preclude the presence of minor or small-scale occurrences that may exist below the detection limit or spatial resolution of orbital instruments. Indeed ferric sulfates have been identified previously at other locations on Mars, both in situ at Meridiani Planum (Christensen et al. 2004; Klingelhöfer et al. 2004) and Gusev Crater (Arvidson et al. 2006; Johnson et al. 2007; Lane et al. 2008) and from orbit at multiple other sites (e.g., Milliken et al. 2008; Bishop et al. 2009; Farrand et al. 2009). Their presence as a minor component of certain bedrock units within Gale Crater has recently been confirmed by the CheMin instrument (Cavanagh et al. 2015; Rampe et al. 2016).

Methodology

The conversion of raw Data Number (DN) values of Mastcam multispectral observations to meaningful radiometric quantities involves the use of both pre-flight calibration measurements as well as near-in-time imaging of the onboard Mastcam calibration target during data acquisition (Bell et al. 2006, 2016 in review1). The calibration pipeline is described in detail by Bell et al. (2016 in review1), but a brief summary is presented here. Raw observations in the form of Experimental Data Records (EDRs) are decompanded from downlinked 8-bit data back to their original 11-bit dynamic range. For most observations, interpolation over unused Bayer pixels in specific shorter wavelength filters is a step handled by the on-board software, as mentioned in the section above; however, for certain observations in which the full pixel array is returned, this step must be done by the user. In these cases we follow the same debayering method as the flight software. Observations are flat-fielded using normalized relative responsivity arrays derived from sky observations acquired on Curiosity mission sols 36–38. Bias and shutter-smear corrections are not yet implemented but are insignificant components of the measured signal for all but the most extreme observational circumstances (not the case for the observations described here). The average dark current, as measured from masked regions of the detector array, was found to be negligible at the operating temperatures and exposure times of the data reported here, and accordingly these pixel columns are excluded by subframing in most observations. Pixels that possess an 8-bit value of 240 or larger in the raw data are considered saturated on the basis of both pre-flight and in-flight observations, and are ignored as “missing data” in subsequent downstream processing. Flatfielded 11-bit DN values are converted to radiance (W/m2/nm/sr) using pre-flight observations of a NIST-calibrated integrating sphere and monochrometer measurements of individual filter bandpasses.

To convert from radiance to radiance factor (unitless I/F, where I is the measured radiance and πF is the incident solar irradiance), multispectral observations are typically immediately followed or preceded by a sequence of images of the on-board calibration target acquired using the same filter set. In cases where they are not, a calibration target observation from another sol imaged within an hour of the appropriate time of day is used (see Table 2). The Mastcam calibration target is located on the right side of the rover deck approximately 1.2 m from the front windows of the cameras, on top of the Rover Pyro Firing Assembly control box. This calibration target or “caltarget,” which is nearly identical to the Pancam caltarget (Bell et al. 2003, 2006), consists of a ball-and-stick central post (gnomon) surrounded by three grayscale rings and four color chips (Fig. 3). Unlike the MER design, the Mastcam calibration target includes six cylindrical magnets embedded just underneath the surface of the color chips and white and gray rings, which keep the center of the magnet regions comparatively clean while attracting a surrounding ring of magnetic dust (cf., Madsen et al. 2003; Goetz et al. 2008). These “clean spots” provide additional information for assessing the performance of the calibration procedure but are not directly utilized in the calibration pipeline.

Radiance values are extracted from region-of-interest (ROI) selections drawn on the grayscale calibration target rings (avoiding the rings of concentrated dust near the magnet locations). These values are plotted against the laboratory-measured directional-hemispherical reflectance values for the caltarget materials, corrected for illumination and viewing geometry by means of a modified He-Torrance model (He et al. 1991; Bell et al. 2003), and corrected for dust deposition (discussed below). The bidirectional reflectance distribution function model, developed on the basis of prior measurements of the caltarget substrate materials at MER/Pancam wavelengths, was judged to be adequate for Mastcam calibration, which uses parameters determined for a nearby Pancam wavelength to model the directional scattering behavior of the calibration target materials. To determine the coefficient for conversion to radiance factor, we assume that the average ROI radiance values for the three grayscale rings, when plotted against their modeled reflectance values, should fall along a straight line passing through the origin (zero radiance at zero reflectance). The coefficient is derived from the slope of the best-fit line.

A significant complication to this procedure is the deposition of martian dust on the calibration target, even early in the mission owing to material raised during the rover's “sky crane” landing, as well as subsequent gradual deposition of airfall dust from the atmosphere. To account for the influence of the dust on the caltarget reflectance values, a two-layer radiative transfer model (Hapke 1993 section 9.D.3) is employed assuming a uniform layer of dust over the selected caltarget ROIs. The model treats single-scattering events in full detail and uses a two-stream formalism (e.g., Zdunkowski et al. 2007) to treat multiple-scattering events. The dust model and procedure follows very closely the one described in full detail in Kinch et al. (2015), which was developed for dust correction on MER as an improvement to the two-stream model (Kinch et al. 2007) presently implemented for the MER Pancam data sets. The utility of this model for analysis of dusty caltarget surfaces was demonstrated by Johnson et al. (2006) in laboratory studies. Preliminary work on employing the dust model for the MSL Mastcam was presented in Kinch et al. (2013), and the model as employed on Mastcam is described in Bell et al. (2016 in review1).

The dust correction procedure fits the scattering model to the three observed caltarget radiances (which are averages of black, gray, and white ring ROI values). The model requires that the dust single-scattering albedo at each wavelength be specified. The spectrum of dust single-scattering albedo values is found from an analysis of all caltarget images over the first 816 sols of the mission. The scattering model, which has two free parameters, is run on all images using many different values for single-scattering albedo. The free parameters are the incoming irradiance and the extinction optical depth of the dust layer on the caltarget. The “correct” albedo value for each filter is the one that results in stable values for the incoming irradiance, given the known variation in Sun-Mars distance and atmospheric dust loading. If the dust was assumed to be too dark, the derived incoming irradiances would drift to higher values as the caltarget gets dustier, and vice versa for dust assumed to be too bright. This procedure as employed on MER Pancam is described in full detail in Kinch et al. (2015). When employed on MSL Mastcam the procedure results in a dust spectrum that is very similar to the spectra derived for the two MERs.

Once the single-scattering albedo spectrum for the dust is determined, the model can be run with those values on every caltarget sequence, returning the best-fit values for the incoming irradiance and dust optical depth. Dusty caltarget reflectance values for each filter can be determined from the model-derived incoming irradiance and radiances. To find the coefficients to convert the images to reflectance, these corrected reflectance values can then be plotted against the measured radiances for the three grayscale rings, and fit with a line passing through the origin, as described above for a clean caltarget.

Observations calibrated to radiance factor (I/F) through the above procedure are divided by the cosine of the solar incidence angle to a quantity referred to as “relative reflectance,” an approximation of the reflectance factor defined in Hapke section 10.B (1993). Relative reflectance spectra of calibrated observations are presented as the mean values derived from manually defined regions-of-interest (ROIs) made in each camera's field of view and plotted against the filter's effective band-center wavelength. The ROIs from which the values are determined are carefully selected by the following criteria: to include only pixels from a spectrally uniform region, to avoid edges, to be as spatially identical as possible between the two cameras, and to avoid sloping, shadowed, or highly textured regions, to the extent feasible. Saturation sometimes occurs in Mastcam observations, especially in the shortest wavelength filters, as a result of the high contrast between strongly absorbing dust and soil, and less dusty disturbed materials. Pixels within an ROI with saturation in any filter of a multispectral sequence are ignored in others from the same camera to prevent biasing some filter values relative to the others. Right-eye filter values are scaled by the ratio of L6/R6 (these filters are less affected by uncertainties in the dust correction than other, shorter-wavelength stereo filters), to remove any offset between the two cameras, and averaged with the left-eye values at overlapping wavelengths to produce one merged spectrum. The vertical error bars at each point represent the standard deviation of the selected pixels within each ROI. Note that, therefore, these error bars do not represent the absolute or filter-to-filter uncertainty of the radiometric calibration, but instead represent the variation of the data values within each ROI, which is largely due to small-scale differences in solar incidence angle as a result of the surface texture. Owing to the smaller spatial resolution of the M-34 camera, the standard deviations may be smaller than the M-100 values. The absolute radiometric accuracy for Mastcam has been estimated at 10–20% (Bell et al. 2013). Relative uncertainties are likely similar to values derived for Pancam, for which a filter-to-filter uncertainty of <5% and pixel-to-pixel variation of <1% were reported (Bell et al. 2006).

Results and analysis

Mastcam has acquired hundreds of multispectral observations of science targets along the rover's traverse. In the sections below, we describe in detail the reflectance spectra from the soil scoop marks, drill tailings piles, and a few additional targets for which inferences on specific spectral features can be made on the basis of comparisons with other rover or orbital data sets. A map showing the position of the rover along the traverse path where each included multispectral observation was acquired is shown in Figure 4, while details of the Mastcam observations are listed in Table 2. These spectra are also available as supplemental material2 published online.

Rocknest

Curiosity's initial traverse took it in the direction of the “Glenelg” site approximately 450 m to the east of the landing site, where three units identified by orbital mapping adjoin (Anderson and Bell 2010; Grotzinger et al. 2014). En route, the rover spent several tens of sols (martian days) studying an aeolian deposit at a site called “Rocknest,” where five scoops of material were collected and processed by the Collection and Handling for Interior Martian Rock Analysis (CHIMRA) unit to remove any residual contamination and to provide initial sand-sized samples for analysis to the internally housed CheMin and SAM instruments (Blake et al. 2013; Anderson et al. 2015b). The soil scooping operation at Rocknest provided a fresh view of the interior of the sand shadow (Figs. 5a–5c). Multispectral observations were acquired on sol 84 in the M-100 filters only (Fig. 5d), owing to operational constraints. Though CheMin analyses of the basaltic soil revealed >40 wt% mafic silicates (olivines and pyroxenes) as crystalline components (Bish et al. 2013), the reflectance spectra of the sub-surface soils exhibit at best only a weak near-infrared absorption feature. Together with the reddish slope in the visible wavelengths, this result suggests that the reflectance is dominated by an amorphous component with strongly wavelength-varying optical properties in the visible, which may include nanophase and/or poorly crystalline ferric oxides. Indeed, the presence of a substantial iron-bearing amorphous component is supported by mass balance considerations from CheMin and APXS data (Blake et al. 2013). The overall reflectance of the material inside the scoop mark is considerably lower than the dusty surface (Fig. 5a), although some vertical heterogeneity exists in the form of a lighter-toned layer visible in the side wall of each scoop (Figs. 5b–5c), perhaps reflecting dust deposition during a past interlude in sand accumulation. Spectra from neighboring scoop marks in Figure 5 are very similar to each other, with minor differences in reflectance likely attributable to a variation in photometric angles from the relative orientation of walls of the troughs.

Yellowknife Bay

Curiosity departed the Rocknest location on sol 102 and continued eastward to Glenelg, where it conducted an extensive scientific campaign in and around a topographic depression known as Yellowknife Bay. The strata of the Yellowknife Bay Formation are interpreted to have been deposited in a fluvial-lacustrine setting (Grotzinger et al. 2014), and it is from the stratigraphically lowest member of this formation, the Sheepbed Mudstone unit, that the first two drill samples were obtained on sols 182 and 279, with accompanying Mastcam multispectral observations [for a detailed stratigraphy of this and other Gale Crater units, see Grotzinger et al. (2015)]. The first of these, “John Klein,” was preceded by a “mini” drill hole on sol 180, a shallow-depth drilling performed to assess the suitability of the target material for the full-depth drill activity. Figure 6a is a reflectance-calibrated color M-100 RGB image from a multispectral sequence acquired on sol 183 and shows both the mini drill hole (nearer the top of the frame) and full drill hole. Reflectance spectra were derived from an ROI positioned on the mini drill tailings pile, which is distributed more evenly on the surface than the full drill tailings pile. In almost all cases (the one exception being Confidence Hills, as mentioned below), spectra drawn from the tailings of adjacent full and mini drill holes reveal no substantial differences in reflectance. Figure 6b shows the Cumberland full drill hole and tailings pile, located just a few meters from John Klein, and the region of interest from which the mean spectrum was taken. Reflectance spectra for these two sets of drill fines are plotted in Figure 6e, together with spectra derived from other drill targets as well as a dust-covered surface near the John Klein target. Differences between the Sheepbed drill fines and other materials are readily apparent from the plot. Compared to the undisturbed, dusty surface, the Sheepbed drill tailings are substantially less “red” (possess a lesser reflectance slope across the visible wavelength range), as indeed are all of the drill fines. The John Klein and Cumberland tailings possess a clear reflectance maximum near 805 nm, with a broad near-infrared absorption band that appears to be centered shortward of 1000 nm. These VNIR characteristics are consistent with a significant contribution to the spectral shape from the high abundance of smectite (~20 wt% of total sample) and/or low-Ca pyroxenes (~10 wt%) detected by the CheMin instrument (Vaniman et al. 2014). The library mineral spectra plotted in Figure 2 include a sample of both trioctahedral saponite and Fe3+-rich, dioctahedral nontronite, in addition to several other reflectance spectra of mineral species similar to those detected by CheMin in these or subsequent drill samples. The shoulder near 650 nm in the two smectite laboratory spectra does not appear in the multispectral data, suggesting it is either too weak to be apparent at the spectral sampling provided by the filters, or else it is not present (perhaps suppressed by other phases). The ChemCam passive data show an absence of such a band, as well as the beginning of a band in the near-infrared [Fig. 7 and Johnson et al. (2015)], consistent with the multispectral data. It can be seen from the convolved spectra in Figure 2 that this near-infrared spectral band, which is centered around 950 nm, could be reproduced with either a smectite spectrum or a pyroxene with features intermediate between the augite and enstatite spectra shown. Unlike the library mineral spectra, the Sheepbed material is a complex mixture of multiple mineral and amorphous phases, and the inherently nonlinear nature of such spectral mixtures makes it not straightforward to evaluate quantitatively the relative contribution of each to the overall spectra shape. Compared to the other drill tailings (which do not possess a band at this position), the Sheepbed samples are distinct in possessing significantly more smectite, which may suggest that it is the clay mineral that is contributing substantially to the unique spectral shape.

The Kimberley

Curiosity did not depart the Yellowknife Bay region until sol 324, and the subsequent emphasis on driving placed nearly 5 km between the Yellowknife Bay drill locations and the next drill hole. The rover arrived at the Kimberley drill location on sol 609 and performed full-depth drilling on sol 621 on a cross-bedded sandstone target in the Dillinger unit (Le Deit et al. 2015; Treiman et al. 2016) named “Windjana,” with an accompanying Mastcam multispectral observation on sol 626. The Windjana drill target produced mini and full drill tailings that are very different from the Sheepbed targets. The observation acquired on sol 626 (Fig. 6c) shows the mini (nearer the bottom of the frame) and full drill holes, the former of which is largely filled in with tailings vibrated back into the hole by the action of the subsequent full drill operation. Spectra derived from the full drill tailings exhibit the lowest overall reflectance of drilled material observed in the mission to date. These low reflectance values are consistent with the relatively higher abundance of strongly absorbing phases, especially higher magnetite (~12 wt% of total sample), reported by the CheMin team (Treiman et al. 2015, 2016). In addition, the reported phyllosilicate abundance in the Windjana drill materials is significantly lower (~10 wt%) than in the Sheepbed drill materials, which may also contribute to the substantially flatter spectral shape near 800 nm as compared to the previous two tailings piles. The general decrease in reflectance at the longest wavelength filters is consistent with absorptions due to mafic minerals.

Hidden Valley

After the Kimberley, the rover continued driving south and southwest, reaching the boundary of the landing ellipse on sol 672. Near an untraversable patch of sand ripples within a topo-graphic low dubbed “Hidden Valley,” a fourth drill attempt was made on sol 724 on what proved to be an unstable rock slab (“Bonanza King”). Images acquired after the failed mini-drill attempt into the target revealed that the rock had shifted during the drill activity (Anderson et al. 2015a), and thus the decision was made not to reattempt drilling on this or a nearby rock. Although no XRD data exist for this target, the mini-drill did produce a small pile of fresh tailings suitable for multispectral analysis (Fig. 6d). The resulting spectrum (Fig. 6e) is distinctly different from both the Sheepbed (John Klein and Cumberland) and Windjana spectra. Compared to those previous tailings piles, Bonanza King is intermediate in overall reflectance and extremely flat over the wavelength range sampled by Mastcam filters. (The unreasonably high value for the blue Bayer filter at 494 nm should be considered an artifact, perhaps resulting from the combination of the broad Bayer bandpass with a spectral radiance curve possessing a different shape than that of most other materials observed.) By analogy with other drill tailings spectra, we can make inferences about the mineralogy of Bonanza King despite the lack of X-ray diffraction data from CheMin. The higher overall reflectance as compared to Windjana implies that Bonanza King either has a smaller percentage of the strongly absorbing species such as magnetite that are present in the Windjana sample, or else also possesses a spectrally neutral species with relatively high reflectance, such as a silica phase. From the nearly featureless spectrum of Bonanza King we can also say, qualitatively, that the percentage of other iron-bearing minerals with strong features in these wavelengths, such as hematite, Fe-saponite, and the Fe-sulfate detected nearby (see the discussion of the “Perdido2” target in Johnson et al. (2016) and below), are likely either absent or lower in abundance compared to targets that show these features more clearly. APXS data show that the Mg/Si and Al/Si ratios are quite low for Bonanza King as compared to most other analyses, which may reflect a history of open-system aqueous alteration (Yen et al. 2015). The lack of distinct reflectance features in the multispectral data is broadly consistent with the breakdown of primary mafic minerals as in such a weathering environment.

Pahrump Hills

After an additional ~600 m of driving, Curiosity began a detailed investigation of basal Mt. Sharp units at Pahrump Hills, where it observed a sequence of predominantly fine-grained mudstones and siltstones that form the lowermost units of the Murray Formation (Stack et al. 2015). Three additional drill samples (“Confidence Hills,” “Mojave 2,” and “Telegraph Peak”) were obtained in these layers on sols 759, 882, and 908 (respectively), with accompanying Mastcam multispectral observations. The first attempt at the Mojave site on sol 867 resulted in only a partial drill hole, which nevertheless produced clean material for multispectral analysis. This target is referred to below as “Mojave”, whereas “Mojave 2” refers to the second, successful drilling that produced a sample for CheMin analysis.

The Confidence Hills drill site was the stratigraphically lowest of the three Pahrump drill locations. Figure 8a shows the mini (lower) and full drill holes and drill tailings. The tailings to the lower left of the drill holes have a slightly higher reflectance (~0.02 higher) than the drill tailings closer to the full drill hole. Unfortunately, the full drill activity significantly disturbed soil from the surrounding surface as well as the tailings pile that previously surrounded the mini drill hole, making it difficult to ascertain the extent to which the tailings piles and soil have been mixed together. The difference in reflectance may be a result of such mixing. An alternative is that differences in composition and/or grain size exist over spatial scales and/or drill depths as small as that sampled here. Besides this small difference in overall reflectance, the tailings are spectrally similar and show a strong absorption in the shorter wavelength filters (around 550 nm) not present in any of the previous drill tailings spectra (cf. Fig. 6). The CheMin X-ray diffraction data revealed that a significant quantity of crystalline hematite (~8 wt% of total sample; Cavanagh et al. 2015) was present in the sieved sample, compared to <1 wt% in prior samples (Bish et al. 2013; Vaniman et al. 2014; Treiman et al. 2015, 2016). Interestingly, the Mastcam spectra do not show any evidence for a band near 860 nm that is typically associated with fine-grained crystalline hematite (e.g., Morris et al. 1985), nor do they show a downturn beyond about 785 nm as is present in the longest wavelengths of ChemCam passive spectra [Fig. 7 and Johnson et al. (2016)], suggesting that such a band, if present, must be very shallow. The 860 nm band may be largely masked by other constituent minerals. Alternatively, the hematite may be present as a combination of nanophase and crystalline phases, mixtures of which have been shown to vary with respect to the strength of either spectral feature mentioned above (Morris et al. 1989).

Two further full drill holes were acquired in strata of the Pahrump Hills. A mini drill attempt into the original Mojave target resulted in a dislodged block and fines (Fig. 8b). The attempt was repeated successfully a few meters away (Fig. 8c). Spectra from both Mojave targets exhibit essentially identical Mastcam spectra (Fig. 8e). Both spectra have a reddish spectral slope with an inflection near the 527 nm filter consistent with the presence of crystalline hematite, although this feature is weaker than in the Confidence Hills spectrum. By contrast, the Telegraph Peak drill tailings (shown in Fig. 8d) possess a spectrum that is flatter and exhibits little evidence of an absorption near 527 nm. We compute a 527 nm band depth using the relative reflectance values at 446 and 676 nm as BD527 = 1 – R527/[0.648R446 + 0.352R676]. The 527 nm band depths for the three Pahrump drill targets are 0.23, 0.13, and 0.04 for Confidence Hills, Mojave 2, and Telegraph Peak, respectively. This parameter correlates well with the measured abundances of crystalline hematite, which were determined by the CheMin team to be ~8, ~4, and ~1 wt%, respectively (Cavanagh et al. 2015; Rampe et al. 2016). Spectral parameters such as these are also influenced by the other constituent minerals; however, within the range of mineralogies sampled by the Pahrump drill holes, such a trend in spectral shape appears to correlate reasonably well with actual hematite abundances.

Observations of select float rocks and distant mound layers

In addition to the previous spectra of drill powders and scoop marks, several other observations display noteworthy and unique features over the spectral range of the Mastcam filter set and are included to further document the spectral diversity observed by the cameras along first 1000 sols of the rover's traverse. Here we show near-infrared absorption features consistent with iron oxides and iron sulfates in Mastcam spectra, the former in observations targeted at a feature approximately 5 km away that was known from orbital spectral data, and the latter as a spatially heterogeneous feature observed over generally small (~millimeters) spatial scales at certain locations in rock units explored by the rover near the basal Mt. Sharp stratigraphy. Each of these observations are discussed in turn below, beginning with a feature observed in several small rock fragments in Hidden Valley.

Fresh surfaces of rocks broken by the rover wheels have been very rewarding targets for multispectral observations. One such observation acquired on sol 721 targeted several broken rocks, including the ChemCam target “Perdido” (Fig. 9). Johnson et al. (2016) report numerous detections of a 430 nm band that suggests the presence of a ferric sulfate mineral such as jarosite, particularly when paired with a near-infrared reflectance downturn at wavelengths greater than about 700 nm (e.g., Rossman 1976; Cloutis et al. 2006), with “Perdido2” (a passive ChemCam raster pointed just below the original LIBS Perdido target) expressing this feature most strongly. The 430 nm band is not detectable with Mastcam's filters, but portions of the rock fragment Per-dido2 and the surrounding rocks do possess a strong, distinctive absorption in the near-infrared at 900 nm, consistent with a ferric mineral such as jarosite. Unfortunately, the rock fragments are so small (centimeter scale and smaller) that the selection of an appropriate ROI is very difficult in the M-34 images. Figure 9c shows a decorrelation stretch of filters R4, R5, and R6 (908, 937, and 1013 nm), in which bluish colors correspond to material possessing the near-infrared feature. Close examination of the Perdido2 target in this false-color image suggests that the rock may show this feature strongly only on a small portion of its surface. For this reason, the red spectrum in Figure 9d is derived from a nearby fragment expressing this feature in more pixels than the Perdido2 target itself and includes only values from the higher spatial resolution right-camera (M-100) filters. The strongest detection of a 430 nm feature in ChemCam was from the final spot of a 5 × 1 dedicated passive observation (Johnson et al. 2016), whose pointing was extremely close to the portion of the Perdido2 fragment that the decorrelation stretch image suggests may have strong near-infrared Mastcam features as well. A jarosite spectrum from the USGS spectral library (Clark et al. 2007) is shown in Figure 2, depicting the long-wavelength absorption that corresponds favorably with the Mastcam data. While this feature is not uniquely indicative of jarosite, the combination of the Mastcam and ChemCam passive data as well as the subsequent identification of jarosite by CheMin in nearby drill samples (Cavanagh et al. 2015; Rampe et al. 2016) strongly suggests that this ferric sulfate is responsible for this spectral feature.

Long-distance multispectral observations targeting the lower layers of the mound were acquired on sols 468 and 475 (the latter is shown in Fig. 10a), in a coordinated effort with Chem-Cam passive observations. The field of view included a ridge referred to as the “hematite ridge” from CRISM observations identifying a layer bearing a spectral signature consistent with that oxide (Milliken et al. 2010; Fraeman et al. 2013). Mastcam spectra show that portions of a layer bear a weak but clearly present absorption feature near the 867 nm filter, consistent with crystalline hematite, as well as a weak inflection near 527 nm (Fig. 10c). These observations are consistent with the enhanced 535 band depths and stronger near-infrared downturns observed by ChemCam (Johnson et al. 2016). The presence of the long-wavelength spectral feature as a characteristic of a spatially coherent region of the ridge is most obvious in the decorrelation stretch image of three near-infrared bands shown in Figure 10b, in which the bluish-purple false color corresponds to the region exhibiting this feature. The other layers of the mound visible from this vantage point do not show any distinguishing spectral features. The dunes show a generally low reflectance and broad long-wavelength decrease consistent with a mafic mineralogy, likely with some influence from a thin dust layer. Brighter material appears to lie in a “fan” of material just above the dunes, but both this material and the topographically higher layers lack any strong or distinctive spectral features aside from the obvious influence of reddish dust at the shorter wavelengths.

Implications

Mastcam has observed an impressive diversity of reflectance spectra within Gale crater. Tailings from each lithologic unit drilled by the rover have distinctive spectral properties, as do certain small-scale materials observed in broken rock fragments and long distance observations of specific layers of Mt. Sharp. The multispectral data, which track changes in Fe mineralogy, are in good general agreement with mineralogical and elemental analyses conducted by other instruments on Curiosity's payload, which corroborate the detection of discrete compositional units. Clear correlations can be drawn between the mineralogy inferred from CheMin X-ray diffraction data and the visible/near-infrared spectral features identified in calibrated Mastcam spectra, further supported by results from geochemical instruments. Mastcam multispectral and ChemCam passive observations provide two independently calibrated measurements of reflectance and display equivalent features within the wavelength range where they overlap.

The drill tailings observed at Yellowknife Bay possess an 800 nm spectral peak and minimum near 930 nm, apparently deriving from Fe-bearing phyllosilicates and pyroxenes. Saponite can accommodate a range of iron content and redox state (Kodama et al. 1988; Treiman et al. 2014), which should be manifested by changes in the visible/near-infrared reflectance. Further study may permit using Mastcam VNIR data to resolve between the ferrian saponites studied by Treiman et al. (2014) and high-Mg ferrosaponites studied by Chemtob et al. (2015) as candidate compositions that fit the 02l diffraction peak measured by CheMin at Yellowknife Bay.

On the basis of data from other instruments, the identity of other minerals producing distinct spectral features in the Mastcam data set can be reasonably ascertained. The presence of magnetite serves to suppress reflectance across the entire VNIR range, which is most apparent in spectra of the Windjana drill tailings. The Pahrump Hills targets of Confidence Hills and Mojave possess a spectral absorption or inflection near the 527 nm filter consistent with the presence of crystalline hematite and stronger than that in surface dust spectra. The strength of this parameter appears to correlate well with hematite abundances measured by X-ray diffraction. A longer wavelength hematite feature near the 867 nm filter is not apparent in the drill fines spectra, but does appear weakly in long-distance observations targeted toward the hematite-bearing ridge, whose hematite signature was detected from orbit (Milliken et al. 2010; Fraeman et al. 2013). The ridge, therefore, may possess a greater relative abundance of crystalline hematite than has yet been sampled. A broad 900 nm near-infrared feature observed in Mastcam spectra is consistent with a ferric sulfate mineral, likely jarosite, found in broken rock near the Perdido target, and is correlated with a ~430 nm feature in ChemCam passive data. Detection of the 900 nm feature to date has been limited to material at small spatial scales close to the rover rather than within an extensive bedrock unit. This limited detection is consistent with the lack of detection in orbital CRISM data as well, which have spatial resolutions of 18 m/pixel. Overall, the Mastcam observations are strongly consistent with the results of ChemCam passive spectra observations on the same targets (Johnson et al. 2016).

Through sol 1000 of the Curiosity mission, the lower mound layers bearing the strongest CRISM signatures of phyllosilicate, sulfate, and iron oxide minerals still lie ahead of the rover. Work by Fraeman et al. (2015) in differentiating CRISM units on the basis of thermophysical and spectral parameters provides a forecast for promising future sites of potential multispectral investigation. Mastcam multispectral observations have demonstrated the ability to identify diversity not readily apparent in RGB color images or orbital spectra. Understanding the mineralogical implications of Mastcam spectral features based on inferences drawn from previous instrument collaborations provides the ability to interpret as well as to distinguish distinct spectral units. The rover's continued ascent up the slopes of Mt. Sharp is expected to bring it into contact with more of the mound's diverse iron mineralogies and associated spectral diversity, including arriving at the hematite ridge itself and the exposures of nontronite expected just beyond (Milliken et al. 2010; Thomson et al. 2011; Fraeman et al. 2013). Further analyses should be able to not only assist reconnaissance spectral imaging for the Curiosity mission, but also provide insight into multispectral analysis of martian materials and imaging strategies relevant for the high-heritage Mastcam-Z instrument (Bell et al. 2014) aboard the upcoming Mars2020 mission, and the similar ExoMars PanCam investigation (Coates et al. 2015) as well.

1 Malin, M.C., Ravine, MA., Caplinger, M.A., Ghaemi, F.T., Schaffner, J.A., Maki, J.N., Bell, J.F. III, Cameron, J.F., Dietrich, W.E., Edgett, K.S., and others, The Mars Science Laboratory (MSL) Mast cameras and Descent imager: I. Investigation and instrument descriptions, submitted to Earth & Space Science, 2016. Bell, J.F. III, Godber, A., McNair, S., Caplinger, M.A., Maki, J.N., Lemmon, M.T., Van Beek, J., Malin, M.C., Wellington, D., Kinch, K.M., and others, The Mars Science Laboratory Curiosity rover Mast Camera (Mastcam) instruments: Pre-flight and in-flight calibration, validation, and data archiving, submitted to Earth and Space Science, 2016.
2Deposit item AM-17-65760, Supplemental Material. Deposit items are free to all readers and found on the MSA web site, via the specific issue's Table of Contents (go to http://www.minsocam.org/MSA/AmMin/TOC/2017/Jun2017_data/Jun2017_data.html).

Acknowledgments

This work was supported by Mars Science Laboratory Project grants and contracts from NASA, JPL/Caltech, and Malin Space Science Systems, Inc. We heartily thank the MSL/Curiosity science and operations team, who made these results possible. The authors also acknowledge specific funding sources, as follows: the MSL Participating Scientist program (J.R. Johnson, JPL contract 1350588, and B.L. Ehlmann), the W.M. Keck Institute for Space Studies and Caltech GPS Division Texaco postdoctoral fellowships (A.A. Fraeman), the NASA Postdoctoral Program (M.S. Rice), and the Danish Council for Independent Research/Natural Sciences (K.M. Kinch, FNU Grant 12-127126). A portion of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. Finally, the authors thank reviewers Edward Guiness and Jack Mustard for greatly appreciated feedback that allowed us to significantly improve upon the initial manuscript.