Regional monitoring strategies frequently employ a nested sampling design where a finite set of study areas from throughout a region are selected within which intensive sub-sampling occurs. This sampling protocol naturally lends itself to a hierarchical analysis to account for dependence among sub-samples. Implementing such an analysis within a classic likelihood framework is computationally prohibitive with species occurrence data when accounting for detection probabilities. Bayesian methods offer an alternative framework to make this analysis feasible. We demonstrate a general approach for estimating occupancy when data come from a nested sampling design. Using data from a regional monitoring program of wood frogs (Lithobates sylvaticus) and spotted salamanders (Ambystoma maculatum) in vernal pools, we analyzed data using static and dynamic occupancy frameworks. We analyzed observations from 2004-2013collected within 14 protected areas located throughout the northeast United States . We use the data set to estimate trends in occupancy at both the regional and individual protected area level. We show that occupancy at the regional level was relatively stable for both species. Much more variation occurred within individual study areas, with some populations declining and some increasing for both species. We found some evidence for a latitudinal gradient in trends among protected areas. However, support for this pattern is overestimated when the hierarchical nature of the data collection is not controlled for in the analysis. For both species, occupancy appeared to be declining in the most southern areas, while occupancy was stable or increasing in more northern areas. These results shed light on the range-level population status of these pond-breeding amphibians and our approach provides a framework that can be used to examine drivers of change including among-year and among-site variation in occurrence dynamics, while properly accounting for nested structure of data collection.