Geological Survey Topographic Maps Research Articles

SCALING-UP DEEP LEARNING PREDICTIONS OF HYDROGRAPHY FROM IFSAR DATA IN ALASKA

Abstract. The United States National Hydrography Dataset (NHD) is a database of vector features representing the surface water features for the country. The NHD was originally compiled from hydrographic content on U.S. Geological Survey topographic maps but is being updated with higher quality feature representations through flow-routing techniques that derive hydrography from high-resolution elevation data. However, deriving hydrography through flow-routing methods is a complex process that needs to be tailored to different geographic conditions, which can lead to varying solutions. To address this problem, this paper evaluates automated deep learning and its transferability to extract hydrography from interferometric synthetic aperture radar (IfSAR) elevation data spanning a range of geographic conditions in Alaska.

Spatial Measurements on USGS Topo Maps

United States Geological Survey topographic maps provide a variety of geospatial measurement tools thanks to their built-in georeferencing capabilities. Georeferencing facilitates display of geographic coordinates and measurement of distances, perimeters, and areas on PDFs using the measure tools built into Adobe Reader and Adobe Acrobat. This 7-page fact sheet written by Hartwig H. Hochmair and published by the UF/IFAS School of Forest Resources and Conservation explains and showcases the use of these geospatial readout functions. It also provides background information about and assistance with estimating point coordinates on the Universal Transverse Mercator grid cast on USGS topographic maps as an optional layer.

How a New Geology and Glacial Paradigm Explains Colorado South Platte-Arkansas River Drainage Divide Topographic Map Evidence, USA

United States Geological Survey (USGS) topographic maps (available at the USGS National Map website) are used to determine development of the asymmetric South Platte River drainage basin (south of Denver) by noting low points (referred to as divide crossings) where south-oriented floodwater channels once crossed the South Platte-Arkansas River drainage divide. Twelve groups of observed divide crossings are described most of which show evidence for flood-formed diverging and converging channels including where divide crossings are cut across Thirtynine Mile volcanic field ejectamenta. A new Cenozoic geologic and glacial history paradigm requires southeast- and south-oriented floods to have flowed across what is now the east-west continental divide to reach the southeast-oriented Arkansas River valley and to have eroded the divide crossings before north-oriented South Platte River valley headward erosion captured the flow. Such floods are not consistent with accepted Cenozoic geologic and glacial history interpretations, but are consistent with new paradigm interpretations (developed to explain Missouri River drainage basin topographic map drainage system evidence) in which a thick continental ice sheet (located where large continental ice sheets are usually reported to have been) deeply eroded the underlying bedrock and caused crustal uplift to create a deep “hole” with a deep “hole” rim segment roughly following today’s Wyoming and northern Colorado east-west continental divide and then continuing eastward along what is now the Missouri-Arkansas River drainage divide. The new paradigm results in a Cenozoic geologic and glacial history in which immense south-oriented continental ice sheet meltwater floods first flowed in a south direction across the rising deep “hole” rim and were then forced by rim uplift to flow along the rim and subsequently in north directions into the developing deep “hole”.

Automated location correction and spot height generation for named summits in the coterminous United States

ABSTRACT Spot elevations published on historical U.S. Geological Survey topographic maps were established as needed to enhance information imparted by the quadrangle’s contours. In addition to other features, labels were routinely placed on mountain summits. While some elevations were established through field survey triangulation, many were computed during photogrammetric stereo-compilation. Today, Global Navigation Satellite System (GNSS) receivers have replaced expensive triangulation methods. However, since GNSS measurements require visiting the feature location, a national dataset containing high-accuracy spot elevations has not yet been created. Consequently, modern U.S. Topo maps are devoid of mountain peak or other spot elevations. Still, topographic map users continue to demand the display of spot heights. Therefore, a pilot study was conducted to evaluate the feasibility of automatically generating elevation values at named U.S. summits using available elevation data. The devised method uses an uphill stepping technique to find the most likely highest point in subsequently higher-resolution elevation models. Resulting elevation values are compared to other published sources. Results from 196 summits indicate that values derived from lidar are generally higher, whereas those populated from the one-third arc-second USGS Seamless 3DEP elevation dataset are generally lower. A thorough understanding of these relationships require the evaluation of more points.

Automated summit spot height generation for modern topographic mapping

Abstract. Spot elevations published on historical U.S. Geological Survey (USGS) topographic maps were established as needed to enhance information imparted by the quadrangle’s contours. In addition to a number of other features like road intersections, section corners, lakes and water wells, labels were routinely placed on mountain peaks, local high points, passes and saddles. While some elevations were established through field survey triangulation, many were “dropped” during photogrammetric stereocompilation. The average accuracy, which varied with method, terrain, and map vintage, was on the order of ± 10 feet. Today, inexpensive consumer devices deliver comparable accuracy using the Global Positioning System (GPS). Professional equipment can achieve even better accuracies. These methods have replaced the field- and funding-intensive traditional triangulation methods in the geographic sciences. However, since GPS measurements require visiting the feature location, a national dataset containing high-accuracy spot elevations has not yet been created. Consequently, modern US Topo maps are devoid of mountain peak or other spot elevations. Additionally, no government agency is currently mandated to create one. Still, US Topo map users continue to demand the display of spot heights, particularly mountain peak elevations. As lidar data are collected through the 3D Elevation Program (3DEP), the source for higher accuracy spot heights increases. Therefore, a pilot study was conducted to evaluate the feasibility of automatically generating elevation values at summits named in the Board on Geographic Name’s Geographic Names Information System (GNIS), using 3DEP data. As the 3DEP incorporates more lidar data, these values should become increasingly more accurate. The first step in the automation process involved “snapping” GNIS summits to the highest and most accurate nearby point (pixel) in the 3DEP one-third arc second seamless dataset (1/3 a-s). Development of this step was completed in a prior study and will be reported elsewhere. In Step 2, a process similar to Step 1 was implemented to identify the highest pixel in the 3DEP one-meter dataset (1-m) from an area defined by a one-pixel buffer around the snapped point pixel. After finding the highest pixel in the 1-m dataset, the same process for Step 2 was repeated on the lidar point cloud dataset (LPC). In the latter case, the highest lidar point within the buffered area was chosen where there were more than one, which was the case in most instances. In the case that the next higher resolution data were unavailable, the summit spot elevation was set to the value obtained from the lower resolution layer. Hence, if 1-m data were unavailable, but LPC data did exist, the program in its current state was unable use the higher resolution data. This is due to the fact that the 1-m data are employed to drive the software to a more precise location of the summit. Without that intermediate step the method fails to identify the correct LPC point. Work in progress is aimed at overcoming this disadvantage, which is particularly important until the time when the 1-m dataset is populated for the entire nation. The National Geospatial Program estimates the completion of lidar acquisition by 2023. Processing from acquisition to update of the 1-m dataset, including quality assessment and revisions, currently requires between two and three years. Resulting elevation values are compared to those published by Peakbagger.com and TOPOZONE.com. Preliminary results from 40 summits indicate that values derived from lidar are generally higher, whereas those populated from the 1/3 arc-second are generally lower. A thorough understanding of these relationships will require the evaluation of more points.

Map Archive Mining: Visual-Analytical Approaches to Explore Large Historical Map Collections.

Historical maps are unique sources of retrospective geographical information. Recently, several map archives containing map series covering large spatial and temporal extents have been systematically scanned and made available to the public. The geographical information contained in such data archives makes it possible to extend geospatial analysis retrospectively beyond the era of digital cartography. However, given the large data volumes of such archives (e.g., more than 200,000 map sheets in the United States Geological Survey topographic map archive) and the low graphical quality of older, manually-produced map sheets, the process to extract geographical information from these map archives needs to be automated to the highest degree possible. To understand the potential challenges (e.g., salient map characteristics and data quality variations) in automating large-scale information extraction tasks for map archives, it is useful to efficiently assess spatio-temporal coverage, approximate map content, and spatial accuracy of georeferenced map sheets at different map scales. Such preliminary analytical steps are often neglected or ignored in the map processing literature but represent critical phases that lay the foundation for any subsequent computational processes including recognition. Exemplified for the United States Geological Survey topographic map and the Sanborn fire insurance map archives, we demonstrate how such preliminary analyses can be systematically conducted using traditional analytical and cartographic techniques, as well as visual-analytical data mining tools originating from machine learning and data science.

Which are the highest peaks in the US Arctic? Fodar settles the debate

Abstract. Though an outstanding achievement for their time, the United States Geological Survey (USGS) topographic maps of the eastern Alaskan Arctic nonetheless contain significant errors, and in this paper we address one of them. Specifically, USGS maps of different scale made in the late 1950s alternate between Mt. Chamberlin and Mt. Isto as the tallest peak in the US Arctic. Given that many of the peaks here are close in height and covered with glaciers, recent climate change may also have changed their height and their order. We resolved these questions using fodar, a new airborne photogrammetric technique that utilizes structure-from-motion (SfM) software and requires no ground control, and validated it using GPS measurements on the peaks as well as airborne lidar. Here we show that Mt. Chamberlin is currently the third tallest peak and that the order and elevations of the five tallest mountains in the US Arctic are Mt. Isto (2735.6 m), Mt. Hubley (2717.6 m), Mt. Chamberlin (2712.3 m), Mt. Michelson (2698.1 m), and an unnamed peak (2694.9 m); these heights are relative to the NAVD88 GEOID12A vertical datum. We find that it is indeed plausible that this ranking has changed over time and may continue to change as summit glaciers continue to shrink, though Mt. Isto will remain the highest under current climate trends. Mt. Isto is also over 100 m taller than the highest peak in Arctic Canada, making it the highest peak in the North American Arctic. Fodar elevations compared to within a few centimeters of our ground-based GPS measurements of the peaks made a few days later and our complete validation assessment indicates a measurement uncertainty of better than ±20 cm (95 % RMSE). By analyzing time series of fodar maps, we were able to detect topographic change on the centimeter level on these steep slopes, indicating that fodar can be used to measure mountain snow packs for water resource availability or avalanche danger, glacier volume change, and slope subsidence, as well as many other applications of benefit to society. Compared to lidar, the current state-of-the-art airborne topographic mapping, we found this SfM technique as accurate, more useful scientifically, and significantly less expensive, suggesting that fodar is a disruptive innovation that will enjoy widespread usage in the future.

Alaska tidewater glacier terminus positions, 1948–2012

A significant portion of the world's glacier ice drains through tidewater outlets, though much remains unknown about the response to recent climate change of tidewater glaciers. We present a 64 year record of length change for 50 Alaska tidewater glaciers. We use U.S. Geological Survey topographic maps to provide a base length for glaciers before 1970. Using all available cloud‐free Landsat images, we manually digitize calving front outlines for each glacier between 1972 and 2012, resulting in a total of more than 10,000 outlines. Tidewater glacier lengths vary seasonally; focusing on the 36 glaciers terminating in tidewater throughout the study period, we find a mean (± standard deviation) seasonal variation of 60± 85 m a−1. We use these oscillations to determine the significance of interannual changes in glacier length. All 36 glaciers underwent at least one period (≥1 year) of significant advance or retreat; 28 glaciers underwent at least one period of both significant advance and retreat. Over the entire period 1948–2012, 24 of these glaciers retreated a total (± uncertainty) of 107.95±0.29 km, 11 advanced a total of 7.71±0.20, and one (Chenega Glacier) did not change significantly. Retreats and advances are highly variable in time; several glaciers underwent rapid, short‐term retreats of a few years duration. These retreats occurred after large changes in summer sea surface temperature anomalies; further study is needed to determine what triggered these retreats. No coherent regional behavior signal is apparent in the length record, although two subregions show a coherence similar to recent observations in Greenland.

21st-century increase in glacier mass loss in the Wrangell Mountains, Alaska, USA, from airborne laser altimetry and satellite stereo imagery

AbstractAlaskan glaciers are among the largest regional contributors to sea-level rise in the latter half of the 20th century. Earlier studies have documented extensive and accelerated ice wastage in most regions of Alaska. Here we study five decades of mass loss on high-elevation, land-terminating glaciers of the Wrangell Mountains (~ 4900 km2) in central Alaska based on airborne center-line laser altimetry data from 2000 and 2007, a digital elevation model (DEM) from ASTER and SPOT5, and US Geological Survey topographic maps from 1957. The regional mass-balance estimates derived from center-line laser altimetry profiles using two regional extrapolation techniques agree well with that from DEM differencing. Repeat altimetry measurements reveal accelerated mass loss over the Wrangell Mountains, with the regional mass-balance rate evolving from –0.07 ± 0.19 m w.e. a–1 during 1957–2000 to –0.24 ± 0.16 m w.e. a–1 during 2000–07. Nabesna, the largest glacier in this region (˜1056 km2), lost mass four times faster during 2000–07 than during 1957–2000. Although accelerated, the mass change over this region is slower than in other glacierized regions of Alaska, particularly those with tidewater glaciers. Together, our laser altimetry and satellite DEM analyses demonstrate increased wastage of these glaciers during the last 50 years.

Geospatial Analysis of Depressional Wetlands near Peace River Watershed Phosphate Mines, Florida, USA

Mining throughout the Floridan aquifer system alters groundwater conditions. Determination of magnitude and extent of groundwater alterations in this regional karst aquifer system is difficult due to preferential flow paths, blasting, and widespread aquifer withdrawals and injections. Additionally, hydrologic models do not reveal subsurface impacts to the biotic environment, such as destruction of vegetative buffers, which are essential habitat for endangered and threatened species and other wildlife, and degradation from invasive species responding to hydroperiod disruptions. Our case study uses remote sensing and a geographic information system (GIS) to evaluate near-infrared (NIR) data for natural herbaceous depressional (wet prairie) wetlands in a central Florida phosphate mining area. These wetlands characteristically are dominated by graminoids, which have lower NIR reflectance and consequent digital numbers (DNs) in remotely sensed imagery than invasive species associated with hydroperiod alterations. Maximum aquifer withdrawals of 76,457 m 3 /d were permitted in November 1977 and remain active for study-area mines, which are surrounded primarily by un-irrigated pasture. Digital color infrared aerial imagery acquired in winter 2003–04 (1-m ground resolution) was used to extract NIR values within depressional wetlands. Digital historic (mid-1950s) U.S. Geological Survey topographic quadrangle maps (1:24,000 scale) were selected to delineate baseline wetland boundaries using ESRI ArcGIS 9.2. The study area included 567 of these wetlands, totaling 1,367 hectares. High mean NIR values (≥100 DN) consistent with hydroperiod-altered wetlands dominated by invasive plant species characterized 284 wetlands (50 and 60 percent total number and area, respectively), including some >5 km from study-area mines. Only 108 wetlands (20 and 10 percent total number and area, respectively) had low mean NIR values (≤80 DN) indicative of natural wet prairies without invasive species. Shallow ditches were not a statistically significant factor influencing NIR signatures in these wetlands. The spatial distribution of wetlands with high NIR DNs was inconsistent with conical groundwater drawdown predicted by groundwater models but suggests more linear fluid movement via subsurface preferential NW-SE flow paths.

Extent and Channel Morphology of Unmapped Headwater Stream Segments of the Quabbin Watershed, Massachusetts1

Brooks, Robert T. and Elizabeth A. Colburn, 2011. Extent and Channel Morphology of Unmapped Headwater Stream Segments of the Quabbin Watershed, Massachusetts. Journal of the American Water Resources Association (JAWRA) 47(1):158-168. DOI: 10.1111/j.1752-1688.2010.00499.x Abstract: Effective regulatory protection and management of headwater resources depend on consistent and accurate identification and delineation of stream occurrence. Published maps and digital resources fail to represent the true occurrence and extent of headwater streams. This study assessed the accuracy of mapped origins of “blue-line” streams depicted on U.S. Geological Survey topographic maps, and, if present, the morphological characteristics of unmapped stream segments. We identified 170 mapped stream origins on the Quabbin Reservoir watershed, Massachusetts. Of 30 mapped stream origins, we identified and examined 26 unmapped stream segments above 25, with an average length of 502 m. Twenty unmapped tributaries occurred on 10 of the 26 unmapped segments, with an average length of 127 m. Wetland reaches occurred more frequently and were larger on unmapped than on mapped stream segments. A significant and complex stream network occurs above most mapped stream origins. For the Quabbin watershed, we estimate that there are 85.8 km of unmapped stream upgradient of 314.5 km of mapped streams. Reliance on mapped stream networks for regulatory standards allows for the potential disturbance or even destruction of the unmapped stream resources. Jurisdictional regulations and guidelines should be revised so that the occurrence of streams should require field validation.

Hyperspectral Analysis of Columbia Spotted Frog Habitat

Abstract: Wildlife managers increasingly are using remotely sensed imagery to improve habitat delineations and sampling strategies. Advances in remote sensing technology, such as hyperspectral imagery, provide more information than previously was available with multispectral sensors. We evaluated accuracy of high‐resolution hyperspectral image classifications to identify wetlands and wetland habitat features important for Columbia spotted frogs (Rana luteiventris) and compared the results to multispectral image classification and United States Geological Survey topographic maps. The study area spanned 3 lake basins in the Salmon River Mountains, Idaho, USA. Hyperspectral data were collected with an airborne sensor on 30 June 2002 and on 8 July 2006. A 12‐year comprehensive ground survey of the study area for Columbia spotted frog reproduction served as validation for image classifications. Hyperspectral image classification accuracy of wetlands was high, with a producer's accuracy of 96% (44 wetlands) correctly classified with the 2002 data and 89% (41 wetlands) correctly classified with the 2006 data. We applied habitat‐based rules to delineate breeding habitat from other wetlands, and successfully predicted 74% (14 wetlands) of known breeding wetlands for the Columbia spotted frog. Emergent sedge microhabitat classification showed promise for directly predicting Columbia spotted frog egg mass locations within a wetland by correctly identifying 72% (23 of 32) of known locations. Our study indicates hyperspectral imagery can be an effective tool for mapping spotted frog breeding habitat in the selected mountain basins. We conclude that this technique has potential for improving site selection for inventory and monitoring programs conducted across similar wetland habitat and can be a useful tool for delineating wildlife habitats.

Hyperspectral Analysis of Columbia Spotted Frog Habitat

Wildlife managers increasingly are using remotely sensed imagery to improve habitat delineations and sampling strategies. Advances in remote sensing technology, such as hyperspectral imagery, provide more information than previously was available with multispectral sensors. We evaluated accuracy of high-resolution hyperspectral image classifications to identify wetlands and wetland habitat features important for Columbia spotted frogs (Rana luteiventris) and compared the results to multispectral image classification and United States Geological Survey topographic maps. The study area spanned 3 lake basins in the Salmon River Mountains, Idaho, USA. Hyperspectral data were collected with an airborne sensor on 30 June 2002 and on 8 July 2006. A 12-year comprehensive ground survey of the study area for Columbia spotted frog reproduction served as validation for image classifications. Hyperspectral image classification accuracy of wetlands was high, with a producer's accuracy of 96% (44 wetlands) correctly classified with the 2002 data and 89% (41 wetlands) correctly classified with the 2006 data. We applied habitat-based rules to delineate breeding habitat from other wetlands, and successfully predicted 74% (14 wetlands) of known breeding wetlands for the Columbia spotted frog. Emergent sedge microhabitat classification showed promise for directly predicting Columbia spotted frog egg mass locations within a wetland by correctly identifying 72% (23 of 32) of known locations. Our study indicates hyperspectral imagery can be an effective tool for mapping spotted frog breeding habitat in the selected mountain basins. We conclude that this technique has potential for improving site selection for inventory and monitoring programs conducted across similar wetland habitat and can be a useful tool for delineating wildlife habitats.

Modern Erosion Rates and Loss of Coastal Features and Sites, Beaufort Sea Coastline, Alaska

This study presents modern erosion rate measurements based upon vertical aerial photography captured in 1955, 1979, and 2002 for a 100 km segment of the Beaufort Sea coastline. Annual erosion rates from 1955 to 2002 averaged 5.6 m a-1. However, mean erosion rates increased from 5.0 m a-1 in 1955–79 to 6.2 m a-1 in 1979–2002. Furthermore, from the first period to the second, erosion rates increased at 60% (598) of the 992 sites analyzed, decreased at 31% (307), and changed less than ± 30 cm at 9% (87). Historical observations and quantitative studies over the past 175 years allowed us to place our erosion rate measurements into a longer-term context. Several of the coastal features along this stretch of coastline received Western place names during the Dease and Simpson expedition in 1837, and the majority of those features had been lost by the early 1900s as a result of coastline erosion, suggesting that erosion has been active over at least the historical record. Incorporation of historical and modern observations also allowed us to detect the loss of both cultural and historical sites and modern infrastructure. U.S. Geological Survey topographic maps reveal a number of known cultural and historical sites, as well as sites with modern infrastructure constructed as recently as the 1950s, that had disappeared by the early 2000s as a result of coastal erosion. We were also able to identify sites that are currently being threatened by an encroaching coastline. Our modern erosion rate measurements can potentially be used to predict when a historical site or modern infrastructure will be affected if such erosion rates persist.

USING GIS TO DETERMINE THE EFFECTS OF TWO COMMON TAPHONOMIC BIASES ON VERTEBRATE FOSSIL ASSEMBLAGES

Abstract Paleoecological statistics calculated for fossil assemblages are often influenced by biases that may act at more than one scale. These biases are difficult to detect without an overarching method that can accommodate the spatial attributes of fossil specimens. Geographic Information Systems (GIS) can be used to manage, analyze, synthesize, and archive large amounts of paleontological data within a geospatial framework. This study uses GIS to test the impact of variation in sampling area and rate of sediment accumulation on two commonly employed paleoecological statistics—species richness and relative body size. Mammals of the Willwood Formation (early Eocene) from the central part of the Bighorn Basin, Wyoming, USA, serve as a test case. Fossil localities outlined on U.S. Geological Survey topographic maps were digitized and GIS was used to estimate their two-dimensional areas. This metric was compared directly with species richness to test sampling area bias and also was used to calculate the de...

Glacial-Geologic Evidence for Decreased Precipitation During The Little Ice Age in The Brooks Range, Alaska

We mapped Little Ice Age (LIA) moraines in the Brooks Range to estimate former equilibrium-line altitudes (ELAs), and combined this information with available proxy temperature estimates to infer precipitation trends because little is known about precipitation changes associated with centennial-scale climate variability in the Arctic during the late Holocene. The Brooks Range, northern Alaska (68°N), hosts hundreds of extant glaciers that exhibit geomorphic evidence for multiple fluctuations in ice extent during the past millennium. Our lichenometric age estimates for LIA moraines in the forefields of five cirque glaciers in the Sagavinerktok River valley and Oolah Valley suggest two intervals of LIA moraine formation centered around a.d. 1250 and 1650. The outermost LIA moraine was mapped on aerial photographs for 114 relatively large (1.2 ± 0.5 km) and geographically simple glaciers along a 700-km-long transect following the range crest. At their maximum extent during the LIA, these glaciers were an average of 0.2 ± 0.1 km longer than the ice margins shown on most recent U.S. Geological Survey topographic maps (1956 and 1972). The ELA was estimated using an accumulation-area-ratio method and GIS analysis. The reconstructed ELAs needed to maintain an equilibrium length for the LIA glaciers were an average of 51 ± 29 m lower than for the smaller glacier sizes of the mid 20th century. This small ELA lowering during the LIA is less than would be expected from available proxy temperature estimates from elsewhere in Alaska that indicate warm-season temperature reductions of about 1 °C. To explain this discrepancy, we suggest that precipitation decreased during the LIA. A prolonged southern displacement of the Arctic front might explain the drier conditions in the Brooks Range during the LIA.

Glacial-Geologic Evidence for Decreased Precipitation During The Little Ice Age in The Brooks Range, Alaska

We mapped Little Ice Age (LIA) moraines in the Brooks Range to estimate former equilibrium-line altitudes (ELAs), and combined this information with available proxy temperature estimates to infer precipitation trends because little is known about precipitation changes associated with centennial-scale climate variability in the Arctic during the late Holocene. The Brooks Range, northern Alaska (68°N), hosts hundreds of extant glaciers that exhibit geomorphic evidence for multiple fluctuations in ice extent during the past millennium. Our lichenometric age estimates for LIA moraines in the forefields of five cirque glaciers in the Sagavinerktok River valley and Oolah Valley suggest two intervals of LIA moraine formation centered around a.d. 1250 and 1650. The outermost LIA moraine was mapped on aerial photographs for 114 relatively large (1.2 ± 0.5 km) and geographically simple glaciers along a 700-km-long transect following the range crest. At their maximum extent during the LIA, these glaciers were an average of 0.2 ± 0.1 km longer than the ice margins shown on most recent U.S. Geological Survey topographic maps (1956 and 1972). The ELA was estimated using an accumulation-area-ratio method and GIS analysis. The reconstructed ELAs needed to maintain an equilibrium length for the LIA glaciers were an average of 51 ± 29 m lower than for the smaller glacier sizes of the mid 20th century. This small ELA lowering during the LIA is less than would be expected from available proxy temperature estimates from elsewhere in Alaska that indicate warm-season temperature reductions of about 1 °C. To explain this discrepancy, we suggest that precipitation decreased during the LIA. A prolonged southern displacement of the Arctic front might explain the drier conditions in the Brooks Range during the LIA.

A Temporal Study of Urban Development for the Municipality of Anchorage, Alaska

A land use/land cover database was produced for a portion of the Municipality of Anchorage, Alaska to document the temporal and spatial extent of urbanization to assist in the analysis of changes in impervious cover and water quality. Data were derived from black and white and color infrared aerial photography, and satellite imagery from the early 1970's to 2000 in roughly ten‐year increments. Aerial photographs and satellite data were manually interpreted to identify and map land use/land cover classes which were then entered into a geographic information system, attributed, and georeferenced to a U.S. Geological Survey topographic map base. The spatial extent of the study was 31,117 hectares. In the early 1970's, approximately 7,356 hectares (24%) of the study area were mapped as urban developed. During the 30‐year analysis period, the largest increase in urban development occurred between the late 1970's and early 1980's when urban area increased to 12,263 hectares (39%). Between 1980 and 1990, and 1990 and 2000, urban area increased to 12,762 hectares (41%), and 13,980 hectares (45%) respectively. Most development occurred in forested or tall shrub areas, although some also occurred in wetlands. Between 1970 and 2000, close to 1,300 hectares of wetlands were lost due to development. Contrary to this, the amount of lake and pond area increased slightly from 261 hectares in 1973 to 334 hectares in 1980, and reduced to 310 hectares by 2000. The increase was primarily due to the filling in of gravel pits with spring melt water.

Extraction of bajadas from digital elevation models and satellite imagery

A methodology was designed for the extraction of bajadas from the 15 min US Geological Survey digital elevation models and Landsat Thematic Mapper imagery. The method was demonstrated for the Death Valley-California where progressive eastward tilting has enabled the west-side fans to coalesce and form bajadas. First, the drainage that crossed the uplands and the bajadas was extracted from the DEM. The drainage pixels were successively grown by checking the surrounding pixels on the basis of their gradient. It was concluded that for gradient in the interval [2°,11°] the upslope bajadas border was segmented. In order to eliminate the drainage pixels that belonged to the uplands, the drainage pixels were subtracted. Then, the isolated small 8-connected foreground pixels were identified and subtracted too. Finally, region growing was performed again to the remaining pixels with the same growing criterion. Isolated 8-connected background pixels, representing almost flat regions inside bajadas, were identified and merged to the segmented pixels. At the end, by taking into account the spectral response in the satellite image, the downslope border of bajadas was segmented. The extracted polygon was in agreement with the information depicted on (a) the US Geological Survey topographic map of scale 1:100,000 and (b) the satellite image and (c) the polygon classified manually by a photo-interpreter.

Exploratory space-time analysis of reported dengue cases during an outbreak in Florida, Puerto Rico, 1991-1992.

The spatial and temporal distributions of dengue cases reported during a 1991-1992 outbreak in Florida, Puerto Rico (population = 8,689), were studied by using a Geographic Information System. A total of 377 dengue cases were identified from a laboratory-based dengue surveillance system and georeferenced by their residential addresses on digital zoning and U.S. Geological Survey topographic maps. Weekly case maps were generated for the period between June and December 1991, when 94.2% of the dengue cases were reported. The temporal evolution of the epidemic was rapid, affecting a wide geographic area within seven weeks of the first reported cases of the season. Dengue cases were reported in 217 houses; of these 56 (25.8%) had between two and six reported cases. K-function analysis was used to characterize the spatial clustering patterns for all reported dengue cases (laboratory-positive and indeterminate) and laboratory-positive cases alone, while the Barton and David and Knox tests were used to characterize spatio-temporal attributes of dengue cases reported during the 1991-1992 outbreak. For both sets of data significant case clustering was identified within individual households over short periods of time (three days or less), but in general, the cases had spatial pattern characteristics much like the population pattern as a whole. The rapid temporal and spatial progress of the disease within the community suggests that control measures should be applied to the entire municipality, rather than to the areas immediately surrounding houses of reported cases. The potential for incorporating Geographic Information System technologies into a dengue surveillance system and the limitations of using surveillance data for spatial studies are discussed.