Description of data and code

csrb_pop_df_JAGS.csv

This is a CSV file of data used to perform final analyses in the population study

Column names (variables):

• Date: date that the lure was retrieved (YYYY-MM-DD)
• Site: site number corresponding to Station
• Station: location at which local springflow (Station_Q) was measured (corresponds to locations listed on Figure S24)
• Sample.period: ordinal sampling period during the population study corresponding to periods listed in Table S1
• Sample: ordinal number assigned to each lure within each sampling period and locality
• Heterelmis.adults: number of Heterelmis comalensis adults on each lure
• Heterelmis.larvae: number of Heterelmis comalensis larvae on each lure
• Microcylloepus.adults: number of Microcylloepus pusillus adults on each lure
• Microcylloepus.larvae: number of Microcylloepus pusillus larvae on each lure
• Biofilm.total: estimated total coverage of biofilm on the lure (percent)
• Station_Q: measured local springflow at each station (cubic feet per second)
• Station_Q_30D: estimated 30-day local springflow average at each Station (cubic feet per second)

cs_fs_df.csv

This is a CSV file of springflow measurements used to generate the local springflow covariate

Column names (variables):

• Date: date that the measurement was taken (YYYY-MM-DD)
• Station: location at which local springflow (Q) was measured (corresponds to locations listed on Figure S24)
• Station_No: number for each Station
• Q: measured springflow at each station (cubic feet per second)
• Unit: device name used to measure springflow

csrb_fall_2024.csv

This is a CSV file of the fall 2024 biomonitoring data used to illustrate the framework for the biological monitoring application in Appendix C.

Column names (variables):

• Date: date that the lure was retrieved (YYYY-MM-DD)
• Site: site number corresponding to Station
• Station: location at which local springflow (Station_Q) was measured (corresponds to locations listed on Figure S24)
• Heterelmis.adults: number of Heterelmis comalensis adults on each lure
• Biofilm.total: estimated total coverage of biofilm on the lure (percent)
• Station_Q: measured local springflow at each station (cubic feet per second)
• Station_Q_30D: estimated 30-day local springflow average at each Station (cubic feet per second)

csrb_custom_functions.R

This is an R file with the code to generate custom functions for the population analysis.

R_Workflow_Local_Springflow_Covariate_Derivation_and_Relative_Abundance_Modeling.R

This is an R file with the code to run the analyses and generate the figures for the population analysis.

Comal-Springs-riffle-beetle-long-term-data_2025-11-05.csv

This is the full long-term dataset CSV file with all samples collected from 2004 to 2025 and associated covariates used in analyses.

Column names (variables):

• N: continuous number (unique identifier) for each sample
• SiteID: identifier for each individual site that was monitored
• Locality: name of the locality in which a lure was set
• Sector: unique identifier for spatial groupings of sites within each locality
• Method: method by which each site was sampled
  • rag: cotton rag used from 2004 to 2013
  • lure: cotton lure used from 2012 to 2025
• Date.start: date that the lure was set in a spring (if known) (YYYY-MM-DD)
• Date.end: date that the lure was removed from the spring (YYYY-MM-DD)
• Day.of.year: ordinal day of the year (1 = 1 January)
• Days.set: number of days that the lure was set in the system
• Days.set.est: number of days that the lure was set in the system, with number of days estimated for sites where Date.start was not known
• Sample.project: name of the project for which lures were set
• Lure.condition: final condition of the lure upon retrieval
  • dry: lure was above water
  • exposed: lure was exposed to light
  • good: lure was in good condition
  • lost: lure could not be found
  • moved: lure was otherwise out of its original location
• Disturbed: whether or not (yes, no) a lure was disturbed in the place it was set
• Heterelmis.larvae: number of Heterelmis comalensis larvae on each lure
• Heterelmis.adults: number of Heterelmis comalensis adults on each lure
• Microcylloepus.larvae: number of Microcylloepus pusillus larvae on each lure
• Microcylloepus.adults: number of Microcylloepus pusillus adults on each lure
• Temp180DayAv: average air temperature (degrees Celsius) over the previous 180 days from the date each lure was retrieved
• Precip30Day: total precipitation (millimeters) over the previous 30 days from the date each lure was retrieved
• Dis180Dav: average total system discharge (cubic feet per second) for Comal Springs over the previous 180 days from the date each lure was retrieved
• year: ordinal year in which each lure was retrieved
• month: ordinal month of the year in which each lure was retrieved
• day: ordinal day of the month in which each lure was retrieved
• Local_30Day_SF_Avg: average local springflow for each sector (cubic feet per second)

Pintar-Sullivan-CSRB-long-term-analysis.R

This is an R file with the code to run the analyses and generate the figures with the long-term dataset.