Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “The 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 1,513 records · Page 84Linked to original sources

Prolonged magmatic activity on Mars inferred from the detection of felsic rocks

Rocks dominated by the silicate minerals quartz and feldspar are abundant in Earth’s upper continental crust. Yet felsic rocks have not been widely identified on Mars, a planet that seems to lack plate tectonics and the associated magmatic processes that can produce evolved siliceous melts on Earth. If Mars once had a feldspar-rich crust that crystallized from an early magma ocean such as that on the Moon, erosion, sedimentation and volcanism have erased any clear surface evidence for widespread felsic materials. Here we report near-infrared spectral evidence from the Compact Reconnaissance Imaging Spectrometer for Mars onboard the Mars Reconnaissance Orbiter for felsic rocks in three geographically disparate locations on Mars. Spectral characteristics resemble those of feldspar-rich lunar anorthosites, but are accompanied by secondary alteration products (clay minerals). Thermodynamic phase equilibrium calculations demonstrate that fractional crystallization of magma compositionally similar to volcanic flows near one of the detection sites can yield residual melts with compositions consistent with our observations. In addition to an origin by significant magma evolution, the presence of felsic materials could also be explained by feldspar enrichment by fluvial weathering processes. Our finding of felsic materials in several locations on Mars suggests that similar observations by the Curiosity rover in Gale crater may be more widely applicable across the planet.

Nature Geoscience↗

Spatial distribution of water level impact to back-barrier bays

Water level in semi-enclosed bays, landward of barrier islands, is mainly driven by offshore sea level fluctuations that are modulated by bay geometry and bathymetry, causing spatial variability in the ensuing response (transfer). Local wind setup can have a secondary role that depends on wind speed, fetch, and relative orientation of the wind direction and the bay. Inlet geometry and bathymetry primarily regulate the magnitude of the transfer between open ocean and bay. Tides and short-period offshore oscillations are more damped in the bays than longer-lasting offshore fluctuations, such as storm surge and sea level rise. We compare observed and modeled water levels at stations in a mid-Atlantic bay (Barnegat Bay) with offshore water level proxies. Observed water levels in Barnegat Bay are compared and combined with model results from the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system to evaluate the spatial structure of the water level transfer. Analytical models based on the dimensional characteristics of the bay are used to combine the observed data and the numerical model results in a physically consistent approach. Model water level transfers match observed values at locations inside the Bay in the storm frequency band (transfers ranging from 70-100%) and tidal frequencies (10-55%). The contribution of frequency-dependent local setup caused by wind acting along the bay is also considered. The approach provides transfer estimates for locations inside the Bay where observations were not available resulting in a complete spatial characterization. The approach allows for the study of the Bay response to alternative forcing scenarios (landscape changes, future storms, and rising sea level). Detailed spatial estimates of water level transfer can inform decisions on inlet management and contribute to the assessment of current and future flooding hazard in back-barrier bays and along mainland shorelines.

Natural Hazards and Earth System Sciences↗

Drought early warning and forecasting

Drought risk management involves three pillars: drought early warning, drought vulnerability and risk assessment, and drought preparedness, mitigation, and response. This book collects in one place a description of all the key components of the first pillar, and describes strategies for fitting these pieces together. The best modern drought early warning systems incorporate and integrate a broad array of environmental information sources: weather station observations, satellite imagery, land surface and crop model simulations, and weather and climate model forecasts, and analyze this information in context-relevant ways that take into account exposure and vulnerability. Drought Early Warning and Forecasting: Theory and Practice assembles a comprehensive overview of these components, providing examples drawn from the Famine Early Warning Systems Network and the United States Drought Monitor. This book simultaneously addresses the physical, social, and information management aspects of drought early warning, and informs readers about the tools, techniques, and conceptual models required to effectively identify, predict, and communicate potential drought-related disasters. This book is a key text for postgraduate scientists and graduate and advanced undergraduate students in hydrology, geography, earth sciences, meteorology, climatology, and environmental sciences programs. Professionals dealing with disaster management and drought forecasting will also find this book beneficial to their work.

Book↗

Rotational motions for teleseismic surface waves

[1] We report the findings for the first teleseismic six degree‐of‐freedom (6‐DOF) measurements including three components of rotational motions recorded by a sensitive rotation‐rate sensor (model R‐1, made by eentec) and three components of translational motions recorded by a traditional seismometer (STS‐2) at the NACB station in Taiwan. The consistent observations in waveforms of rotational motions and translational motions in sections of Rayleigh and Love waves are presented in reference to the analytical solution for these waves in a half space of Poisson solid. We show that additional information (e.g., Rayleigh wave phase velocity, shear wave velocity of the surface layer) might be exploited from six degree‐of‐freedom recordings of teleseismic events at only one station. We also find significant errors in the translational records of these teleseismic surface waves due to the sensitivity of inertial translation sensors (seismometers) to rotational motions. The result suggests that the effects of such errors need to be counted in surface wave inversions commonly used to derive earthquake source parameters and Earth structure.

Geophysical Research Letters↗

Viscosity and finite strength of the mantle as determined from water and ice loads

Some recent examples of transient Earth loads (Lake Bonneville, Utah; Glacier Bay, Alaska; northeast Greenland) indicate that both the viscosity and finite strength of the mantle are lower than is commonly presumed. A time constant (1/ e ) of 4000 years is estimated for Lake Bonneville, and of 1000 years for northeast Greenland. A strain rate of 10 −14 is typical. These figures imply viscosities in a homogeneous half space ranging from 10 20 to 10 21 poises. An upper limit of finite strength is set by Lake Bonneville at a few times 10 6 dyn/cm 2 . If mountain ranges like the Sierra Nevada or Himalaya are regarded as dynamically supported rather than static systems, this low value is not incompatible with other geologic observations.

Alaska, Utah↗

Knowledge, transparency, and refutability in groundwater models, an example from the Death Valley regional groundwater flow system

This work demonstrates how available knowledge can be used to build more transparent and refutable computer models of groundwater systems. The Death Valley regional groundwater flow system, which surrounds a proposed site for a high level nuclear waste repository of the United States of America, and the Nevada National Security Site (NNSS), where nuclear weapons were tested, is used to explore model adequacy, identify parameters important to (and informed by) observations, and identify existing old and potential new observations important to predictions. Model development is pursued using a set of fundamental questions addressed with carefully designed metrics. Critical methods include using a hydrogeologic model, managing model nonlinearity by designing models that are robust while maintaining realism, using error-based weighting to combine disparate types of data, and identifying important and unimportant parameters and observations and optimizing parameter values with computationally frugal schemes. The frugal schemes employed in this study require relatively few (10–1000 s), parallelizable model runs. This is beneficial because models able to approximate the complex site geology defensibly tend to have high computational cost. The issue of model defensibility is particularly important given the contentious political issues involved.

California, Nevada↗

Four major Holocene earthquakes on the Reelfoot fault recorded by sackungen in the New Madrid seismic zone, USA

Three sequences of well-documented, major ~M7+ earthquakes (1811-1812 CE, ~1450 CE, and ~900 CE) in the New Madrid seismic zone, USA, contribute significantly to seismic hazard in the region. However, it is unknown whether this <550 yr recurrence interval has been constant throughout the Holocene given limited geomorphic evidence of prior earthquakes. We extend the record of paleoearthquakes along the Reelfoot fault via investigation of ridge-top gravitational failure features, interpreted as sackungen. The sackungen occur in bluffs along the eastern margin of the Mississippi River floodplain and are concentrated near (<15 km) the southwest-dipping Reelfoot reverse fault. A paleoseismic trench excavated across sackungen at the Paw Paw site exposed four packages of colluvial sediment that postdate 30-11 ka Peoria loess. We interpret the colluvial packages to have been deposited following episodic failure of the sackungen as a result of strong ground motions from the following sequence of earthquakes: event 4, 1640 ± 1730 BCE; event 3, 340 ± 670 CE; event 2, 1430 ± 380 CE; and event 1, 1810 ± 50 CE (2-sigma). Event timing corresponds to previously documented earthquakes and represents the longest archive of paleoearthquakes on the Reelfoot fault. If the trenched sackungen record all major Reelfoot fault earthquakes, our observations in combination with prior investigations indicate a period of quiescence from at least 11 – 4.7 ka, followed by four major seismic events culminating in the 1811-1812 CE sequence. This clustered earthquake recurrence helps place bounds on seismic-hazard and geodynamic models in the New Madrid seismic zone.

Arkansas, Missouri↗

Aqueously altered igneous rocks sampled on the floor of Jezero crater, Mars

The Perseverance rover landed in Jezero crater, Mars, to investigate ancient lake and river deposits. We report observations of the crater floor, below the crater’s sedimentary delta, finding the floor consists of igneous rocks altered by water. The lowest exposed unit, informally named Séítah, is a coarsely crystalline olivine-rich rock, which accumulated at the base of a magma body. Fe-Mg carbonates along grain boundaries indicate reactions with CO 2 -rich water, under water-poor conditions. Overlying Séítah is a unit informally named Máaz, which we interpret as lava flows or the chemical complement to Séítah in a layered igneous body. Voids in these rocks contain sulfates and perchlorates, likely introduced by later near-surface brine evaporation. Core samples of these rocks were stored aboard Perseverance for potential return to Earth.

Science↗

Aftershocks of the 2010 M w 7.2 El Mayor-Cucapah earthquake revealcomplex faulting in the Yuha Desert, California

We detect and precisely locate over 9500 aftershocks that occurred in the Yuha Desert region during a 2 month period following the 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake. Events are relocated using a series of absolute and relative relocation procedures that include Hypoinverse, Velest, and hypoDD. Location errors are reduced to ~40 m horizontally and ~120 m vertically.Aftershock locations reveal a complex pattern of faulting with en echelon fault segments trending toward the northwest, approximately parallel to the North American-Pacific plate boundary and en echelon, conjugate features trending to the northeast. The relocated seismicity is highly correlated with published surface mapping of faults that experienced triggered surface slip in response to the EMC main shock. Aftershocks occurred between 2 km and 11 km depths, consistent with previous studies of seismogenic thickness in the region. Three-dimensional analysis reveals individual and intersecting fault planes that are limited in their along-strike length. These fault planes remain distinct structures at depth, indicative of conjugate faulting, and do not appear to coalesce onto a throughgoing fault segment. We observe a complex spatiotemporal migration of aftershocks, with seismicity that jumps between individual fault segments that are active for only a few days to weeks. Aftershock rates are roughly consistent with the expected earthquake production rates of Dieterich (1994). The conjugate pattern of faulting and nonuniform aftershock migration patterns suggest that strain in the Yuha Desert is being accommodated in a complex manner.

California↗

Frequency-dependent effects of rupture for the 2004 Parkfield mainshock, results from UPSAR

The frequency-dependent effects of rupture propagation of the Parkfield, California earthquake (Sept. 28, 2004, M6) to the northwest along the San Andreas fault can be seen in acceleration records at UPSAR (USGS Parkfield Seismic Array) in at least two ways. First, we can see the effects of directivity in the acceleration traces at UPSAR, which is about 11.5 km from the epicenter. Directivity or the seismic equivalent of a Doppler shift has been documented in many cases by comparing short duration, high-amplitude pulses (P or S) in the forward direction with longer duration body waves in the backward direction. In this case we detect a change from a relatively large amplitude, coherent, high-frequency signal at the start of rupture to a low-amplitude, low-coherent, low-frequency signal at about the time the rupture front transfers from the forward azimuth to the back azimuth at about 34-36 s (time is UTC and are the seconds after day 272 and 17 hours and 15 minutes. S arrival is just after 30s) for rays leaving the fault and propagating to UPSAR. The frequency change is obvious in the band about 5 to 30 Hz, which is significantly above the corner frequency of the earthquake (about 0.11Hz). From kinematic source models, the duration of faulting is about 9.2 s and the change in frequency is during faulting as the rupture extends to the northwest. Understanding the systematic change in frequency and amplitude of seismic waves in relation to the propagation of the rupture front is important for predicting strong ground motion. Second, we can filter the acceleration records from the array to determine if the low frequency energy emerges from the same part of the fault as the high frequency signal (e.g. has the same back azimuth and apparent velocity at UPSAR) an important clue to the dynamics of rupture. Analysis of sources of strong motion (characterized by relatively high frequencies) compared to kinematic slip models (relatively low frequency) for the March 11, 2011 Tohoku earthquake as well as Maule (Feb. 27, 2010) and Chi-Chi (Sept. 20, 1999) earthquakes show that high- and low-frequency sources do not have the same locations on the fault. In this paper we filter the accelerograms from UPSAR for the 2004 mainshock in various passbands and then re-compute the cross correlations to determine the vector slowness of the incoming waves. At Parkfield, it appears that for seismic waves with frequencies above 1 Hz there is no discernible frequency-dependent difference in source position (up to 8 Hz) based on estimates of back azimuth and apparent velocity. However at lower frequencies, sources appear to be from shallower depths and trail the high frequencies as the rupture proceeds down the fault. This result is greater than one standard deviation of an estimate of error, based on a new method of estimating error that is a measure of how broad the peak in correlation is and an estimate of the variance of the correlation values. These observations can be understood in terms of a rupture front that is more energetic and coherent near the front of rupture (radiating higher frequencies) and less coherent and less energetic (radiating in a lower frequency band) behind the initial rupture front. This result is a qualitative assessment of changes in azimuth and apparent velocity with frequency and time and does not include corrections to find the source location on the fault.

California↗

In situ measurements of post-fire debris flows in southern California: Comparisons of the timing and magnitude of 24 debris-flow events with rainfall and soil moisture conditions

Debris flows often occur in burned steeplands of southern California, sometimes causing property damage and loss of life. In an effort to better understand the hydrologic controls on post-fire debris-flow initiation, timing and magnitude, we measured the flow stage, rainfall, channel bed pore fluid pressure and hillslope soil-moisture accompanying 24 debris flows recorded in five different watersheds burned in the 2009 Station and Jesusita Fires (San Gabriel and Santa Ynez Mountains). The measurements show substantial differences in debris-flow dynamics between sites and between sequential events at the same site. Despite these differences, the timing and magnitude of all events were consistently associated with local peaks in short duration (< = 30 min) rainfall intensity. Overall, debris-flow stage was best cross-correlated with time series of 5-min rainfall intensity, and lagged the rainfall by an average of just 5 min. An index of debris-flow volume was also best correlated with short-duration rainfall intensity, but found to be poorly correlated with storm cumulative rainfall and hillslope soil water content. Post-event observations of erosion and slope stability modeling suggest that the debris flows initiated primarily by processes related to surface water runoff, rather than shallow landslides. By identifying the storm characteristics most closely associated with post-fire debris flows, these measurements provide valuable guidance for warning operations and important constraints for developing and testing models of post-fire debris flows. copyright. 2011 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface↗

Observations of an extreme atmospheric river storm with a diverse sensor network

Observational networks enhance real‐time situational awareness for emergency and water resource management during extreme weather events. We present examples of how a diverse, multitiered observational network in California provided insights into hydrometeorological processes and impacts during a 3‐day atmospheric river storm centered on 14 February 2019. This network, which has been developed over the past two decades, aims to improve understanding and mitigation of effects from extreme storms influencing water resources and natural hazards. We combine atmospheric reanalysis output and additional observations to show how the network allows: (1) the validation of record cool season precipitable water observations over southern California; (2) the identification of phenomena that produce natural hazards and present difficulties for short‐term weather forecast models, such as extreme precipitation amounts and snow level variability; (3) the use of soil moisture data to improve hydrologic model forecast skill in northern California's Russian River basin; and (4) the combination of meteorological data with seismic observations to identify when a large avalanche occurred on Mount Shasta. This case study highlights the value of investments in diverse observational assets and the importance of continued support and synthesis of these networks to characterize climatological context and advance understanding of processes modulating extreme weather.

California↗

Landscape-scale modeling to forecast fluvial-aeolian sediment connectivity in river valleys

Sedimentary landforms on Earth and other planetary bodies are built through scour, transport, and deposition of sediment. Sediment connectivity refers to the hypothesis that pathways of sediment transport do not occur in isolation, but rather are mechanistically linked. In dryland river systems, one such example of sediment connectivity is the transport of fluvially deposited sediment by wind. However, predictive tools that can forecast fluvial-aeolian sediment connectivity at meaningful scales are rare. Here we develop a suite of models for quantifying the availability of river-sourced sediment for aeolian transport as a function of river flow, wind regime, and land cover across 168 km of the Colorado River in Grand Canyon, USA. We compare and validate these models using topographic changes observed over 10 years in a coupled river sandbar-aeolian dunefield setting. The models provide a path forward for directly linking fluvial hydrology with the management and understanding of aeolian landscapes.

Arizona↗

Multi-scale intense geoelectric and geomagnetic field perturbations observed after an interplanetary magnetic field turning

Intense geoelectric fields during geomagnetic storms generate geomagnetically induced currents in power grids and other infrastructure, necessitating an understanding of their causes, for example, through coordinated space and ground observations. This study investigates localized intense geoelectric ( E ) and geomagnetic ( B ) field perturbations following an Interplanetary Magnetic Field (IMF) turning during a geomagnetic storm on 25 October 2011. Observations from EarthScope magnetotelluric sites in the upper Midwest United States revealed shorter period ( ∼ 1 min) ultra-low-frequency (ULF) waves superimposed on longer period ( ∼ 10 min) perturbations in both E and B fields. These sites, located at ∼ 19 hr magnetic local time and 56 −57 ° magnetic latitude, recorded large amplitude E and B perturbations. Ground-based all-sky imagers showed auroral brightening with sunward and poleward propagation, while upstream spacecraft linked the perturbations to an IMF turning and solar wind dynamic pressure impulse. The longer-period E and B field perturbations likely stem from localized ionospheric currents tied to substorm auroral activity post-IMF turning. The combination of ionospheric currents, ULF waves, and the Earth's varying conductivity produces intense geoelectric fields of ≥ 2 V/km in the upper Midwest. A comparison using input data and software compatible with the NOAA/USGS geoelectric field nowcast model revealed its limitations in capturing such events due to the temporal and spatial resolution of the underlying data. Using 1-s geomagnetic field data can improve geoelectric field models by capturing short-period and large spatial scale waves, although localized magnetic perturbations remain underestimated due to insufficient ground magnetometer density.

Space Weather↗

A conserved interdomain microbial network underpins cadaver decomposition despite environmental variables

Microbial breakdown of organic matter is one of the most important processes on Earth, yet the controls of decomposition are poorly understood. Here we track 36 terrestrial human cadavers in three locations and show that a phylogenetically distinct, interdomain microbial network assembles during decomposition despite selection effects of location, climate and season. We generated a metagenome-assembled genome library from cadaver-associated soils and integrated it with metabolomics data to identify links between taxonomy and function. This universal network of microbial decomposers is characterized by cross-feeding to metabolize labile decomposition products. The key bacterial and fungal decomposers are rare across non-decomposition environments and appear unique to the breakdown of terrestrial decaying flesh, including humans, swine, mice and cattle, with insects as likely important vectors for dispersal. The observed lockstep of microbial interactions further underlies a robust microbial forensic tool with the potential to aid predictions of the time since death.

Nature Microbiology↗

Using mineral geochemistry to decipher slab, mantle, and crustal input in the generation of high-Mg andesites and basaltic andesites from the northern Cascade Arc

To better understand the role of slab melt in the petrogenesis of North Cascades magmas, this study focuses on petrogenesis of high-Mg lavas from the two northernmost active volcanoes in Washington. High-Mg andesites (HMA) and basaltic andesites (HMBA) in the Cascade Arc have high Mg# [molar Mg/(Mg+Fe 2+ )] relative to their SiO 2 contents, elevated Nd/Yb, and are Ni- and Cr-enriched. The rock units examined here include the Tarn Plateau HMBA (51.8–54.0 wt% SiO 2 , Mg# 68–70) and Glacier Creek HMA (58.3–58.7 wt% SiO 2 , Mg# 63–64) from the Mount Baker Volcanic Field, and the Lightning Creek HMBA (54.8–54.6 SiO 2 , Mg# 69–73) from Glacier Peak. This study combines major and trace element compositions of minerals and whole rocks to test several petrogenetic hypotheses and to determine which, if any, are applicable to North Cascades HMA and HMBA. In the Tarn Plateau HMBA, rare earth element (REE) equilibrium liquids calculated from clinopyroxene compositions have high Nd/Yb that positively correlates with Mg#. This correlation suggests an origin similar to that proposed for Aleutian adakites, where intermediate, high Nd/Yb slab-derived melts interact with the overlying mantle to become Mg-rich, and subsequently mix with low Nd/Yb, mantle-derived mafic magmas with lower Mg#. In the Glacier Creek HMA, elevated whole-rock MgO and SiO 2 contents resulted from accumulation of xenocrystic olivine and differentiation processes, respectively, but the cause of high Nd/Yb is less clear. However, high whole-rock Sr/P (fluid mobile/fluid immobile) values indicate a mantle source that was fluxed by an enriched, hydrous slab component, likely producing the observed high Nd/Yb REE signature. The Lightning Creek HMBA is a hybridized rock unit with at least three identifiable magmatic components, but only one of which has HMA characteristics. Cr and Mg contents in Cr-spinel and olivine pairs in this HMA component suggest that its source is a strongly depleted mantle, and high whole-rock Sr/P values indicate mantle melting that was induced through hydration, likely adding the component responsible for the observed high Nd/Yb REE pattern. The elevated SiO 2 contents (54.6 wt%) of the HMA component resulted from differentiation or high degrees of partial melting of ultramafic material through the addition of H 2 O. Therefore the Lightning Creek HMBA is interpreted to have originated from a refractory mantle source that underwent melting through interaction with an enriched slab component. Our results indicate that in addition to slab-derived fluids, slab-derived melts also have an important role in the production of HMA and HMBA in the north Cascade Arc.

American Mineralogist↗

The seasonal cycle of diabatic heat storage in the Pacific Ocean

This study quantifies uncertainties in closing the seasonal cycle of diabatic heat storage (DHS) over the Pacific Ocean from 20??S to 60??N through the synthesis of World Ocean Circulation Experiment (WOCE) reanalysis products from 1993 to 1999. These products are DHS from Scripps Institution of Oceanography (SIO); near-surface geostrophic and Ekman currents from Earth and Space Research (ESR); and air-sea heat fluxes from Comprehensive Ocean-Atmosphere Data Set (COADS), National Centers for Environmental Prediction (NCEP), and European Center for Mid-Range Weather Forecasts (ECMWF). With these products, we compute residual heat budget components by differencing long-term monthly means from the long-term annual mean. This allows the seasonal cycle of the DHS tendency to be modeled. Everywhere latent heat flux residuals dominate sensible heat flux residuals, shortwave heat flux residuals dominate longwave heat flux residuals, and residual Ekman heat advection dominates residual geostrophic heat advection, with residual dissipation significant only in the Kuroshio-Oyashio current extension. The root-mean-square (RMS) of the differences between observed and model residual DHS tendencies (averaged over 10??latitude-by-20??longitude boxes) is <20 W m-2 in the interior ocean and <100 W m-2 in the Kuroshio-Oyashio current extension. This reveals that the residual DHS tendency is driven everywhere by some mix of residual latent heat flux, shortwave heat flux, and Ekman heat advection. Suppressing bias errors in residual air-sea turbulent heat fluxes and Ekman heat advection through minimization of the RMS differences reduces the latter to <10 W m-2 over the interior ocean and <25 W m -2 in the Kuroshio-Oyashio current extension. This reveals air-sea temperature and specific humidity differences from in situ surface marine weather observations to be a principal source of bias error, overestimated over most of ocean but underestimated near the Intertropical Convergence Zone. ?? 2005 Elsevier Ltd. All rights reserved.

Progress in Oceanography↗

A 1.2 billion pixel human-labeled dataset for data-driven classification of coastal environments

The world’s coastlines are spatially highly variable, coupled-human-natural systems that comprise a nested hierarchy of component landforms, ecosystems, and human interventions, each interacting over a range of space and time scales. Understanding and predicting coastline dynamics necessitates frequent observation from imaging sensors on remote sensing platforms. Machine Learning models that carry out supervised (i.e., human-guided) pixel-based classification, or image segmentation, have transformative applications in spatio-temporal mapping of dynamic environments, including transient coastal landforms, sediments, habitats, waterbodies, and water flows. However, these models require large and well-documented training and testing datasets consisting of labeled imagery. We describe “Coast Train,” a multi-labeler dataset of orthomosaic and satellite images of coastal environments and corresponding labels. These data include imagery that are diverse in space and time, and contain 1.2 billion labeled pixels, representing over 3.6 million hectares. We use a human-in-the-loop tool especially designed for rapid and reproducible Earth surface image segmentation. Our approach permits image labeling by multiple labelers, in turn enabling quantification of pixel-level agreement over individual and collections of images.

Scientific Data↗