Search NASA⌕ Search

SEARCH · Search NASA

Results for “mixed effect model”

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 19 records

HighDimMixedModels.jl: Robust high-dimensional mixed-effects models across omics data

High-dimensional mixed-effects models are an increasingly important form of regression in which the number of covariates rivals or exceeds the number of samples, which are collected in groups or clusters. The penalized likelihood approach to fitting these models relies on a coordinate descent algorithm that lacks guarantees of convergence to a global optimum. Here, we empirically study the behavior of this algorithm on simulated and real examples of three types of data that are common in modern biology: transcriptome, genome-wide association, and microbiome data. Our simulations provide new insights into the algorithm’s behavior in these settings, and, comparing the performance of two popular penalties, we demonstrate that the smoothly clipped absolute deviation (SCAD) penalty consistently outperforms the least absolute shrinkage and selection operator (LASSO) penalty in terms of both variable selection and estimation accuracy across omics data. To empower researchers in biology and other fields to fit models with the SCAD penalty, we implement the algorithm in a Julia package, HighDimMixedModels.jl .

Gorstein, Evan↗

Evaluation of partial volume correction and analysis of longitudinal [ 18 F]GTP1 tau PET imaging in Alzheimer's disease using linear mixed-effects models

Purpose: We evaluated the impact of partial volume correction (PVC) methods on the quantification of longitudinal [ 18 F]GTP1 tau positron-emission tomography (PET) in Alzheimer's disease and the suitability of describing the tau pathology burden temporal trajectories using linear mixed-effects models (LMEM). Methods: We applied van Cittert iterative deconvolution (VC), 2-compartment, and 3-compartment, and the geometric transfer matrix plus region-based voxelwise methods to data acquired in an Alzheimer's disease natural history study over 18 months at a single imaging site. We determined the optimal PVC method by comparing the standardized uptake value ratio change (%ΔSUVR) between diagnostic and tau burden–level groups and the longitudinal repeatability derived from the LMEM. The performance of LMEM analysis for calculating %ΔSUVR was evaluated in a natural history study and in a multisite clinical trial of semorinemab in prodromal to mild Alzheimer's disease by comparing results to traditional per-visit estimates. Results: The VC, 2-compartment, and 3-compartment PVC methods had similar performance, whereas region-based voxelwise overcorrected regions with a higher tau burden. The lowest within-subject variability and acceptable group separation scores were observed without PVC. The LMEM-derived %ΔSUVR values were similar to the per-visit estimates with lower variability. Conclusion: The results indicate that the tested PVC methods do not offer a clear advantage or improvement over non-PVC images for the quantification of longitudinal [ 18 F]GTP1 PET data. LMEM offers a robust framework for the longitudinal tau PET quantification with low longitudinal test–retest variability. Clinical trial registration: NCT02640092 and NCT03289143.

60 APPLIED LIFE SCIENCES↗

Estimating basis functions in massive fields under the spatial mixed effects model

Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.

Pazdernik, Karl↗

Mixed Effect Modeling of Dose and Linear Energy Transfer Correlations With Brain Image Changes After Intensity Modulated Proton Therapy for Skull Base Head and Neck Cancer

Intensity modulated proton therapy (IMPT) could yield high linear energy transfer (LET) in critical structures and increased biological effect. For head and neck cancers at the skull base this could potentially result in radiation-associated brain image change (RAIC). The purpose of the current study was to investigate voxel-wise dose and LET correlations with RAIC after IMPT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Evaluation of Plasma Phosphorylated Tau217 for Differentiation Between Alzheimer Disease and Frontotemporal Lobar Degeneration Subtypes Among Patients With Corticobasal Syndrome

Plasma phosphorylated tau217 (p-tau217), a biomarker of Alzheimer disease (AD), is of special interest in corticobasal syndrome (CBS) because autopsy studies have revealed AD is the driving neuropathology in up to 40% of cases. This differentiates CBS from other 4-repeat tauopathy (4RT)–associated syndromes, such as progressive supranuclear palsy Richardson syndrome (PSP-RS) and nonfluent primary progressive aphasia (nfvPPA), where underlying frontotemporal lobar degeneration (FTLD) is typically the primary neuropathology. To validate plasma p-tau217 against positron emission tomography (PET) in 4RT-associated syndromes, especially CBS. This multicohort study with 6, 12, and 24-month follow-up recruited adult participants between January 2011 and September 2020 from 8 tertiary care centers in the 4RT Neuroimaging Initiative (4RTNI). All participants with CBS (n = 113), PSP-RS (n = 121), and nfvPPA (n = 39) were included; other diagnoses were excluded due to rarity (n = 29). Individuals with PET-confirmed AD (n = 54) and PET-negative cognitively normal control individuals (n = 59) were evaluated at University of California San Francisco. Operators were blinded to the cohort. Plasma p-tau217, measured by Meso Scale Discovery electrochemiluminescence, was validated against amyloid-β (Aβ) and flortaucipir (FTP) PET. Imaging analyses used voxel-based morphometry and bayesian linear mixed-effects modeling. Clinical biomarker associations were evaluated using longitudinal mixed-effect modeling. Of 386 participants, 199 (52%) were female, and the mean (SD) age was 68 (8) years. Plasma p-tau217 was elevated in patients with CBS with positive Aβ PET results (mean [SD], 0.57 [0.43] pg/mL) or FTP PET (mean [SD], 0.75 [0.30] pg/mL) to concentrations comparable to control individuals with AD (mean [SD], 0.72 [0.37]), whereas PSP-RS and nfvPPA showed no increase relative to control. Within CBS, p-tau217 had excellent diagnostic performance with area under the receiver operating characteristic curve (AUC) for Aβ PET of 0.87 (95% CI, 0.76-0.98; P < .001) and FTP PET of 0.93 (95% CI, 0.83-1.00; P < .001). At baseline, individuals with CBS-AD (n = 12), defined by a PET-validated plasma p-tau217 cutoff 0.25 pg/mL or greater, had increased temporoparietal atrophy at baseline compared to individuals with CBS-FTLD (n = 39), whereas longitudinally, individuals with CBS-FTLD had faster brainstem atrophy rates. Individuals with CBS-FTLD also progressed more rapidly on a modified version of the PSP Rating Scale than those with CBS-AD (mean [SD], 3.5 [0.5] vs 0.8 [0.8] points/year; P = .005). In this cohort study, plasma p-tau217 had excellent diagnostic performance for identifying Aβ or FTP PET positivity within CBS with likely underlying AD pathology. Plasma P-tau217 may be a useful and inexpensive biomarker to select patients for CBS clinical trials.

59 BASIC BIOLOGICAL SCIENCES↗

Single reader between-cases AUC estimator with nested data

The area under the receiver operating characteristic curve (AUC) is widely used in evaluating diagnostic performance for many clinical tasks. It is still challenging to evaluate the reading performance of distinguishing between positive and negative regions of interest (ROIs) in the nested-data problem, where multiple ROIs are nested within the cases. To address this issue, we identify two kinds of AUC estimators, within-cases AUC and between-cases AUC. We focus on the between-cases AUC estimator, since our main research interest is in patient-level diagnostic performance rather than location-level performance (the ability to separate ROIs with and without disease within each patient). Another reason is that as the case number increases, the number of between-cases paired ROIs is much larger than the number of within-cases ROIs. We provide estimators for the variance of the between-cases AUC and for the covariance when there are two readers. We derive and prove the above estimators’ theoretical values based on a simulation model and characterize their behavior using Monte Carlo simulation results. We also provide a real-data example. Moreover, we connect the distribution-based simulation model with the simulation model based on the linear mixed-effect model, which helps better understand the sources of variation in the simulated dataset.

Health Care Sciences & Services↗

A Robust Segmented Mixed Effect Regression Model for Baseline Electricity Consumption Forecasting

Renewable energy production has been surging around the world in recent years. To mitigate the increasing uncertainty and intermittency of the renewable generation, proactive demand response algorithms and programs are proposed and developed to further improve the utilization of load flexibility and increase the efficiency of power system operation. One of the biggest challenges to efficient control and operation of demand response resources is how to forecast the baseline electricity consumption and estimate the load impact from demand response resources accurately. In this paper, we propose a mixed effect segmented regression model and a new robust estimate for forecasting the baseline electricity consumption in Southern California, USA, by combining the ideas of random effect regression model, segmented regression model, and the least trimmed squares estimate. Since the log-likelihood of the considered model is not differentiable at breakpoints, we propose a new backfitting algorithm to estimate the unknown parameters. The estimation performance of the new estimation procedure has been demonstrated with both simulation studies and the real data application for the electric load baseline forecasting in Southern California.

42 ENGINEERING↗

The temporal onset of associations of cortical proteins with cognitive resilience vary during late life

Background: Cortical proteins associated with cognitive resilience have been identified but their temporal onset in older adults is unknown. We present a multistage approach to first identify cortical proteins associated with cognitive resilience and then examine their associated temporal onset. Methods: We used data from a subset of 1088 decedents from two cohort-studies who had selected reaction monitoring proteomics from the dorsolateral prefrontal cortex, and at least 3 cognitive assessments. Cognition was assessed using a composite derived from 19 tests. We first used linear mixed-effects models to identify cortical proteins associated with cognitive resilience. We then used functional mixed-effects models to examine non-linear associations between proteins and cognitive resilience to identify their temporal onset. Results: Mean age at death was 90 years (SD = 6.4); 69 % were female. On average, cognition started to decline at around 15 years before death, with accelerated decline in the last 7 years. We identified 40 proteins associated with cognitive resilience, of which 17 proteins also showed non-linear associations. Non-linear associations indicated that higher levels of 10 proteins were associated with slower cognitive decline between 23 and 4 years before death. In contrast, higher levels of 7 proteins were associated with faster decline only within the last 7 years before death. Conclusions: Cognitive resilience proteins are differentially related to late-life cognitive aging; the onset of proteins that maintain cognition may begin many years before the onset of proteins that hasten cognitive decline. The temporal onset of cognitive resilience proteins may be crucial for timing efficacious interventions.

Zammit, Andrea↗

Do Remote Camera Arrangements and Image Capture Settings Improve Individual Identification of Golden Eagles?

Individual identification of animals from camera traps has become an important task in wildlife research, but camera deployment methods often do not facilitate this important undertaking. Identification of individual golden eagles (Aquila chrysaetos) is possible using uniquely marked rectrices, but no studies have explored methods to maximize the rate of individual identification from camera images. Furthermore, our objectives were to assess whether different camera heights (1 m vs. 3 m), image capture settings (one image after a 1-min delay vs. burst of 5 images after a 30 sec delay), and arrangements relative to bait (dorsally vs. ventrally aimed) affected views of rectrices on golden eagles and our ability to identify individuals. We conducted our study from 15 December 2016 to 3 March 2017 on the Savannah River Site, South Carolina. First, we developed a scoring system based on views of rectrices and used a linear mixed-effects model to compare image scores among different camera arrangements and image settings. Next, after identifying individual eagles, we used generalized linear mixed-effects models to compare total individual eagle detections, total days an individual was detected, and probability of obtaining an unknown individual identification among camera arrangements and settings. Overall, we scored a total of 27,499 images, with 8,083 providing views of marked rectrices that allowed identification of 18 individual eagles. Average image scores and proportion of images suitable for individual identification were higher from elevated (3 m) camera arrangements than standard arrangements (1 m) across sites. Regardless of camera height, faster frequency of image capture provided more images that could be used to identify individuals and the most trap days per individual. Researchers and managers should consider deploying elevated cameras traps with faster frequency of image capture to improve data quality and potential for analysis of golden eagle populations and trends across the species’ range.

60 APPLIED LIFE SCIENCES↗

Comparative Analyses of Bioequivalence Assessment Methods for In Vitro Permeation Test Data

ABSTRACT For topical, dermatological drug products, an in vitro option to determine bioequivalence (BE) between test and reference products is recommended. In particular, in vitro permeation test (IVPT) data analysis uses a reference‐scaled approach for two primary endpoints, cumulative penetration amount (AMT) and maximum flux ( J max ), which takes the within donor variability into consideration. In 2022, the Food and Drug Administration (FDA) published a draft IVPT guidance that includes statistical analysis methods for both balanced and unbalanced cases of IVPT study data. This work presents a comprehensive evaluation of various methodologies used to estimate critical parameters essential in assessing BE. Specifically, we investigate the performance of the FDA draft IVPT guidance approach alongside alternative empirical and model‐based methods utilizing mixed‐effects models. Our analyses include both simulated scenarios and real‐world studies. In simulated scenarios, empirical formulas consistently demonstrate robustness in approximating the true model, particularly in effectively addressing treatment–donor interactions. Conversely, the effectiveness of model‐based approaches heavily relies on precise model selection, which significantly influences their results. The research emphasizes the importance of accurate model selection in model‐based BE assessment methodologies. It sheds light on the advantages of empirical formulas, highlighting their reliability compared to model‐based approaches and offers valuable implications for BE assessments. Our findings underscore the significance of robust methodologies and provide essential insights to advance their understanding and application in the assessment of BE, employed in IVPT data analysis.

Leon, Sami↗

Arm and shoulder muscle segmentation in axial MRI with UNet deep learning model

Quantifying individual upper-limb muscle volumes from MRI provides key insight into muscle-specific strength, deficits, and adaptations. Manual delineation is the gold standard but time‑intensive, and the performance of current deep learning approaches, particularly for small or anatomically complex muscles, remains incompletely characterized. We evaluated a state‑of‑the‑art deep learning framework across the entire upper limb and analyzed factors governing segmentation performance, with attention to the forearm. Three previously published MRI datasets (1.5 T, 3D GRE T1‑weighted; total n = 39) spanning young, middle‑aged, and older adults were curated and quality‑checked, including expert manual segmentations for 31 muscles. Following multiclass mask reconstruction, we trained three 3D nnU‑Net multiclass models matched to the muscle subsets present across datasets, using five‑fold cross‑validation and a composite Dice Similarity Coefficient (DSC) + cross entropy loss. Segmentation accuracy was assessed with DSC. Performance varied across muscles (mean DSC = 0.806 ± 0.098), ranging from 0.920 (Deltoid) to 0.461 (Extensor pollicis brevis). In uncertainty‑weighted regressions, muscle volume was positively associated with DSC (R2 = 0.36, p < 0.001), whereas training segmentation count and muscle orientation showed negligible associations (R2 ≤ 0.06). A weighted mixed‑effects model identified volume as the strongest evaluated predictor, explaining 23.9% of variance in DSC; orientation and training count each contributed <1%, leaving 61.5% unexplained. These results indicate that deep learning–based segmentation can accurately quantify muscle volume for many upper‑limb muscles but remains constrained for small, low‑contrast forearm muscles.

Gillespie, Samuel↗

Quantitative methods and modeling to assess COVID–19–interrupted in vivo pharmacokinetic bioequivalence studies with two reference batches

The coronavirus disease 2019 (COVID-19) has presented unprecedented challenges to the generic drug development, including interruptions in bioequivalence (BE) studies. Per guidance published by the US Food and Drug Administration (FDA) during the COVID-19 public health emergency, any protocol changes or alternative statistical analysis plan for COVID-19-interrupted BE study should be accompanied with adequate justifications and not lead to biased equivalence determination. In this study, we used a modeling and simulation approach to assess the potential impact of study outcomes when two different batches of a Reference Standard (RS) were to be used in an in vivo pharmacokinetic BE study due to the RS expiration during the COVID-19 pandemic. Simulations were performed with hypothetical drugs under two scenarios: (1) uninterrupted study using a single batch of an RS, and (2) interrupted study using two batches of an RS. The acceptability of BE outcomes was evaluated by comparing the results obtained from interrupted studies with those from uninterrupted studies. The simulation results demonstrated that using a conventional statistical approach to evaluate BE for COVID-19-interrupted studies may be acceptable based on the pooled data from two batches. An alternative statistical method which includes a “batch” effect to the mixed effects model may be used when a significant “batch” effect was found in interrupted four-way crossover studies. However, such alternative method is not applicable for interrupted two-way crossover studies. Overall, the simulated scenarios are only for demonstration purpose, the acceptability of BE outcomes for the COVID19-interrupted studies could be case-specific.

60 APPLIED LIFE SCIENCES↗

Environmental controls on the light use efficiency of terrestrial gross primary production

Abstract Gross primary production (GPP) by terrestrial ecosystems is a key quantity in the global carbon cycle. The instantaneous controls of leaf‐level photosynthesis are well established, but there is still no consensus on the mechanisms by which canopy‐level GPP depends on spatial and temporal variation in the environment. The standard model of photosynthesis provides a robust mechanistic representation for C 3 species; however, additional assumptions are required to “scale up” from leaf to canopy. As a consequence, competing models make inconsistent predictions about how GPP will respond to continuing environmental change. This problem is addressed here by means of an empirical analysis of the light use efficiency (LUE) of GPP inferred from eddy covariance carbon dioxide flux measurements, in situ measurements of photosynthetically active radiation (PAR), and remotely sensed estimates of the fraction of PAR (fAPAR) absorbed by the vegetation canopy. Focusing on LUE allows potential drivers of GPP to be separated from its overriding dependence on light. GPP data from over 100 sites, collated over 20 years and located in a range of biomes and climate zones, were extracted from the FLUXNET2015 database and combined with remotely sensed fAPAR data to estimate daily LUE. Daytime air temperature, vapor pressure deficit, diffuse fraction of solar radiation, and soil moisture were shown to be salient predictors of LUE in a generalized linear mixed‐effects model. The same model design was fitted to site‐based LUE estimates generated by 16 terrestrial ecosystem models. The published models showed wide variation in the shape, the strength, and even the sign of the environmental effects on modeled LUE. These findings highlight important model deficiencies and suggest a need to progress beyond simple “goodness of fit” comparisons of inferred and predicted carbon fluxes toward an approach focused on the functional responses of the underlying dependencies.

54 ENVIRONMENTAL SCIENCES↗

Accuracy of Tau Positron Emission Tomography as a Prognostic Marker in Preclinical and Prodromal Alzheimer Disease: A Head-to-Head Comparison Against Amyloid Positron Emission Tomography and Magnetic Resonance Imaging

Importance: Tau positron emission tomography (PET) tracers have proven useful for the differential diagnosis of dementia, but their utility for predicting cognitive change is unclear. Objective: To examine the prognostic accuracy of baseline fluorine 18 ( 18 F)-flortaucipir and [18F]RO948 (tau) PET in individuals across the Alzheimer disease (AD) clinical spectrum and to perform a head-to-head comparison against established magnetic resonance imaging (MRI) and amyloid PET markers. Design, setting, and participants: This prognostic study collected data from 8 cohorts in South Korea, Sweden, and the US from June 1, 2014, to February 28, 2021, with a mean (SD) follow-up of 1.9 (0.8) years. A total of 1431 participants were recruited from memory clinics, clinical trials, or cohort studies; 673 were cognitively unimpaired (CU group; 253 [37.6%] positive for amyloid-β [Aβ]), 443 had mild cognitive impairment (MCI group; 271 [61.2%] positive for Aβ), and 315 had a clinical diagnosis of AD dementia (315 [100%] positive for Aβ). Exposures: [ 18 F]Flortaucipir PET in the discovery cohort (n = 1135) or [ 18 F]RO948 PET in the replication cohort (n = 296), T1-weighted MRI (n = 1431), and amyloid PET (n = 1329) at baseline and repeated Mini-Mental State Examination (MMSE) evaluation. Main outcomes and measures: Baseline [ 18 F]flortaucipir/[ 18 F]RO948 PET retention within a temporal region of interest, MRI-based AD-signature cortical thickness, and amyloid PET Centiloids were used to predict changes in MMSE using linear mixed-effects models adjusted for age, sex, education, and cohort. Mediation/interaction analyses tested whether associations between baseline tau PET and cognitive change were mediated by baseline MRI measures and whether age, sex, and APOE genotype modified these associations. Results: Among 1431 participants, the mean (SD) age was 71.2 (8.8) years; 751 (52.5%) were male. Findings for [ 18 F]flortaucipir PET predicted longitudinal changes in MMSE, and effect sizes were stronger than for AD-signature cortical thickness and amyloid PET across all participants (R 2 , 0.35 [tau PET] vs 0.24 [MRI] vs 0.17 [amyloid PET]; P < .001, bootstrapped for difference) in the Aβ-positive MCI group (R2, 0.25 [tau PET] vs 0.15 [MRI] vs 0.07 [amyloid PET]; P < .001, bootstrapped for difference) and in the Aβ-positive CU group (R 2 , 0.16 [tau PET] vs 0.08 [MRI] vs 0.08 [amyloid PET]; P < .001, bootstrapped for difference). These findings were replicated in the [ 18 F]RO948 PET cohort. MRI mediated the association between [ 18 F]flortaucipir PET and MMSE in the groups with AD dementia (33.4% [95% CI, 15.5%-60.0%] of the total effect) and Aβ-positive MCI (13.6% [95% CI, 0.0%-28.0%] of the total effect), but not the Aβ-positive CU group (3.7% [95% CI, -17.5% to 39.0%]; P = .71). Age (t = -2.28; P = .02), but not sex (t = 0.92; P = .36) or APOE genotype (t = 1.06; P = .29) modified the association between baseline [ 18 F]flortaucipir PET and cognitive change, such that older individuals showed faster cognitive decline at similar tau PET levels. Conclusions and relevance: The findings of this prognostic study suggest that tau PET is a promising tool for predicting cognitive change that is superior to amyloid PET and MRI and may support the prognostic process in preclinical and prodromal stages of AD.

59 BASIC BIOLOGICAL SCIENCES↗

Herpes simplex virus and rates of cognitive decline or whole brain atrophy in the Dominantly Inherited Alzheimer Network

Abstract Objective To investigate whether herpes simplex virus type 1 (HSV‐1) infection was associated with rates of cognitive decline or whole brain atrophy among individuals from the Dominantly Inherited Alzheimer Network (DIAN). Methods Among two subsets of the DIAN cohort (age range 19.6–66.6 years; median follow‐up 3.0 years) we examined (i) rate of cognitive decline ( N = 164) using change in mini‐mental state examination (MMSE) score, (ii) rate of whole brain atrophy ( N = 149), derived from serial MR imaging, calculated using the boundary shift integral (BSI) method. HSV‐1 antibodies were assayed in baseline sera collected from 2009–2015. Linear mixed‐effects models were used to compare outcomes by HSV‐1 seropositivity and high HSV‐1 IgG titres/IgM status. Results There was no association between baseline HSV‐1 seropositivity and rates of cognitive decline or whole brain atrophy. Having high HSV‐1 IgG titres/IgM was associated with a slightly greater decline in MMSE points per year (difference in slope − 0.365, 95% CI: −0.958 to −0.072), but not with rate of whole brain atrophy. Symptomatic mutation carriers declined fastest on both MMSE and BSI measures, however, this was not influenced by HSV‐1. Among asymptomatic mutation carriers, rates of decline on MMSE and BSI were slightly greater among those who were HSV‐1 seronegative. Among mutation‐negative individuals, no differences were seen by HSV‐1. Stratifying by APOE4 status yielded inconsistent results. Interpretation We found no evidence for a major role of HSV‐1, measured by serum antibodies, in cognitive decline or whole brain atrophy among individuals at high risk of early‐onset AD.

Warren‐Gash, Charlotte↗

Fungal Spore Seasons Advanced Across the US Over Two Decades of Climate Change

Abstract Phenological shifts due to climate change have been extensively studied in plants and animals. Yet, the responses of fungal spores—organisms important to ecosystems and major airborne allergens—remain understudied. This knowledge gap limits our understanding of their ecological and public health implications. To address this, we analyzed a long‐term (2003–2022), large‐scale (the continental US) data set of airborne fungal spores collected by the US National Allergy Bureau. We first pre‐processed the spore data by gap‐filling and smoothing. Afterward, we extracted 10 metrics describing the phenology (e.g., start and end of season) and intensity (e.g., peak concentration and integral) of fungal spore seasons. These metrics were derived using two complementary but not mutually exclusive approaches—ecological and public health approaches, defined as percentiles of total spore concentration and allergenic thresholds of spore concentration, respectively. Using linear mixed‐effects models, we quantified annual shifts in these metrics across the continental US. We revealed a significant advancement in the onset of the spore seasons defined in both ecological (11 days, 95% confidence interval: 0.4–23 days) and public health (22 days, 6–38 days) approaches over two decades. Meanwhile, total spore concentrations in an annual cycle and in a spore allergy season tended to decrease over time. The earlier start of the spore season was significantly correlated with climatic variables, such as warmer temperatures and altered precipitations. Overall, our findings suggest possible climate‐driven advanced fungal spore seasons, highlighting the importance of climate change mitigation and adaptation in public health decision‐making.

Environmental Sciences & Ecology↗

Continental-scale patterns of extracellular enzyme activity in the subsoil: an overlooked reservoir of microbial activity

Abstract Chemical stabilization of microbial-derived products such as extracellular enzymes (EE) onto mineral surfaces has gained attention as a possibly important mechanism leading to the persistence of soil organic carbon (SOC). While the controls on EE activities and their stabilization in the surface soil are reasonably well-understood, how these activities change with soil depth and possibly diverge from those at the soil surface due to distinct physical, chemical, and biotic conditions remains unclear. We assessed EE activity to a depth of 1 m (10 cm increments) in 19 soil profiles across the Critical Zone Observatory Network, which represents a wide range of climates, soil orders, and vegetation types. For all EEs, activities per mass of soil correlated positively with microbial biomass (MB) and SOC, and all three of these variables decreased logarithmically with depth ( p < 0.05). Across all sites, over half of the potential EE activities per mass soil consistently occurred below 20 cm for all measured EEs. Activities per unit MB or SOC were substantially higher at depth (soils below 20 cm accounted for 80% of whole-profile EE activity), suggesting an accumulation of stabilized (i.e. mineral sorbed) EEs in subsoil horizons. The pronounced enzyme stabilization in subsurface horizons was corroborated by mixed-effects models that showed a significant, positive relationship between clay concentration and MB-normalized EE activities in the subsoil. Furthermore, the negative relationships between soil C, N, and P and C-, N-, and P-acquiring EEs found in the surface soil decoupled below 20 cm, which could have also been caused by EE stabilization. This finding suggests that EEs may not reflect soil nutrient availabilities deeper in the soil profile. Taken together, our results suggest that deeper soil horizons hold a significant reservoir of EEs, and that the controls of subsoil EEs differ from their surface soil counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗