Patchy distribution fields: sampling distance unit of a zigzag survey and reconstruction adequacy

A mathematical model was used to examine the effects of choosingvarious units of sampling distance of a zigzag survey on the adequacy of reconstructing patchy distribution fields. The modelsimulates fish or plankton patches (or gaps) of different shapesand spatial orientations, and an acoustic survey by zigzag or parallel transects along which a unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). A priori information on the autocorrelation radii for the field in the directions of the survey (R_s) and perpendiculardirection (R_p) allows optimisation of the survey design andthe algorithm of data analysis. A field can be reconstructed properly (r² > 0.70) if the distance between transects D < (1.0–1.5)R_s and the unit of sampling distance d < (1.0–1.5)R_p. A posteriori determination of patchorientation allows reconstruction of the best field attainable onthe basis of the survey data. In cases of field movement, if thedimension of patches in the direction of movement exceeds that ofa surveyed area, a survey in the opposite direction gives bestresults; in contrast, if the dimension of moving patches is smaller than that of a surveyed area, it is reasonable to carryout a survey in the same direction. The criterion remains validwhen a survey is carried out by zigzag transects and a unit ofsampling distance is set along them. The results obtained indicate that, for a fixed transect spacing and a given number of sampling points on each full transect, zigzag pattern allowsless adequate reconstruction of an original distribution field (in cases of both immovable and movable fields) than corresponding parallel pattern.     

Patchy distribution fields: A spiral survey design and reconstruction adequacy

A mathematical model was used to examine the effects of a spiral survey design on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and acoustic surveys by a spiral of Archimedes; for comparison, surveys by parallel or zigzag transects are imitated. Adequacy of the reconstructed fields to those originally generated was evaluated by calculating their correlations (r). The mathematical experiments conducted showed that spiral surveys ensure, practically speaking, the same adequacy of field reconstruction (both in cases of immovable or movable patches) as do surveys by parallel or zigzag transects with greater sampling effort (overall path). In the case of a spiral survey, a patchy field can be reconstructed properly (r² = 0.70) if the overall survey path is not less than S/R_av = 20–30, where R_av is the autocorrelation radius averaged for various directions. Thus, the results obtained allow us to conclude that a spiral survey design is expedient in cases that the minimal duration of a survey is a decisive factor for its conduction and there is a priori information that no onshore-offshore gradients of fish density exist in a region under study.     

Patchy Distribution Fields: Sampling Distance Unit of a Spiral Survey and Reconstruction Adequacy

A mathematical model is used to examine the effects of choosing various units of sampling distance of a spiral acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and an acoustic survey by a spiral of Archimedes along which a unit of sampling distance is set. For comparison, surveys are imitated by parallel transects along which the same unit of sampling distance is set. Adequacy of the reconstructed fields to those originally generated is evaluated by calculating their correlations (r). In the case of a spiral survey, the mathematical experiments conducted show that an immovable field can be reconstructed properly (r ² > 0.70) if the ratio of the units of sampling distance to the autocorrelation radius for the field averaged in various directions d/R _av < 1.0–1.5. Regarding immovable fields, the spiral surveys ensure, practically speaking, the same adequacy of the field reconstruction as do surveys by parallel transects with the same unit of sampling distance. In regard to movable fields, a comparison of the results of spiral surveys with those of surveys by parallel transects indicates that the former may ensure even higher adequacy of the field reconstruction than do the latter, provided that the units of sampling distances in these surveys are equal to each other.     

Patchy Distribution Fields: Sampling Distance Unit of an Interleaved Survey and Reconstruction Adequacy

A mathematical model was used to examine the effects of choosing various units of sampling distance of an interleaved (two-pass) acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches (or gaps) of different shapes and spatial orientations, and an interleaved survey by parallel or zigzag transects, along which a unit of sampling distance is set. The efficiency of a survey is determined by the adequacy of a reconstructed field to that originally generated, which is evaluated by calculating their correlations (r). Regarding immovable fields, the experiments conducted indicate that a patchy field can be reconstructed properly (r ² > 0.70) if the distance between transects D < (1.5–2.0)R and the unit of sampling distance d < (1.5– 2.0)R _p. As, for regular surveys d < (1.0–1.5)R _p, it may be concluded that interleaved surveys are more efficient than regular ones because of the factor studied. In regard to movable fields, a comparison of the results of interleaved surveys with those of regular surveys directly indicates that the former may ensure a more adequate field reconstruction than the latter do. This fact confirms the previous conclusion that an interleaved survey is expedient in cases where there is no preference regarding the position of a vessel for further work.     

Patchy Distribution Fields: An Interleaved Survey Design and Reconstruction Adequacy

A mathematical model was used to compare the effects of a regular (one-pass) or interleaved (two-pass) acoustic survey on the adequacy of reconstructing patchy distribution fields. The model simulates fish or plankton patches of different shapes and spatial orientations, and a set of parallel or zigzag transects forming a regular or interleaved acoustic survey. The efficiency of a survey is determined by the adequacy of a reconstructed field to that originally generated, which is evaluated by calculating their correlations. Regarding the immovable fields, the efficiency of a regular or interleaved acoustic survey was tested with the following two alternative assumptions: (1) the entire survey was completed; (2) the survey was interrupted for some reason at the moment when one transect remained non-accomplished. In the former case, the efficiencies of both acoustic surveys were nearly the same; in the latter case, the efficiency of an interleaved survey was superior to that of a regular one. With respect to movable fields, the efficiency of the completed interleaved surveys was even higher than that of the regular ones. Thus, the results obtained allow us to conclude that an interleaved survey is expedient in cases where there is no preference regarding the position of a vessel for further work.     

Spatial Scale of Autocorrelation of Assemblages of Benthic Invertebrates in Two Upland Rivers in South-Eastern Australia and Its Implications for Biomonitoring and Impact Assessment in Streams

Spatial autocorrelation in ecological systems is a critical issue for monitoring (and a general understanding of ecological dynamics) yet there are very few data available, especially for riverine systems. Here, we report here on assemblage-level autocorrelation in the benthic-invertebrate assemblages of riffles in two adjacent, relatively pristine rivers in south-eastern Victoria, Australia (40-km reaches of the Wellington [surveys in summers of 1996 and 1997] and Wonnangatta Rivers [survey in summer of 1996 only], with 16 sites in each river). We found that analyses were similar if the data were resolved to family or to species level. Spatial autocorrelation was assessed by using Mantel-tests for the data partitioned into different sets of spatial separations of survey sites (e.g. 0–6 km, 6–12 km, etc.). We found strong small-scale (≤6 km) autocorrelation in the Wellington River, which is consistent with known dispersal abilities of many aquatic invertebrates. Surprisingly, there were strong negative correlations at longer distance classes for the Wellington River in one of the two summers (20–40 km) and the Wonnangatta River (12–20 km). That two largely unimpacted, adjacent rivers should have such different autocorrelation patterns suggests that impact assessment cannot assume dependence or independence of sites a priori. We discuss the implications of these results for use of “reference” sites to assess impacts at nominally affected sites.     

Modeling and Measurement of Structural Changes at a Landscape Scale in Dryland Areas

We present a technique to quantify and model the intensity of structural changes produced by management of dry grazing lands at a landscape scale. The technique is illustrated with the analysis of digitized black–white (b/w) imagery and an application to the study of changes induced by grazing gradients. Structural changes in patchy vegetation canopies were studied in the Patagonian Monte (Chubut, Argentina) at two resolution scales by means of linear transects in the field (50 m) and others drawn on aerial b/w photographs (2–5 km) of grazed paddocks. Spatial series of plant cover values along transects in the field and on photographs were analyzed with standard techniques of spectral analyses, including auto-correlation spectra and Fourier transforms. In order to test the internal consistency of the techniques used, synthetic plant canopies with patches of varying cover and size were generated by means of a stochastic model of plant growth under different stocking rates or after varying periods of recovery. The behavior of the simulation model is consistent with the observed dynamics of plant canopies in semiarid environments. There is a consistent relation between the number and geometric properties of plant patches (patch number, patch size, patch connectivity) and the signal/noise ratios of the Fourier decomposition describing plant density data. Signal/noise ratios corresponding to plant cover data in paddocks with different grazing treatment are consistent with the assumptions derived from modeled canopies, as well as those estimated from optical density of b/w aerial photographs of paddocks. We tested the hypotheses that patch arrangements as quantified by the signal/noise ratios vary in accordance with grazing gradients in paddocks with a permanent corner-located watering point. The use of digitized b/w images allows inspecting permanent changes over time periods when other types of images were not yet available.     

Abundance, behavior, and movement patterns of western gray whales in relation to a 3-D seismic survey, Northeast Sakhalin Island, Russia

A geophysical seismic survey was conducted in the summer of 2001 off the northeastern coast of Sakhalin Island, Russia. The area of seismic exploration was immediately adjacent to the Piltun feeding grounds of the endangered western gray whale (Eschrichtius robustus). This study investigates relative abundance, behavior, and movement patterns of gray whales in relation to occurrence and proximity to the seismic survey by employing scan sampling, focal follow, and theodolite tracking methodologies. These data were analyzed in relation to temporal, environmental, and seismic related variables to evaluate potential disturbance reactions of gray whales to the seismic survey. The relative numbers of whales and pods recorded from five shore-based stations were not significantly different during periods when seismic surveys were occurring compared to periods when no seismic surveys were occurring and to the post-seismic period. Univariate analyses indicated no significant statistical correlation between seismic survey variables and any of the eleven movement and behavior variables. Multiple regression analyses indicated that, after accounting for temporal and environmental variables, 6 of 11 movement and behavior variables (linearity, acceleration, mean direction, blows per surfacing, and surface-dive blow rate) were not significantly associated with seismic survey variables, and 5 of 11 variables (leg speed, reorientation rate, distance-from-shore, blow interval, and dive time) were significantly associated with seismic survey variables. In summary, after accounting for environmental variables, no correlation was found between seismic survey variables and the linearity of whale movements, changes in whale swimming speed between theodolite fixes, mean direction of whale movement, mean number of whale exhalations per minute at the surface, mean time at the surface, and mean number of exhalations per minute during a whales surface-to-dive cycle. In contrast, at higher received sound energy exposure levels, whales traveled faster, changed directions of movement less, were recorded further from shore, and stayed under water longer between respirations.     

Heterogeneous Traffic Induced Effects on Vertical Dispersion Parameter in the Near Field of Roadways - a Wind Tunnel Study

The varying traffic parameters such as traffic volume, speed, shape and size, and terrain roughness conditions play a vital role on dispersion of pollutants in the near field of roadways. Simulation experiments were carried out in the Environmental Wind Tunnel (EWT) to evaluate the traffic induced effects on vertical dispersion parameter (σ _z ) for heterogeneous traffic conditions in the near field of roadways for evaluating the effect of variations in traffic volume, terrain roughness condition and approaching wind direction. The model vehicle movement system was fabricated and made operational in the EWT, which allowed the variation in traffic volume, speed and wind road inclination. Sixty-six hydrocarbon tracer experiments were performed to evaluate σ _z in the near field of roadways for variable traffic volume, three terrain roughness conditions and two approaching wind directions (i.e., 90° and 60°). The values of σ _z for heterogeneous traffic conditions were found to be higher for low roughness conditions in comparison to other two higher roughness conditions for various traffic volumes and approaching wind directions. This revised version was published online in July 2006 with corrections to the Cover Date.     

Evaluating the ecological integrity of Atlantic forest remnants by using rapid ecological assessment

The need for quick identification of priority areas for biodiversity protection makes rapid assessment methods important management tools for defining conservation strategies. An increasingly used rapid assessment method is rapid ecological assessment (REA), a fast and flexible survey directed toward selected indicator species and vegetation forms. The purpose of this study was to propose and test REA based on plant community features of the semideciduous Atlantic forest (SAF). Correlation tests were performed between data collected by REA and plant species diversity, richness, and abundance collected by conventional woody plant inventory methods. The study was conducted in 21 SAF patches in Northern Paraná State, Brazil. The phytosociological inventory was conducted along a single transect and required 2 days to complete (excluding time spent for herbarium identification), whereas REA was conducted along three to four transects per working day. REA results correlated positively with woody plant diversity, proving REA to be an efficient method for defining the conservation status of SAF fragments, but accuracy of evaluations of threats to biological integrity are relatively low. Both the selection of appropriate variables and the skill level of field staff are critical and can strongly influence REA results.     

Spatial factors of white-tailed deer herbivory assessment in the central Appalachian Mountains

Because moderate to over-abundant white-tailed deer (Odocoileus virginianus) herbivory impacts biodiversity and can alter community function, ecological benchmarks of herbivory impact are needed to assess deer impacts. We evaluated spatial patterns of deer herbivory and their relation to herbivory assessment by evaluating woody vegetation along 20 transects at each of 30 sites spread across a wide range of deer herd densities and vegetative condition throughout the biodiverse Appalachian Mountains of Virginia, USA. Surprisingly, herbivory patterns and the availability of woody forage generally were unchanged among physiographic regions and land use diversity classes. However, some relationships between browsing pattern and vegetation varied with scale. The total quantity of vegetation browsed on a given site and at the transect scale were related positively to the availability of forage, as the proportion of stems browsed decreased as stem density increased. However, this was only true when all stems were considered equally. When stem densities by species were weighted for deer preference, the proportion of stems browsed had no relationship or increased with stem density. Compared to the value from all transects sampled, on average, the mean of ≥?3 transects within a site was within 0.1 of the browsing ratio and stem densities were within 0.5 stems m^?2. Our results suggest that one transect per square kilometer with a minimum of three transects may be sufficient for most browsing intensity survey requirements to assess herbivory impacts in the Appalachian region of Virginia. Still, inclusion of spatial factors to help partition variation of deer herbivory potentially may allow for improved precision and accuracy in the design of field herbivory impact assessment methods and improve their application across various landscape contexts.     

Relationship Between Bird Abundances and Landscape Characteristics: The Influence of Scale

Scale is important to consider when investigating effects of the environment on a species. Breeding Bird Survey (BBS) data and landscape metrics derived from aerial photographs were evaluated to determine how relationships of bird abundances with landscape variables changed over a continuous range of 16 spatial scales. We analyzed the average number of birds per stop (1985–1994) for five songbird species (family Cardinalidae) for each of 50 stops on 198 BBS transects throughout six states in the Central Plains, USA. Land along each transect was categorized into six cover types, and landscape metrics of fractal dimension (a measure of shape complexity of habitat patches), edge density, patch density, and percent area were calculated, with principal components used to construct composite environmental variables. Associations of bird abundances and landscape variables changed in accordance with small scale changes. Abundances of three species were correlated with edge density and one with component I, which subsumes initial variables of patch density for urban, closed forest, open forest, and open country. Fractal dimension and component II (summarizing amount of closed forest versus open country) were associated with the most species. Correlation patterns of fractal dimension with northern cardinal (Cardinalis cardinalis) and painted bunting (Passerina ciris) abundances were similar, with highest correlations at intermediate to small scales, suggesting indirectly that these species thrive in areas where local habitat conditions are most important. Multiscale analysis can provide insight into the spatial scale(s) at which species respond, a topic of intrinsic scientific interest with applied implications for researchers establishing protocols to assess and monitor avian populations.     

Accuracy assessment of a small-area method for estimating the spatial distribution of the degree of tree damage

Tree damage, gauged by the amount of defoliation, is one of the basic criteria used to determine treatments for protected and economic forests. Monitoring should include an assessment of the degree of tree damage in different spatial scales. Therefore, in addition to the commonly applied large-area methods, small-area methods should be used. The aim of the paper is to present the results of the accuracy assessment of a small-area method, proposed by Podlaski (2005) [Podlaski, R. (2005). Inventory of the degree of tree defoliation in small areas. Forest Ecology and Management, 215, 361–377], for monitoring the degree of tree damage. The degree of tree damage was shown in sub-blocks P₃ of the system of information on natural environment (SINUS). To estimate the spatial distribution of the degree of tree defoliation, survey sampling, based on simple random sampling with replacement (SRSWR), was used. The degree of damage to fir (Abies alba Mill.) and beech (Fagus sylvatica L.) was analysed in the Święty Krzyż forest section in the Świętokrzyski National Park. The maximum total estimation errors for the proportion of trees with a degree zero of damage, and with second and third degrees of damage together (for α = 0.05) were at most 30.8% for fir and 24.3% for beech trees. For standard, small-area evaluations, these are satisfactory values. In the Święty Krzyż forest section, the number of P₃ sub-blocks with 0.00–5.00% of undamaged trees and with 80.01–100.00% of moderately- or severely-damaged trees was significantly greater for fir than for beech. These results indicate that the fir population was unhealthier than the beech group in the study area. P₃ sub-blocks of the SINUS system, in which the proportion of the healthiest trees was highest, were situated at the forest margin, bordering on meadows and arable fields (in the case of fir) and forming dense patches consisting of several sub-blocks, or occurring singly in the whole study area (in the case of beech). The results show the significant differentiation of forest tree health in small areas.     

Evaluation of Terrestrial and Streamside Salamander Monitoring Techniques at Shenandoah National Park

In response to concerns about amphibian declines, a study evaluating and validating amphibian monitoring techniques was initiated in Shenandoah and Big Bend National Parks in the spring of 1998. We evaluate precision, bias, and efficiency of several sampling methods for terrestrial and streamside salamanders in Shenandoah National Park and assess salamander abundance in relation to environmental variables, notably soil and water pH. Terrestrial salamanders, primarily redback salamanders (Plethodon cinereus), were sampled by searching under cover objects during the day in square plots (10 to 35 m²). We compared population indices (mean daily and total counts) with adjusted population estimates from capture-recapture. Analyses suggested that the proportion of salamanders detected (p) during sampling varied among plots, necessitating the use of adjusted population estimates. However, adjusted population estimates were less precise than population indices, and may not be efficient in relating salamander populations to environmental variables. In future sampling, strategic use of capture-recapture to verify consistency of p's among sites may be a reasonable compromise between the possibility of bias in estimation of population size and deficiencies due to inefficiency associated with the estimation of p. The streamside two-lined salamander (Eurycea bislineata) was surveyed using four methods: leaf litter refugia bags, 1 m² quadrats, 50 × 1 m visual encounter transects, and electric shocking. Comparison of survey methods at nine streams revealed congruent patterns of abundance among sites, suggesting that relative bias among the methods is similar, and that choice of survey method should be based on precision and logistical efficiency. Redback and two-lined salamander abundance were not significantly related to soil or water pH, respectively.     

An assessment and monitoring program for birds in the Lower Urubamba region,Peru

We developed an assessment and monitoring plan for birds in connection with the exploration and potential development of a large natural gas field in the Lower Urubamba drainage of Peru, a project of Shell Prospecting and Development Peru (SPDP). Our objectives were to: (1) inventory the birds in the area, including information on habitat use and abundance, and (2) devise long-term monitoring protocols for birds. We sampled birds through a combination of visual and auditory surveys and mist-netting at 4 well sites and 3 sites along the Urubamba and Camisea rivers. We recorded 420 speciesduring 135 days of field work. We consider the highest prioritiesfor a future monitoring program to be: (1) establish whether edge effects are occurring at well sites, along roads and alongthe planned pipeline route and determine the significance and extent of these effects and (2) assess the impact of increased human access to the area on game and other exploited species. The remoteness of the area, its rugged terrain and dense vegetation and the lack of trained personnel limit the choice of survey and monitoring methods. We recommend use of mist-netting and transects for monitoring edge effects and useof transects for monitoring game and other exploited species.     

Measuring differences between overlapping but unpaired spatial surveys using a geographic information system

To compare the results from two overlapping but unpaired synoptic surveys, a method using a geographic information system is proposed. The new method uses a GIS to create (i) minimum convex polygons (MCP) enclosing each set of survey data and (ii) Voronoi tesselations assigning area weights to each datum. The pairs of MCP and Voronoi maps are overlaid to produce a mosaic of polygons each with one datum assigned from each survey. The differences between the pairs weighted by the polygon areas provides the basis for statistical testing. Area-weighted means and variances of paired differences are computed and a z-statistic measures the significance of differences for the whole intersection area. A cross-product autocorrelation statistic provides an assessment of the spatial distribution of differences. Alternate, conventional methods are compared with the new method: analysis of variance (ANOVA), analysis of covariance (ANCOVA), and contour overlays. The four methods were applied to the comparison of two macrobenthic surveys conducted in Lake Erie. Measures for five different taxa were examined. The ANOVA and ANCOVA methods found many significant differences between surveys. Both methods were judged inappropriate as sampling data are not expected to be drawn from normally-distributed populations in spatial surveys. Differences between surveys were detected but were difficult to assess using contour overlays. The new method only found a significant difference in one taxon, an expected difference as members of that taxon were absent from most sites in one survey.     

A Spatial Model to Estimate Habitat Fragmentation and Its Consequences on Long-Term Persistence of Animal Populations

The increasing use of the landscape by humans has led to important diminutions of natural surfaces. The remaining patches of wild habitat are small and isolated from each other among a matrix of inhospitable land-uses. This habitat fragmentation, by disabling population movements and stopping their spread to new habitats, is a major threat to the survival of numerous plant and animal species. We developed a general model, adaptable for specific species, capable of identifying suitable habitat patches within fragmented landscapes and investigating the capacity of populations to move between these patches. This approach combines GIS analysis of a landscape, with spatial dynamic modeling. Suitable habitat is identified using a threshold area to perimeter ratio. Potential movement pathways of species between habitat patches are modeled using a cellular automaton. Habitat connectivity is estimated by overlaying habitat patches with movement pathways. The maximum potential population is calculated within and between connected habitat patches and potential risk of inbreeding within meta-populations is considered. The model was tested on a sample map and applied to scenario maps of predicted land-use change in the Peoria Tri-county region (IL). It (1) showed area of natural area alone was insufficient to estimate the consequences on animal populations; (2) underscored the necessity to use approaches investigating the effect of land-use change spatially through the landscape and the importance of considering species-specific life history characteristics; and (3) highlighted the model's potential utility as an indicator of species likelihood to be affected negatively by land-use scenarios and therefore requiring detailed investigation.     

Degradability and Sediment Sorption of an Alcohol Polyglycol Ether Surfactant Putatively Useful for the Control of Red Swamp Crayfish in Rice Fields

This work reports studies of the degradation rates of a fattyalcohol polyglycol ether non-ionic surfactant, Genapol OXD-080, putatively useful for the control of red swamp crayfish (Procambarus clarkii Girard) in rice fields under laboratory and field conditions. The influence of temperature,sediment site specificity and sorption were taken into account.The degradation kinetics of the surfactant depends on the experimental conditions: type of inocula and temperature. Thedistribution of this chemical in aquatic systems was also examined. Genapol OXD-080 was removed into the sedimentsreadily after application, and sorption was considered the majorpath of removal from the water phase. Data suggest that furtherstudies are required regarding the effects of Genapol OXD-080 in aquatic organisms resident in rice fields, in parallelwith the development of technologies related with the use ofsurfactants to control P. clarkii populations.     

Spatial Analysis of Volatile Organic Compounds from a Community-Based Air Toxics Monitoring Network in Deer Park,Texas, USA

In the summer of 2003, ambient air concentrations of volatile organic compounds (VOCs) were measured at 12 sites within a 3-km radius in Deer Park, Texas near Houston. The purpose of the study was to assess local spatial influence of traffic and other urban sources and was part of a larger investigation of VOC spatial and temporal heterogeneity influences in selected areas of Houston. Seventy 2-h samples were collected using passive organic vapor monitors. Most measurements of 13 VOC species were greater than the method detection limits. Samplers were located at 10 residential sites, a regulatory air monitoring station, and a site located at the centroid of the census tract in which the regulatory station was located. For residential sites, sampler placement locations (e. g., covered porch vs. house eaves) had no effect on concentration with the exception of methyl tertiary-butyl ether (MTBE). Relatively high correlations (Pearson r > 0.8) were found between toluene, ethylbenzene, and o,m,p-xylenes suggesting petroleum-related influence. Chloroform was not correlated with these species or benzene (Pearson r < 0.35) suggesting a different source influence, possibly from process-related activities. As shown in other spatial studies, wind direction relative to source location had an effect on VOC concentrations.     

Precision Measurements from Very-Large Scale Aerial Digital Imagery

Managers need measurements and resource managers need the length/width of a variety of items including that of animals, logs, streams, plant canopies, man-made objects, riparian habitat, vegetation patches and other things important in resource monitoring and land inspection. These types of measurements can now be easily and accurately obtained from very large scale aerial (VLSA) imagery having spatial resolutions as fine as 1 millimeter per pixel by using the three new software programs described here. VLSA images have small fields of view and are used for intermittent sampling across extensive landscapes. Pixel-coverage among images is influenced by small changes in airplane altitude above ground level (AGL) and orientation relative to the ground, as well as by changes in topography. These factors affect the object-to-camera distance used for image-resolution calculations. ‘ImageMeasurement’ offers a user-friendly interface for accounting for pixel-coverage variation among images by utilizing a database. ‘LaserLOG’ records and displays airplane altitude AGL measured from a high frequency laser rangefinder, and displays the vertical velocity. ‘Merge’ sorts through large amounts of data generated by LaserLOG and matches precise airplane altitudes with camera trigger times for input to the ImageMeasurement database. We discuss application of these tools, including error estimates. We found measurements from aerial images (collection resolution: 5–26 mm/pixel as projected on the ground) using ImageMeasurement, LaserLOG, and Merge, were accurate to centimeters with an error less than 10%. We recommend these software packages as a means for expanding the utility of aerial image data.