mass wasting on the moon: implications for seismicity · mass wasting on the moon: implications for...

1
Mass wasting on the Moon: Implications for seismicity RENEE WEBER - NASA MARSHALL SPACE FLIGHT CENTER AMANDA NAHM - DLR INSTITUTE FOR PLANETARY RESEARCH NICK SCHMERR - UNIV. OF MARYLAND BRIAN YANITES - UNIV. OF IDAHO Overview Seismicity estimates play an important role in creating regional geological characterizations, which are useful for understanding a planet’s formation and evolution, and are of key importance to site selection for landed missions. Here we investigate the regional effects of seismicity in planetary environments with the goal of determining whether such surface features on the Moon, could be triggered by fault motion (Fig. 1). Fig. 1: (left) Landslide deposits (granular flow) on an interior slope of Marius crater (11.9° N, -50.8° E). (right) Boulder track emanating from the central peak complex of Schiller crater (-51.8° N, -40.0° E). Lobate scarps, the typical surface expressions of thrust faults resulting from tectonic compression, are widely observed on the Moon (Figs. 2&3). Compared to other types of faults, surface-cutting thrust faults require the largest amount of stress to form and/or slip, so they could possibly generate large quakes. While normal faults, graben, and wrinkle ridges may be more abundant on Mars, the Moon, and Mercury respectively, these struc- tures would create smaller theoretical maximum quakes than lobate scarp thrust faults. Thus, we optimize our chances of finding mass wasting associated with faults by studying lobate scarps. Lobate scarps Fig. 2: Examples of lobate scarps Evershed S1 center lat/lon 33°N/197.1°E Utopia Planitia #s 1801, 1802, 1804 center lat/lon 52.9°N/119.2°E Beagle Rupes center lat/lon -3.5°N/100.7°E Fig. 3: Map of digitized locations of lobate scarps on the Moon. Over 3200 lunar scarps (red) have been mapped. Mare basalt units are shown in tan. From Watters et al., 2015. (Geology Vol. 43 No. 10) Wavefield modeling 0 50 100 0 50 100 0 1 2 3 4 5 6 0 50 100 0 50 100 0 1 2 LOG 10 [Amplitude] (cm/s) Elevation (km) Distance (km) Distance (km) Distance (km) Topography Maximum Vertical Velocity Fig. 5: (left) LOLA surface topography input into the wave propagation simulation. The Evershed scarp is centered in the image (see Fig. 2). (right) Predicted maximum vertical ground motion for a M W =7.8 quake on a subjacent reverse fault, with a 2.25 km depth of faulting. The surface trace of the scarp is marked with a red line. A random distribution of heterogeneity of 25% in S- and P-wave velocity with 100 km scale length scatterers is placed in the megaregolith to simulate the scattering observed on Apollo seismograms. The expected damage area indi- cated by seismic wavefield modeling is compared to mapped imagery to determine the likeli- hood of a quake having triggered mass wasting. Future work Map instances of mass wasting around faults Rule out impact- generated shaking 100m Methodology Following the method outlined in Nahm & Velasco, 2013 (LPSC 44th Abstract #1422), we derive a theoretical quake magnitude from basic fault properties. These are esti- mated either from imagery, laboratory rock experiments, or elastic dislocation models, and include the length (L), dip angle (δ), depth of faulting (T), displacement (D), and fault width (w). Fault displacement is calculated using displacement-length scaling such that D=γL, where γ is determined by rock type and tectonic setting. To determine the dimensions of an area affected by seismic shaking, we model the ground motion resulting from the theoretical maximum quake along a given fault (Figs. 4&5). We use a numerical code for simulating seismic wave propagation through a 3-D structure model including topography. Peak vertical ground motion typi- cally occurs within a few kilometers of the main shock and drops off rapidly from there. This implies that we should expect most of the mass wasting phenomena to occur in the immediate vicinity of the fault. However, this result may depend on regional effects like surface slope and megaregolith thickness. A thicker megaregolith (as might be expected in the vicinity of craters) would tend to focus shaking in some of the crater basins. Fig. 4: Seismic structure of the lunar crust and megaregolith. (left) The crustal velocity model after Weber et al., 2011 (Science Vol. 331), showing the geology of the subsurface. The megaregolith extends from the surface to just above the structurally disturbed crust. (right) Wave propagation for a model with a 1 km thick layer of low-velocity megaregolith showing the development of surface waves and reverberations trapped in the layer. From Weber et al., 2015 (Extraterrestrial Seismology, Cambridge Univ. Press). 75 150 225 300 450 600 750 900 1000 Epicentral Distance (km) Time(s) 0 500 1000 1500 2000 2500 3000 3500 REGOLITH EJECTA FRACTURED CRUST STRUCTURALLY DISTURBED CRUST LUNAR CRUST LUNAR MANTLE V p (km/s) 0 2 15 38 1 3 5 7 Depth (km) V s (km/s) Geomorphological analysis: Establish a method to translate quake parameters into mass wasting estimates e.g. quake magnitude ~ area of damage zone? e.g. fault length ~ boulder trail density? https://ntrs.nasa.gov/search.jsp?R=20160006721 2018-08-27T06:34:47+00:00Z

Upload: trantruc

Post on 27-Aug-2018

212 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Mass wasting on the Moon: Implications for seismicity · Mass wasting on the Moon: Implications for seismicity RENEE WEBER - NASA MARSHALL SPACE FLIGHT CENTER AMANDA NAHM - …

Mass wasting on the Moon: Implications for seismicityRENEE WEBER - NASA MARSHALL SPACE FLIGHT CENTER AMANDA NAHM - DLR INSTITUTE FOR PLANETARY RESEARCH

NICK SCHMERR - UNIV. OF MARYLAND BRIAN YANITES - UNIV. OF IDAHO

Overview

Seismicity estimates play an important role in creating regional geological characterizations, which are useful for understanding a planet’s formation and evolution, and are of key importance to site selection for landed missions. Here we investigate the regional e�ects of seismicity in planetary environments with the goal of determining whether such surface features on the Moon, could be triggered by fault motion (Fig. 1).

Fig. 1: (left) Landslide deposits (granular �ow) on an interior slope of Marius crater (11.9° N, -50.8° E).

(right) Boulder track emanating from the central peak complex of Schiller crater (-51.8° N, -40.0° E).

Lobate scarps, the typical surface expressions of thrust faults resulting from tectonic compression, are widely observed on the Moon (Figs. 2&3). Compared to other types of faults, surface-cutting thrust faults require the largest amount of stress to form and/or slip, so they could possibly generate large quakes. While normal faults, graben, and wrinkle ridges may be more abundant on Mars, the Moon, and Mercury respectively, these struc-tures would create smaller theoretical maximum quakes than lobate scarp thrust faults. Thus, we optimize our chances of �nding mass wasting associated with faults by studying lobate scarps.

Lobate scarps

Fig. 2: Examples of lobate scarps

Evershed S1center lat/lon 33°N/197.1°E

Utopia Planitia #s 1801, 1802, 1804 center lat/lon52.9°N/119.2°E

Beagle Rupescenter lat/lon

-3.5°N/100.7°E

Fig. 3: Map of digitized locations of lobate scarps on the Moon. Over 3200 lunar scarps (red) have been mapped. Mare basalt units are shown in tan. From Watters et al., 2015. (Geology Vol. 43 No. 10)

Wave�eld modeling

0

50

100

0 50 100

0123456

0

50

100

0 50 100

0

1

2

LOG

10 [Am

plitude] (cm/s)

Elevation (km

)Dis

tanc

e (k

m)

Distance (km) Distance (km)

Topography Maximum Vertical Velocity

Fig. 5: (left) LOLA surface topography input into the wave propagation simulation. The Evershed scarp is centered in the image (see Fig. 2). (right) Predicted maximum vertical ground motion for a MW=7.8 quake on a subjacent reverse fault, with a 2.25 km depth of faulting. The surface trace of the scarp is marked with a red line. A random distribution of heterogeneity of 25% in S- and P-wave velocity with 100 km scale length scatterers is placed in the megaregolith to simulate the scattering observed on Apollo seismograms. The expected damage area indi-cated by seismic wave�eld modeling is compared to mapped imagery to determine the likeli-hood of a quake having triggered mass wasting.

Future work

Map instances of masswasting around faults

Rule outimpact-

generatedshaking

100m

Methodology

Following the method outlined in Nahm & Velasco, 2013 (LPSC 44th Abstract #1422), we derive a theoretical quake magnitude from basic fault properties. These are esti-mated either from imagery, laboratory rock experiments, or elastic dislocation models, and include the length (L), dip angle (δ), depth of faulting (T), displacement (D), and fault width (w). Fault displacement is calculated using displacement-length scaling such that D=γL, where γ is determined by rock type and tectonic setting.

To determine the dimensions of an area a�ected by seismic shaking, we model the ground motion resulting from the theoretical maximum quake along a given fault (Figs. 4&5). We use a numerical code for simulating seismic wave propagation through a 3-D structure model including topography. Peak vertical ground motion typi-cally occurs within a few kilometers of the main shock and drops o� rapidly from there. This implies that we should expect most of the mass wasting phenomena to occur in the immediate vicinity of the fault. However, this result may depend on regional e�ects like surface slope and megaregolith thickness. A thicker megaregolith (as might be expected in the vicinity of craters) would tend to focus shaking in some of the crater basins.

Fig. 4: Seismic structure of the lunar crust and megaregolith. (left) The crustal velocity model after Weber et al., 2011 (Science Vol. 331), showing the geology of the subsurface. The megaregolith extends from the surface to just above the structurally disturbed crust. (right) Wave propagation for a model with a 1 km thick layer of low-velocity megaregolith showing the development of surface waves and reverberations trapped in the layer. From Weber et al., 2015 (Extraterrestrial Seismology, Cambridge Univ. Press).

75

150

225

300

450

600

750

900

1000

Epic

entra

l Dis

tanc

e (k

m)

Time(s)0 500 1000 1500 2000 2500 3000 3500

REGOLITH

EJECTA

FRACTURED CRUST

STRUCTURALLY DISTURBED CRUST

LUNAR CRUST

LUNAR MANTLE

Vp (km/s)

0

2

15

38

1 3 5 7

Dep

th (k

m)

Vs (km/s)

Geomorphological analysis:Establish a method to translate

quake parameters into mass wasting estimates

e.g. quake magnitude

~area of damage

zone?

e.g.fault length

~boulder trail

density?

https://ntrs.nasa.gov/search.jsp?R=20160006721 2018-08-27T06:34:47+00:00Z