Results for “spatial mapping” — Page 5

NASA NTRS2017

Using High Spatial Resolution Satellite Imagery to Map Forest Burn Severity Across Spatial Scales in a Pine Barrens Ecosystem

As a primary disturbance agent, fire significantly influences local processes and services of forest ecosystems. Although a variety of remote sensing based approaches have been developed and applied to Landsat mission imagery to infer burn severity at 30 m spatial resolution, forest burn severity have still been seldom assessed at fine spatial scales (less than or equal to 5 m) from very-high-resolution (VHR) data. We assessed a 432 ha forest fire that occurred in April 2012 on Long Island, New York, within the Pine Barrens region, a unique but imperiled fire-dependent ecosystem in the northeastern United States. The mapping of forest burn severity was explored here at fine spatial scales, for the first time using remotely sensed spectral indices and a set of Multiple Endmember Spectral Mixture Analysis (MESMA) fraction images from bi-temporal - pre- and post-fire event - WorldView-2 (WV-2) imagery at 2 m spatial resolution. We first evaluated our approach using 1 m by 1 m validation points at the sub-crown scale per severity class (i.e. unburned, low, moderate, and high severity) from the post-fire 0.10 m color aerial ortho-photos; then, we validated the burn severity mapping of geo-referenced dominant tree crowns (crown scale) and 15 m by 15 m fixed-area plots (inter-crown scale) with the post-fire 0.10 m aerial ortho-photos and measured crown information of twenty forest inventory plots. Our approach can accurately assess forest burn severity at the sub-crown (overall accuracy is 84% with a Kappa value of 0.77), crown (overall accuracy is 82% with a Kappa value of 0.76), and inter-crown scales (89% of the variation in estimated burn severity ratings (i.e. Geo-Composite Burn Index (CBI)). This work highlights that forest burn severity mapping from VHR data can capture heterogeneous fire patterns at fine spatial scales over the large spatial extents. This is important since most ecological processes associated with fire effects vary at the less than 30 m scale and VHR approaches could significantly advance our ability to characterize fire effects on forest ecosystems.

Meng, Ran↗

NASA NTRS2019

A Long-Term Vision for Space-Based Interferometry

The processes leading to the formation of planets; the extreme physics occurring near the event horizon of black holes; detailed studies of exoplanets through spectral-spatial mapping: new and unique insights into the physical processes involved across nearly the whole gamut of astrophysics await discovery at small angular scales. The fine spatial resolution needed to explore these processes, however, lies beyond the capabilities of current astronomical facilities and nearly all proposed future facilities. Interferometers can crack this angular resolution problem, and space-based interferometry missions promise to explore entirely new regions of scientific phase space, providing unique new insights into the physical processes lurking at small angular scales.

Rinehart, S. A.↗

NASA NTRS1982

High spatial resolution Mg/Al maps of the western Crisium and Sulpicius Gallus regions

High spatial resolution Mg/Al ratio maps of the western Crisium and Sulpicius Gallus regions of the moon are presented. The data is from the X-ray fluorescence experiment and the image enhancement technique in the Laplacian subtraction method using a special least-squares version of the Laplacian to reduce noise amplification. In the highlands region west of Mare Crisium several relatively small patches of smooth material have high local Mg/Al ratio similar to values found in mare sites, suggesting volcanism in the highlands. In the same highland region there were other smooth areas with no high Mg/Al local values and they are probably Cayley Formation material produced by impact mass wasting. The Sulpicius Gallus region has variable Mg/Al ratios. In this region there are several high Mg/Al ratio spots, two of which occur at the highland-mare interface. Another high Mg/Al ratio area corresponds to the Sulpicius Gallus Rima I region. The high Mg/Al ratio material in the Sulpicius Gallus region is probably pyroclastic.

Schonfeld, E.↗

DOE OSTI2020

Visualizing Charge Transport and Nanoscale Electrochemistry by Hyperspectral Kelvin Probe Force Microscopy

Charge-transport and electrochemical processes are heavily influenced by the local microstructure. Kelvin probe force microscopy (KPFM) is a widely used technique to map electrochemical potentials at the nanometer scale; however, it offers little information on local charge dynamics. In this study, we implement a hyperspectral KPFM approach for spatially mapping bias-dependent charge dynamics in timescales ranging from the sub-millisecond to the second regime. As a proof of principle, we investigate the role mobile surface charges play in a three-unit-cell LaAlO 3 /SrTiO 3 oxide heterostructure. We explore machine learning approaches to assist with visualization, pattern recognition, and interpretation of the information-rich data sets. Linear unmixing methods reveal hidden bias-dependent interfacial processes, most likely water splitting, which are essentially unnoticed by functional fitting of the dynamic response alone. Hyperspectral KPFM will be beneficial for investigating nanoscale charge transport and local reactivity in systems involving a possible combination of electronic, ionic, and electrochemical phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Stanford Thermal Earth Model for the Conterminous United States

Provided here are various forms of the Stanford Thermal Earth Model, as well as the data and methods used for its creation. The predictions produced by this model were visualized in two-dimensional spatial maps across the modeled depths (0-7 km) for the conterminous United States. The thermal earth model is made available as an application programming interface (API) and as feature layers on ArcGIS, which are both provided via links below. A data-driven spatial interpolation algorithm based on physics-informed graph neural networks was used to develop these national temperature-at-depth maps. The model satisfied the three-dimensional heat conduction law by predicting subsurface temperature, surface heat flow, and rock thermal conductivity. Many physical quantities, including bottomhole temperature, depth, geographic coordinates, elevation, sediment thickness, magnetic anomaly, gravity anomaly, gamma-ray flux of radioactive elements, seismicity, and electric conductivity were used as model inputs. Surface heat flow, temperature, and thermal conductivity predictions were constructed for depths of 0-7 km at an interval of 1 km with spatial resolution of 18 km2 per grid cell. The model showed superior temperature, surface heat flow and thermal conductivity mean absolute errors of 4.8C, 8.1 mW/m2 and 0.07 W/(C-m), respectively..

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Spatiotemporally Registered In-Situ and Ex-Situ Datasets for Laser-based Blown Powder Directed Energy Deposition

This dataset is comprised of in situ sensing data collected during laser-based, blown powder directed energy deposition (DED) of Inconel 718 representing eight different printing conditions: (1) nominal, (2) +15% scan speed, (3) +12% laser power, (4) +42% powder feed rate, (5) +100% jerk limit, (6) +10% layer height, (7) +20% hatch spacing, (8) +20% carrier gas flow. All eight DED builds constructed an identical test coupon geometry consisting of geometric features representative of industrial print requirements (e.g., bulk deposition, thin walls, overhangs). In situ data consists of xyz-coordinates (100 Hz) and on-axis melt pool camera video (60 Hz), both of which have been temporally synchronized to spatially map the melt pool camera data. In addition, post-build X-ray computed tomography (XCT) data for each of the eight test geometries have been spatially registered to the recorded xyz-coordinates, allowing for comparisons between melt pool camera data and flaws identified in the XCT data.

additive manufacturing↗

NASA NTRS1991

Digital array scanned interferometers for astronomy

Investigations are reported of digital array scanned interferometers (DASI) with silicon CCD array detectors to define the operational capabilities of a mapping (polarimetric) spectrometer for astronomical applications based on these instruments. For spectral mapping, spatially resolved spectra using a cylindrical lens to image in the interferometers's redundant coordinate are given. The signal-to-noise characteristics of the Fourier transformed data are demonstrated with regard to the effects of a rectangular sampling function, spectral multiplexing and the pixel-to-pixel variation of the CCD array. These data indicate that DASIs can offer simple, versatile (polarization) mapping spectrometers suitable for spectral mapping observations from the ultraviolet to the infrared of extended sources at variable spatial resolution, particularly where long term stable operation is essential, as for spacecraft instruments.

Smith, Wm. Hayden↗

NASA NTRS2015

Mission Adaptive UAS Platform for Earth Science Resource Assessment

NASA Ames Research Center has led a number of important Earth science remote sensing missions including several directed at the assessment of natural resources. A key asset for accessing high risk airspace has been the 180 kg class SIERRA UAS platform, providing mission durations of up to 8 hrs at altitudes up to 3 km. Recent improvements to this mission capability are embodied in the incipient SIERRA-B variant. Two resource mapping problems having unusual mission characteristics requiring a mission adaptive capability are explored here. One example involves the requirement for careful control over solar angle geometry for passive reflectance measurements. This challenges the management of resources in the coastal ocean where solar angle combines with sea state to produce surface glint that can obscure the ocean color signal. Furthermore, as for all scanning imager applications, the primary flight control priority to fly the UAS directly to the next waypoint should compromise with the requirement to minimize roll and crab effects in the imagery. A second example involves the mapping of natural resources in the Earth's crust using precision magnetometry. In this case the vehicle flight path must be oriented to optimize magnetic flux gradients over a spatial domain having continually emerging features, while optimizing the efficiency of the spatial mapping task. These requirements were highlighted in several recent Earth Science missions including the October 2013 OCEANIA mission directed at improving the capability for hyperspectral reflectance measurements in the coastal ocean, and the Surprise Valley Mission directed at mapping sub-surface mineral composition and faults, using high-sensitivity magentometry. This paper reports the development of specific aircraft control approaches to incorporate the unusual and demanding requirements to manage solar angle, aircraft attitude and flight path orientation, and efficient (directly geo-rectified) surface and sub-surface mapping, including the near-time optimization of these sometimes conflicting requirements. *

Flight Control↗

DOE OSTI2023

Spatially Resolved Top-Down Proteomics of Tissue Sections Based on a Microfluidic Nanodroplet Sample Preparation Platform

Conventional proteomics measures the averaged signal from mixed cell populations or bulk tissues, leading to the dilution of significant changes in subpopulations of cells that might serve as important biomarkers. Recent developments in bottom-up proteomics have enabled spatial mapping of cellular heterogeneity in tissue microenvironments. However, bottom-up proteomics cannot precisely infer the abundance changes of intact proteins, which are presented as proteoforms. Herein, we described a spatially resolved top-down proteomics (TDP) platform for proteoform identification and quantification directly on thin tissue sections. The spatial TDP platform consisted of a nanoPOTS (nanodroplet Processing in One pot for Trace Samples)-based sample preparation system and an LCM (laser capture microdissection)-based cell isolation system. We improved the nanoPOTS sample preparation by adding benzonase in the extraction buffer to enhance the coverage of nucleus proteins. Using ~200 cultured cells as model samples, the improved approach increased proteoform identifications from 493 to 700; newly identified proteoforms primarily corresponded to nuclear proteins. To demonstrate the spatial TDP platform in tissue samples, we analyzed LCM-isolated tissue voxels from rat brain cortex and hypothalamus regions. We quantified 426 proteoforms by combining identifications from TopPIC and TDPortal with the quantitation from ProMex. Several proteoforms corresponding to the same gene exhibited mixed abundance profiles between two tissue regions, suggesting potential PTM-specific spatial distributions. The spatial TDP workflow has prospects for biomarker discovery at proteoform level from small tissue sections.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2003

Jupiter Atmospheric Science with JIMO: Linking Science Objectives and Measurement Goals

Although the primary focus of the Jupiter Icy Moons Orbiter (JIMO) mission will be the characterization and study of Jupiter's icy moons, there will be opportunities throughout the mission for unprecendented observations of Jupiter. With an adaptable suite of payload instruments, the atmospheric data collected by JIMO can help to answer fundamental questions about the largest planet in our solar system that remain after (or were generated by) previous spacecraft reconnaissance (e.g. Voyager, Galileo, and Cassini). Near-IR (0.7-4 micron) spectral imaging will most likely be used to identify mineralogies and ices on the Jovian satellites by virtue of their spectral signatures. This same capability is very well tailored for studies of Jovian atmospheric dynamics and structure. Near-IR methane absorption bands allow 2-D mapping of the horizontal wind field at size scales to tens of kms, as well as the height dependence of this field above the ammonia cloud deck (700 to a few mbar), constraining current models of atmospheric vertical structure. Likewise, atmospheric ice aerosols with unique spectroscopic signatures (ammonia ice near 1.5, 2.0, and 2.8 microns and water ice between 3.0 - 3.5 microns) can be detected and mapped using spectral difference imaging or spectrally inclusive principal-component methods. Spectral imaging of the Jovian aurora via (3)H(+) emission lines between 3 - 4 microns can be used to spatially map the interplay between the satellites) Jupiter's magnetosphere, and Jupiter's atmosphere. Each of these measurements addresses one or more fundamental questions related to the energy balance in Jupiter's atmosphere. All of these tunable imaging objectives can be achieved using acousto-optic tunable filters (AOTF's), which have been used for years in ground-based observing instruments and which have been proposed for numerous planetary missions. The application of this technology to the science objectives of both the icy satellites and Jovian atmospheric components of the JIMO mission will be discussed.

Chanover, N. J.↗

NASA NTRS1994

Spatial patterns in the length of the sea ice season in the Southern Ocean, 1979-1986

The length of the sea ice season summarizes in one number the ice coverage conditions for an individual location for an entire year. It becomes a particularly valuable variable when mapped spatially over a large area and examined for regional and interannual differences, as is done here for the Southern Ocean over the years 1979-1986, using the satellite passive microwave data of the Nimbus 7 scanning multichannel microwave radiometer. Three prominent geographic anomalies in ice season lengths occur consistently in each year of the data set, countering the general tendency toward shorter ice seasons from south to north: (1) in the Weddell Sea the tendency is toward shorter ice seasons from southwest to northeast, reflective of the cyclonic ice/atmosphere/ocean circulations in the Weddell Sea region. (2) Directly north of the Ross Ice Shelf anomalously short ice seasons occur, lasting only 245-270 days, in contrast to the perennial ice coverage at comparable latitudes in the southern Bellingshausen and Amundsen Seas and in the western Weddell Sea. The short ice season off the Ross Ice Shelf reflects the consistently early opening of the ice cover each spring, under the influence of upwelling along the continental slope and shelf and atmospheric forcing from winds blowing off the Antarctic continent. (3) In the southern Amundsen Sea, anomalously short ice seasons occur adjacent to the coast, owing to the frequent existence of coastal polynyas off the many small ice shelves bordering the sea. Least squares trends in the ice season lengths over the 1979-1986 period are highly coherent spatially, with overall trends toward shorter ice seasons in the northern Weddell and Bellingshausen seas and toward longer ice seasons in the Ross Sea, around much of East Antarctica, and in a portion of the south central Weddell Sea.

Parkinson, Claire L.↗

NASA NTRS2016

Titan Polar Landscape Evolution

With the ongoing Cassini-era observations and studies of Titan it is clear that the intensity and distribution of surface processes (particularly fluvial erosion by methane and Aeolian transport) has changed through time. Currently however, alternate hypotheses substantially differ among specific scenarios with respect to the effects of atmospheric evolution, seasonal changes, and endogenic processes. We have studied the evolution of Titan's polar region through a combination of analysis of imaging, elevation data, and geomorphic mapping, spatially explicit simulations of landform evolution, and quantitative comparison of the simulated landscapes with corresponding Titan morphology. We have quantitatively evaluated alternate scenarios for the landform evolution of Titan's polar terrain. The investigations have been guided by recent geomorphic mapping and topographic characterization of the polar regions that are used to frame hypotheses of process interactions, which have been evaluated using simulation modeling. Topographic information about Titan's polar region is be based on SAR-Topography and altimetry archived on PDS, SAR-based stereo radar-grammetry, radar-sounding lake depth measurements, and superposition relationships between geomorphologic map units, which we will use to create a generalized topographic map.

Cassini↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2023

Equivalence between Fermion-to-Qubit Mappings in two Spatial Dimensions

We argue that all locality-preserving mappings between fermionic observables and Pauli matrices on a two-dimensional lattice can be generated from the exact bosonization in Chen et al., whose gauge constraints project onto the subspace of the toric code with emergent fermions. Starting from the exact bosonization and applying Clifford finite-depth generalized local unitary transformation, we can achieve all possible fermion-to-qubit mappings (up to the re-pairing of Majorana fermions). In particular, we discover a new supercompact encoding using 1.25 qubits per fermion on the square lattice. We prove the existence of finite-depth quantum circuits to obtain fermion-to-qubit mappings with qubit-fermion ratios r = 1 + 1/2k for positive integers k, utilizing the trivialness of quantum cellular automata in two spatial dimensions. Also, we provide direct constructions of fermion-to-qubit mappings with ratios arbitrarily close to 1. When the ratio reaches 1, the fermion-to-qubit mapping reduces to the one-dimensional Jordan-Wigner transformation along a certain path in the two-dimensional lattice. Finally, we explicitly demonstrate that the Bravyi-Kitaev superfast simulation, the Verstraete-Cirac auxiliary method, Kitaev’s exactly solved model, the Majorana loop stabilizer codes, and the compact fermion-to-qubit mapping can all be obtained from the exact bosonization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Enzymatically hydrolyzed fluorescence-based chemical probe enables in situ mapping of chitinase activity in the rhizosphere

Rhizosphere microbiomes and root exudates play a pivotal biochemical role in helping to catalyze chitin catabolism. Chitin is a recalcitrant and ubiquitous soil biopolymer, estimated to be the second most abundant organic soil biopolymer on Earth. Despite its abundance, role as a source of C and N in soil, and importance to ecosystem function, the biochemical mechanisms controlling chitin fate in the rhizosphere are elusive and poorly understood. To enable spatial mapping of chitinase activity in the rhizosphere, we designed and synthesized an enzymatically activated fluorogenic substrate, chitotriose-TokyoGreen (chitotriose-TG), by incorporating a fluorescein derivative (TG) onto the trimeric unit of chitin. This non-fluorescent substrate is selectively hydrolyzed by chitinase to release TG and yield a high fluorescence signal, which can be used to spatially image and measure chitinase activity in the rhizosphere. To demonstrate the application of this technique, we grew switchgrass (Panicum virgatum) in rhizoboxes amended with a horizontal layer of chitin. We extracted mobile proteins from the rhizobox using a nitrocellulose membrane blotting technique which offer non-destructive enzyme extraction while preserving the 2D spatial position of the enzymes. We then subjected these membranes to the synthesized chitotriose-TG stain to spatially visualize the distribution of chitinase activity within the rhizosphere. Furthermore, we observed increased chitinase activity near plant roots and higher activity within the soil zone enriched in chitin, showing an adaptive response of chitinase activity with spatial focusing in areas of higher chitin abundance. Thus, the enzyme extraction and visualization strategy we describe here can help enlighten efforts to better understand spatial controls on chitin breakdown in rhizosphere, further elucidating the role of chitin as a C and N source in these systems.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability: Supporting Data and Code

This repository contains R code and associated datasets for reproducing the analysis described in the manuscript titled “Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability” (DOI: 10.1029/2024JG008604). The provided scripts facilitate a comprehensive analysis of snow depth variability influenced by microtopography and vegetation distribution in a subarctic watershed. Included datasets are high-resolution spatial maps of snow depth, terrain elevation, vegetation height, and distance from shrubs taller than 1 meter, all formatted as text files (.txt). These data are fully describe in doi:10.15485/2316038. Users can adapt the provided R scripts to accommodate different data formats or larger spatial domains, noting that some output files may require modification due to their size.The code includes implementations for boosted regression tree analysis adapted from methods outlined in Elith et al. (2008). Users interested in understanding or modeling landscape-scale snow distribution patterns, particularly in Arctic or subarctic ecosystems, will find this package useful. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗