Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Earth Observer”

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.

At least 703 records · Page 39Linked to original sources

Polar paleoenvironmental perspectives on modern climate change

In today’s rapidly changing climate, society needs a better understanding of climate impacts on sea level, ice sheets and glaciers, sea ice, ocean circulation, ecosystems, biodiversity, and other aspects of planet Earth. Paleoenvironmental records provide a unique and invaluable source of insight into these complex issues, and place recent observations into a broader historical context. This essay discusses why paleoclimate reconstructions from polar regions provide critical information to help anticipate possible future climate impacts. By highlighting some key research examples, this essay explains the value of expanding proxy-based research in Arctic/Antarctic regions, and makes a case for paying greater attention to the lessons already distilled from it.

PLOS Climate↗

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗

Rare-earth element geochemistry of the island-arc volcanic rocks of Rabaul and Talasea, New Britain

The island-arc volcanic rocks of Rabaul and Talasea, New Britain, range in composition from basalt through rhyolite. Rare-earth elements have been determined by mass-spectrometric isotope dilution in 16 samples. Chondrite-normalized rare-earth element patterns are distinct for each volcanic center, but all are relatively flat (Ce/Yb E.F. = 1.1 to 3.0). Within each center, rare-earth element concentrations increase from basalt to dacite, and Eu anomalies become progressively more negative from basalt to dacite (Eu/Eu* = 1.0 to 0.8). Lavas inferred to contain cumulate phenocrysts show positive Eu anomalies (Eu/Eu* of up to 1.2). The observed variations are consistent with fractional crystallization of basalt or basaltic andesite by removal or accumulation of the observed phenocrysts, including olivine, plagioclase, pyroxene, and opaque minerals. Quantitative trace-element models for Talasea lavas indicate precipitation of 50 wt percent of phenocrysts from basalt to produce andesite and an additional 22 wt percent to produce dacite. A total-equilibrium model provides the best approximation to observed concentrations, implying that the process took place at a depth sufficient to allow slow cooling of the magma as precipitation proceeded. The “least fractionated” basalts and basaltic andesites are similar in rare-earth element pattern to those of corresponding type suites in other oceanic-island arcs, implying a widespread uniformity of source and process in the production of each magma type.

Geological Society of America Bulletin↗

Constraints from GPS measurements on plate coupling within the Makran subduction zone and tsunami scenarios in the western Indian Ocean

Plate-coupling estimates and previous seismicity indicate that portions of the Makran megathrust of southern Pakistan and Iran are partially coupled and have the potential to produce future magnitude 7+ earthquakes. However, the GPS observations needed to constrain coupling models are sparse and lead to an incomplete understanding of regional earthquake and tsunami hazard. In this study, we assess GPS velocities for plate coupling of the Makran subduction zone with specific attention to model resolution and the accretionary prism rheology. We use finite element model-derived Green's functions to invert for the interseismic slip deficit under both elastic and viscoelastic Earth assumptions. We use the model resolution matrix to characterize plate-coupling scenarios that are consistent with the limited spatial resolution afforded by GPS observations. We then forward model the corresponding tsunami responses at major coastal cities within the western Indian Ocean basin. Our plate-coupling results show potential segmentation of the megathrust with varying coupling from west to east, but do not rule out a scenario where the entire length of the megathrust could rupture in a single earthquake. The full subduction zone rupture scenarios suggest that the Makran may be able to produce earthquakes up to M w 9.2. The corresponding tsunami model from the largest earthquake event ( M w 9.2) estimates maximum wave heights reaching 2–5 m at major port cities in the northern Arabian Sea region. Cities on the west coast of India are less affected (1–2 m). Coastlines bounding eastern Africa, and the Strait of Hormuz, are the least affected (<1 m).

Geophysical Journal International↗

Variability in geyser eruptive timing and its causes: Yellowstone National Park

Field data from Upper Geyser Basin, Yellowstone, indicate that geyser frequency is less sensitive to elastic deformation than might be surmised from a review of the literature. Earth-tide influences are not identifiable in any of the geysers we monitored. Though atmospheric-pressure influences are observed, only long-period variations on the order of 5 mBars or greater seem to influence geyser frequency. Long-distance interconnections between geysers are common and add to the difficulty of identifying strain influences. Additional variations in geyser periodicity may be governed by the internal dynamics of the geysers rather than external influences.

Geophysical Research Letters↗

Hotspot of accelerated sea-level rise on the Atlantic coast of North America

Climate warming does not force sea-level rise (SLR) at the same rate everywhere. Rather, there are spatial variations of SLR superimposed on a global average rise. These variations are forced by dynamic processes, arising from circulation and variations in temperature and/or salinity, and by static equilibrium processes, arising from mass redistributions changing gravity and the Earth's rotation and shape. These sea-level variations form unique spatial patterns, yet there are very few observations verifying predicted patterns or fingerprints. Here, we present evidence of recently accelerated SLR in a unique 1,000-km-long hotspot on the highly populated North American Atlantic coast north of Cape Hatteras and show that it is consistent with a modelled fingerprint of dynamic SLR. Between 1950&ndash;1979 and 1980&ndash;2009, SLR rate increases in this northeast hotspot were ~ 3&ndash;4 times higher than the global average. Modelled dynamic plus steric SLR by 2100 at New York City ranges with Intergovernmental Panel on Climate Change scenario from 36 to 51 cm (ref. 3); lower emission scenarios project 24&ndash;36 cm (ref. 7). Extrapolations from data herein range from 20 to 29 cm. SLR superimposed on storm surge, wave run-up and set-up will increase the vulnerability of coastal cities to flooding, and beaches and wetlands to deterioration.

Atlantic Coast, Cape Cod, Cape Hatteras, North Ame↗

Rapid cooling rates at an active mid-ocean ridge from zircon thermochronology

Oceanic spreading ridges are Earth's most productive crust generating environment, but mechanisms and rates of crustal accretion and heat loss are debated. Existing observations on cooling rates are ambiguous regarding the prevalence of conductive vs. convective cooling of lower oceanic crust. Here, we report the discovery and dating of zircon in mid-ocean ridge dacite lavas that constrain magmatic differentiation and cooling rates at an active spreading center. Dacitic lavas erupted on the southern Cleft segment of the Juan de Fuca ridge, an intermediate-rate spreading center, near the intersection with the Blanco transform fault. Their U–Th zircon crystallization ages (29.3 − 4.6 + 4.8 ka; 1σ standard error s.e.) overlap with the (U–Th)/He zircon eruption age (32.7 ± 1.6 ka) within uncertainty. Based on similar 238 U− 230 Th disequilibria between southern Cleft dacite glass separates and young mid-ocean ridge basalt (MORB) erupted nearby, differentiation must have occurred rapidly, within ~ 10–20 ka at most. Ti-in-zircon thermometry indicates crystallization at 850–900 °C and pressures > 70–150 MPa are calculated from H 2 O solubility models. These time-temperature constraints translate into a magma cooling rate of ~ 2 × 10 − 2 °C/a. This rate is at least one order-of-magnitude faster than those calculated for zircon-bearing plutonic rocks from slow spreading ridges. Such short intervals for differentiation and cooling can only be resolved through uranium-series ( 238 U– 230 Th) decay in young lavas, and are best explained by dissipating heat convectively at high crustal permeability.

Juan de Fuca ridge and the Blanco transform fault ↗

Mesoscale raised rim depressions (MRRDs) on Earth: A review of the characteristics, processes, and spatial distributions of analogs for Mars

Fields of mesoscale raised rim depressions (MRRDs) of various origins are found on Earth and Mars. Examples include rootless cones, mud volcanoes, collapsed pingos, rimmed kettle holes, and basaltic ring structures. Correct identification of MRRDs on Mars is valuable because different MRRD types have different geologic and/or climatic implications and are often associated with volcanism and/or water, which may provide locales for biotic or prebiotic activity. In order to facilitate correct identification of fields of MRRDs on Mars and their implications, this work provides a review of common terrestrial MRRD types that occur in fields. In this review, MRRDs by formation mechanism, including hydrovolcanic (phreatomagmatic cones, basaltic ring structures), sedimentological (mud volcanoes), and ice-related (pingos, volatile ice-block forms) mechanisms. For each broad mechanism, we present a comparative synopsis of (i) morphology and observations, (ii) physical formation processes, and (iii) published hypothesized locations on Mars. Because the morphology for MRRDs may be ambiguous, an additional tool is provided for distinguishing fields of MRRDs by origin on Mars, namely, spatial distribution analyses for MRRDs within fields on Earth. We find that MRRDs have both distinguishing and similar characteristics, and observation that applies both to their mesoscale morphology and to their spatial distribution statistics. Thus, this review provides tools for distinguishing between various MRRDs, while highlighting the utility of the multiple working hypotheses approach. ?? 2008 Elsevier Ltd.

Planetary and Space Science↗

Challenges in observational seismology

Earthquake seismology became a quantitative scientific discipline after instruments were developed to record seismic waves in the late 19th century ( Dewey and Byerly, 1969 ; Chapter 1 by Agnew). Earthquake seismology is essentially based on field observations. The great progress made in the past several decades was primarily due to increasingly plentiful and high-quality data that are readily distributed. Our ability to collect, process, and analyze earthquake data has been accelerated by advances in electronics, communications, computers, and software (see Chapter 85 edited by Snoke and Garcia-Fernandez). Instrumental observation of earthquakes has been carried out for a little over 100 years by seismic stations and networks of various sizes, from local to global scales (see Chapter 87 edited by Lahr and van Eck). The observed data have been used, for example, (1) to compute the source parameters of earthquakes, (2) to determine the physical properties of the Earth's interior, (3) to test the theory of plate tectonics , (4) to map active faults, (5) to infer the nature of damaging ground shaking, and (6) to carry out seismic hazard analysis. Construction of a satisfactory theory of the earthquake process has not yet been achieved within the context of physical laws. Good progress, however, has been made in building a physical foundation of the earthquake source process, partly as a result of research directed toward earthquake prediction. This chapter is intended for a general audience. Technical details are not given, but relevant references and chapters in this Handbook are referred to. The first part of this chapter presents a brief overview of the observational aspects of earthquake seismology, concentrating on instrumental observations of seismic waves generated by earthquakes (i.e., seismic monitoring), and readers are referred to Chapter 49 by Musson and Cecic for noninstrumental observations. A few key developments and practices are summarized by taking a general view, since many national and regional developments have been chronicled in national and institutional reports (see Chapter 79 edited by Kisslinger). In the latter part of this chapter, the nature of seismic monitoring and some challenges in observational seismology are discussed from a personal perspective. Comments of a technical or philosophical nature are given in the Notes at the end of the chapter, and they are referenced by superscript numbers in the text.

International Geophysics↗

Reconstructing the dynamics of the highly similar May 2016 and June 2019 Iliamna Volcano, Alaska ice–rock avalanches from seismoacoustic data

Surficial mass wasting events are a hazard worldwide. Seismic and acoustic signals from these often remote processes, combined with other geophysical observations, can provide key information for monitoring and rapid response efforts and enhance our understanding of event dynamics. Here, we present seismoacoustic data and analyses for two very large ice–rock avalanches occurring on Iliamna Volcano, Alaska (USA), on 22 May 2016 and 21 June 2019. Iliamna is a glacier-mantled stratovolcano located in the Cook Inlet, ∼200 km from Anchorage, Alaska. The volcano experiences massive, quasi-annual slope failures due to glacial instabilities and hydrothermal alteration of volcanic rocks near its summit. The May 2016 and June 2019 avalanches were particularly large and generated energetic seismic and infrasound signals which were recorded at numerous stations at ranges from ∼9 to over 600 km. Both avalanches initiated in the same location near the head of Iliamna's east-facing Red Glacier, and their ∼8 km long runout shapes are nearly identical. This repeatability – which is rare for large and rapid mass movements – provides an excellent opportunity for comparison and validation of seismoacoustic source characteristics. For both events, we invert long-period (15–80 s) seismic signals to obtain a force-time representation of the source. We model the avalanche as a sliding block which exerts a spatially static point force on the Earth. We use this force-time function to derive constraints on avalanche acceleration, velocity, and directionality, which are compatible with satellite imagery and observed terrain features. Our inversion results suggest that the avalanches reached speeds exceeding 70 m s−1, consistent with numerical modeling from previous Iliamna studies. We lack sufficient local infrasound data to test an acoustic source model for these processes. However, the acoustic data suggest that infrasound from these avalanches is produced after the mass movement regime transitions from cohesive block-type failure to granular and turbulent flow – little to no infrasound is generated by the initial failure. At Iliamna, synthesis of advanced numerical flow models and more detailed ground observations combined with increased geophysical station coverage could yield significant gains in our understanding of these events.

Alaska↗

A plate flexure approximation to postseismic and interseismic deformation

The rather large postseismic deformation that is associated with two‐dimensional dip‐slip faulting in the lithosphere is related to the bending of a free plate generated by dip‐slip faulting. In the absence of gravity, asthenosphere relaxation eventually permits the faulted lithosphere to assume the dihedral configuration of a faulted free plate. For thrust faulting, the faulted area is depressed into the asthenosphere, and the flanks of the plate slope uniformly upward. In the presence of gravity, buoyancy forces act upon the plate, and the ultimate ( t →∞) postseismic uplift is approximated ( y /α > 0.2) by w 1 exp (− y /α) (cos y /α ‐ sin y /α), where w 1 is an explicit function of the fault parameters and lithosphere properties, y is the horizontal distance from the downdip end of the fault, and α is the flexural parameter for the lithosphere plate. The relaxed‐asthenosphere response is the sum of the coseismic deformation (unrelaxed‐astheriosphere response) and this ultimate postseismic deformation. The annual, steady deformation associated with strain accumulation at a subduction zone is simply the relaxed‐asthenosphere response to virtual, normal slip on the main thrust zone of an amount equal to the annual plate convergence. From these relations we have estimated the deformation (strain and uplift) expected along the Pacific coast of northern Honshu. The estimated deformation is about twice the observed deformation, but the predicted geographical distribution of the deformation is similar to the observed deformation.

Journal of Geophysical Research B: Solid Earth↗

Apophis specific action team report

This report about Asteroid (99942) Apophis's Earth close approach on April 13, 2029 was generated by a Specific Action Team (SAT) formed by the Small Body Assessment Group (SBAG) at the request of NASAs Planetary Science Division (PSD). The SAT assessed the current predictions for the effects that may occur due to the close encounter, evaluated observing capabilities, and identified possible investigations then sorted them into priority categories. In addition, the SAT evaluated whether or not a spacecraft sent to Apophis could increase the risk of a future Earth impact.

Report↗

Plume-driven plumbing and crustal formation in Iceland

Through combination of surface wave and body wave constraints we derive a three-dimensional (3-D) crustal S velocity model and Moho map for Iceland. It reveals a vast plumbing system feeding mantle plume melt into upper crustal magma chambers where crustal formation takes place. The method is based on the partitioned waveform inversion to which we add additional observations. Love waves from six local events recorded on the HOTSPOT- SIL networks are fitted, S n travel times from the same events measured, previous observations of crustal thickness are added, and all three sets of constraints simultaneously inverted for our 3-D model. In the upper crust (0–15 km) an elongated low-velocity region extends along the length of the Northern, Eastern and Western Neovolcanic Zones. The lowest velocities (−7%) are found at 5–10 km below the two most active volcanic complexes: Hekla and Bárdarbunga-Grímsvötn. In the lower crust (>15 km) the low-velocity region can be represented as a vertical cylinder beneath central Iceland. The low-velocity structure is interpreted as the thermal halo of pipe work which connects the region of melt generation in the uppermost mantle beneath central Iceland to active volcanoes along the neovolcanic zones. Crustal thickness in Iceland varies from 15–20 km beneath the Reykjanes Peninsula, Krafla and the extinct Snæfellsnes rift zone, to 46 km beneath central Iceland. The average crustal thickness is 29 km. The variations in thickness can be explained in terms of the temporal variation in plume productivity over the last ∼20 Myr, the Snæfellsnes rift zone being active during a minimum in plume productivity. Variations in crustal thickness do not depart significantly from an isostatically predicted crustal thickness. The best fit linear isostatic relation implies an average density jump of 4% across the Moho. Rare earth element inversions of basalt compositions on Iceland suggest a melt thickness (i.e., crustal thickness) of 15–20 km, given passive upwelling. The observed crustal thickness of up to 46 km implies active fluxing of source material through the melt zone by the mantle plume at up to 3 times the passive rate.

Journal of Geophysical Research B: Solid Earth↗

Meteorites on Mars observed with Mars Exploration Rovers

Reduced weathering rates due to the lack of liquid water and significantly greater typical surface ages should result in a higher density of meteorites on the surface of Mars compared to Earth. Several meteorites were identified among the rocks investigated during Opportunity's traverse across the sandy Meridiani plains. Heat Shield Rock is a IAB iron meteorite and has been officially recognized as 'Meridiani Planum.' Barberton is olivine-rich and contains metallic Fe in the form of kamacite, suggesting a meteoritic origin. It is chemically most consistent with a mesosiderite silicate clast. Santa Catarina is a brecciated rock with a chemical and mineralogical composition similar to Barberton. Barberton, Santa Catarina, and cobbles adjacent to Santa Catarina may be part of a strewn field. Spirit observed two probable iron meteorites from its Winter Haven location in the Columbia Hills in Gusev Crater. Chondrites have not been identified to date, which may be a result of their lower strengths and probability to survive impact at current atmospheric pressures. Impact craters directly associated with Heat Shield Rock, Barberton, or Santa Catarina have not been observed, but such craters could have been erased by eolian-driven erosion.

Journal of Geophysical Research E: Planets↗

Ultraviolet and visible remote sensing of volcanic gases

As magma rises in volcanic systems, volatile species exsolve from the silicate melt and are emitted as gases into the atmosphere. Measuring the magnitude and composition of gas emissions from volcanoes provides insights into processes occurring deep within the Earth and helps constrain the impact of volcanic degassing on atmospheric chemistry. Optical remote sensing techniques allow volcanic gas emissions to be characterized without the need to access hazardous areas near active volcanic vents. This paper reviews the state of the art in ultraviolet and visible volcanic gas remote sensing from the ground, air, and space. Special attention is given to discussing the physics of atmospheric radiative transfer on which these techniques are based. The functionality and limitations of different remote sensing instruments are examined, making clear that the ideal choice of instrumentation will depend on the volcanic system to which it is applied and the sought measurement parameters. Common algorithms for determining trace gas column densities, gas burdens, and volcanic emission rates from measurements of spectral radiance are outlined and compared, showing how some algorithms attempt to model the physics of the measurement while others maximize sensitivity. Several examples demonstrate how remote sensing measurements continue to advance our understanding of volcanic systems and their impact on the atmosphere. Finally, a few promising directions of inquiry are suggested that could lead to improvements in remote sensing instrumentation and analysis techniques. By combining spectroscopic and imaging techniques, improving our understanding of atmospheric radiative transfer, expanding the suite of target gases, and increasing the coverage and frequency of observations, we stand to significantly improve our ability to detect and quantify volcanic gas emissions and gain new insights into important Earth-system processes.

Journal of Volcanology and Geothermal Research↗

Coupling ice-sheet and climate models for simulation of former ice sheets

This chapter explores the development of coupled climate and ice-sheet models over the past two decades, discusses the current technical and physical capabilities of models, and identifies future work for developing a better understanding of ice-climate events that have punctuated Earth history. The chapter also illustrates the complex behavior of the climate system and the modeling challenges posed by the observations. Climate and ice-sheet models continue to improve, both in terms of model physics and technical capabilities. Regional climate model simulations should be carried out for each temporal snapshot or matrix element for improved mass–balance calculation over the ice sheets. Fully coupled atmosphere-ocean-cryosphere-land surface models are required for addressing a number of paleoclimatic puzzles, particularly with respect to millennial climate variability. The concerted effort to better understand West Antarctic Ice Sheet dynamics, and the development of subglacial hydrological models should lead to improvements in the next generation of ice-sheet models and help to address millennial-scale variability in ice-sheet/climate models.

Developments in Quaternary Sciences↗

Rapid shallow megathrust afterslip from the 2021 M8.2 Chignik, Alaska earthquake revealed by seafloor geodesy

The shallower portions of subduction zone megathrust faults host Earth’s most hazardous tsunamigenic earthquakes, yet understanding how and when they slip remains elusive because of challenges making seafloor observations. We performed Global Navigation Satellite System Acoustic seafloor geodetic surveys before and ~2.5 months after the 29 July 2021 M w (moment magnitude) 8.2 Chignik, Alaska, earthquake and determine ~1.4 meters cumulative co- and post-seismic horizontal displacement ~60 kilometers from the megathrust front. Only for the 2011 M w 9 Tohoku event have closer subduction zone earthquake displacements been observed. We estimate ~2 to 3 meters of megathrust afterslip shallower than 20 kilometers, a portion of the megathrust on which both inter- and co-seismic slip likely had occurred previously. Our analysis demonstrates that by 2.5 months, shallower and deeper moment had effectively equilibrated on the megathrust, suggesting that its tsunamigenic potential remains no more elevated than before the earthquake.

Alaska↗