Research Program · Quantitative, Population & Community Ecology

Environmental history, persistence, and species interactions

The distributions and communities we observe today reflect both contemporary demography and evolutionary history: related species may retain similar performance surfaces, while populations within a species carry genetic differences reflecting past environments. My research develops quantitative approaches grounded in ecological and evolutionary theory to connect genomic, environmental, and phylogenetic data — translating heterogeneous biological records into testable predictions with quantified uncertainty.

−22% / −45%
Mean / max projected range reduction when stochastic growth replaces climate means (XSDM)
>100,000
Verbatim sentence-anchored ecological interactions in EcoTagger with blind human review
8,014 tips
Dated mammalian phylogeny in PhyloSeek (replicate-tree uncertainty in development)
11 on CRAN
Maintained open-source R & C++ packages (~34,400 cumulative downloads)

Program synthesis

Three coupled ecological axes + open ecoinformatics
Integrative quantitative ecology research program diagram Three interconnected axes — Axis 1 Niche Evolution, Axis 2 Stochastic Demography and Coexistence, and Axis 3 Ecological Interaction Networks — linked horizontally by mechanistic arrows and supported beneath by an auditable biodiversity informatics foundation. UNIFYING QUESTION Given fluctuating environmental sequences rather than static means, which populations and species combinations persist? AXIS 1 · NICHE EVOLUTION How evolutionary history shapes environmental niches BM / OU evolving niche ellipsoids • Conserved response vs. accessible environment (M) & sampling bias • Withheld-record & clade simulations • Quaternary stability → pop-gen (SPA) arfun · xnicher · Ecography 2025 AXIS 2 · STOCHASTIC DEMOGRAPHY Performance curves, persistence & multi-species coexistence pt+1 = λtpt • Nonlinear averaging across concave curves lowers long-run growth • Experimental vital-rate priors • Next: ordering, stage & invasion rate xsdm · xsdmbayes · bioRxiv AXIS 3 · INTERACTION NETWORKS Ancestry vs. encounter in community networks observed (+) vs. uncertain (U) • Inherited compatibility vs. shared climate, geography & sampling • Parasitic, mutualistic & trophic webs • Positive–unlabeled (PU) inference phylopred · interseekr · PNAS 2022 dated clades & pop-gen annual climate & PAMs verbatim-verified edges AUDITABLE BIODIVERSITY INFORMATICS & HIGH-PERFORMANCE COMPUTING FOUNDATION PhyloSeek & GenoSeek 8,014-tip mammal trees, GPU MSA (genoaligner) & SRA mining xbioclim & xpam Year-resolved ERA5-Land series & C++17 streaming PAM/Jaccard EcoTagger & WildSeek >100k sentence-traced edges, NCBI/GBIF taxonomy & QA portal 11 CRAN + Stan / C++ Uncertainty quantification, held-out & simulation benchmarks
How the three directions connect: Ancestral evolution shapes the position, breadth, and asymmetry of species' environmental response curves (Axis 1); those curves govern long-term stochastic growth rates and multi-species coexistence under fluctuating climate sequences (Axis 2); and divergent or conserved demographic responses rewire encounter opportunities and specialization across ecological interaction networks (Axis 3). Every test rests on open, versioned data and software where primary records remain auditable.

Core research directions

Direction 01 · Eco-Evolutionary Macroecology & Population Genetics

How much does evolutionary history tell us about a species' environmental niche?

Related species often occupy similar environments, yet similarity can reflect conserved physiological responses, shared accessible geography (M), or spatially uneven sampling. With Jorge Soberón, I co-develop arfun to model how ellipsoidal niche position and covariance evolve under Brownian Motion and Ornstein–Uhlenbeck processes along dated phylogenies. Using controlled simulations and withheld-record experiments on well-sampled clades, we quantify when relatives' information genuinely improves predictions for sparsely sampled species and when it misleads.

Within species, I developed the Suitability Prevalence Area (SPA) theoretical framework (Ecography, 2025; co-first author) to reconstruct where environmental suitability persisted across the Late Quaternary and link that temporal stability to standing population-genetic diversity in fragmented montane populations (Sciurus aberti). Developing infrastructure in ancemb (with Henry Secaira-Morocho) and my independent phyloseek / genoaligner phylogenomic pipeline extends these comparative tests across deep clades.

Direction 02 · Stochastic Population Dynamics & Coexistence

What do differences in environmental response mean for persistence and coexistence?

Standard distribution models rely on long-term climate averages, yet population trajectories respond nonlinearly to year-to-year environmental sequences. With Daniel Reuman and collaborators, I co-developed XSDM (bioRxiv), embedding Lewontin–Cohen stochastic demography (pt+1 = λtpt) and long-term stochastic growth rates into geographic range modeling. Because nonlinear averaging across concave performance curves depresses long-run growth even where mean climate appears benign, accounting for interannual variability reduced projected distributions by 22% on average and up to 45% across focal species.

Ongoing work develops Bayesian XSDM (xsdmbayes) informed by experimental thermal and physiological performance curves to quantify parameter and forecast uncertainty. Because the baseline model estimates single-species density-independent growth from temporal means of log-multipliers, my next theoretical and empirical steps extend the framework to temporal ordering and autocorrelation, stage-structured demography, and density-dependent low-density invasion rates — testing when response-breadth and asymmetry differences promote or erode multi-species coexistence under shared climate fluctuations.

Tools & data: xsdm (CRAN) xsdmbayes (in dev) xbioclim.org maxentcpp
Direction 03 · Community Assembly, Ecological Networks & Ecoinformatics

When does evolutionary ancestry explain specialization in ecological networks?

In ecological networks, narrow partner breadth can arise from conserved trait compatibility or simply from which species co-occur and are sampled together. Across wildlife host–pathogen systems (PNAS, 2022), ectoparasite assemblages (Parasitology, 2025), and parasitic-plant networks (Phoradendron; Am. J. Bot., 2025), I showed that phylogenetic, environmental, and geographic distances carry system-specific explanatory weight — treating relatedness as a testable mechanistic hypothesis rather than a universal proxy.

Evaluating network assembly at macroecological scales requires traceable primary data and inference methods that treat unobserved pairs as uncertain rather than true negatives. I designed and maintain EcoTagger, which normalizes taxa against NCBI Taxonomy and the GBIF Backbone and accepts an interaction only when anchored to a verbatim sentence or table row (>100,000 records across host–parasite, mutualistic, and trophic systems, audited via a blind multi-reviewer portal). On this foundation, I develop positive–unlabeled learning methods (phylopred, xplus, xpam, interseekr) and prospective spillover-surveillance prototypes (WildSeek, in development) to test how climate variability rewires encounter opportunities and context-dependent interactions.

Representative publications & preprints

Full publication list in CV →
  • ECOGRAPHY · 2025 · Co-first author

    Hernández, N.A.H.*, Robles-Fernández, Á.L.*, & Upham, N. (2025). Ecography, 2025(3), e07202.

    Introduces the Suitability Prevalence Area (SPA) theoretical framework, demonstrating that Late Quaternary persistence of suitable habitat predicts contemporary genetic diversity across fragmented montane populations of Sciurus aberti.

  • PNAS · 2022 · First author

    Robles-Fernández, Á.L., Santiago-Alarcón, D., & Lira-Noriega, A. (2022). Proceedings of the National Academy of Sciences, 119(35), e2122851119.

    Integrates phylogenetic, environmental, and geographic distances across global host–pathogen networks (coronaviruses, West Nile virus, avian malaria), showing that the relative weight of evolutionary history vs. ecological overlap varies systematically among disease systems.

  • BIORXIV PREPRINT · v2026 · Stochastic demography & range forecasting

    Berti, E., Robles-Fernández, Á.L., Rosenbaum, B., Peterson, T.A., Soberón, J., & Reuman, D.C. (2026). bioRxiv, 2024.10.30.621023.

    Couples stochastic population growth with interannual climate time series in XSDM (implemented in the CRAN package xsdm); accounting for temporal variability reduces projected species distributions by 22% on average and up to 45%.

  • AMERICAN JOURNAL OF BOTANY & PARASITOLOGY · 2025 · Comparative interaction networks

    Tinoco-Domínguez et al. (2025, Am. J. Bot. 112:e70025) · Krasnov et al. (2025, Parasitology 152:745–756).

    Extends comparative phylogenetic–environmental network inference to hemiparasitic plants and continental mammal–ectoparasite assemblages.

Open science & teaching: As an RStudio Certified Trainer and maintainer of 11 CRAN packages, I pair quantitative research with open educational tools such as EcoSeek Learning.

Download current CV (PDF) →