## H O 2 2 WITHIN CHAOS TERRAIN ON EUROPA'S LEADING HEMISPHERE

Samantha K. Trumbo, Michael E. Brown

Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125, USA

## Kevin P. Hand

Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA

## ABSTRACT

Hydrogen peroxide is part of an important radiolytic cycle on Europa and may be a critical source of oxidants to the putative subsurface ocean. The surface geographic distribution of hydrogen peroxide may constrain the processes governing its abundance as well as its potential relevance to the subsurface chemistry. However, maps of Europa's hydrogen peroxide beyond hemispherical averages have never been published. Here, we present spatially resolved L-band (3.16-4 µ m) observations of Europa's 3.5 µ m hydrogen peroxide absorption, which we obtained using the near-infrared spectrometer NIRSPEC and the adaptive optics system on the Keck II telescope. Using these data, we map the strength of the 3.5 µ mabsorption across the surface at a nominal spatial resolution of ∼ 300 km. Though previous diskintegrated data seemed consistent with the laboratory expectation that Europa's hydrogen peroxide exists primarily in its coldest and iciest regions, we find nearly the exact opposite at this finer spatial scale. Instead, we observe the largest hydrogen peroxide absorptions at low latitudes on the leading and anti-Jovian hemispheres, correlated with chaos terrain, and relative depletions toward the cold, icy high latitudes. This distribution may reflect the effects of decreased hydrogen peroxide destruction due to efficient electron scavenging by CO 2 within chaos terrain.

Keywords: planets and satellites: composition - planets and satellites: individual (Europa) - planets and satellites: surfaces

## 1. INTRODUCTION

The continuous bombardment of Europa's surface by energetic particles within Jupiter's magnetosphere results in extensive radiolytic processing of the surface ice (e.g. Johnson &amp; Quickenden 1997). The impinging electrons, protons, and ions lead to the dissociation of water molecules, producing OH radicals, which can recombine to produce hydrogen peroxide (H O ) and other oxi2 2 dants (e.g Johnson et al. 2003; Cooper et al. 2003; Loeffler et al. 2006). H 2 is lost in the process, creating an increasingly oxidizing surface. Understanding this radiolytic cycle on Europa is not only important in terms of Europa's surface composition and the general study of surface-magnetosphere interactions on icy bodies, but it is also critical for assessing the potential chemical energy sources within Europa's subsurface ocean. Hydrothermal processes at the seafloor may be a source of reductants, but the potential habitability of Europa's ocean will likely depend on a complementary supply of oxidants, such as H O , from the radiolytically processed 2 2 surface (Chyba 2000; Hand et al. 2009).

The geographic distribution of Europa's H O 2 2 may hold clues to the processes controlling its presence on the surface as well as to its potential delivery to the subsurface ocean. The Galileo Near Infrared Mapping Spectrometer (NIMS) definitively detected H O 2 2 on the surface via a prominent 3.5 µ mabsorption in an average spectrum of the leading/anti-Jovian quadrant (Carlson et al. 1999). However, the intense radiation encountered during closer flybys of Europa resulted in poor-quality spectra beyond 3 µ m, and no spatially resolved NIMS maps of H O 2 2 were ever published. To date, groundbased spectra, in which the 3.5 µ m absorption of H O 2 2 was rotationally resolved across four nights of observation, have provided the best constraints on its geography (Hand &amp; Brown 2013). These spectra show the largest absorptions on the leading and anti-Jovian hemispheres, in agreement with the NIMS observation, and almost no detectable absorption on the trailing hemisphere. At this scale, Europa's H O 2 2 distribution was easily attributed to the availability of surface water ice and the low surface temperatures associated with Europa's iciest regions, which lie at the mid to high latitudes of the leading and anti-Jovian hemispheres (Brown

&amp; Hand 2013). Water molecules are a necessary precursor for H O 2 2 production, and laboratory experiments studying the production of H O 2 2 via water ice radiolysis (Moore &amp; Hudson 2000; Loeffler et al. 2006; Zheng et al. 2006; Hand &amp; Carlson 2011) consistently observe lower equilibrium abundances with increasing temperature, in some cases by nearly a factor of seven across the ∼ 80-120 K temperature range relevant to Europa (Hand &amp; Carlson 2011).

Considering all of the published observations and laboratory data, one might reasonably expect that higher spatial resolution observations show an abundance of H O 2 2 at the icy high latitudes of Europa's leading and anti-Jovian hemispheres, and relative depletions at the warm equator, which likely contains a significant nonice fraction in the form of endogenous salt deposits associated with largescale regions of geologically young and highly disrupted chaos terrain (Fischer et al. 2015, 2017; Trumbo et al. 2019). In order to address this hypothesis, we obtained spatially resolved L-band spectra of Europa's surface, with the goal of mapping the geographic distribution of its 3.5 µ m H O 2 2 band.

## 2. OBSERVATIONS AND DATA REDUCTION

We observed Europa on Feb. 24th and 25th of 2016 and again on June 6th of 2018 using the near-infrared spectrograph NIRSPEC and the adaptive optics system on the KECK II telescope. We used the 3.92 ′′ x 0.072 ′′ slit in low-resolution mode (R ∼ 2000). Our data cover the wavelength range of ∼ 3.16-4 µ min the L-band. As a telluric calibrator for the 2016 observations, we observed HD 98947, a V = 6.9 G5 star that was ∼ 1.2 ◦ from Europa on the sky. In 2018, we used HD 128596, a V = 7.5 G2 star with a 3.9 ◦ separation from Europa. We observed all targets in an ABBA nodding pattern. For our 2016 data, each Europa pointing consisted of 30 20second coadds, and each calibrator pointing consisted of 15 2-second coadds. In 2018, we reduced the integration times to 30 10-second co-adds for each Europa pointing and 10 2-second coadds for each calibrator pointing.

During both sets of observations, Europa had an angular diameter of nearly 1 ′′ , corresponding to ∼ 10 300km resolution elements at the diffraction limit of Keck at 3.5 µ m. For each Europa pointing, we aligned the slit in either an E/W or N/S orientation with respect to Europa's north pole. We obtained a corresponding SCAM guide camera image for each exposure as a record of the exact slit positions on the disk. The telescope pointing tended to drift significantly during nods for the 2018 observations. Thus, we used these images in real time to manually maintain a consistent slit position across each ABBA set.

We reduce all data in Python, following the standard methodology of image rectification, image pair subtrac- tion, residual sky subtraction, and telluric and wavelength calibration, utilizing both the Astropy (Robitaille et al. 2013) and skimage.transform (van der Walt et al. 2014) packages. We use an ATRAN atmospheric transmission spectrum (Lord 1992) for wavelength calibration. Some of the raw spectral files display readout artifacts, which artificially brighten every 8th pixel in the dispersion direction. Thus, we replace each affected pixel by the average of the two immediately adjacent pixels.

Following this initial data reduction, we extract spectra for individual spatial resolution elements within each slit position. Using the Python Basemap package and the calculated size of Europa in SCAM pixels (0.0168 ′′ /pixel), we align the SCAM image corresponding to each NIRSPEC exposure with an orthographic projection and estimate the coordinates of the slit on the disk. Though we made real-time manual corrections to the inaccurate telescope nodding experienced during the 2018 observations, some ABBA sets still exhibited significant position discrepancies between nods. Thus, we align SCAM images and estimate the corresponding slit coordinates for each nod independently. We then locate Europa in the 2D spectral images using its calculated size in NIRSPEC pixels (0.013 ′′ /pixel) and extract individual spectra corresponding to an 8-pixel spatial resolution element, stepping by 4 spatial pixels between extractions. We determine the geographic coordinates of each extracted spectrum using the previously obtained slit positions and the aforementioned NIRSPEC pixel scale.

To minimize geographic smearing, while still maximizing the signal-to-noise, we average spectra from overlapping slit positions, as determined by manually inspecting the SCAM images and the extracted coordinates corresponding to each nod. In the case of the trailing/subJovian slits, we tolerate moderate spatial smearing in favor of the enhanced signal-to-noise needed to quantify weak H O 2 2 absorptions. In all cases, we account for any resulting adjustments to the effective slit width and to the geographic area covered by each spectrum. Table 1 summarizes our observations and the geometries of geographically unique slit positions.

## 3. H O 2 2 MAPPING

In order to determine the geographic distribution of Europa's 3.5 µ mH O absorption, we calculate its band 2 2 area in each individual extracted spectrum. For most spectra, we fit a second-order polynomial from 3.37 to 3.715 µ m, excluding the portion corresponding to the H O 2 2 absorption ( ∼ 3.4-3.65 µ m). For the trailing hemisphere observations, we use a third-order polynomial, which better matches the continuum shape. We inspect all fits by eye and, if necessary to achieve a satis-

Table 1 . Table of Observations

| Date (UT)   | Hemisphere          | Slit Orientation   | Airmass Range   |   Total Int. Time (min) | Central Longitude   | Central Latitude   | Telluric Calibrator   |
|-------------|---------------------|--------------------|-----------------|-------------------------|---------------------|--------------------|-----------------------|
| 2016 Feb 24 | trailing/sub-Jovian | E/W                | 1.10-1.75       |                     120 | 338 W               | 1 S                | HD 98947              |
| 2016 Feb 24 | sub-Jovian          | N/S                | 1.03-1.07       |                     100 | 348 W               | 2 S                | HD 98947              |
| 2016 Feb 25 | leading             | E/W                | 1.24-1.38       |                      40 | 79 W                | 0 N                | HD 98947              |
| 2016 Feb 25 | leading             | E/W                | 1.11-1.18       |                      40 | 82 W                | 12 S               | HD 98947              |
| 2016 Feb 25 | leading             | N/S                | 1.04-1.08       |                      40 | 88 W                | 2 S                | HD 98947              |
| 2016 Feb 25 | leading             | E/W                | 1.08-1.14       |                      40 | 95 W                | 0 N                | HD 98947              |
| 2018 Jun 06 | leading             | N/S                | 1.40-1.42       |                      10 | 119 W               | 4 S                | HD 128596             |
| 2018 Jun 06 | leading             | N/S                | 1.36-1.38       |                      10 | 124 W               | 4 S                | HD 128596             |
| 2018 Jun 06 | leading/anti-Jovian | N/S                | 1.24-1.30       |                      20 | 166 W               | 3 S                | HD 128596             |
| 2018 Jun 06 | leading/anti-Jovian | N/S                | 1.22-1.26       |                      20 | 146 W               | 3 S                | HD 128596             |
| 2018 Jun 06 | leading/anti-Jovian | N/S                | 1.25-1.25       |                      10 | 162 W               | 3 S                | HD 128596             |
| 2018 Jun 06 | leading             | N/S                | 1.22-1.22       |                      15 | 115 W               | 4 S                | HD 128596             |
| 2018 Jun 06 | leading             | N/S                | 1.23-1.25       |                      20 | 104 W               | 3 S                | HD 128596             |
| 2018 Jun 06 | leading             | N/S                | 1.26-1.30       |                      20 | 87 W                | 3 S                | HD 128596             |
| 2018 Jun 06 | leading             | N/S                | 1.30-1.32       |                      10 | 120 W               | 4 S                | HD 128596             |

factory continuum fit, make small adjustments to these parameters. We then remove the calculated continuum from each spectrum and integrate the residual absorption to calculate the equivalent width (i.e. the width of a 100% absorption of the same integrated area). We map the resulting band areas using the slit positions and widths obtained as in Section 2. Our 2016 data show generally stronger absorptions than do our 2018 data, with maximum band areas ∼ 25% larger than those observed in 2018. This is perhaps unsurprising given that H O 2 2 concentrations on Europa reflect a dynamic equilibrium between constant formation and decay (Hand &amp; Carlson 2011) that may be influenced by the temporal variability of the radiation environment or of the local surface temperature. Thus, as we are concerned primarily with geographic correlations, we map the 2016 and 2018 slits separately.

K (Spencer et al. 1999), and visible and infrared spectra suggest that it is likely salty in composition (Fischer et al. 2015, 2017; Trumbo et al. 2019). Yet, its 3.5 µ m H O 2 2 band area reaches more than twice that of the most spectrally icy region. Figure 1 shows spectra representative of both locations.

Figure 1 shows a single mapped slit from our 2016 observations. This slit position crosses both the most spectrally icy location on the surface at ∼ 30 ◦ N and 90 ◦ Won the leading hemisphere (Brown &amp; Hand 2013) and the comparatively warm and ice-poor chaos region Tara Regio, making it ideal for testing the hypothesis that Europa's H O 2 2 prefers cold, icy terrain. Puzzlingly, contrary to this hypothesis, the 3.5 µ m H O 2 2 absorptions are strongest in exactly the warmest and least icy portion of the slit, nearly perfectly corresponding to the bounds of Tara Regio (mapped approximately from Doggett et al. (2009)). Tara Regio is located at low latitudes, where daytime temperatures approach ∼ 130

The apparent preference of H O 2 2 for low-latitude chaos terrain persists across all of the mapped 2016 slits (Figure 2A). Like the N/S slit of Figure 1, the leading-hemisphere E/W slit positions also show strong absorptions across Tara Regio. In addition, they demonstrate a marked difference between low-latitude chaos and low-latitude plains terrain, displaying enhanced absorption across the other large-scale chaos regions of the leading/anti-Jovian hemisphere (also mapped approximately from Doggett et al. (2009)), but much weaker band strengths across the ridged plains immediately east of Tara Regio ( ∼ 30 ◦ W). In fact, even the resolution elements that include the narrow space between Tara Regio ( ∼ 85 ◦ W) and eastern Powys Regio ( ∼ 125 ◦ W) appear to exhibit somewhat weakened absorptions on average. The trailing/sub-Jovian regions, by comparison, show only very weak 3.5 µ m bands, consistent with the full-disk data of Hand &amp; Brown (2013). However, what little H 2 O 2 is present still seems correlated with chaos, as suggested by the slight enhancement near the outlined western-most portion of Annwn Regio ( ∼ 0 ◦ N, 360 ◦ W). We note, however, that accurately quantifying such weak absorptions in this region is difficult, especially given the substantially different continuum shape

Figure 1 . A single N/S slit from our 2016 observations, which falls on the leading hemisphere and crosses both the most spectrally icy location on the surface ( ∼ 30 ◦ N, 90 ◦ W, Brown &amp; Hand 2013) and the salty, low-latitude chaos region Tara Regio. Contrary to the hypothesis that Europa's H O 2 2 should follow the cold, icy terrain of the upper latitudes, the 3.5 µ m band appears anti-correlated with both temperature and ice availability. Instead, the strongest absorptions fall nearly perfectly within the outlined bounds of Tara Regio, with band areas more than twice that of the aforementioned ice-rich region. We show representative spectra from both locations to the right of the map, where the dashed red lines outline the H O 2 2 band. Second-order polynomial continua are indicated by the dashed black curves. We also include the continuum-removed absorptions to ease comparison of the band strengths, although the differences are readily apparent in the spectra themselves. Both spectra are normalized to their individual peaks in the 3.6-3.7 µ m region.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01093v1-figure-fa050094-15f7-43bf-8fc2-defa11cc3496.jpg)

of spectra from the low-latitudes of the trailing/subJovian hemisphere.

While the 2016 map certainly suggests that Europa's H O 2 2 favors chaos terrain, it is also potentially consistent with a simple preference for low latitudes combined with the known hemispherical trend observed previously in full-disk spectra (Hand &amp; Brown 2013). In addition, the E/W slits, though approximately constant in latitude, still each cover a range of surface temperatures reflecting the early morning through late afternoon, which may impart unknown diurnal effects onto the longitudinal patterns we observe. To account for these caveats, our 2018 map (Figure 2B) is composed of only N/S slit positions. This orientation more concretely differentiates between a preference for low latitudes and a specific correlation with the chaos regions, which are asymmetric about the equator. N/S slits also have the advantage of capturing the latitudes within a given slit position at similar times of the local day, thus minimizing diurnal complications to the observed H O 2 2 pattern within a slit.

Our 2018 dataset confirms the trends suggested by the 2016 map-Europa's H O 2 2 prefers warm, low-latitude, ice-poor chaos terrain to the comparatively cold and icerich regions of the upper latitudes. The 2018 N/S slit positions expand our coverage of the large-scale chaos regions of the leading/anti-Jovian hemisphere, and all individual slit positions consistently show the largest 3.5 µ m absorptions within the outlined chaos terrain. In most cases, the observed band areas increase sharply near the largescale chaos boundaries, reaching values up to three times as large as those immediately above or below the chaos terrain. In all cases, especially those of Tara Regio and eastern Powys Regio, this requires asymmetric displacement south of the equator, which is inconsistent with a simple preference for low latitudes. In fact, the complete 2018 map indicates that most of Europa's H O 2 2 falls between ∼ 40 ◦ S and 20 ◦ N, in a latitudinal pattern that possesses similar N/S asymmetries to those of the overall pattern of leading/anti-Jovian chaos terrain. This observed distribution requires dominating influences on Europa's H O 2 2 other than temperature or water ice abundance.

## 4. DISCUSSION

Outside of ice availability and surface temperature, perhaps the most obvious potential factor governing the geographic distribution of H O 2 2 on Europa is the local radiation environment. Though the entirety of Europa's surface is irradiated by energetic particles, the energy flux and particle composition is geographically nonuniform (Paranicas et al. 2001, 2009; Dalton et al. 2013; Nordheim et al. 2018). In fact, the observed H O 2 2 enhancement at low latitudes on the leading hemi-

Figure 2 . All mapped slits from our 2016 (A) and 2018 (B) observations, which demonstrate the spatial distribution of Europa's 3.5 µ m H O 2 2 absorption. The deepest absorptions map to low latitudes of the leading and anti-Jovian hemispheres and appear correlated with the geologically young chaos regions (outlined in black). The trailing/sub-Jovian slits, in comparison, display very weak absorptions, consistent with prior full-disk spectra (Hand &amp; Brown 2013).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01093v1-figure-938cd8a1-e481-4caa-9ddd-036fb25488a9.jpg)

sphere bears some resemblance to the leading hemisphere 'lens' pattern of bombardment by very highenergy ( glyph[greaterorsimilar] 20 MeV) electrons, which move against the co-rotation direction of Jupiter's magnetic field (Paranicas et al. 2009; Nordheim et al. 2018). This lens centers on the leading hemisphere apex, weakens toward the mid latitudes, and tapers toward the sub- and antiJovian points. However, laboratory irradiation experiments have yielded more efficient H 2 O production and 2 higher H O 2 2 concentrations using energetic protons and ions (Loeffler et al. 2006) than they have using electrons (Hand &amp; Carlson 2011). This is due to the higher linear energy transfer (LET) of ions (Loeffler et al. 2006; Hand &amp; Carlson 2011), which reflects a higher density of resulting excitations and ionizations in the ice. Although H O 2 2 production by 20+ MeV electrons has never been directly studied in the lab, the LET of electrons at these energies is still more than an order of magnitude less than that of the majority of energetic protons and ions striking the surface (Hand &amp; Carlson 2011). In contrast to the leading hemisphere electrons, energetic protons and ions bombard Europa much more uniformly (Paranicas 2002; Paranicas et al. 2009) and potentially even favor the high latitudes (Dalton et al. 2013). Indeed, the combined bombardment patterns of energetic ions and glyph[lessorsimilar] 20 MeV electrons, which primarily impact the trailing hemisphere (Paranicas et al. 2001), result in a relative sheltering of the leading-hemisphere equator in terms of total flux (Paranicas et al. 2009; Dalton et al. 2013). Thus, while the leading hemisphere electron lens may help produce the H O 2 2 at low latitudes, it is unlikely to be the sole explanation. Indeed, because the lens is expected to be symmetric across the equator, it certainly cannot account for the apparent correlation of H O 2 2 with chaos terrain.

The geographic association of H 2 O with chaos terrain 2

instead may imply that compositional differences across the surface are a dominating factor in controlling H 2 O 2 abundance. Chaos regions are geologically young and extensively disrupted areas of the surface interpreted to reflect recent interaction with subsurface material (e.g Collins &amp; Nimmo 2009; Schmidt et al. 2011). In fact, leading-hemisphere chaos regions have previously been shown to be compositionally distinct from their surroundings (Fischer et al. 2015, 2017; Trumbo et al. 2019), likely reflecting the contributions of endogenous material rich in chloride salts. To our knowledge, the effects of salty ice or frozen brines on the radiolytic formation of H 2 O 2 have not been studied, leaving us with no reason to expect H O 2 2 enhancements due to salt. However, the addition of electron-accepting contaminants, such as O 2 and CO , to the ice has been shown to en2 hance H O 2 2 yields (Moore &amp; Hudson 2000). In these cases, the CO 2 and O 2 molecules slow the destruction of newly formed H O 2 2 by scavenging destructive electrons created during the continued irradiation of the ice. In fact, it has been noted that very limited NIMS data of portions of the leading hemisphere show a correlation between Europa's surface CO 2 and its H O 2 2 (Carlson 2001; Carlson et al. 2009), although the spectra and corresponding maps were never published. Critically, these data cover portions of Tara Regio and Powys Regio, and demonstrate a clear enhancement of CO 2 in these chaos regions relative to the intervening plains (R. W. Carlson personal communication). Our NIRSPEC datasets map the same regions at the large scale. If we extrapolate the NIMS correlations, then the presence of CO 2 at low latitudes and in chaos terrain may help explain our observed H O 2 2 distribution. Perhaps more excitingly, if CO 2 is truly a dominant factor in controlling the H O 2 2 distribution on Europa, then the widespread coverage of our H 2 O data may hint at a broad associa2

tion of Europa's CO 2 with its geologically young chaos terrain, and therefore with an interior carbon source. However, as no comprehensive map of Europa's CO 2 was ever made, and as the corresponding spectra were never published, more complete investigations of the distribution and abundance of Europa's CO 2 are certainly needed to address this hypothesis.

While we suggest that the observed H O 2 2 distribution may reflect that of Europa's CO , it is worth ad2 dressing some caveats to this hypothesis. One potential complication is the apparent lack of significant H 2 O on 2 the trailing hemisphere, which contains extensive chaos terrain (Doggett et al. 2009) and potential spectral evidence for CO 2 (Hansen &amp; McCord 2008), although the spectra are poor quality. While water ice abundance does not appear to play a significant role in the latitudinal distribution of H O 2 2 on the leading and antiJovian hemispheres, it is possible that the more dramatic contrast in ice-fraction between the leading and trailing hemispheres accounts for their observed differences in H O 2 2 abundance. The widespread sulfur radiolysis experienced on the trailing hemisphere may also contribute to these differences. SO , a minor product 2 of the radiolytic sulfur cycle on the trailing hemisphere (Carlson et al. 2002), reacts efficiently with H 2 O in wa2 ter ice at temperatures relevant to Europa, rapidly producing sulfate (Loeffler &amp; Hudson 2013). Thus, while the composition of the leading hemisphere may enhance local H 2 O concentrations, the composition of the trail2 ing hemisphere may act to deplete them.

## 5. CONCLUSIONS

Using ground-based L-band spectroscopy, we have mapped Europa's 3.5 µ m H O 2 2 absorption across its surface. Although laboratory data would seem to suggest that Europa's H O 2 2 should be most abundant at the cold, icy high latitudes of the leading hemisphere, we find the strongest absorptions at low latitudes, preferentially within geologically young chaos terrain. We suggest that this distribution may largely reflect efficient electron scavenging by CO 2 within the chaos regions, which can enhance H O 2 2 concentrations. This hypothesis would benefit from continued mapping of Europa's H O 2 2 and CO 2 and from further laboratory experiments involving ion and electron irradiation of ice mixtures potentially reflective of the composition of leading-hemisphere chaos terrain.

Lastly, we would like to acknowledge the possibility that the 3.5 µ m H O 2 2 band areas may not directly reflect relative H 2 O 2 abundances, because they may also include the effects of geographically varying grain sizes. With only a single absorption band, it is impossible to disentangle the two effects. However, laboratory data suggest that H O 2 2 on Europa likely exists as isolated trimers (H O 2 2 · 2H O) within the surface ice (Loeffler 2 &amp; Baragiola 2005). It has been suggested that the ice grain size on Europa increases with latitude (Carlson et al. 2009; Cassidy et al. 2013), which should produce the opposite latitudinal trend in band area from what we observe. Thus, we feel that the 3.5 µ mband of H O 2 2 likely represents a good proxy for its abundance.

This research was supported by Grant 1313461 from the National Science Foundation. This work was also supported by NASA Headquarters under the NASA Earth and Space Science Fellowship Program Grant 80NSSC17K0478. K. P. H. acknowledges support from the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration and funded in part through the internal Research and Technology Development program. The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. We thank Robert W. Carlson for sharing enlightening details concerning Galileo NIMS observations of CO 2 on Europa.

Software: Astropy (Robitaille et al. 2013), skimage.transform (van der Walt et al. 2014)

## REFERENCES

Brown, M. E., &amp; Hand, K. P. 2013, The Astronomical Journal, 145, 110

Carlson, R. W., Calvine, W. M., Dalton, J. B., et al. 2009,

Carlson, R. W. 2001, in Bulletin of the American Astronomical

Society, Vol. 33, 1125

Carlson, R. W., Anderson, M. S., Johnson, R. E., Schulman,

M. B., &amp; Yavrouian, A. H. 2002, Icarus, 157, 456

```
<_Python_>,
```

- Cassidy, T., Paranicas, C., Shirley, J., et al. 2013, Planetary and Space Science, 77, 64 , surfaces, atmospheres and magnetospheres of the outer planets and their satellites and ring systems: Part VIII
- Chyba, C. F. 2000, Nature, 403, 381
- Collins, G., &amp; Nimmo, F. 2009, Europa, ed. R. T. Pappalardo, W. B. McKinnon, &amp; K. Khurana (Tucson, AZ: The University of Arizona Press), 259-282
- Cooper, P. D., Johnson, R. E., &amp; Quickenden, T. I. 2003, Icarus, 166, 444
- Dalton, J. B., Cassidy, T., Paranicas, C., et al. 2013, Planetary and Space Science, 77, 45
- Doggett, T., Greeley, R., Figueredo, P., &amp; Tanaka, K. 2009, Europa, ed. R. T. Pappalardo, W. B. McKinnon, &amp; K. Khurana (Tucson, AZ: University of Arizona Press), 137-160
- Fischer, P. D., Brown, M. E., &amp; Hand, K. P. 2015, The Astronomical Journal, 150, 164
- Fischer, P. D., Brown, M. E., Trumbo, S. K., &amp; Hand, K. P. 2017, The Astronomical Journal, 153, 13
- Hand, K. P., &amp; Brown, M. E. 2013, The Astrophysical Journal, 766, L21
- Hand, K. P., &amp; Carlson, R. W. 2011, Icarus, 215, 226 Hand, K. P., Chyba, C. F., Priscu, J. C., Carlson, R. W., &amp; Nealson, K. H. 2009, Europa, ed. R. T. Pappalardo, W. B. McKinnon, &amp; K. Khurana (Tucson, AZ: The University of Arizona Press), 589-630
- Hansen, G. B., &amp; McCord, T. B. 2008, Geophysical Research Letters, 35, doi:10.1029/2007gl031748
- Johnson, R. E., &amp; Quickenden, T. I. 1997, Journal of Geophysical Research: Planets, 102, 10985
- Johnson, R. E., Quickenden, T. I., Cooper, P. D., McKinley, A. J., &amp; Freeman, C. G. 2003, Astrobiology, 3, 823, pMID: 14987485
- Loeffler, M., Raut, U., Vidal, R., Baragiola, R., &amp; Carlson, R. 2006, Icarus, 180, 265
- Loeffler, M. J., &amp; Baragiola, R. A. 2005, Geophysical Research Letters, 32, L17202
- Loeffler, M. J., &amp; Hudson, R. L. 2013, Icarus, 224, 257 Lord, S. D. 1992, NASA Technical Memorandum 103957, Ames Research Center, Moffett Field, CA
- Moore, M., &amp; Hudson, R. 2000, Icarus, 145, 282
- Nordheim, T. A., Hand, K. P., &amp; Paranicas, C. 2018, Nature Astronomy, 2, 673
- Paranicas, C. 2002, Geophysical Research Letters, 29, doi:10.1029/2001gl014127
- Paranicas, C., Carlson, R. W., &amp; Johnson, R. E. 2001, Geophysical Research Letters, 28, 673
- Paranicas, C., Cooper, J. F., Garrett, H. B., Johnson, R. E., &amp; Sturner, S. J. 2009, Europa, ed. R. T. Pappalardo, W. B. McKinnon, &amp; K. Khurana (Tucson, AZ: The University of Arizona Press), 529-544
- Robitaille, T. P., Tollerud, E. J., Greenfield, P., et al. 2013, Astronomy and Astrophysics, 558, 9
- Schmidt, B. E., Blankenship, D. D., Patterson, G. W., &amp; Schenk, P. M. 2011, Nature, 479, 502
- Spencer, J. R., Tamppari, L. K., Martin, T. Z., &amp; Travis, L. D. 1999, Science, 284, 1514
- Trumbo, S. K., Brown, M. E., &amp; Hand, K. P. 2019, Science Advances, 5, doi:10.1126/sciadv.aaw7123
- van der Walt, S., Sch¨nberger, J. L., Nunez-Iglesias, J., et al. o 2014, PeerJ, 2
- Zheng, W., Jewitt, D., &amp; Kaiser, R. I. 2006, The Astrophysical Journal, 648, 753
Samantha K. Trumbo Michael E. Brown Kevin P. Hand
H$_2$O$_2$ within chaos terrain on Europa's leading hemisphere
Hydrogen peroxide is part of an important radiolytic cycle on Europa and may
be a critical source of oxidants to the putative subsurface ocean. The surface
geographic distribution of hydrogen peroxide may constrain the processes
governing its abundance as well as its potential relevance to the subsurface
chemistry. However, maps of Europa's hydrogen peroxide beyond hemispherical
averages have never been published. Here, we present spatially resolved L-band
(3.16 - 4 $\mu$m) observations of Europa's 3.5 $\mu$m hydrogen peroxide
absorption, which we obtained using the near-infrared spectrometer NIRSPEC and
the adaptive optics system on the Keck II telescope. Using these data, we map
the strength of the 3.5 $\mu$m absorption across the surface at a nominal
spatial resolution of $\sim$300 km. Though previous disk-integrated data seemed
consistent with the laboratory expectation that Europa's hydrogen peroxide
exists primarily in its coldest and iciest regions, we find nearly the exact
opposite at this finer spatial scale. Instead, we observe the largest hydrogen
peroxide absorptions at low latitudes on the leading and anti-Jovian
hemispheres, correlated with chaos terrain, and relative depletions toward the
cold, icy high latitudes. This distribution may reflect the effects of
decreased hydrogen peroxide destruction due to efficient electron scavenging by
CO$_2$ within chaos terrain.