SEARCH · Search NASA
Results for “disequilibrium”
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.
The Influence of Air‐Sea CO 2 Disequilibrium on Carbon Sequestration by the Ocean's Biological Pump
Abstract The ocean's biological carbon pump (BCP) affects the Earth's climate by sequestering CO 2 away from the atmosphere for decades to millennia. One primary control on the amount of carbon sequestered by the biological pump is air‐sea CO 2 disequilibrium, which is controlled by the rate of air‐sea CO 2 exchange and the residence time of CO 2 in surface waters. Here, we use a data‐assimilated model of the soft tissue BCP to quantify carbon sequestration inventories and time scales from remineralization in the water column to equilibration with the atmosphere. We find that air‐sea CO 2 disequilibrium enhances the global biogenic carbon inventory by ∼35% and its sequestration time by ∼70 years compared to identical calculations made assuming instantaneous air‐sea CO 2 exchange. Locally, the greatest enhancement occurs in the subpolar Southern Ocean, where air‐sea disequilibrium increases sequestration times by up to 600 years and the biogenic dissolved inorganic carbon inventory by >100% in the upper ocean. Contrastingly, in deep‐water formation regions of the North Atlantic and Antarctic margins, where biological production creates undersaturated surface waters which are subducted before fully equilibrating with the atmosphere, air‐sea CO 2 disequilibrium decreases the depth‐integrated sequestration inventory by up to ∼150%. The global enhancement of carbon sequestration by air‐sea disequilibrium is particularly important for carbon respired in deep waters that upwell in the Southern Ocean. These results highlight the importance of accounting for air‐sea CO 2 disequilibrium when evaluating carbon sequestration by the biological pump and for assessing the efficacy of ocean‐based CO 2 removal methods.
Linking Community‐Climate Disequilibrium to Ecosystem Function
Turnover in species composition often lags behind the pace of climate change, resulting in mismatches between climate and communities. However, the impact of these community‐climate disequilibria on ecosystem functions is rarely considered, and current methods for measuring disequilibria assume that species ranges were, until recently, in equilibrium with climate. Here, in this work, we develop a simple theoretical model to address both of these problems by linking community‐climate disequilibrium with ecosystem functioning. We show how disequilibrium can impair functioning in the near‐term even when climate change is expected to enhance functioning in the long‐term. Responses are most likely to change over time in communities where turnover is slow, the impact of disequilibrium counteracts the direct effects of climate on ecosystem function, and pre‐existing disequilibrium is large. These findings emphasise the importance of precise and unbiased estimates of community‐climate disequilibria for improving ecological forecasts. By fitting our model to time series of both climate and ecosystem function from a metacommunity simulation, we show the potential for community‐climate disequilibrium to be inferred without direct knowledge about species' distributions or climatic tolerances. We end by outlining a research agenda to apply dynamic disequilibrium concepts and test novel hypotheses across diverse ecosystems.
Disequilibrium Chemistry in Exoplanet Atmospheres Observed with the Hubble Space Telescope
Literature on the theory of exoplanet atmospheric disequilibrium chemistry is rich, although its observational counterpart has yet to emerge beyond the hints provided by a few targets in dedicated studies. We report results from an uniform data reduction and analysis for a catalog of 62 Hubble Space Telescope exoplanet transit spectra where we assess the atmospheric model preference for disequilibrium chemistry (i.e., water vapor is not the dominant absorption spectral signature) over thermal equilibrium chemistry in a comparative planetology context. Where model preference assessment is possible, we find that disequilibrium occurs in about half of the atmospheres, indicating that disequilibrium processes play an important role in the composition of exoplanet atmospheres. While very hot atmospheres, over 1800 K, prefer equilibrium chemistry, we find a clustering of preference for disequilibrium in the 1200–1800 K temperature range. We suggest that UV-augmented thermochemistry may play a significant role for those atmospheres.
The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era
Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope (JWST) era, high-resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constraints on the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models. In this paper, we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide-separation giant planets with key molecular species such as CH 4 , H 2 O, CO, and NH 3 in disequilibrium. Our grid covers atmospheres with T eff ∈ [500 K, 1300 K], log g ∈ [3.0, 5.5] (cgs) and an eddy diffusion parameter of logK zz = 2,4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH 4 , H 2 O, CO, and NH 3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.
Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs
Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects to be discovered by future direct imaging surveys. We present a low-resolution, R ∼ 300, M-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K {sub zz}) to estimate the strength of vertical mixing in each object. The K {sub zz} values of cooler gaseous objects are close to their theoretical maximum, and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K {sub zz} values imply that disequilibrium phosphine (PH{sub 3}) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH{sub 3} absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's M-band spectrum by including both CO and water clouds in the atmospheric model.
Ionization disequilibrium in K- and L-shell ions
Time-gated Sc K-shell and Ge L-shell spectra are presented from a range of characterized thermodynamic states spanning ion densities of 1019–1020cm−3 and plasma temperatures around 2000 eV. For the higher densities studied and temperatures from 1000 to 3000 eV, the Sc and Ge x-ray emission spectra are consistent with steady-state calculations from the modern atomic kinetics model SCRAM. At the lower ion densities achieved through plasma expansion, however, the model calculations require a higher plasma temperature to reproduce the observed Ge spectrum. We attribute this to ionization disequilibrium of the Sc because the ionization time scales exceed the hydrodynamic timescale when the inferred temperatures diverge.
Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes
Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.
Global Ocean Cooling of 2.3°C During the Last Glacial Maximum
Abstract Quantitative constraints on past mean ocean temperature (MOT) critically inform our historical understanding of Earth's energy balance. A recently developed MOT proxy based on paleoatmospheric Xe, Kr, and N 2 ratios in ice core air bubbles is a promising tool rooted in the temperature dependences of gas solubilities. However, these inert gases are systematically undersaturated in the modern ocean interior, and it remains unclear how air‐sea disequilibrium may have changed in the past. Here, we carry out 30 tracer‐enabled model simulations under varying circulation, sea ice cover, and wind stress regimes to evaluate air‐sea disequilibrium in the Last Glacial Maximum (LGM) ocean. We find that undersaturation of all three gases was likely reduced, primarily due to strengthened high‐latitude winds, biasing reconstructed MOT by −0.38 ± 0.37°C (1σ). Accounting for air‐sea disequilibrium, paleoatmospheric inert gases indicate that LGM MOT was 2.27 ± 0.46°C (1σ) colder than the pre‐industrial era.
Dynamic Surface Incorporation of Pb 2+ Ions at the Actively Dissolving Calcite (104) Surface
The reaction of dissolved Pb 2+ with calcite surfaces at near-equilibrium conditions involves adsorption of Pb 2+ and precipitation of secondary heteroepitaxial Pb-carbonate minerals. A more complex behavior is observed under far-from-equilibrium conditions, including strong inhibition of calcite dissolution, development of microtopography, and near-surface incorporation of multiple monolayers (ML) of Pb 2+ without precipitation of secondary phases [where 1 ML ≡ 1 Ca/20.2 Å 2 , the crystallographic site density of the calcite (104) lattice plane]. However, the mechanistic controls governing far-from-equilibrium reactivity are not well understood. Here, in this study, we observe the interfacial incorporation of dissolved Pb 2+ during the dissolution of calcite (104) surfaces at pH ~3.7 in a flow-through reaction cell, revealing the formation of a ~1 nm thick Pb-rich calcite layer with a total Pb coverage of ~1.4 ML. These observations of the sorbed Pb distribution used resonant anomalous X-ray reflectivity, X-ray fluorescence, and nanoinfrared atomic force microscopy. We propose that this altered surface layer represents a novel sorption mode that is stabilized by conditions of sustained disequilibrium. This behavior may significantly impact the transport of dissolved metals during disequilibrium processes occurring in acid mine drainage and subsurface CO 2 injection and, if appropriately accounted for, could improve the predictive capability of geochemical reactive-transport models.
Structure-informed clustering for population stratification in association studies
Background: Identifying variants associated with complex traits is a challenging task in genetic association studies due to linkage disequilibrium (LD) between genetic variants and population stratification, unrelated to the disease risk. Existing methods of population structure correction use principal component analysis or linear mixed models with a random effect when modeling associations between a trait of interest and genetic markers. However, due to stringent significance thresholds and latent interactions between the markers, these methods often fail to detect genuinely associated variants. Results: To overcome this, we propose CluStrat, which corrects for complex arbitrarily structured populations while leveraging the linkage disequilibrium induced distances between genetic markers. It performs an agglomerative hierarchical clustering using the Mahalanobis distance covariance matrix of the markers. In simulation studies, we show that our method outperforms existing methods in detecting true causal variants. Applying CluStrat on WTCCC2 and UK Biobank cohorts, we found biologically relevant associations in Schizophrenia and Myocardial Infarction. CluStrat was also able to correct for population structure in polygenic adaptation of height in Europeans. Conclusions: CluStrat highlights the advantages of biologically relevant distance metrics, such as the Mahalanobis distance, which captures the cryptic interactions within populations in the presence of LD better than the Euclidean distance.
pyFLANK, a graph neural network based null distribution inference model for F ST outlier detection
Detecting genomic regions under selection is essential for understanding how populations adapt to different environments, yet it remains challenging due to the confounding effects of demographic history and linkage disequilibrium (LD). Fixation index (F ST ) is a widely used statistic to identify genomic regions under adaptation. However, identifying genes under selection by defining F ST outliers often remains challenging, owing to confounding effects of underlying demographic history. Traditional methods assume independence among loci and rely on simple demographic models, while newer models perform much better but are computationally expensive and not easily scalable. Here, we present pyFLANK, an open-source and automated Python implementation which detects F ST outliers using a null distribution inferred from quasi-independent loci. Our tool integrates three approaches to identify loci obeying a null distribution: graph neural network (GNN) inference, linkage disequilibrium (LD)-based inference, and user-defined input. Because pyFLANK uses GNN-based inference of quasi-independent loci, it yields a more accurate null model with less need for user parameter input. In simulation experiments, pyFLANK achieved lower false positive rates than current methods while maintaining comparable detection power, indicating that its refined null model better distinguishes true adaptive loci from background variation. The GNN-based model, in particular, detected additional loci associated with phenotypic variance that were not identified by existing methods. Assessments of simulation and real data from different species demonstrate that pyFLANK achieves lower false positive rates compared with other commonly used F ST outlier detectors, while maintaining comparable detection power and excellent computational performance, providing a robust and user-friendly tool for identifying loci under divergent selection. It extends existing F ST outlier frameworks by incorporating explicit LD-aware strategies for null model calibration. The method is intended as a practical and scalable complement to existing genome scan approaches.
PyPop: a mature open-source software pipeline for population genomics
Python for Population Genomics (PyPop) is a software package that processes genotype and allele data and performs large-scale population genetic analyses on highly polymorphic multi-locus genotype data. In particular, PyPop tests data conformity to Hardy-Weinberg equilibrium expectations, performs Ewens-Watterson tests for selection, estimates haplotype frequencies, measures linkage disequilibrium, and tests significance. Standardized means of performing these tests is key for contemporary studies of evolutionary biology and population genetics, and these tests are central to genetic studies of disease association as well. Here, we present PyPop 1.0.0, a new major release of the package, which implements new features using the more robust infrastructure of GitHub, and is distributed via the industry-standard Python Package Index. New features include implementation of the asymmetric linkage disequilibrium measures and, of particular interest to the immunogenetics research communities, support for modern nomenclature, including colon-delimited allele names, and improvements to meta-analysis features for aggregating outputs for multiple populations.
Simplification of the Grid Model and its Impact on the Analysis of Electrical Power Systems
Here, this paper analyzes the impact of the use of Kron reduction on the state variables of a three-phase electrical system, even when it does not meet the necessary conditions for its application. Reduction is applied to a power line model to eliminate the equation corresponding to the neutral conductor of the line. The ATP program is used to model and simulate the behavior of an electrical system considering different degrees of disequilibrium as a reference for the comparison of results. The results show that under certain conditions of disequilibrium the Kron reduction can lead to significant errors in the state variables of the system.
Apatite geochemistry as a tool for understanding the petrogenesis of layered mafic-ultramafic rocks in the Bushveld Complex, South Africa
The sources of the magmas that formed the Rustenburg Layered Suite of the Bushveld Complex in South Africa remain debated, despite decades of research. Vertical and lateral variation in bulk rock and mineral separate Sr-Nd isotopic compositions, which generally indicate enriched sources, demonstrate that the layered sequence was formed by the emplacement of multiple batches of magma, crucially resulting in episodes of PGE-Cr-V mineralisation. The Lu-Hf isotope compositions of zircon are, however, at odds with the bulk rock Sr-Nd isotopic heterogeneity as they show near homogeneous compositions throughout the layered sequence (εHf (2.06 Ga) =−8). This lack of variation in Hf isotope composition has been attributed to deep, continental lithospheric mantle-related and/or crustal contamination of plume-derived Bushveld magmas. In this study, we analysed the major, trace element and Sr-Nd isotope geochemistry of apatite in the Rustenburg Layered Suite. Apatite occurs as an intercumulus mineral in the lowermost regions and a cumulus mineral in the uppermost regions of the layered sequence and can therefore be used to test existing models for the isotopic disequilibrium between bulk rock Sr-Nd and zircon Hf isotopic compositions. Apatite is largely chlorapatite in the lowermost regions and fluorapatite in the uppermost regions of the layered sequence. The Merensky Reef is unusual in that it comprises both chlorapatite and fluorapatite. Apatite throughout the layered sequence is generally unzoned and shows no evidence of late-stage alteration. Trace element data show that apatite is enriched in L/HREE, with common negative Eu-Sr anomalies. These trace element signatures are consistent with a magmatic origin for the apatite grains, with prior, or concurrent, plagioclase crystallization from the same melt. Variability in in situ Sr and Nd isotope compositions of apatite is recorded throughout the layered sequence with εNd (2.06 Ga) compositions varying between −2.5 and − 10.2 and initial 87 Sr/ 86 Sr compositions varying between 0.7079 and 0.7103 (for the Marikana dikes only). The variability in Sr-Nd isotope compositions of apatite is consistent with the bulk rock (and mineral separate) variation in Sr-Nd isotope compositions, suggesting apatite preserves primary magmatic compositions in the Rustenburg Layered Suite.
Ecological acclimation: A framework to integrate fast and slow responses to climate change
Ecological responses to climate change occur across vastly different time-scales, from minutes for physiological plasticity to decades or centuries for community turnover and evolutionary adaptation. Accurately predicting the range of ecosystem trajectories will require models that incorporate both fast processes that may keep pace with climate change and slower ones likely to lag behind and generate disequilibrium dynamics. However, the knowledge necessary for this integration is currently fragmented across disciplines. We develop ‘ecological acclimation’ as a unifying framework to emphasize the similarity of dynamics driven by processes operating on dramatically different time-scales and levels of biological organization. The framework focuses on ecoclimate sensitivities, measured as the change in an ecological response variable per unit of climate change. Acclimation processes acting at different time-scales cause these sensitivities to shift in magnitude and even direction over time. We highlight shifting ecoclimate sensitivities in case studies from diverse ecosystems, including terrestrial plant communities, coral reefs and soil microbiomes. Models predicting future ecosystem states inevitably make assumptions about acclimation processes; these assumptions must be explicit for users to evaluate whether a model is appropriate for a given forecast horizon. Similarly, decision frameworks that clearly account for multiple acclimation processes and their distinct time-scales will help natural resource managers plan for ecological impacts of climate change from years to many decades into the future. We outline a synthetic research programme focused on the time-scales of ecological acclimation to reduce uncertainty in ecological forecasts.
A preference for cold dark matter over Superfluid Dark Matter in local Milky Way data
There are many well-known correlations between dark matter and baryons that exist on galactic scales. These correlations can essentially be encompassed by a simple scaling relation between observed and baryonic accelerations, historically known as the Mass Discrepancy Acceleration Relation (MDAR). The existence of such a relation has prompted many theories that attempt to explain the correlations by invoking additional fundamental forces on baryons. The standard lore has been that a theory that reduces to the MDAR on galaxy scales but behaves like cold dark matter (CDM) on larger scales provides an excellent fit to data, since CDM is desirable on scales of clusters and above. However, this statement should be revised in light of recent results showing that a fundamental force that reproduces the MDAR is challenged by local Milky Way dynamics and rotation curve data between 5–18 kpc. In this study, we test this claim on the example of Superfluid Dark Matter. We find that a standard CDM model is preferred over a static superfluid profile assuming a steady-state Galactic disk and discuss the robustness of this conclusion to disequilibrium effects. This preference is due to the fact that the superfluid model over-predicts vertical accelerations, even while reproducing galactic rotation curves. Our results establish an important criterion that any dark matter model must satisfy within the Milky Way.
Formation of zinc carbonate phases on dissolving calcite, aragonite, and vaterite in acidic aqueous solutions
Calcium carbonate (CaCO 3 ) minerals serve as a major sink to retain metal ions through mineral-water interfacial reactions in which the capacity and long-term stability of contaminant uptake are influenced by coupled dissolution and precipitation reactions of both primary and secondary carbonate minerals (i.e., mineral replacement). Notably, recent studies of calcite reactivity in acidic solutions containing high levels of metal ions demonstrated complex behavior under conditions of sustained disequilibrium. Here, we explored the reactivity of three CaCO 3 polymorphs (calcite, aragonite, and vaterite) with acidic Zn 2+ -containing aqueous solutions using a suite of imaging techniques including optical and scanning electron microscopies, synchrotron-based micro X-ray fluorescence, and transmission X-ray microscopy. Zn uptake by calcite occurred through a two-step process: the formation of a thin layer of the zinc precipitate on the substrate surface followed by the growth of fibrous and radiating hydrozincite particles from the layer. In contrast, Zn uptake by aragonite occurred by mineral replacement where the secondary Zn carbonate phase preserved the external morphology of the original crystal (i.e., a pseudomorph). The replacement of vaterite by hydrozincite occurred within the confined space beneath the porous shell of vaterite, signifying that the primary mechanism driving Zn carbonate precipitation was chemical exchange through the pores. When multiple CaCO 3 polymorphs coexisted, the replacement of aragonite and vaterite occurred preferentially over that of calcite. These results demonstrate the distinct morphological and mineralogical controls over the reactivities of calcium carbonate minerals with Zn 2+ under acidic conditions.