Occupancy with multi-season, state, and scale?

[cnagy]

Hi everyone,
Occupancy modeling question here.

I am wondering if there is a way in Mark or Presence, or any other software, to run all of the multi’s: multi-season, multi-state, multi-scale. I know that (but haven’t dived into it), Mark does multi-season multi-scale, but I have not found anything that seems to do all three.

If you’d like more detail:
We’re examining 8 years of camera trap data from parks and other greenspaces in a large city, looking at patterns of site colonization and extinction for coyotes and possibly other species.
That leads me to a multi-season design, where we can estimate col and ext rates across many years.

I’m also interested, re: coyotes, as to which sites have a single coyote, a pair, and those where pups are photographed, and possibly what covariates would make a park suitable for a lone coyote vs allowing a pair to consistently raise pups.

This leads me to look at multi-state models, with, eg, 0, 1 (a single coyote), 2 (a pair), 3 (pups observed).
Also, we have varying numbers of cameras in each greenspace, primarily based on the size of the park. It ranges from 1 – 6 cameras per space. It would be nice to examine patterns of what areas or covariates are associated with coyote activity, or the presence of pups, for example, within occupied greenspaces.
So this makes me think of multi-state models.

If all three are not possible, I thought perhaps doing a multi season + state analysis at the park level (pooling all cameras in each park), followed by a multi-season + scale analysis using detection histories for all cameras that only looks at breeding state, since what makes a place suitable for breeding (versus not) is the main question we’d be interested in within parks.

Any advice is greatly appreciated.
Thanks!