| iNat page showing representative reptile species included in this analysis. |
Changes 2018-2024
| iNat page showing representative reptile species included in this analysis. |
In my previous post "Sonoran Desert Tortoise Population Status", I argued that a decline in the ratio of tortoise observations to total iNaturalist observations indicated a possible decline in the actual population of tortoises.
However, that conclusion rests on the assumption that opportunistic "citizen science" observers have not changed their preferences in photographing other animals. This could happen if people became less interested in tortoises, or if they became more interested in other taxa.
Part 1: Trends in Other Taxa
I analyzed various other well-represented taxa in the iNat data to look for possible trends in observer preferences. Specifically, I looked at all Research Grade (RG) iNat observations of amphibians, reptiles, insects, birds, and plants in the study area:
| Data table with total RG amphibian, reptile, insect, bird, and plant observations in the study area each year. |
| Count of observations by taxa 2013-2024 |
Observations of all taxa have increased over the last 10 years, but some appear to have increased more than others. I divided the count of each taxa's observations by total observations to investigate proportional changes in the ratio of various taxa:
| Ratio of various taxa to total observations 2013-2024 |
Instead, I think this must be an indication of a changed bias toward plants, perhaps due to specific iNat users who focus on plants, or a general trend toward more plant-focused observers on iNat. Plants, because they don't run away, are arguably the easiest organisms to photograph; it could be that as iNat has grown there are now more casual users biased toward photographing easier to observe organisms.
Regardless of why iNat users are observing proportionally more plants, this appears to be a source of bias that should be corrected in my analysis.
To correct for the large increase in plant observations, I compared various taxa to total non-plant observations:
| Ratio of various taxa to total non-plant ("total2") observations 2013-2024 |
Except for reptiles, all taxa ratios stayed within approximately +/- 15%:
| Table: percent change in ratios of various taxa since 2013 and 2017. |
Part 2: Another Look at the Tortoise Trend
The analysis in Part 1 led to a refinement in the total observation count used to create ratios, and helped set a baseline for expected change in ratios over the last 10 years. I used this information to reassess the observed decline in ratio of tortoise observations.
Tortoise observations were compared to various taxa and to total non-plant observations:
| Percent change in the ratio of tortoises to various taxa and to total non-plant ("total2") observations 2013-2024 |
While most animal taxa ratios show less than 15% change from baseline (Part 1), tortoises are observed 66-84% (mean 72%) in 2024 compared to 2017, and 35-57% (mean 43%) in 2024 compared to 2013:
| Table: percent change of tortoise observations to various taxa and to total non-plant observations. |
Conclusions
Tortoises show large declines compared to various representative animal taxa and compared to all non-plant observations. This is the same result I found in my original blog post when I compared tortoises to all observations.
The result was not affected by removing a potential source of user bias, showing that the original result is robust to some observational biases. Of course, there will always be more sources of bias that could be analyzed and corrected for. However, the fact that correcting for one source of bias didn't change the result makes me somewhat more confident that this result is directionally correct.
Conversely, the fact that I did find a large source of observer bias makes me wonder whether there are other large biases in iNat observation trends. Without analyzing all sources of bias (why did reptiles change in 2017?) these results must remain clouded by potential uncertainties.
The result was also not affected by which taxa I compared tortoises with, showing that the original result is robust to choice of comparison. Tortoises appear to be declining, whether they are compared to all observations, all animals, amphibians, reptiles, insects, or birds. While it is possible that one or more of these taxa are affected by observational bias, the fact that they all point in the same direction makes me more confident that this result is directionally correct and reflects an actual downward population trend.
Abstract
Population trend data for species under consideration for federal protection is often limited. This study evaluated Sonoran desert tortoise (Gopherus morafkai) population trends in Arizona using iNaturalist citizen science data. We analyzed 1,402 research-grade tortoise observations and normalized them against total observations in the area to account for sampling bias. While both tortoise observations and total observations increased since the early 2010s, the ratio of tortoise to total observations declined consistently from 0.12% to 0.08% between 2017-2025. This decreasing ratio suggests either shifting observer preferences or declining tortoise encounter rates, potentially indicating population decline. These results provide concerning evidence of negative population trends for this species of conservation concern and demonstrate the utility of citizen science data for monitoring species lacking formal survey programs.
Introduction
Sonoran desert tortoise (Gopherus morafkai) live in the Sonoran desert of Arizona and northern Mexico. A closely related species in California and Nevada, the Mojave desert tortoise (Gopherus agassizii), is Federally protected as a Threatened species under the Endangered Species Act (ESA). The Sonoran desert tortoise is listed by AZ as a Species of Greatest Conservation Need (SGCN) and was proposed for protection under the ESA but determined "Not Warranted" in 2022. On June 4, 2025 wildlife groups filed suit in federal court challenging that determination.
It is difficult to find current population trend data for species that may warrant listing under the ESA, so I attempted to do so using publicly available iNaturalist data. iNaturalist is a public repository of photographic observations of species.
Methods
I filtered iNat for Research Grade observation of Sonoran desert tortoises in AZ and found 1,402 observations. I considered limiting the analysis to only live tortoises, but eventually decided to use all research grade observations due to the difficulty filtering. Only ~3% (44 observations) were annotated as dead , and only ~1% (12 observations) were annotated as scat.
I created an Area of Interest (AOI) bounding box around the tortoise data to search for total RG observations. There were a total of 1.2 million RG observations in the Sonoran desert of southern AZ.
| RG tortoise observations in AZ with bounding box showing AOI. |
| Data table of AOI with Total RG Observations, Total Observations, RG:Total Observations Ratio, RG Tortoise Observations, and RG Tortoise: RG Total observations. |
Total RG observations in the AZ Sonoran desert have been increasingly consistently since the early 2010s.
| Count of total RG observations in the AOI per year. |
| Count of tortoise observations per year. |
By comparing tortoise observations to total observations, a normalized ratio can be derived. This ratio decreases over time. There is a step change decrease in 2017, which could be a real population decrease or could be due to the huge increase in popularity of iNat in 2017. However, even if only looking at 2017 -2025, there is still a long term decrease in the ratio of tortoise observations, from around 0.12% to 0.08% of total observations.
The decreasing ratio of tortoise observations indicates that either people are trending to preferentially observe other species more than tortoises, or that tortoises are making up a smaller percentage of the animals people encounter.
Despite limitations of citizen science data, these results provide concerning evidence of negative population trends for this species of conservation concern and highlight the utility of iNaturalist data for monitoring species that lack formal survey programs.
The iNaturalist website collects species observations from people all over the world. It started in 2008 and grew slowly at first and then entered a period of rapid growth in 2017. As a consequence, the number of species recorded on the website is constantly increasing, passing 300,000 in 2020. The website is currently adding more than 50 million observations a year. This raises an interesting biodiversity question: how long can the number of species keep increasing? Another way of stating the question: how many species are there?
Biodiversity scientists use species accumulation curves to estimate the total number of species in a given area. As they investigate a new study site, they record new species and the date/time the species was observed. For most sites, the number of new species increases rapidly as scientists describe common species; the number of new species slows as scientists search for more and more rare species. Graphing the number of species over time should reveal a logarithmic curve. Based on the equation for that curve, scientists can estimate the asymptote - the number of species the curve will eventually reach given enough time. This allows scientists to estimate the total number even if they don't finish counting all of the species.
| From iNat 2024 Year in Review |
| From 2024 Year in Review |
Nutritional yeast is tasty and nutritious! How nutritious is it? Well, that depends on whether it is fortified or not, and how much information you find on the nutrition label.
After looking at the labels of about a dozen brands, I found that nutritional yeast is consistently a good source of B2, B3, and B5, and with fortification is a great source of all major B vitamins except probably B7.
Fortified nutritional yeast contains extra B vitamins and sometimes iron. These extra nutrients are added during the processing of the product to make it more nutritional. The vitamin content of fortified nutritional yeast depends on how much it is fortified.
For example, here are two popular brands, Braggs and Bob's Red Mill:
| % Daily Value (DV) per 2 tablespoon serving. Data reported from nutritional labels. |
Both brands add supplement B1, B2, B3, B6, B9, and B12, but they add different amounts. In general, both seem to try to add more than 200% and less than 500% of the daily value. Bragg adds more of most B vitamins, but Bob's Red Mill adds more Folate (B9). Neither brand reports any B5 or B7, but that doesn't mean they don't have any. According to my analysis of unfortified nutritional yeast, they probably just didn't bother to measure and/or report it.
Non Fortified or Unfortified Nutritional Yeast contains only the vitamins and minerals found in the yeast cells and has no added vitamins and minerals. It contains variable amounts of B vitamins due to differences in yeast strain, growth substrate, growing conditions, and deactivation temperature.
While some brands (like Anthony's) don't report any nutritional values, I was able to find 7 brands on Amazon that reported nutritional testing results for their products:
| % Daily Value (DV) per 2 tablespoon serving, data reported from nutritional labels. |
Another way to visualize this information is by Brand:
| % Daily Value (DV) per 2 tablespoon serving, data reported from nutritional labels. |
The differences between brands are larger than the similarities. Some brands reported zero or almost zero B1, B2, B6, and B7. Conversely, other brands reported over 100% DVs of B1, B2, B5, and B6. Only B3 and B5 were consistently reported with decent amounts across all brands.
If we assume unreported values are not zero but were just not measured and look only at the vitamins reported by each brand, the mean vitamin content of unfortified nutritional yeast looks pretty good. However, the standard deviations are as large as the reported values for 4 of the 7 B vitamins. While some brands are good sources of B1 and B6, some were quite low for these nutrients. It seems that unfortified nutritional yeast can only be relied on to supply more than 20% daily value per serving for B2, B3, and B5:
Four of the seven unfortified brands reported B7, so it is likely that the other unfortified and fortified brands also contain this nutrient, but fail to test and/or report it. However, even the brands that report B7 report such a small and variable amount, nutritional yeast probably cannot be relied on as a good source of this important vitamin.
** Update **
I reached out to Sari Food to ask if their package listed nutritional values for all B vitamins, or if some were omitted. They confirmed that some were omitted and provided me with a new nutritional analysis that differed significantly from the nutrition reported on the package.
The new analysis reported 7 B vitamins, compared to the package that only lists 3. Their analysis showed significant amounts of B1 and B7 that were not listed on the package. Interestingly, the new analysis also differed from the package for the 3 listed on the package. I think this shows the variability between different batches of the same product.
| % DV of B vitamins for Sari Nutritional Yeast |
| Updated with new analysis from Sari Foods. |
How
likely a given area is to encounter wildfire is important for planning and
wildfire mitigation. Historic or expected fire return statistics are often cited
for ecosystems in Arizona, but I was curious how often wildfire actually burns across
different Arizona ecosystems.
Figure 1 Example wildfire polygons around Clints Well, AZ
showing overlapping fires from newer (blue, labelled), to older (shades of
brown, unlabelled). Data sources include
WFIGS and GeoMAC.
I used the WFIGS Interagency Fire Perimeter GIS data, which has good data on wildfires from 2000-2023. I limited this analysis to USFS land in Arizona.
Out
of a total of 11.168 million acres of USFS land in AZ, wildfire has burned 4.8
million cumulative acres in the last 24 years.
This counts areas that burned more than once as additional acres. It includes natural and human-ignitions, as
well as wildfire managed for resource benefit.
Figure 2 Example WFIGS Interagency Fire Perimeters in AZ
Figure 3 Wildfire acreage over time in AZ. 2011 was the Wallow fire.
Vegetation Types
To
evaluate wildfire probabilities in Arizona ecosystems, I looked at the 15 most
common ecosystem types, as defined by the USFS Ecosystem Response Unit (ERU)
vegetation type GIS layer. Together,
these 15 ecosystems account for 9.1 million out of the 11.1 million acres of
USFS land in Arizona.
Figure 6 Example ERU polygons showing aspen (pink) and
mixed conifer around the San Franscisco Peaks, AZ.
To calculate the percent of each ERU burned per year, I divided total acres burned by total acres of ERU and divided that by 24 years. Spruce-Fire forest and Mixed Conifer is most likely to burn, whereas Mixed Conifer with Aspen is least likely. Ponderosa pine ecosystems rank in the middle, at around 3% chance.
This analysis counts acres more than once if they burned more than once in the 24 year time period. For example, Spruce Fir Forest ERU has more acres of wildfire than there are total acres of ERU. This does not mean that every acre burned, but some acres burned more than once.
|
ERU |
ERU Acres |
Wildfire Acres |
% burned in 24 years |
% burned per year |
|
Ponderosa
Pine Forest |
1,966,603 |
1,431,424 |
72.79% |
3.03% |
|
PJ
Woodland |
1,175,545 |
208,685 |
17.75% |
0.74% |
|
PJ
Evergreen Shrub |
1,136,221 |
311,254 |
27.39% |
1.14% |
|
Mojave-Sonoran
Desert Scrub |
779,939 |
386,363 |
49.54% |
2.06% |
|
Semi-Desert
Grassland |
730,015 |
300,189 |
41.12% |
1.71% |
|
Interior
Chaparral |
713,754 |
533,678 |
74.77% |
3.12% |
|
Juniper
Grass |
539,830 |
299,074 |
55.40% |
2.31% |
|
Colorado
Plateau / Great Basin Grassland |
367,114 |
41,812 |
11.39% |
0.47% |
|
Ponderosa
Pine – Evergreen Oak |
362,838 |
238,365 |
65.69% |
2.74% |
|
Madrean
Pinyon-Oak Woodland |
354,836 |
92,160 |
25.97% |
1.08% |
|
Mixed
Conifer - Frequent Fire |
349,006 |
304,104 |
87.13% |
3.63% |
|
Mixed
Conifer w/ Aspen |
242,169 |
9,782 |
4.04% |
0.17% |
|
Montane
/ Subalpine Grassland |
157,163 |
92,461 |
58.83% |
2.45% |
|
Spruce-Fir
Forest |
112,827 |
124,593 |
110.43% |
4.60% |
|
PJ
Grass |
96,016 |
8,995 |
9.37% |
0.39% |
|
Madrean
Encinal Woodland |
93,939 |
23,092 |
24.58% |
1.02% |
Figure 7 ERU acres, wildfire acres, percent burned in 24
years, and percent burned per year.
Table ranked from most to least common ERU.
Fire
return interval is the average length of time until fire returns at a given
point in the landscape. The chance that
any given acre burns depends on a large number of complex factors, including
when it last burned, the topography, fuel reduction treatments, proximity to
WUI and/or human use. Still, percent
burned per year in the table above (Wildfire/Year, W) can be used to calculate
expected return intervals of fire, all else being equal.
To calculate expected return intervals, first calculate the probability (P) that fire will not occur in a given span of time (X).
P = (1-W)^X
For
example, for Ponderosa Pine Forest over 10 years:
P = (1-0.0303)^10
P = (0.9697)^10
P = 73.5% chance that fire will not occur, or 26.5% chance that fire will occur in 10 years.
20
years:
(0.9697)^20=54% chance that fire will not occur, or 46% chance that fire will occur.
Figure 8 Cumulative probability of wildfire in AZ
Ponderosa Pine ERU
If we determine a Probability, but need to know the span of time until fire occurs, we can solve for X:
P = (1-W)^X
P = log x / log (1-W)
For
example, if we determine "expected return interval" to be the length
of time necessary for 50% chance of fire:
0.5 = (0.9697)^X
X = log (0.5) / log (0.9697)
X = 22 years until there is a 50% chance of fire in Ponderosa Pine Forest.
However, if we interpret "expected return interval" to be the length of time necessary for 90% chance of fire:
0.1 = (0.9697)^x
X = log (0.1) / log (0.9697)
X
= 75 years
Over time, the probability approaches, but never actually reaches, 100% that a wildfire will occur:
Figure 9 Cumulative Probability of Fire in Ponderosa Pine
ERU
The length of time until fire returns at a given point in the landscape depends on how certain we want to be of the chance of fire. If we want to be very certain (90% probability), then we would expect to wait 75 years on average. If we are OK taking the flip of a coin (50% probability), than we would expect fire to return at any given point in 22 years. If we are risk adverse, and can only tolerate a 10% chance of fire visiting our chosen point, we should expect fire every 3.5 years, on average.