Search NASA⌕ Search

SEARCH · Search NASA

Results for “projected gradient”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

Novel strategies for modal-based structural material identification

Here, we present modal-based methods for model calibration in structural dynamics, and address several key challenges in the solution of gradient-based optimization problems with eigenvalues and eigenvectors, including the solution of singular Helmholtz problems encountered in sensitivity calculations, non-differentiable objective functions caused by mode swapping during optimization, and cases with repeated eigenvalues. Unlike previous literature that relied on direct solution of the eigenvector adjoint equations, we present a parallel iterative domain decomposition strategy (Adjoint Computation via Modal Superposition with Truncation Augmentation) for the solution of the singular Helmholtz problems. For problems with repeated eigenvalues we present a novel Mode Separation via Projection algorithm, and in order to address mode swapping between inverse iterations we present a novel Injective mode ordering metric. We present the implementation of these methods in a massively parallel finite element framework with the ability to use measured modal data to extract unknown structural model parameters from large complex problems. A series of increasingly complex numerical examples are presented that demonstrate the implementation and performance of the methods in a massively parallel finite element framework [7], [5], using gradient-based optimization techniques in the Rapid Optimization Library (ROL) [21].

36 MATERIALS SCIENCE↗

Extended MHD modeling of nonlinear instabilities in fusion and space plasmas

A number of different sub-projects where pursued within this DOE early career project. The primary focus was on using fully nonlinear, curvilinear, extended MHD simulations of instabilities with applications to fusion and space plasmas. In particular, we performed comprehensive studies of the dynamics of the double tearing mode in different regimes and confi gurations, using Cartesian and cyclindrical geometry and investigating both linear and non-linear dynamics. In addition to traditional extended MHD involving Hall term and electron pressure gradient, we also employed a new multi-fluid moment model, which shows great promise to incorporate kinetic effects, in particular off-diagonal elements of the pressure tensor, in a fluid model, which is naturally computationally much cheaper than fully kinetic particle or Vlasov simulations. We used our Vlasov code for detailed studies of how weak collisions effect plasma echos. In addition, we have played an important supporting role working with the PPPL theory group around Will Fox and Amitava Bhattacharjee on providing simulation support for HED plasma experiments performed at high-powered laser facilities like OMEGA-EP in Rochester, NY. This project has support a great number of computational advances in our fluid and kinetic plasma models, and has been crucial to winning multiple INCITE computer time awards that supported our computational modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of boundary-layer wind predictions during nocturnal low-level jet events using the Weather Research and Forecasting model

Simulations of two periods featuring three consecutive low level jet (LLJ) events in the US Upper Great Plains during the autumn of 2011 were conducted to explore the impacts of various setup configurations and physical process models on simulated flow parameters within the lowest 200 m above the surface, using the Weather Research and Forecasting (WRF) model. Sensitivities of simulated flow parameters to the horizontal and vertical grid spacing, planetary boundary layer (PBL) and land surface model (LSM) physics options, were assessed. Data from a Light Detection and Ranging (lidar) system, deployed to the Weather Forecast Improvement Project (WFIP; Finley et al. 2013) were used to evaluate the accuracy of simulated wind speed and direction at 80 m above the surface, as well as their vertical distributions between 120 and 40 m, covering the typical span of contemporary tall wind turbines. All of the simulations qualitatively captured the overall diurnal cycle of wind speed and stratification, producing LLJs during each overnight period, however large discrepancies occurred at certain times for each simulation in relation to the observations. 54-member ensembles encompassing changes of the above discussed configuration parameters displayed a wide range of simulated vertical distributions of wind speed and direction, and potential temperature, reflecting highly variable representations of stratification during the weakly stable overnight conditions. Root mean square error (RMSE) statistics show that different ensemble members performed better and worse in various simulated parameters at different times, with no clearly superior configuration . Simulations using a PBL parameterization designed specifically for the stable conditions investigated herein provided superior overall simulations of wind speed at 80 m, demonstrating the efficacy of targeting improvements of physical process models in areas of known deficiencies. However, the considerable magnitudes of the RMSE values of even the best performing simulations indicate ample opportunities for further improvements.

Mirocha, Jeff D.↗

NSA Site Science: Use of ARM Observations from Northern Alaska to Evaluate and Improve Prediction Capabilities

The Arctic is warming at a rate nearly double that of the rest of the planet, leading to profound changes in atmospheric, oceanic, and ice processes. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility has played a significant role in Arctic research, operating observatories in Alaska's North Slope for over 25 years. These observatories provide a rich, and wide-reaching dataset that offers insight into atmospheric processes in northern Alaska. This report details the results of a nine-year research project (2015–2024) supported by the DOE Atmospheric Systems Research (ASR) program, that leverages data from ARM’s deployment of observing facilities at Utqiaġvik (known as the North Slope of Alaska, or NSA, site) and Oliktok Point, Alaska. The project was conducted in two phases: - Phase 1 (2015–2019): Focused on understanding key atmospheric processes at Oliktok Point, including cloud formation, high-latitude precipitation, aerosol-cloud interactions, and cloud properties. - Phase 2 (2019–2024): Extended the research to the broader North Slope region, using data from both Oliktok Point and Utqiaġvik. Topics explored included surface energy budgets, atmospheric stability, ice nucleation processes, and microphysics in Arctic clouds. The project resulted in numerous research products, including 50 peer-reviewed publications and dissertations, 169 presentations, and 10 data products. These products cover a variety of topics, including: - Cloud Macro- and Microphysical Properties: Arctic clouds play a crucial role in energy transfer, and accurate representation in models is critical. The study explored cloud transitions, ice crystal shapes, and dual-wavelength radar data to understand ice crystal habits and size distributions. - Aerosol Properties and Processes: The team examined aerosol sources in the Arctic, including industrial emissions and natural sources. Observations showed significant spatial gradients in aerosol concentrations due to human activities and wildfire smoke. The influence of aerosols on cloud formation and the surface energy budget was also assessed. - Aerosol-Cloud Interactions: Research revealed that aerosols might suppress cloud ice production, affecting cloud radiative forcing and precipitation. The impact of local industrial emissions on cloud properties was also investigated. - Contextualizing the North Slope of Alaska in the context of the broader Arctic: To understand broader trends, the project evaluated large-scale circulation patterns and the influence of weather systems on the Arctic. Studies indicated that large-scale processes play a significant role in temperature patterns and the timing of snowmelt. - Advancing ARM Observational and Modeling Capabilities: The project developed new radar data products and advanced measurement techniques, including clutter mitigation and drizzle detection. Uncrewed aerial systems (UAS) and tethered balloon systems (TBS) were deployed to gather detailed atmospheric data. Additionally, the project supported 10 early career scientists, providing training and mentorship to undergraduate interns, graduate students, postdoctoral researchers, and early career researchers. These efforts contributed to the advancement of ARM research capabilities and fostered a new generation of scientists skilled in Arctic atmospheric research. Ultimately, this ASR-supported project has provided valuable insights into Arctic atmospheric processes and their broader climate implications. Recommendations for future work include continuing support for long-term observing at Arctic locations to foster additional research, further exploration of aerosol-cloud interactions and the potential impacts of enhanced industrialization of the Arctic, and expanded use of uncrewed systems to gather data in this remote and harsh environment. Additionally, the data products developed by this work, and the data products developed through the ARM infrastructure, leave a treasure-trove of additional information that should be explored for many years to come to gain additional insight into physical processes in the Arctic atmosphere that drive the rapid changes occurring in at high latitudes and their global impact.

58 GEOSCIENCES↗

Age gradients in the stellar populations of massive star forming regions based on a new stellar chronometer

A major impediment to understanding star formation in massive star-forming regions (MSFRs) is the absence of a reliable stellar chronometer to unravel their complex star formation histories. We present a new estimation of stellar ages using a new method that employs near-infrared (NIR) and X-ray photometry, Age {sub JX} . Stellar masses are derived from X-ray luminosities using the L{sub X} -M relation from the Taurus cloud. J-band luminosities are compared to mass-dependent pre-main-sequence (PMS) evolutionary models to estimate ages. Age {sub JX} is sensitive to a wide range of evolutionary stages, from disk-bearing stars embedded in a cloud to widely dispersed older PMS stars. The Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) project characterizes 20 OB-dominated MSFRs using X-ray, mid-infrared, and NIR catalogs. The Age {sub JX} method has been applied to 5525 out of 31,784 MYStIX Probable Complex Members. We provide a homogeneous set of median ages for over 100 subclusters in 15 MSFRs; median subcluster ages range between 0.5 Myr and 5 Myr. The important science result is the discovery of age gradients across MYStIX regions. The wide MSFR age distribution appears as spatially segregated structures with different ages. The Age {sub JX} ages are youngest in obscured locations in molecular clouds, intermediate in revealed stellar clusters, and oldest in distributed populations. The NIR color index J – H, a surrogate measure of extinction, can serve as an approximate age predictor for young embedded clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Enhanced jet stream waviness induced by suppressed tropical Pacific convection during boreal summer

Abstract Consensus on the cause of recent midlatitude circulation changes toward a wavier manner in the Northern Hemisphere has not been reached, albeit a number of studies collectively suggest that this phenomenon is driven by global warming and associated Arctic amplification. Here, through a fingerprint analysis of various global simulations and a tropical heating-imposed experiment, we suggest that the suppression of tropical convection along the Inter Tropical Convergence Zone induced by sea surface temperature (SST) cooling trends over the tropical Eastern Pacific contributed to the increased summertime midlatitude waviness in the past 40 years through the generation of a Rossby-wave-train propagating within the jet waveguide and the reduced north-south temperature gradient. This perspective indicates less of an influence from the Arctic amplification on the observed mid-latitude wave amplification than what was previously estimated. This study also emphasizes the need to better predict the tropical Pacific SST variability in order to project the summer jet waviness and consequent weather extremes.

54 ENVIRONMENTAL SCIENCES↗

Scaling FFAG accelerator for muon acceleration

Recent developments in scaling fixed field alternating gradient (FFAG) accelerators have opened new ways for lattice design, with straight sections, and insertions like dispersion suppressors. Such principles and matching issues are detailed in this paper. An application of these new concepts is presented to overcome problems in the PRISM project.

43 PARTICLE ACCELERATORS↗

Support Structure Design of the $\hbox{Nb}_{3}\hbox{Sn}$ Quadrupole for the High Luminosity LHC

New low-β quadrupole magnets are being developed within the scope of the High Luminosity LHC (HL-LHC) project in collaboration with the US LARP program. The aim of the HLLHC project is to study and implement machine upgrades necessary for increasing the luminosity of the LHC. The new quadrupoles, which are based on the Nb₃Sn superconducting technology, will be installed in the LHC Interaction Regions and will have to generate a gradient of 140 T/m in a coil aperture of 150 mm. In this paper, we describe the design of the short model magnet support structure and discuss results of the detailed 3D numerical analysis performed in preparation for the first short model test.

43 PARTICLE ACCELERATORS↗

Cyclotrons and FFAG Accelerators as Drivers for ADS

Our review summarizes projects and studies on circular accelerators proposed for driving subcritical reactors. The early isochronous cyclotron cascades, proposed about 20 years ago, and the evolution of these layouts up to the most recent solutions or designs based on cyclotrons and fixed field alternating gradient accelerators, are reported. Additionally, the newest ideas and their prospects for development are discussed.

43 PARTICLE ACCELERATORS↗

Final Report - Impact of Microstructure on the Containment and Migration of CO 2 in Fractured Basalts

The overall objective of the project was to advance the scientific and technical understanding of fracture microstructure on the flow and mineralization of CO 2 injected into fractured basalt. The project combined bench-scale CO 2 -water-rock testing, geochemical characterization of rock cores, and advanced characterization of the evolution of fracture structure and carbon trapping mechanisms. A series of natural and artificial basalts were acquired and characterized for use in experiments. High pressure and high temperature experiments were performed in both batch modes and in flow-through systems. The dissolution of silicate minerals in the basalts was tracked, and the timing, identity and location of carbonate minerals that trap that added CO 2 were determined. A reactive transport model was developed that generated simulations that were consistent with the major trends in the experimental data. Novel X-ray computed tomography and nuclear magnetic resonance tools were developed as part of the project. The project found that the formation of carbonate minerals in basalts subjected to injection of CO 2 -rich aqueous solutions can be rapid and extensive. The mineral trapping rate and capacity of porous basalt is much greater than that of sandstone reservoirs. Geochemical gradients that develop upon the exposure of basalt to CO 2 -rich aqueous solutions lead to spatial localization of carbonate minerals in confined environments (e.g., narrow fractures, dead-end fractures, and vesicles). While carbonate minerals form in dead-end fractures, they do not form in ways that inhibit subsequent carbonation of the basalt.

58 GEOSCIENCES↗

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗

Synthesis of finite displacements and displacements in continental margins

The scope of the project is the analysis of displacement-rate fields in the transitional regions between cratonal and oceanic lithospheres over Phanerozoic time (last 700 ma). Associated goals are an improved understanding of range of widths of major displacement zones; the partition of displacement gradients and rotations with position and depth in such zones; the temporal characteristics of such zones-the steadiness, episodicity, and duration of uniform versus nonunifrom fields; and the mechanisms and controls of the establishment and kinematics of displacement zones. The objective is to provide a context of time-averaged kinematics of displacement zones. The initial phase is divided topically among the methodology of measurement and reduction of displacements in the lithosphere and the preliminary analysis from geologic and other data of actual displacement histories from the Cordillera, Appalachians, and southern North America.

Speed, R. C.↗

The Changing California Coast: The Effect of a Variable Water Budget on Coastal Vegetation Succession

The land-ocean interface along the central coast of California is one of the most diverse biogeographic regions of the state. This area is composed of a species-rich mosaic of coastal grassland, shrubland, and forest vegetation types. An acceleration of conifer encroachment into shrublands and shrub encroachment into grasslands along the coast has been recently documented. These vegetation changes are believed to be driven primarily by fire suppression and changing grazing patterns. Climatic variables such as precipitation, fog, cloud cover, temperature, slope, and elevation also play an important role in vegetation succession. Our study area is located along the central California coast, which is characterized by a precipitation gradient from the relatively wetter and cooler north to the drier and warmer south. Some studies indicate changing fog patterns along this coast, which may greatly impact vegetation. A decrease in water availability could slow succession processes. The primary objective of this project is to determine if vegetation succession rates are changing for the study area and to identify climate and ecosystem variables which contribute to succession, specifically the transition among grassland, shrubland, and forest. To identify vegetation types and rates of succession, we classified two Landsat TM 5 scenes from 1985 to 2010 with a resulting overall accuracy of 82.4%. Vegetation succession was correlated to changes in maximum and minimum temperatures, precipitation, and elevation for each sub-region of the study area. Fog frequency was then compared between the northern and southern regions of the study area for determining the spatial relation between fog frequency and the percent of vegetation change.

Hsu, Wei-Chen↗

Use of Physics Based Management Optimization for Groundwater Remedy Design - 18671

Numerical modeling coupled with formal mathematical optimization techniques was used to design a groundwater remedy strategy to address a carbon tetrachloride (CCl{sub 4}) plume covering several square miles. The optimization approach, methodology, and project costs are described. The site consists of a large grain elevator in Nebraska. CCl{sub 4} was used as a fumigant and subsequently impacted soils and groundwater. Remedy design focused on the more diffuse plume, which extended approximately 4.5 miles down gradient of the source area. An eleven-layer model was developed and calibrated to observed groundwater conditions. This model was linked with formal optimization techniques to identify the best design capable of meeting client specified requirements and well location limitations. The planned system included requirements on the location of extraction and injection wells, total system flow rate, number of management periods, and specific constraints limiting further plume expansion. A separate trinitrotoluene (TNT) contaminant plume from a third party was located vertically above the leading edge of the CCl{sub 4} plume within the same aquifer. A further constraint imposed was that the design could not affect or cause additional migration of the TNT plume. Optimal system design was determined using HydroGeoLogic's Physics Based Management Optimization (PBMO{sup TM}) system. Over 38,000 candidate solutions consisting of different combinations of pumping and injection wells locations and rates, were evaluated through flow modeling coupled with either particle tracking or fully convective-dispersive transport simulations to determine the optimal design. Planning for installation of the optimal design is underway. Simulations were completed on a distributed computing platform consisting of 3 Quad Core PCs. Flow and particle tracking simulations required less than 2 minutes per simulation. Coupled flow and transport simulation run times were generally between 30-40 minutes per simulation. The total costs for developing the groundwater flow model and completing the subsequent optimal design were approximately $400,000. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Implementation of Enhanced Attenuation at the DOE Mound Site OU-1 Landfill: Accelerating Progress and Reducing Costs - 16270

A field demonstration of enhanced attenuation (EA) was initiated in 2014 to expedite the remediation of industrial solvents in the groundwater impacted by the former landfill Operable Unit 1 (OU-1) at the US Department of Energy (DOE) Mound, Ohio, Site. EA uses active engineering solutions to alter the target site so that the contaminant plume will passively stabilize and shrink. For OU-1, the EA was implemented by targeted injection of an electron donor (vegetable oil) to create 'structured geochemical zones'. Volatile organic compounds (VOCs) in the altered subsurface system encounter alternating anaerobic and aerobic environments along the plume flow path. The effect of these structured geochemical zones is effective and rapid degradation of parent compounds such as trichloroethene (TCE) and tetrachloroethene (PCE), and daughter compounds such as cis-1,2-dichloroethene and vinyl chloride. EA at Mound OU-1 provides a transition to natural attenuation and is an alternative to the baseline pump-and-treat (P and T) system. Industrial solvents (primarily TCE) and other VOCs that originated from the former solid-waste landfill contaminated the groundwater in the Buried Valley Aquifer beneath the Mound site. A groundwater P and T system was initiated in 1996 to control the plume, and a soil vapor extraction system was installed and operated from 1997 to 2003 to remove the bulk of the VOC sources from the vadose zone. Physical removal of the landfill waste and contaminated soil was performed between 2006 and 2010. Approximately 76,100 cubic meters (99,500 cubic yards) of material were removed from the OU-1 landfill area. Operation of the groundwater P and T system continued throughout this period. The principles and technical basis of the ongoing EA field demonstration were developed by the Interstate Technology and Regulatory Council. The goal of the EA implementation is to provide a bridge between active groundwater treatments and monitored natural attenuation. EA relies on understanding and manipulating the contaminant mass balance. Technologies are applied to sustainably reduce the contaminant loading from the source and/or increase the rate of natural attenuation processes within the plume. In 2014, a core team that included DOE, the US Environmental Protection Agency (EPA), the Ohio EPA, and technical support organizations evaluated remedial options for OU-1. During the 2014 evaluation, the team projected a cleanup time frame for the baseline P and T system of 26 additional years and a cleanup time frame for the EA remedy of 13 years. On the basis of that evaluation, EPA approved the implementation of the EA remedy that required the shutdown of the P and T system. The combination of technologies implemented for OU-1 EA included (a) neat vegetable oil in the vadose zone in the former landfill and (b) emulsified vegetable oil within the footprint of the groundwater plume. In the first part of the deployment, neat oil spreads laterally and forms a thin layer on the water table beneath residual soil sources to intercept and reduce future VOC loading and to reduce oxygen inputs to the local groundwater. In the second part of the deployment, emulsified oil forms active bioremediation anaerobic treatment zones within the plume footprint to degrade existing groundwater contaminants (via reductive dechlorination, co-metabolism, or both) and stimulates long-term attenuation capacity in the down-gradient aerobic plume (via co-metabolism). The results from the first year of the field demonstration indicate that the dissolved TCE and PCE plumes have decreased in size and mass. Statistical tests indicate that the VOC concentrations in the majority of the source-area wells are decreasing. The VOC concentrations throughout the plume are projected to be below the target maximum contaminant level of 2 to 5 μg/L within 5 to 10 years-sooner than the design projection of 13 years. The biogeochemistry and microbial community within the structured geochemical zones are consistent with the design and appear to be stable. At Mound OU-1, EA has accelerated progress toward the remedial action objectives and reduced costs. Ongoing monitoring will determine whether the EA remedy has transitioned the site into monitored natural attenuation status and whether the EA will continue to be effective, timely, and sustainable. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Extending Parsimonious Bayesian Inference

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

97 MATHEMATICS AND COMPUTING↗

Gyrokinetic projection of the divertor heat-flux width from present tokamaks to ITER

Here, the XGC1 edge gyrokinetic code is used to study the width of the heat-flux to divertor plates in attached plasma condition. The flux-driven simulation is performed until an approximate power balance is achieved between the heat-flux across the steep pedestal pressure gradient and the heat-flux on the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nitrogen fixing shrubs advance the pace of tall-shrub expansion in low-Arctic tundra

Abstract Tall deciduous shrubs are critically important to carbon and nutrient cycling in high-latitude ecosystems. As Arctic regions warm, shrubs expand heterogeneously across their ranges, including within unburned terrain experiencing isometric gradients of warming. To constrain the effects of widespread shrub expansion in terrestrial and Earth System Models, improved knowledge of local-to-regional scale patterns, rates, and controls on decadal shrub expansion is required. Using fine-scale remote sensing, we modeled the drivers of patch-scale tall-shrub expansion over 68 years across the central Seward Peninsula of Alaska. Models show the heterogeneous patterns of tall-shrub expansion are not only predictable but have an upper limit defined by permafrost, climate, and edaphic gradients, two-thirds of which have yet to be colonized. These observations suggest that increased nitrogen inputs from nitrogen-fixing alders contributed to a positive feedback that advanced overall tall-shrub expansion. These findings will be useful for constraining and projecting vegetation-climate feedbacks in the Arctic.

54 ENVIRONMENTAL SCIENCES↗