Editorial: Exploring new frontiers with JGR–Solid Earth
No abstract available.
Journal of Geophysical Research B: Solid Earth↗
SEARCH · Geology Reports
Results for “Geophysical Journal”
Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Anisotropy, repeating earthquakes, and seismicity associated with the 2008 eruption of Okmok Volcano, Alaska
We use shear wave splitting (SWS) analysis and double-difference relocation to examine temporal variations in seismic properties prior to and accompanying magmatic activity associated with the 2008 eruption of Okmok volcano, Alaska. Using bispectrum cross-correlation, a multiplet of 25 earthquakes is identified spanning five years leading up to the eruption, each event having first motions compatible with a normal fault striking NE–SW. Cross-correlation differential times are used to relocate earthquakes occurring between January 2003 and February 2009. The bulk of the seismicity prior to the onset of the eruption on 12 July 2008 occurred southwest of the caldera beneath a geothermal field. Earthquakes associated with the onset of the eruption occurred beneath the northern portion of the caldera and started as deep as 13 km. Subsequent earthquakes occurred predominantly at 3 km depth, coinciding with the depth at which the magma body has been modeled using geodetic data. Automated SWS analysis of the Okmok catalog reveals radial polarization outside the caldera and a northwest-southeast polarization within. We interpret these polarizations in terms of a magma reservoir near the center of the caldera, which we model with a Mogi point source. SWS analysis using the same input processing parameters for each event in the multiplet reveals no temporal changes in anisotropy over the duration of the multiplet, suggesting either a short-term or small increase in stress just before the eruption that was not detected by GPS, or eruption triggering by a mechanism other than a change of stress in the system.
Alaska↗
Lessons from (triggered) tremor
I test a “clock-advance” model that implies triggered tremor is ambient tremor that occurs at a sped-up rate as a result of loading from passing seismic waves. This proposed model predicts that triggering probability is proportional to the product of the ambient tremor rate and a function describing the efficacy of the triggering wave to initiate a tremor event. Using data mostly from Cascadia, I have compared qualitatively a suite of teleseismic waves that did and did not trigger tremor with ambient tremor rates. Many of the observations are consistent with the model if the efficacy of the triggering wave depends on wave amplitude. One triggered tremor observation clearly violates the clock-advance model. The model prediction that larger triggering waves result in larger triggered tremor signals also appears inconsistent with the measurements. I conclude that the tremor source process is a more complex system than that described by the clock-advance model predictions tested. Results of this and previous studies also demonstrate that (1) conditions suitable for tremor generation exist in many tectonic environments, but, within each, only occur at particular spots whose locations change with time; (2) any fluid flow must be restricted to less than a meter; (3) the degree to which delayed failure and secondary triggering occurs is likely insignificant; and 4) both shear and dilatational deformations may trigger tremor. Triggered and ambient tremor rates correlate more strongly with stress than stressing rate, suggesting tremor sources result from time-dependent weakening processes rather than simple Coulomb failure.
Cascadia↗
Volcano monitoring using GPS: Developing data analysis strategies based on the June 2007 Kīlauea Volcano intrusion and eruption
The global positioning system (GPS) is one of the most common techniques, and the current state of the art, used to monitor volcano deformation. In addition to slow (several centimeters per year) displacement rates, GPS can be used to study eruptions and intrusions that result in much larger (tens of centimeters over hours-days) displacements. It is challenging to resolve precise positions using GPS at subdaily time intervals because of error sources such as multipath and atmospheric refraction. In this paper, the impact of errors due to multipath and atmospheric refraction at subdaily periods is examined using data from the GPS network on Kīlauea Volcano, Hawai'i. Methods for filtering position estimates to enhance precision are both simulated and tested on data collected during the June 2007 intrusion and eruption. Comparisons with tiltmeter records show that GPS instruments can precisely recover the timing of the activity.
Hawai'i↗
Cyclic ground tilt associated with the 2004–2008 eruption of Mount St. Helens
The long‐term behavior of the 2004–2008 effusive eruption of Mount St. Helens was characterized by a gradual decline in the rates of seismicity, dome growth, and broad‐scale ground deformation, but shallow near‐periodic “drumbeat” earthquakes over timescales of minutes indicated episodic short‐term behavior. In part to better characterize this behavior and any associated ground deformation, a network of tiltmeters was installed and recorded thousands of cyclic tilt events within the crater. The duration of these events, from minutes to hours, was too long for them to be directly caused by the drumbeat seismicity. Tilt events were characterized by a recoverable, asymmetric pattern of rapid tilt away from the vent followed by a more gradual reversal, were highly correlated between different tiltmeters in the crater, sometimes occurred in association with volcanic seismicity, and ceased at the end of the eruption. Tilt vectors converged on a point just south of the center of the preexisting 1980s lava dome, and the absence of detectable tilt outside the crater suggests a shallow source (<1 km). We examine several models, including cycles of conduit pressurization and plug slip or gas loss and stick‐slip behavior on the interface between the ascending plug and the 1980s lava dome. The small number of stations within the crater prevents a unique determination of source type or geometry, but results are consistent with a mechanism involving extrusion of the semisolid dacite plug and/or cycles of conduit pressurization.
Washington↗
Coulomb stress interactions among M≥5.9 earthquakes in the Gorda deformation zone and on the Mendocino Fracture Zone, Cascadia megathrust, and northern San Andreas fault
The Gorda deformation zone, a 50,000 km 2 area of diffuse shear and rotation offshore northernmost California, has been the site of 20 M ≥ 5.9 earthquakes on four different fault orientations since 1976, including four M ≥ 7 shocks. This is the highest rate of large earthquakes in the contiguous United States. We calculate that the source faults of six recent M ≥ 5.9 earthquakes had experienced ≥0.6 bar Coulomb stress increases imparted by earthquakes that struck less than 9 months beforehand. Control tests indicate that ≥0.6 bar Coulomb stress interactions between M ≥ 5.9 earthquakes separated by <9 months are unlikely to occur by random chance, suggesting that the multiple short-term stress interactions observed among the recent Gorda zone earthquakes are not an apparent effect. In all well-constrained ≥0.2 bar Coulomb stress interactions between earthquakes that occurred within 4 years of each other, the second earthquake is promoted. On longer timescales, calculated stress changes imparted by the 1980 M w = 7.3 Trinidad earthquake are consistent with the locations of M ≥ 5.9 earthquakes in the Gorda zone until at least 1995, as well as earthquakes on the Mendocino Fault Zone in 1994 and 2000. Coulomb stress changes imparted by the 1980 earthquake are also consistent with its distinct elbow-shaped aftershock pattern. From these observations, we derive generalized static stress interactions among right-lateral, left-lateral and thrust faults near triple junctions.
California↗
Origins of large-volume, compositionally zoned volcanic eruptions: New constraints from U-series isotopes and numerical thermal modeling for the 1912 Katmai-Novarupta eruption
We present the results of a combined U-series isotope and numerical modeling study of the 1912 Katmai-Novarupta eruption in Alaska. A stratigraphically constrained set of samples have compositions that range from basalt through basaltic andesite, andesite, dacite, and rhyolite. The major and trace element range can be modeled by 80–90% closed-system crystal fractionation over a temperature interval from 1279°C to 719°C at 100 MPa, with an implied volume of parental basalt of ∼65 km 3 . Numerical models suggest, for wall rock temperatures appropriate to this depth, that 90% of this volume of magma would cool and crystallize over this temperature interval within a few tens of kiloyears. However, the range in 87 Sr/ 86 Sr, ( 230 Th/ 238 U), and ( 226 Ra/ 230 Th) requires open-system processes. Assimilation of the host sediments can replicate the range of Sr isotopes. The variation of ( 226 Ra/ 230 Th) ratios in the basalt to andesite compositional range requires that these were generated less than several thousand years before eruption. Residence times for dacites are close to 8000 years, whereas the rhyolites appear to be 50–200 kyr old. Thus, the magmas that erupted within only 60 h had a wide range of crustal residence times. Nevertheless, they were emplaced in the same thermal regime and evolved along similar liquid lines of descent from parental magmas with similar compositions. The system was built progressively with multiple inputs providing both mass and heat, some of which led to thawing of older silicic material that provided much of the rhyolite.
Alaska↗
Variability and trends in dry day frequency and dry event length in the southwestern United States
Daily precipitation from 22 National Weather Service first-order weather stations in the southwestern United States for water years 1951 through 2006 are used to examine variability and trends in the frequency of dry days and dry event length. Dry events with minimum thresholds of 10 and 20 consecutive days of precipitation with less than 2.54 mm are analyzed. For water years and cool seasons (October through March), most sites indicate negative trends in dry event length (i.e., dry event durations are becoming shorter). For the warm season (April through September), most sites also indicate negative trends; however, more sites indicate positive trends in dry event length for the warm season than for water years or cool seasons. The larger number of sites indicating positive trends in dry event length during the warm season is due to a series of dry warm seasons near the end of the 20th century and the beginning of the 21st century. Overall, a large portion of the variability in dry event length is attributable to variability of the El Niño–Southern Oscillation, especially for water years and cool seasons. Our results are consistent with analyses of trends in discharge for sites in the southwestern United States, an increased frequency in El Niño events, and positive trends in precipitation in the southwestern United States.
Journal of Geophysical Research D: Atmospheres↗
Identification of hydrated silicate minerals on Mars using MRO-CRISM: Geologic context near Nili Fossae and implications for aqueous alteration
[1] The Noachian terrain west of the Isidis basin hosts a diverse collection of alteration minerals in rocks comprising varied geomorphic units within a 100,000 km 2 region in and near the Nili Fossae. Prior investigations in this region by the Observatoire pour l'Minéralogie, l'Eau, les Glaces, et l'Activité (OMEGA) instrument on Mars Express revealed large exposures of both mafic minerals and iron magnesium phyllosilicates in stratigraphic context. Expanding on the discoveries of OMEGA, the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) aboard the Mars Reconnaissance Orbiter (MRO) has found more spatially widespread and mineralogically diverse alteration minerals than previously realized, which represent multiple aqueous environments. Using CRISM near-infrared spectral data, we detail the basis for identification of iron and magnesium smectites (including both nontronite and more Mg-rich varieties), chlorite, prehnite, serpentine, kaolinite, potassium mica (illite or muscovite), hydrated (opaline) silica, the sodium zeolite analcime, and magnesium carbonate. The detection of serpentine and analcime on Mars is reported here for the first time. We detail the geomorphic context of these minerals using data from high-resolution imagers onboard MRO in conjunction with CRISM. We find that the distribution of alteration minerals is not homogeneous; rather, they occur in provinces with distinctive assemblages of alteration minerals. Key findings are (1) a distinctive stratigraphy, in and around the Nili Fossae, of kaolinite and magnesium carbonate in bedrock units always overlying Fe/Mg smectites and (2) evidence for mineral phases and assemblages indicative of low-grade metamorphic or hydrothermal aqueous alteration in cratered terrains. The alteration minerals around the Nili Fossae are more typical of those resulting from neutral to alkaline conditions rather than acidic conditions, which appear to have dominated much of Mars. Moreover, the mineralogic diversity and geologic context of alteration minerals found in the region around the Nili Fossae indicates several episodes of aqueous activity in multiple distinct environments.
Journal of Geophysical Research E: Planets↗
Seasonal H2O and CO2 ice cycles at the Mars Phoenix landing site: 1. Prelanding CRISM and HiRISE observations
The condensation, evolution, and sublimation of seasonal water and carbon dioxide ices were characterized at the Mars Phoenix landing site from Martian northern midsummer to midspring (L s ∼ 142° – L s ∼ 60°) for the year prior to the Phoenix landing on 25 May 2008. Ice relative abundances and grain sizes were estimated using data from the Compact Reconnaissance Imaging Spectrometer for Mars and High Resolution Imaging Science Experiment aboard Mars Reconnaissance Orbiter and a nonlinear mixing model. Water ice first appeared at the Phoenix landing site during the afternoon in late summer (L s ∼ 167°) as an optically thin layer on top of soil. CO 2 ice appeared after the fall equinox. By late winter (L s ∼ 344°), the site was covered by relatively pure CO 2 ice (∼30 cm thick), with a small amount of ∼100 μ m diameter water ice and soil. As spring progressed, CO 2 ice grain sizes gradually decreased, a change interpreted to result from granulation during sublimation losses. The combined effect of CO 2 sublimation and decreasing H 2 O ice grain sizes allowed H 2 O ice to dominate spectra during the spring and significantly brightened the surface. CO 2 ice disappeared by early spring (L s ∼ 34°) and H 2 O ice by midspring (L s ∼ 59°). Spring defrosting was not uniform and occurred more rapidly over the centers of polygons and geomorphic units with relatively higher thermal inertia values.
Journal of Geophysical Research E: Planets↗
Mineralogy and stratigraphy of phyllosilicate-bearing and dark mantling units in the greater Mawrth Vallis/west Arabia Terra area: Constraints on geological origin
Analyses of MRO/CRISM images of the greater Mawrth Vallis region of Mars affirm the presence of two primary phyllosilicate assemblages throughout a region ∼1000 × 1000 km. These two units consist of an Fe/Mg-phyllosilicate assemblage overlain by an Al-phyllosilicate and hydrated silica assemblage. The lower unit contains Fe/Mg-smectites, sometimes combined with one or more of these other Fe/Mg-phyllosilicates: serpentine, chlorite, biotite, and/or vermiculite. It is more than 100 m thick and finely layered at meter scales. The upper unit includes Al-smectite, kaolin group minerals, and hydrated silica. It is tens of meters thick and finely layered as well. A common phyllosilicate stratigraphy and morphology is observed throughout the greater region wherever erosional windows are present. This suggests that the geologic processes forming these units must have occurred on at least a regional scale. Sinuous ridges (interpreted to be inverted channels) and narrow channels cut into the upper clay-bearing unit suggesting that aqueous processes were prevalent after, and possibly during, the deposition of the layered units. We propose that layered units may have been deposited at Mawrth Vallis and then subsequently altered to form the hydrated units. The Fe/Mg-phyllosilicate assemblage is consistent with hydrothermal alteration or pedogenesis of mafic to ultramafic rocks. The Al-phyllosilicate/hydrated silica unit may have formed through alteration of felsic material or via leaching of basaltic material through pedogenic alteration or a mildly acidic environment. These phyllosilicate-bearing units are overlain by a darker, relatively unaltered, and indurated material that has probably experienced a complex geological history.
Journal of Geophysical Research↗
Detection and mapping of hydrocarbon deposits on Titan
We report the identification of compounds on Titan's surface by spatially resolved imaging spectroscopy methods through Titan's atmosphere, and set upper limits to other organic compounds. We present evidence for surface deposits of solid benzene (C 6 H 6 ), solid and/or liquid ethane (C 2 H 6 ), or methane (CH 4 ), and clouds of hydrogen cyanide (HCN) aerosols using diagnostic spectral features in data from the Cassini Visual and Infrared Mapping Spectrometer (VIMS). Cyanoacetylene (2-propynenitrile, IUPAC nomenclature, HC 3 N) is indicated in spectra of some bright regions, but the spectral resolution of VIMS is insufficient to make a unique identification although it is a closer match to the feature previously attributed to CO 2 . We identify benzene, an aromatic hydrocarbon, in larger abundances than expected by some models. Acetylene (C 2 H 2 ), expected to be more abundant on Titan according to some models than benzene, is not detected. Solid acetonitrile (CH 3 CN) or other nitriles might be candidates for matching other spectral features in some Titan spectra. An as yet unidentified absorption at 5.01- μ m indicates that yet another compound exists on Titan's surface. We place upper limits for liquid methane and ethane in some locations on Titan and find local areas consistent with millimeter path lengths. Except for potential lakes in the southern and northern polar regions, most of Titan appears “dry.” Finally, we find there is little evidence for exposed water ice on the surface. Water ice, if present, must be covered with organic compounds to the depth probed by 1–5- μ m photons: a few millimeters to centimeters.
Journal of Geophysical Research - Planets↗
A comparison of Martian north and south polar cold spots and the long‐term effects of the 2001 global dust storm
In the 1970s, Viking and Mariner observed areas in the polar regions of Mars with winter brightness temperatures below the expected kinetic temperatures for CO 2 ice sublimation. These areas have since been termed “cold spots” and have been identified as surface deposits of CO 2 atmospheric condensates and, occasionally, active CO 2 storms. Three Mars years of data from the Mars Global Surveyor Thermal Emission Spectrometer were used to observe autumn and winter cold spot activity. In this study, cold spots that occur near and on the southern perennial cap were compared to those found near or on the northern perennial cap. On the southern perennial cap, cold spots associated with topographic features (induced by orographic lifting) were less common than cold spots independent of topography, similar to the north. However, the cold spots in the south lasted longer than those observed in the north. There is also evidence that cold spot formation in the south was affected by the global dust storm of 2001, even though the dust storm occurred during the southern spring and summer seasons. Prior to the dust storm, the amount of overall cold spot activity closer to the perennial cap increased and the average CO 2 grain size for most of the cold spots increased as well. Following the dust storm, the majority of cold spots in the south increased in size and duration but they did not form north of 62°S latitude, whereas, in other years, cold spots formed as far north as 48°S.
Journal of Geophysical Research E: Planets↗
Crater population and resurfacing of the Martian north polar layered deposits
Present‐day accumulation in the north polar layered deposits (NPLD) is thought to occur via deposition on the north polar residual cap. Understanding current mass balance in relation to current climate would provide insight into the climatic record of the NPLD. To constrain processes and rates of NPLD resurfacing, a search for craters was conducted using images from the Mars Reconnaissance Orbiter Context Camera. One hundred thirty craters have been identified on the NPLD, 95 of which are located within a region defined to represent recent accumulation. High Resolution Imaging Science Experiment images of craters in this region reveal a morphological sequence of crater degradation that provides a qualitative understanding of processes involved in crater removal. A classification system for these craters was developed based on the amount of apparent degradation and infilling and where possible depth/diameter ratios were determined. The temporal and spatial distribution of crater degradation is interpreted to be close to uniform. Through comparison of the size‐frequency distribution of these craters with the expected production function, the craters are interpreted to be an equilibrium population with a crater of diameter D meters having a lifetime of ∼30.75 D 1.14 years. Accumulation rates within these craters are estimated at 7.2 D −0.14 mm/yr, which corresponds to values of ∼3–4 mm/yr and are much higher than rates thought to apply to the surrounding flat terrain. The current crater population is estimated to have accumulated in the last ∼20 kyr or less.
Journal of Geophysical Research E: Planets↗
Bathymetry and absorptivity of Titan's Ontario Lacus
Ontario Lacus is the largest and best characterized lake in Titan's south polar region. In June and July 2009, the Cassini RADAR acquired its first Synthetic Aperture Radar (SAR) images of the area. Together with closest approach altimetry acquired in December 2008, these observations provide a unique opportunity to study the lake's nearshore bathymetry and complex refractive properties. Average radar backscatter is observed to decrease exponentially with distance from the local shoreline. This behavior is consistent with attenuation through a deepening layer of liquid and, if local topography is known, can be used to derive absorptive dielectric properties. Accordingly, we estimate nearshore topography from a radar altimetry profile that intersects the shoreline on the East and West sides of the lake. We then analyze SAR backscatter in these regions to determine the imaginary component of the liquid's complex index of refraction (κ). The derived value, κ = (6.1 −1.3 +1.7 ) × 10 −4 , corresponds to a loss tangent of tan Δ = (9.2 −2.0 +2.5 ) × 10 −4 and is consistent with a composition dominated by liquid hydrocarbons. This value can be used to test compositional models once the microwave optical properties of candidate materials have been measured. In areas that do not intersect altimetry profiles, relative slopes can be calculated assuming the index of refraction is constant throughout the liquid. Accordingly, we construct a coarse bathymetry map for the nearshore region by measuring bathymetric slopes for eleven additional areas around the lake. These slopes vary by a factor of ∼5 and correlate well with observed shoreline morphologies.
Journal of Geophysical Research E: Planets↗
Complexities in barrier island response to sea level rise: Insights from numerical model experiments, North Carolina Outer Banks
Using a morphological-behavior model to conduct sensitivity experiments, we investigate the sea level rise response of a complex coastal environment to changes in a variety of factors. Experiments reveal that substrate composition, followed in rank order by substrate slope, sea level rise rate, and sediment supply rate, are the most important factors in determining barrier island response to sea level rise. We find that geomorphic threshold crossing, defined as a change in state (e.g., from landward migrating to drowning) that is irreversible over decadal to millennial time scales, is most likely to occur in muddy coastal systems where the combination of substrate composition, depth-dependent limitations on shoreface response rates, and substrate erodibility may prevent sand from being liberated rapidly enough, or in sufficient quantity, to maintain a subaerial barrier. Analyses indicate that factors affecting sediment availability such as low substrate sand proportions and high sediment loss rates cause a barrier to migrate landward along a trajectory having a lower slope than average barrier island slope, thereby defining an “effective” barrier island slope. Other factors being equal, such barriers will tend to be smaller and associated with a more deeply incised shoreface, thereby requiring less migration per sea level rise increment to liberate sufficient sand to maintain subaerial exposure than larger, less incised barriers. As a result, the evolution of larger/less incised barriers is more likely to be limited by shoreface erosion rates or substrate erodibility making them more prone to disintegration related to increasing sea level rise rates than smaller/more incised barriers. Thus, the small/deeply incised North Carolina barriers are likely to persist in the near term (although their long-term fate is less certain because of the low substrate slopes that will soon be encountered). In aggregate, results point to the importance of system history (e.g., previous slopes, sediment budgets, etc.) in determining migration trajectories and therefore how a barrier island will respond to sea level rise. Although simple analytical calculations may predict barrier response in simplified coastal environments (e.g., constant slope, constant sea level rise rate, etc.), our model experiments demonstrate that morphological-behavior modeling is necessary to provide critical insights regarding changes that may occur in environments having complex geometries, especially when multiple parameters change simultaneously.
North Carolina↗
Estimating the timing and location of shallow rainfall-induced landslides using a model for transient, unsaturated infiltration
Shallow rainfall-induced landslides commonly occur under conditions of transient infiltration into initially unsaturated soils. In an effort to predict the timing and location of such landslides, we developed a model of the infiltration process using a two-layer system that consists of an unsaturated zone above a saturated zone and implemented this model in a geographic information system (GIS) framework. The model links analytical solutions for transient, unsaturated, vertical infiltration above the water table to pressure-diffusion solutions for pressure changes below the water table. The solutions are coupled through a transient water table that rises as water accumulates at the base of the unsaturated zone. This scheme, though limited to simplified soil-water characteristics and moist initial conditions, greatly improves computational efficiency over numerical models in spatially distributed modeling applications. Pore pressures computed by these coupled models are subsequently used in one-dimensional slope-stability computations to estimate the timing and locations of slope failures. Applied over a digital landscape near Seattle, Washington, for an hourly rainfall history known to trigger shallow landslides, the model computes a factor of safety for each grid cell at any time during a rainstorm. The unsaturated layer attenuates and delays the rainfall-induced pore-pressure response of the model at depth, consistent with observations at an instrumented hillside near Edmonds, Washington. This attenuation results in realistic estimates of timing for the onset of slope instability (7 h earlier than observed landslides, on average). By considering the spatial distribution of physical properties, the model predicts the primary source areas of landslides.
Journal of Geophysical Research F: Earth Surface↗
Flow structures and sandbar dynamics in a canyon river during a controlled flood, Colorado River, Arizona
In canyon rivers, debris fan constrictions create rapids and downstream pools characterized by secondary flow structures that are closely linked to channel morphology. In this paper we describe detailed measurements of the three‐dimensional flow structure and sandbar dynamics of two pools along the Colorado River in the Grand Canyon during a controlled flood release from Glen Canyon Dam. Results indicate that the pools are characterized by large lateral recirculation zones (eddies) resulting from flow separation downstream from the channel constrictions, as well as helical flow structures in the main channel and eddy. The lateral recirculation zones are low‐velocity areas conducive to fine sediment deposition, particularly in the vicinity of the separation and reattachment points and are thus the dominant flow structures controlling sandbar dynamics. The helical flow structures also affect morphology but appear secondary in importance to the lateral eddies. During the controlled flood, sandbars in the separation and reattachment zones at both sites tended to build gradually during the rising limb and peak flow. Deposition in shallow water on the sandbars was accompanied by erosion in deeper water along the sandbar slope at the interface with the main channel. Erosion occurred via rapid mass failures as well as by gradual boundary shear stress driven processes. The flow structures and morphologic links at our study sites are similar to those identified in other river environments, in particular sharply curved meanders and channel confluences where the coexistence of lateral recirculation and helical flows has been documented.
Journal of Geophysical Research F: Earth Surface↗