Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 433 records · Page 24

A twist grain boundary phase in aqueous solutions of the nucleic acid tetramer GTAC

At high concentration, long Watson/Crick (WC) double-helixed DNA forms columnar crystal or liquid crystal phases of linear, parallel duplex chains packed on periodic lattices. This can also be a structural motif of short NA oligomers such as the 5’-GTAC-3’ studied here, which makes four-base WC duplexes having hydrophobic blunt ends. End-to-end aggregation then assembles these duplexes into columns and columnar phases are stabilized, in spite of breaks in the double helix every four bases. But the new degrees of freedom introduced by such breaks also enable opportunities for a more diverse palette of self-assembly modes, producing striking self-assemblies of DNA that would not be achievable with contiguous polymers. These include recently reported three-dimensional (3D) periodic low-density nanoscale networks of GCCG, and the twist grain boundary (TGB) phase presented here. In the TGB, columns of GTAC pairs assemble into monolayer sheets in which the duplex columns are mutually parallel. However, unlike in the columnar crystals, these sheets stack in helical fashion into lamellar arrays in which the column axis of each layer is rotated through a 60° angle with respect to the columns in neighboring layers. This assembly of DNA is unique in that it the fills a 3D volume wherein the major grooves of columns in each layer mutually enter and interlock with the major grooves of columns in neighboring layers. This locking is optimized by small adjustments in structure enabled by the breaks in the duplex backbones.

DNA↗

Unveiling and Mapping Polymorphs in Fluorite Y2TiO5 Using 4D-STEM and Unsupervised Machine Learning

Y2TiO5 belongs to the Ln2TiO5 (Ln = lanthanide or Y) family of ceramic materials and exhibits a range of desirable material properties such as radiation tolerance, frustrated magnetism, and large dielectric constant. However, understanding the complex crystal structure of Y2TiO5 remains elusive, given that Y2TiO5 can adopt multiple polymorphs such as cubic, orthorhombic, and hexagonal phases within the lattice. In this work, we report a detailed structural analysis of Y2TiO5 using four-dimensional scanning transmission electron microscopy coupled with unsupervised machine learning. The pyrochlore nanodomains, characterized by the ordered arrangement of yttrium cations on the A site of their A2BO5 structure, are present within the matrix of a predominantly fluorite-structured Y2TiO5 along with a third polymorph, the hexagonal phase. The pyrochlore phase is found to form 2 nm boundary regions around hexagonal phase stacking faults, highlighting the potential influence of the hexagonal phase on the occurrence and distribution of the pyrochlore phase. Lastly, we identify a unique pyrochlore phase with asymmetric arrangement of cation ordering along a single planar direction. Our findings provide invaluable insights into the possible mechanisms stabilizing pyrochlore nanodomains within the fluorite lattice of Y2TiO5.

36 MATERIALS SCIENCE↗

Exploring Local Ionic Motion in Perovskite Solar Cells at Nanoscale through Scanning Thermo-Ionic Microscopy

We investigate the local ionic motion in the absorber layer of perovskite solar cells using scanning thermo-ionic microscopy (STIM). STIM images of perovskite films show higher STIM amplitude at most of the grain boundaries due to higher ionic motion at those regions. The perovskite absorber layer of devices aged under blue (450 nm) light soaking shows higher amplitude in STIM signal compared to controls kept in dark storage. Such enhancement of STIM amplitude upon aging under light implies an increased concentration or diffusivity (or both) of mobile ionic species in stressed devices. Our results demonstrate the utilization of the STIM technique to assess and understand the intrinsic local ionic motion in perovskite absorbers for further advancement in device stability and functionality.

aging↗

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating Global Terrestrial Carbon and Nitrogen Biogeochemical Cycles With Implicit and Explicit Representations of Soil Microbial Activity

Abstract Nutrient limitation is widespread in terrestrial ecosystems. Accordingly, representations of nitrogen (N) limitation in land models typically dampen rates of terrestrial carbon (C) accrual, compared with C‐only simulations. These previous findings, however, rely on soil biogeochemical models that implicitly represent microbial activity and physiology. Here we present results from a biogeochemical model testbed that allows us to investigate how an explicit versus implicit representation of soil microbial activity, as represented in the MIcrobial‐MIneral Carbon Stabilization (MIMICS) and Carnegie‐Ames‐Stanford Approach (CASA) soil biogeochemical models, respectively, influence plant productivity, and terrestrial C and N fluxes at initialization and over the historical period. When forced with common boundary conditions, larger soil C pools simulated by the MIMICS model reflect longer inferred soil organic matter (SOM) turnover times than those simulated by CASA. At steady state, terrestrial ecosystems experience greater N limitation when using the MIMICS‐CN model, which also increases the inferred SOM turnover time. Over the historical period, however, warming‐induced acceleration of SOM decomposition over high latitude ecosystems increases rates of N mineralization in MIMICS‐CN. This reduces N limitation and results in faster rates of vegetation C accrual. Moreover, as SOM stoichiometry is an emergent property of MIMICS‐CN, we highlight opportunities to deepen understanding of sources of persistent SOM and explore its potential sensitivity to environmental change. Our findings underscore the need to improve understanding and representation of plant and microbial resource allocation and competition in land models that represent coupled biogeochemical cycles under global change scenarios.

54 ENVIRONMENTAL SCIENCES↗

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Characterization of Flashback and Flame-Holding in a Jet-In-Crossflow Mixing Configuration with Methane-Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence (CL) diagnostics techniques. The JICF burner has an optically accessible pre-mixing tube allowing the optical diagnostics. The testing was conducted for varied pre-mixer velocities (V) and equivalence ratio (ϕ) for 90%-100% (H2, by mole) H2/CH4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashbacks were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower flux ratio which bends the jet sharply, closer toward the injection plane. The flashback was found to occur along the boundary layer, assisted by the inner wall of the pre-mixer tube. The mean OH-PLIF images characterized the flame-holding behavior where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port.

flame holding↗

Pit Stability Predictions of Additively Manufactured SS316 Surfaces Using Finite Element Analysis

Stainless steels are susceptible to localized forms of corrosion attack, such as pitting. The size and lifetime of a nucleated pit can vary, depending on a critical potential or current density criterion, which determines if the pit repassivates or continues growing. This work uses finite element method (FEM) modeling to compare the critical pit radii predicted by thermodynamic and kinetic repassivation criteria. Experimental electrochemical boundary conditions are used to capture the active pit kinetics. Geometric and environmental parameters, such as the pit shape and size (analogous to additively manufactured lack-of-fusion pores), solution concentration, and water layer thickness were considered to assess their impact on the pit repassivation criterion. The critical pit radius (the transition point from stable growth to repassivation) predicted for a hemispherical pit was larger when using the repassivation potential (E rp ) criteria, as opposed to the current density criteria (pit stability product). Including both the pit stability product and E rp into its calculations, the analytical maximum pit model predicted a critical radius two times more conservative than the FEA approach, under the conditions studied herein. The complex pits representing lack-of-fusion pores were shown to have minimal impact on the critical radius in atmospheric conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Absorbing boundary conditions in material point method adopting perfectly matched layer theory

This study focuses on solving the numerical challenges of imposing absorbing boundary conditions for dynamic simulations in the material point method (MPM). To attenuate elastic waves leaving the computational domain, the current work integrates the Perfectly Matched Layer (PML) theory into the implicit MPM framework. The proposed approach introduces absorbing particles surrounding the computational domain that efficiently absorb outgoing waves and reduce reflections, allowing for accurate modeling of wave propagation and its further impact on geotechnical slope stability analysis. The study also includes several benchmark tests to validate the effectiveness of the proposed method, such as several types of impulse loading and symmetric and asymmetric base shaking. The conducted numerical tests also demonstrate the ability to handle large deformation problems, including the failure of elasto-plastic soils under gravity and dynamic excitations. The findings extend the capability of MPM in simulating continuous analysis of earthquake-induced landslides, from shaking to failure.

58 GEOSCIENCES↗

Validation of NSFsim as a Grad-Shafranov equilibrium solver at DIII-D

Plasma shape is a significant factor that must be considered for any Fusion Pilot Plant (FPP) as it has significant consequences for plasma stability and core confinement. A new simulator, NSFsim, has been developed based on a historically successful code, DINA [1], offering tools to simulate both transport and plasma shape. Specifically, NSFsim is a free boundary equilibrium and transport solver and has been configured to match the properties of the DIII-D tokamak. This paper is focused on validating the Grad-Shafranov (GS) solver of NSFsim by analyzing its ability to recreate the plasma shape, the poloidal flux distribution, and the measurements of the simulated diagnostic signals originating from flux loops and magnetic probes in DIII-D. Five different plasma shapes are simulated to show the robustness of NSFsim to different plasma conditions; these shapes are Lower Single Null (LSN), Upper Single Null (USN), Double Null (DN), Inner Wall Limited (IWL), and Negative Triangularity (NT). The NSFsim results are compared against real measured signals, magnetic profile fits from EFIT [2], and another plasma equilibrium simulator, GSevolve [3]. EFIT reconstructions of shots are readily available at DIII-D, but GSevolve was manually ran by us to provide simulation data to compare against.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase Switching as the Origin of Large Piezoelectric Response in Organic-Inorganic Perovskites: A First-Principles Study

Piezoelectrics are critical functional components of many practical applications such as sensors, ultrasonic transducers, actuators, medical imaging, telecommunications. So far, the best performing piezoelectrics are ferroelectric ceramics, many of which are toxic, heavy, hard, and cost-ineffective. Recently, a groundbreaking discovery of extraordinarily large piezoelectric coefficients in the family of organic-inorganic perovskites gave a hope for a cheaper, environmentally friendly, inexpensive, light-weight and flexible alternative. However, the origin of such response in organic-inorganic ferroelectrics whose spontaneous polarization is an order of magnitude smaller than for inorganic counterparts remains unclear. In this study we employ first-principles simulations to predict that the mechanism associated with large piezoelectric constants is of extrinsic origin and associated with switching between the stable phase and previously overlooked energetically competitive metastable phase that can be stabilized by the external stress. Here, the phase switching changes the polarization direction and, therefore, gives origin to large piezoelectric response, similar to PbZr$_{1-x}$Ti$_x$O$_3$ near morphotropic phase boundary. Existence of such metastable phases is likely to manifest as the dynamical molecular disorder above the Curie temperature and, therefore, could be intrinsic to the entire family of organic-inorganic ferroelectrics with such disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electropolishing-induced topographic defects in niobium: insights and implications for superconducting radio frequency applications

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. Here, in this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Hryhorenko, Oleksandr [Thomas Jefferson National A↗

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗

Degradation science: Integrating modeling and experiments to predict localized corrosion processes (Annual Progress Report)

Additively manufactured (AM) eutectic high-entropy alloys (EHEAs), such as nano-lamellar AlCoCrFeNi 2.1 , have excellent strength, ductility, and wear resistance even at elevated temperatures, but their corrosion behavior in aggressive acids at different length scales remain poorly understood. This work investigates the corrosion behavior of laser powder bed–fused (L-PBF) AlCoCrFeNi 2.1 as a function of annealing temperatures, probing degradation mechanisms from nanoscopic to macroscopic length scales. The alloy is dual phase consisting of a ductile FCC L1 2 phase and a high-strength BCC B2 phase. Rapid solidification during L-PBF produces a far-from-equilibrium nano-lamellar structure with nearly homogeneous elemental distribution, which tends to evolve upon annealing toward Cr/Co/Fe-enriched FCC and Al/Ni-enriched B2. Three conditions were studied: as-printed, 600 °C/5 h, and 1000 °C/1 h, over which B2 lamellae coarsen, lamellar spacing increases, and elemental segregation becomes more prominent. Microstructure and chemistry were characterized by scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS), while in-situ electrochemical atomic force microscopy (EC-AFM) was used to link early (<5 h) local dissolution to microstructure after exposure in sulfuric acid. EC-AFM highlights preferential dissolution of the BCC/B2 phase where the surrounding matrix is Cr-depleted and directly quantifies the dissolution rates within individual phases, tracks the transition from early nano-scale attack to partial repassivation, to correlate height differences with current and impedance responses. To monitor longer-term behavior (up to 96 h), ex-situ AFM, SEM, and confocal imaging were combined with conventional bulk electrochemical tests, bridging nanoscale observations to micro/meso-scale damage morphologies. At the meso-scale, the deepest dissolution channels align with the build-direction lamellae and melt-pool boundaries, indicating that printing directionality guides the propagation of these localized corrosion sites. Annealing modifies corrosion by restructuring BCC/FCC phase fractions, lamellar spacing, and Cr/Al segregation, thereby changing the cathode/anode ratio and passive film stability. The results clarify how as-printed and annealed nano-lamellar architectures differ in their susceptibility to selective dissolution; how elemental segregation competes with residual stresses along the build direction. With these insights, future work will use CALPHAD-guided alloy modification to stabilize higher Cr contents in the B2 phase while retaining a dominant FCC+B2/BCC microstructure, with the goal of designing mechanically robust, corrosion-resistant EHEAs for safety-critical applications to leverage the LLNL’s broader national and global security mission.

36 MATERIALS SCIENCE↗

Electropolishing-induced topographic defects in niobium: Insights and Implications for SRF

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. In this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Lechner, Eric [Thomas Jefferson National Accelerat↗

Elemental segregation to lattice defects in the CrMnFeCoNi high-entropy alloy during high temperature exposures

The influence of small plastic pre-strains on the elevated-temperature stability and microstructure of the equiatomic CrMnFeCoNi FCC solid solution is investigated. Particular attention is given to whether any of the alloy elements segregate to individual dislocations. To that end, CrMnFeCoNi samples were first deformed in tension at room temperature to plastic strains of 0.2 and 2.3%, and subsequently annealed at 973 K for 800 hours. The pre-strains activated planar slip of 1/2<110>-type dislocations on {111}-type glide planes. Interactions of this planar slip with special Σ3 grain boundaries formed a large number of dislocation segments with a <110>-type crystallographic orientation suitable for a credible end-on analysis of dislocation cores in HR-STEM. The cores of the 1/2<110> dislocations pushed up against the investigated grain boundaries were found to be close to the compact configuration. Within the sensitivity of the Super-X EDS mapping, no concentration gradient was detected near dislocations that would indicate enrichment at dislocation cores of any of the elemental constituents of the alloy after the pre-deformation and annealing. However, a Cr-rich tetragonal sigma phase nucleated and grew at grain boundary triple junctions during this anneal, processes that were not accelerated by the enhanced dislocation density present after pre-strain. A clear chromium gradient was observed in the Cr-depleted zones near grain boundaries suggesting that Cr transport occurred by relatively slow diffusion from the bulk to the grain boundaries and then by relatively fast diffusion along the grain boundaries to the precipitates. Accompanying the Cr depletion near grain boundaries is a simultaneous Ni and Mn enrichment, which promotes formation of the L1 0 NiMn phase that is observed on the grain boundaries after prolonged annealing.

36 MATERIALS SCIENCE↗

The persistence and conversion of coastal foredune and swale vegetation community distributions 63 years later

Vegetation shifts can directly alter habitat dynamics and indirectly impact habitat stability relative to disturbance response. Barrier island dune habitats exhibit spatiotemporally dynamic topography that is affected by vegetation. However, vegetation distribution data can be rare and vegetation persistence is largely unknown such that species turnover over in communities can occur unnoticed. This is true despite concerns and documented cases of woody encroachment related to climate change in these ecogeomorphic habitats where physical stability to resist storm erosion varies with vegetation distribution and density. In 1956 and 1957, the vegetation of Island Beach State Park, NJ, was mapped as a permanent record for subsequent ecological study. In June 2020, we remapped the vegetation of 5.4 of 17 km north to south, seaward of the thicket community boundary. We maintained the same classification system as the historic record, physically mapping vegetation patches via GPS. We quantified changes in thicket, heather, and grass community distribution between the two time periods. Habitat persistence and conversion varied in the 63 years. Heath communities saw 80%–97% habitat loss in conversion to woody thicket. Conversely, thicket community distribution drastically increased, replacing heath where it was previously prevalent. This represents the first known documented instances of woody encroachment for Morella pensylvanica. When they did not expand, thicket communities receded landward or underwent turnover to an invasive species. Woody species of interest for dune stabilization occupied areas of similar habitat characteristics. Dune vegetation distribution persistence from 1957 to 2020 was relatively consistent and stable with the exception of heath habitat conversion. Barrier island stability is directly related to vegetation stability such that understanding where community shifts might occur over time can aid in managing and modeling efforts surrounding these dynamic ecogeomorphic habitats.

54 ENVIRONMENTAL SCIENCES↗