Critical U.S.- Mexico Borderland Watershed Analysis, Twin Cities Area Of Nogales, Arizona And Nogales, Sonora

By Laura Margaret Brady, Floyd Gray, Mario Castaneda, Mark Bultman, and Karen Sue Bolm

US Geological Survey and Arizona Department of Environmental Quality

Abstract
The United States - Mexico border area faces the challenge of integrating aspects of its binational physical boundaries to form a unified or, at least, compatible natural resource management plan. Specified geospatial components such as stream drainages, mineral occurrences, vegetation, wildlife, and land-use can be analyzed in terms of their overlapping impacts upon one another. Watersheds have been utilized as a basic unit in resource analysis because they contain components that are interrelated and can be viewed as a single interactive ecological system. In developing and analyzing critical regional natural resource databases, the Environmental Protection Agency (EPA) and other federal and non-governmental agencies have adopted a “watershed by watershed” approach to dealing with such complicated issues as ecosystem health, natural resource use, urban growth, and pollutant transport within hydrologic systems.
These watersheds can facilitate the delineation of both large scale and locally important hydrologic systems and urban management parameters necessary for sustainable, diversified land-use. The twin border cities area of Nogales, Sonora and Nogales, Arizona, provide the ideal setting to demonstrate the utility and application of a complete, cross-border, GIS-based, watershed analysis in the characterization of a wide range of natural resource as well as urban features and their interactions. In addition to the delineation of a unified, cross-border watershed, the database contains sewer/water line locations and status, well locations, geology, hydrology, topography, soils, geomorphology, and vegetation data, as well as remotely sensed imagery. The project goal is to develop a GIS database that is widely accessible to the general public, researchers, and the local land management community with a broad range of application and utility.


Introduction

The twin border cities area of Nogales, Sonora and Nogales, Arizona, provide the ideal setting to demonstrate the utility and application of a complete, cross-border, GIS-based watershed analysis of the characterization of a wide range of natural resource and urban-industrial components and their common interactions. The twin border cities area, referred to as “Ambos Nogales”, typifies the condition of a number of border areas with its assortment of water issues, including surface- and groundwater contamination, inadequate water supplies inequitably distributed, flash floods, and endangered riparian habitat (Ingram and others, 1994). Population growth, nearing 250,000 (Sanchez and Lara, 1992), exacerbates already depleted ground water sources. Poor waste containment infrastructures threaten existing and future water resources as well as create health-related concerns based on environmental exposures encompassing physical, chemical, and biological agents. With these and other socio-economic issues occurring across the United States - Mexico border area, the region faces the challenge of integrating aspects of its binational physical and political boundaries to form a unified or, at least, compatible natural resource database upon which to execute a viable international management plan. The objective of this paper is to provide an overview of the newly generated watershed-based framework and key natural resource components, attributed features, and analytical tools database developed for a GIS-platform. In some cases a transnational, compatible database needed to be created, e.g. topographical base, hydrological and geological maps. Each section focuses on issues and solutions associated with the gathering and integrating of transnational geospatial data and preliminary results from the databases are presented.
The construction of the database is accomplished in three phases: (1) GIS database compilation and development; (2) creation of an integrative platform, available in Arc/Info and ArcView, that provides a means for new input of scientific research and enables integrated watershed management; and (3) presentation of the database as a means of continuing education, community outreach, and distribution to decision makers and the general public. This study is the result of an informal data-sharing partnership between the U.S. Geological Survey, Arizona Department of Environmental Quality (ADEQ), The University of Arizona’s Udall Center and others. The goals of the partnership are to combine and focus the resources of our investigative teams, thereby to eliminate unnecessary overlap in the execution of this study, and to make the database as widely available as possible through the combined publication and web-accessibility of the government and university community. Some possible release mechanisms for the data are a deliverable spatial data archive CD-ROM or an interactive data-warehouse Internet site.
Major issues in the Nogales wash, including natural resource degradation, technical data acquisition and manipulation, and logistics of publication and presentation, can be addressed with an integrated watershed GIS:
1) Pollution leaks; inefficient water usage and inadequate sewage containment characterize the urban infrastructure of this semi-arid region. Bacteriological contamination, seepage of nitrates, heavy metals and volatile organic compounds (VOCs) from industrial sites into soils and underground water supply have been identified in the waters of both cities (Varady and Mack, 1995). Surface and groundwater hydrological models to determine runoff and sub-surface flow paths could help to pinpoint sources of pollution.
2) Lack of digital data in a unified framework; the data concerning watershed, demographic, environmental, and infrastructure characteristics are not currently available in a contained geospatial database.
3) Current border maps are fragmented; most of the existing digital and analog copy datasets originating from either country terminate at the border, and resource policy decisions are implemented without much regard for impact on adjacent populations. Geospatial inconsistencies occur due to differences in mapping, cartography, resolution and time periods of creation. The international nature of the project increases our challenge of a regular GIS project involved with integrating databases due to national, linguistic, and physical barriers (Wright and Winckell, 1999).
4) Lack of communication; transborder data sharing between research and educational facilities and policy makers is lacking. This situation is partly due to the binational/bilingual nature of the region. A graphic digital representation of the factors influencing the transborder watershed could aid in enhancing understanding by sidestepping communication barriers and could be the first step in water quality conflict resolution.


Location of the Study Area

Figure 1: Location map depicting the study area of the twin cities of Nogales, Arizona and Nogales, Sonora.


The twin cities of Nogales lie astride the Arizona-Sonora section of the U.S. – Mexico Border (fig.1), at Latitude 31.34, Longitude –110.94, approximately 110 km south of the metropolitan area of Tucson, Arizona (fig. 1). The area lies within the Sonoran Desert, locally occurring in a dissected, foothills-like setting that contains the headwaters and main reach of Nogales Wash. Several northeast flowing tributaries, Las Canoas, Mariposa Canyon, Potrero Creek, Alamo Canyon and Pesquiera Canyon flow directly into the main reach of the Nogales Wash system. Within the study area, specified geospatial components such as stream drainage, mineral occurrences, vegetation, wildlife corridors, and land-use can be analyzed.
The average high and low maximum temperature range is 93.9 - 27.4 *F (34.4 - -2.5 *C) with only an average of 17.66 inches (44.9 centimeters) of annual precipitation, almost half of which occurs in the summer months of July and August (http://www.wrcc.dri.edu/summary/climsmaz.html). These semi-arid conditions indicate a potential water shortage as well as limit pollution dilution potential (Liverman and others, 1999).


GIS Database Development

The purpose of this section is to describe the production of unified, consistent, digital, geographically registered cross-border watershed data to accurately portray the transboundary nature of the selected hydrologic sub-region. Binational water resource priorities were identified and systematically mapped to characterize transboundary relationships. Basic geologic and hydrologic data, in digital and analog format, that would be compatible with a GIS were initial research for the Nogales study area. The physical setting of the region is herein described using a watershed approach to aid in dealing with such complicated issues as ecosystem health, natural resource use, urban growth, and pollutant transport within surface and ground-water hydrologic systems. The procedures for delineating watersheds and creating stream networks are described in the following section, along with descriptions of the created transboundary watershed.


Hypsography

One of the projects preliminary accomplishments has been to generate integrated elevation data from the U.S. and Mexican sources in order to create a watershed wide topographic surface and derived themes such as slope, aspect, shaded relief, hydrologic features, and hypsographic contours (fig. 2). Due to projection, resolution and contour interval differences between datasets, the integration across the border was difficult. Digital Elevation Models (DEM’s) were acquired and concatenated for the areas of interest; resolution of the grids is 90 meters. The USGS National Elevation Data set (NED) (http://edcnts12.cr.usgs.gov/ned/) supplied 30-meter resolution DEM’s that were resampled to a 90-meter grid cell size for that portion of the analysis that occurred on the U.S. side of the border. The Mexico DEM’s were purchased through Resource Science, Inc. (http://www.resourcescience.com). These data came packaged in 3-arc-second projection using decimal degree units and were re-projected to Universal Transverse Mercator (UTM) and also resampled to a 90-meter resolution to achieve a uniform projection and grid cell size. Topographic contour lines were generated within Arc/Info from this newly created DEM, at a contour interval of every 100 meters to portray the elevation (fig. 2). The lowest elevation in the Ambos Nogales is 1,060 meters above sea level, and the highest peak is 1,710 meters. The average elevation is approximately 1,278 meters. The topography reveals the typical mountain and intervening valley structure of the Basin and Range Province. There is a regional North-Northwest structure grain in the area, typical of this portion of SE Arizona. The basins have varying amounts of alluvial sediments transected by stream channels.

Figure 2: Elevation contours of the Nogales area at a 30-m. resolution.


Watershed Delineation

A watershed is defined as a catchment basin, which is delineated topographically and drained by a stream system; it is a hydrological unit used for the planning and management of natural resources (Brooks and others, 1997). Components of a watershed include water, topography, soils, rock, vegetation, wildlife, climate and land-use. Stream networks were created from the newly generated DEM, using the GRID module of Arc/Info and compared to lower resolution vector datasets in the area. The point at which the Nogales wash enters the Santa Cruz River was determined to be the outlet of the Nogales wash tributary system and that point was used as a pour point in the watershed’s delineation. This was done using a script developed in ArcView (see Appendix), which utilizes the 90- meter resolution DEM and the previously created flow direction and flow accumulation grids to calculate the catchment basin of that specified point. The result of this digital construction of transborder watershed delineation identifies the areas of land within Mexico that are immediately contributing to the pour point or headwaters located in the United States. Figure 3 shows the newly delineated watershed area overlain with generated stream networks and a 1:250,000-scale digital raster graphic (DRG) of the area. The Nogales Wash watershed covers 235 square kilometers, just under half of the area is in Mexico. The topography is the major factor that identifies this basin boundary. Drainage in the basin is asymmetric and clearly focused along the eastern edge.

Figure 3: Nogales Wash watershed delineation portrayed in a hill shade relief of a DEM.


Hydrography

The stream flow emanating from a watershed is the result of many integrated influences of the watersheds physical composition. The amount of water received though precipitation and how this water is distributed though space and time defines the watershed itself. The Nogales wash is a gaining stream, which receives water from the local groundwater table in addition to discharge from the border area upstream. The Santa Cruz River, which is located in southern Arizona and northern Sonora, Mexico, flows into the Gila River and then into the Lower Colorado River. There are numerous problems associated with the understanding of the Santa Cruz River system’s surface flows and groundwater/surface water relationship (US-EPA, 1996).
The Nogales Wash is a nested tributary within the Santa Cruz River drainage basin and has been identified as major sources of possible contaminant contribution to water in the downstream reach of the Santa Cruz River. Utilizing the predefined Upper Santa Cruz watershed boundary in the United States, identified by the U.S. Environmental Protection Agency in the “Surf Your Watershed” Website (US-EPA, 2000), pour points were identified and Upper Santa Cruz transborder watershed was created utilizing the Script in ArcView (see Appendix). Figure 4 depicts those stream networks generated within the GRID module of Arc/Info from the newly generated DEM, associated with the Upper Santa Cruz transborder watershed and the placement of the Nogales Wash watershed within.

Figure 4: Map showing the relationship of the Nogales Watershed to the Upper Santa Cruz River Transborder Watershed.


A Digital Elevation Model was used to model the effects of land topography in the watershed delineation by identifying the streams, flow direction and a catchment perimeter or that portion of the land that effects a given outlet (Maidmont, 1993). Surface runoff flow paths are from ridge tops. Runoff flows from watershed boundaries toward the dissecting channels in which the flow accumulates toward the downstream discharge point. The surface runoff from each segment of land surface has its own unique flow path and velocities by which it reaches any designated point. In the Nogales area, during periods of surface runoff, the two sources of stream flow are surface runoff and groundwater discharge. When precipitation stops, the surface runoff rapidly ceases.
A surface runoff flow direction map (fig. 5) was generated in Arc/Info’s GRID module, from a 90-meter resolution transborder Digital Elevation Model (DEM), for purposes of study and analysis of water quality. As depicted, surface flow runs predominantly north and northeast.

Figure 5: Flow direction map of the Nogales Watershed as depicted by a DEM in Arc/ Info’s GRID module.


Within the Nogales wash, as mentioned in the site description, are many contributing tributaries. A sub-watershed was delineated at each junction to the Nogales wash, using the script developed in ArcView (see Appendix) for further investigations of the possible contaminant sources within the area (fig. 6).

Figure 6: Sub-watersheds of the Nogales Watershed defined by streams and DEM input.


Geologic Setting

Known geologic maps produced in this portion of the US-Mexico border are country-specific. At least one effort to merge known geologic databases (Gray and others, in prep.) reveals the problems of scale and nomenclature that are initially encountered at the frontier. For a watershed study of this dimension, the scales used (1:1,000,000 and 1:2,000,000) are somewhat inadequate for a portrayal of the local geology (fig. 7). A preliminary geologic map derived from larger scale base maps and satellite imagery was prepared for this report (fig. 8). This cross-border geology was created using Thematic Mapper (Landsat 5) in a similar process as is used with aerial photography. Known points were examined for color, texture, and slope styles; extension of those areas was approximated where a level of certainty could be maintained.
The sedimentary rocks in the Nogales Wash- Upper Santa Cruz Valley are Miocene to Holocene, chiefly alluvial sand and gravel deposits of fans, valley centers, terraces, and channels. On the basis of age, stratigraphic position, and degree of consolidation, these rocks can be separated into two basin-fill1 units overlain by surficial deposits as follows; (1) lower basin-fill unit or Nogales Formation, probably lower and middle Miocene, and poorly to moderately well consolidated, (2) upper basin-fill unit, upper Miocene to lower Pleistocene, and unconsolidated to poorly consolidated, and (3) Pleistocene and Holocene surficial deposits including alluvium of stream channels, flood plains, and terraces, unconsolidated overall but locally well indurated (Houser and others, 1985; Dickinson, 1991). Significant differences in porosity indurated (Nelson, 1963; Drewes, 1971, 1972,1980; Cooper, 1973; Simons, 1974); and saturated bulk density exist between the lower and upper basin –fill units with the Nogales formation having the greater induration and the lower yields of water wells as compared to the upper basin fill (Gettings and Houser, 1997).
The lower slopes of the Atascosa Mountains on the western side of the valley are composed of Tertiary volcanic rocks. The Pajarito (US) and Cerro Pedregoso Mountains (MX), located along the international border west of Nogales, consist of Cretaceous volcanics. The San Cayetano and Patagonia Mountains are on the eastern side of the Nogales Wash and upper Santa Cruz Valley. These mountains are made up of a variety of rocks including igneous, metamorphic, volcanic, and sedimentary rocks ranging in age from Precambrian to Miocene (Drewes, 1980; Simons, 1974). The northern slopes of Cerro El Ocotilloso, lying on the northwestern flanks of Sierra El Pinito, is made up of Cretaceous granitic rocks cropping out from beneath a cover of Tertiary sediments equivalent to the Nogales Formation.
Two roughly parallel sub basins are defined in the immediate area of the Nogales twin cities by a complex interpolation of Bouguer gravity data and stratigraphic-structural geologic information: 1) the Nogales sub basin, and 2) the upper Santa Cruz sub basin (Gettings and Houser, 1997). The Nogales sub basin forms a northwest-trending approximately 5 km wide trough that begins near Pesquiera Canyon at its northwest parameter, narrows through the Nogales urban corridor and terminates in the southeast near La Calera, Sonora. The basin is estimated to reach a depth of about 700 m in the northwest and shallows to some 250 m in its southeastern edge. The upper Santa Cruz sub basin and the course of the river appear to be controlled along the Mt. Benedict fault. The sub basin underlies one of the narrower valleys in southern Arizona (8 to 12 km wide) suggesting that the basin it overlies has undergone only minor to moderate lateral extension during the late Cenozoic (Gettings and Houser, 1997). Because there is a general positive correlation between basin width, amount of extension, and basin depth, the narrow span of the valley suggests that the basin is shallow. Other indications of a shallow depth to the basin include: (1) weak Bouguer gravity anomalies of the basin (Gettings and Houser, 1997), (2) the absence of lacustrian or playa sediments in the basin fill implying that the basin was never closed, and (3) ubiquitous outcrops of the Nogales Formation in the area suggesting that basin subsidence stopped after the Formation was deposited.
Geophysical data (aeromagnetic and gravity) indicate the possibility of at least two other sub basins within the Ambos Nogales area (Consejo de Recursos Minerales, 2000) that may merit further delineation. The Mariposas sub basin is a northeast-trending anomaly that roughly follows along Las Mariposas tributary. The approximately 5 km-wide area is connected to the Nogales sub basin in the northeast and extends to the la Bellotosa -El Pajarito pueblos region in to the southwest in Mexico. The southernmost sub basin intersects the southern tip of the watershed and is a roughly rectangular feature located due east of Cerro El Ocotilloso. The tributary La Cason and Mexico highway 15 dissects the area.

Figure 7: This image depicts the discrepancy in resolution and nomenclature encountered at the international border by combining geologic maps of State of Arizona at 1:1,000,000 (Reynolds, 1988) and the country of Mexico at 1:2,000,000 (Ortega-Gutierrez and others, 1992).

Figure 8: Preliminary geologic map derived from larger –scale base maps (Consejo de Recursos Naturales no Renovables, 1967; Gettings and Houser, 1997) and Landsat 7 TM imagery in the Nogales Wash watershed.


Vegetation

According to the vegetation regions of North America (FAO-UNESCO, 1975), the Nogales area lies within number 20: Southwestern Desert Shrubs region. This is the driest range; average annual precipitation is 6-12 inches/year. It is dominated by sage, creosote bush, sagebrush and greasewood shrubs. Short grasses occur with desert cacti and some Juniper and Pinyon Pine exist in the Northern areas. The shrub land/ rangeland is hydrologically very similar to historical accounts although some cattle trampling and vegetation removal increases runoff and erosion unnaturally. ADEQ has acquired aerial photographs of the study area, that are less than 1-meter resolution, which are being used for re-vegetation analyses by the Bureau of Applied Research and Anthropology (BARA) group at The University of Arizona (fig. 9).

Figure 9: Aerial photograph taken of the central commercial downtown Ambos Nogales, portraying the border crossing which overlies the tunnel containing the Nogales Wash.


Soils

The soil acts like a large sponge to hold infiltrated water and percolation increases as soils get wetter until the point of saturation, which is rare in dry areas like Nogales, where the soil mantle has the first opportunity to intercept the precipitation and little to no groundwater recharge occurs (Saxton and Shiau, 1990). According to the FAO-UNESCO of 1975, the Nogales area lies within the characterized 12b soil region. This is dominated by Vermosols, hot and dry, dominated by plains and alluvial fans with only occasional mountain ranges. Surface runoff is generally limited to local floods caused by intense rain and is often absorbed by drier regions downstream. No digital soil maps was found at this time for the Mexican portion of the watershed area. The USDA, Natural Resource Conservation Service (NRCS), offers a couple of digital soils databases; the State Soil Geographic (STATSGO) Data Base is the only digital publication of an area within Santa Cruz County (USDA-NRCS, 1995). These data were downloaded as Arc Info coverages and unzipped and untarred from their compressed deliverable. The projection was converted from Albers, NAD27 to UTM, Zone 12, NAD83 (fig.10). .

Figure 10: Soil types of the Nogales Watershed as defined by the STATSGO database (USDA-NRCS, 1995).


Geomorphology

The basin area is approximately 24,435 hectares (60,380 acres). In arid lands, where channels absorb stream flow, runoff decreases with increasing basin size. This was described in watershed studies in nearby Tombstone, Arizona, according to a relationship described by Glymph and Horton (1969): due to the 244 square kilometer basin size, Nogales Watershed yields only 8 millimeters mean annual runoff to its outlet from precipitation. Compared to the 448.5 millimeters of precipitation received by the watershed, only 1.7% of what hits the ground ends up in the creek.
Slope of the ground surface is a factor in the overland flow process and very important hydrologically in small basins where it may be the dominant factor in determining hydrograph shape. The influence of slope is most relevant to indices of peak flow and to hydrograph shape. Average slope in the Nogales Wash watershed is approximately 3 degrees, the standard deviation from the average is approximately 3 while minimum slope = 0 and maximum slope = 29 degrees (fig.11).

Figure 11: Degrees of slope depicted in the Nogales Watershed.


The shape of the basin also affects the hydrograph characteristics of lag time, the time of rise, and the peak flow rate. Given the pear like shape of the watershed, the time of rise in a hydrograph will be later than normal allowing for a longer lag time.
Channel network patterns express the shape of a hydrograph as well. The pattern displayed by the Nogales Wash watershed depicts a slower rise in the hydrograph, but a higher peak compared to a more dispersed network. As area of watershed increases, so does the order of each watershed according to the Strahler stream ordering method (Strahler, 1957, Strahler, 1964). The Strahler stream order number increases when streams of the same order intersect. This is useful when identifying hierarchy within a stream system (fig. 12). Characteristics of streams are inherent under this type of ordering system. For example: first order streams are headwaters, which have no streams running into them. Therefore, these are most susceptible to non-point source pollution problems (Strahler, 1957). Further, 3rd and 5th order streams transport the majority of sediment that they receive because of the inherent steeper gradients and contact with faster moving streams (Naiman and others, 1992).

Figure 12: Strahler stream order depicted for the Nogales Watershed.


Drainage density is a measurement of channel length over area of basin. All of the arcs describing stream channels were listed, exported and summed in Microsoft Excel to get a total stream length for all the drainage within the basin. The total length of all streams in the watershed is 173 kilometers, which when divided by the basin area of 244 square kilometers, gives a drainage density of 0.7 km/km^2. Very low drainage density occurs where soil materials are resistant to erosion or, conversely, are very permeable and the relief is small. The hydrologic significance of drainage density results from the fact that water and sediment yields are very much influenced by the length of the watercourses per unit area. The channel networks are developed by long-term hydrology and erosion but, in turn, affect the amount and rate of watershed discharge. Melton (1957) discusses the relationship between low runoff rates and low drainage density both true in the Nogales watershed.


Land Use

Land use is described by the Land-use regions of North America (FAO-UNESCO, 1975) as #6, which contains rapidly developing residential and commercial cities, where water is available. Figure 13 is an image of Band 7 (red), Band 4 (green), and Band 2 (blue) color composite of two Landsat 7 scenes (path 35 row 38 and path 35 row 39 both acquired on 13 Nov 1999). Data is UTM projected using NAD83 and the resolution is 30 m pixels. This was combined with Landsat Band 8, 15-m. panchromatic data. Resolution is 15m pixels. All satellite imagery was purchased from USGS-Eros Data Center, Sioux Falls, SD. This image can be used to identify geologic features, such as faults, land use, and vegetation coverage.

Figure 13: A satellite image of the twin cities of Nogales.


Infrastructure

Water delivery lines, well locality, sewer lines, roads and city boundaries maps are examples of infrastructure that describe an area (see Figures 14, 15, & 16). The ADEQ and the Arizona Department of Health Services (ADHS) tested the groundwater and surface water in the 1980’s, detecting high levels of perchloroethylene (PCE) and tetrachloraethylene (TCE) in the groundwater (Varady and Mack, 1995). As a result, ADEQ began monitoring in 1990, finding higher levels of these substances on the Mexican side of the border, also detecting fecal coliform and volatile organic compounds (VOC’s). A detailed comparison of water quality data from specific wells with the potential source maps would be useful to guide future monitoring efforts, as attempts are made to link known ground water contamination to specific potential sources of that contaminant (ADEQ, 1997; IBWC, 1998).

Figure 14: Well locations and water delivery lines within the Nogales watershed(modified after Varady and Mack, 1995).



Figure 15: Sewer lines in the cities of Nogales, Arizona and Nogales, Sonora (modified after Varady and Mack, 1995).


Similar to the water delivery lines, sewer lines in the city of Nogales, Sonora are distributed according to monetary supply vs. human need; the colonias are not equipped with sewer lines at all (Varady and Mack, 1995).

Figure 16: Image depicting those colonias in Nogales, Sonora that have no access to sewer and water delivery lines (modified after Varady and others, 1995).


Due to the erosion rates incurred by the increasing construction of the colonias, the water catchment capacity has decreased at a drastic rate. This explains why pipes are collapsing and why the water quality is so poor (Varady and Mack, 1995).


Potential Sources of Contamination

In the summer of 2000, Chris Donovan, a science teacher at Desert View High School in Tucson, approached the U.S. Geological Survey’s Southwest Field Office with a request for assistance with a GIS project for which she had received a grant. She wanted the students to be involved in a real project, one for which they could see the results of their class work being used. The Southwest Mineral and Environmental Investigations project was working with contamination of water in the Nogales area. The project staff indicated a willingness to work with the class and allow them to map possible sources of contamination in Nogales. Equipped with geographic positioning systems, the students located possible sites in the area and took and analyzed water samples (fig. 17). This work combined with numerous sources previously identified by ADEQ provided a layer of data that can be analyzed to see if relationships exist between any of the possible contamination sources and sites of known water contamination.

Figure 17: Michael and Christina, students from Desert View High School, collect water samples in Nogales.


Broken sewer lines, unmonitored industrial sources, outhouses, storm runoff mixed with street refuse and septic tank fields are sources of untreated sewage (Varady and others, 1995). High levels of lead and arsenic have been found in the Potrero Creek surface water but not the groundwater.
A plume is a concentrated ‘pool’ of contaminant in water. It can migrate as a result of direction and velocity of water movement in the aquifer and due to characteristics of the contaminant itself. Trichloroethylene (TCE) had been dumped at the United Musical Instrument factory near Meadow Hills in the 1960’s. Since 1988, there has been a plan to remove a 13-acre TCE plume that was threatening the local wells and a wetlands area (fig. 18). The people of Nogales fear carcinogens in the plume that may cause lupus and cancer (ADEQ, 1997).

Figure 18: Photo taken by students from Desert View High School of the Nogales wetlands area..


Potential sources of contamination were identified in the Wellhead Protection Area of Valle Verde and were identified by the ADEQ; some sites include industrial and warehouse areas, which handle hazardous materials that could be spilled, automotive service and fueling stations and laundromats (see fig.19). Another area of concern is the Union Pacific Railroad corridor (formerly known as Southern Pacific), which links to the Mexican National Railroad at the border. This corridor was identified due to its history of spills of unknown substances and the fact that railroad operations and maintenance require use of potentially toxic materials (i.e.: diesel fuel, etc.) (ADEQ, 1997).

Figure 19: Sites of possible contamination identified by the ADEQ and by Desert View High School.


Surface and Ground Water Evaluation

One of the more significant aspects of the database collection project involves the compilation and generation of water quality information within the Nogales watershed area. Initial compilations will summarize existing data from ADEQ and successful cooperative efforts undertaken by the US and Mexico pursuant to the International Boundary and Water Commission (IBWC) Joint Report or Principal Engineers Relative to the Joint Monitoring of the Quality of the Groundwaters in the Ambos Nogales area, dated January 25th, 1996. This cooperative effort serves as a model for successful binational investigations of borderland area (IBWC, 1998). Future cooperative efforts between ADEQ, the USGS, and Mexican partners will continue to investigate aspects of organic and inorganic contaminant chemistry of the area. The geospatial database will be available to GIS software as well as spreadsheet or database software.


ElectroMagnetic Maps

An airborne electromagnetic (EM) survey was conducted in the Upper Santa Cruz River area in an attempt to define the geometry of the basin, the composition of sediments in the basin, and to map conductivity of basin sediments (including faults, water-saturated sediments); parameters were optimized to obtain maximum depth of penetration. The Nogales area of the Upper Santa Cruz EM survey was carried out by Geoterrex, Ltd., under the supervision of the USGS in March, 1998. The preliminary interpretation presented here is based on partial data released to the USGS as of 1998. The available data include Conductivity Depth Transforms (“CDTs”), of the 20-channel airborne EM data; a merged aeromagnetic map acquired during the survey; and a graphical representation of the flight line locations. Detailed descriptions of airborne EM surveys of this type are given in Wynn and Gettings, (1997) and Bultman and others, (1997). Examples of CDTs can be seen in Figure 20. The color scale indicates relative conductivity. Hotter colors are more conductive and cooler colors are less conductive or more resistive. Crystalline rock with low conductivity (for instance a granite basement rock) will generally be blue to black, while sedimentary rock saturated with slightly saline water typical of southwestern U.S. aquifers will generally show up as red to purple colors (Wynn and Gettings, 1997, Bultman and others, 1997).
Figure 20 depicts the EM interpretation along flight line 504 of the 1998 Santa Cruz EM survey. The red line approximates the flight line. The colored cross section located below the red line is the interpretation of the EM data along flight line 504. The top of the cross section displays the surface topography along the flight line. The maximum depth of the data shown on this cross section is about 700 meters. Figure 20 may display the ground water associated with the Santa Cruz River Valley as the conductive area on the eastern side of the cross section. The water table on the western side of flight line 504 may be displayed as the bright red band near the surface. An area of high conductivity exists just to the west of an abandoned instrument factory. This area of high conductivity may be associated with known TCE ground water contamination in this area. These are very preliminary interpretations and much more work needs to be done in order to access their veracity.

Figure 20: This Digital Raster Graphic (DRG) shows the EM flight line # 504.


This dataset contains 20 flight lines in the Nogales, Arizona area. Future work will be concentrated on detailed interpretation of this information in conjunction with well log records, water quality and other stratigraphic details. It would be extremely beneficial if similar flight lines in the Nogales, Sonora area could be acquired.


Conclusion

The Ambos Nogales area, spanning the International border between US and Mexico, can be characterized within a GIS as being physically situated within a 94 square mile sub-watershed of the Upper Santa Cruz River system. Approximately half of this watershed is within Sonora, Mexico. A digitally accurate Nogales area watershed delineation provides the base unit within to view existing natural resource information and urban infrastructure data as well as from a context database for future data collection. This boundary-less watershed unit is a tool for analysis and presentation of various databases that are needed to respond to the water resource situation occurring in Ambos Nogales.
This database, when completed along the guidelines indicated in this text, may enable the modeling of potential environmental changes that can be accurately portrayed and perhaps predicted using an integrative and dynamic GIS system. With the data layers that have been created as part of this project, staff is better able to study the dynamics of the Nogales Watershed and the effects of its various characteristics on the flow of water and the quality of that water. This database will be analyzed for the exposure of contaminants in disease demographic context and merged with medical data to look for correlations. The Borderlands area is a proposed setting for this type of medical examination & this Nogales Wash watershed geospatial analysis serves as a preliminary investigation. Such studies open the door to development of methodologies that can be used to address similar problems in other watersheds or to better understand how varying factors can change the dynamics within an area.


Future Research

In the course of this project, those involved have seen many directions for future efforts. One goal is the incorporation of new soils, vegetation, land-use, economic, and demographic data that are more robust, current, and of higher resolution than what now exists into an ongoing database that can be easily accessed by researchers for a minimal cost. Lines of communication with government agencies, educational institutions, and other interested parties across the border must be opened and improved. A more thorough analysis of all of the electromagnetic data available in the Nogales area may provide a better understanding of the aquifer condition and impacts upon it. Ongoing collection and collaboration are vital to the discovery of the dynamic relationships in the Ambos Nogales area.


References Cited

ADEQ, 1997, Draft final project report: Nogales/Santa Cruz County Wellhead Protection Program, 103 p.

Brooks, K.N., Ffolliett, P. F., Gregerson, H. M., and Thames, J. L., 1997, Hydrology and the management of watersheds, (2nd Edition): Ames, Iowa, Iowa State University Press; p. 438- 439.

Bultman, M. W., Gettings, M. E. and Wynn, Jeff, 1999, An interpretation of the 1997 Airborne ElectroMagnetic (AEM) survey, Fort Huachuca vicinity, Cochise County, Arizona, Open-File Report 99-7-B.

Consejo de Recursos Naturales no Renovables, 1967, MAPA GEOLOGICO de la Parte Septentrional del Estado de Sonora, sheet 7: Exploracion de Yacimientos de minerals metalicos con la Coloboracion de las NACIONES UNIDAS; scale: 1:100,000.

Cooper, J.R., 1973, Geologic map of the Twin Buttes quadrangle, southwest of Tucson, Pima County, Arizona: U.S. Geological Survey Miscellaneous Investigations Series Map I-745, scale 1:48,000.

Dickinson, W.R., 1991, Tectonic setting of faulted Tertiary strata associated with the Catalina core complex in southern Arizona: Geological Society of America Special Paper 264, 106 p.

Drewes, Harald, 1971, Geologic map of the Mount Wrightson quadrangle, southeast of Tucson, Santa Cruz and Pima Counties, Arizona: U.S. Geological Survey Miscellaneous Investigations Series Map I-614, scale 1:48,000.

-----------, 1972, Cenozoic rocks of the Santa Rita Mountains, Southeast of Tucson: U.S. Geological Survey Professional Paper 746, 66p.

------------, 1980, Tectonic map of southeast Arizona: U.S. Geological Survey Miscellaneous Investigations Series Map I-1109, scale 1:125,000.

FAO-UNESCO, 1975, Soil map of the world; Volume II, North America: Paris, United Nations Educational, Scientific, and Cultural Organization, 220 p.

Gettings, M.E., and Houser, B.B., 1997, Basin Geology of the Upper Santa Cruz Valley, Pima and Santa Cruz Counties, Southeastern Arizona: U.S. Geological Survey Open-File Report 97-676.

Gylmph, L.M., and Horton, H.N., 1969, Land treatment in agricultural watershed hydrology research, in Moore, W.L., and Morgan, C.W., eds., Effects of watershed change on streamflow: American Society of Civil Engineers, pp. 44-68.

Houser, B.B., Richter, D.H., and Shafiqullah, M., 1985, Geologic map of the Safford quadrangle, Graham County, Arizona: U.S. Geological Survey, Miscellaneous Investigations Series Map I-617, scale 1:48,000.

IBWC (International Boundary and Water Commission), 1998, Binational Nogales Wash United States/Mexico Groundwater Monitoring Program Interim Report, May 1998, 61 p.

Ingram, Helen, Milich, Lenard, and Varady, R.G., 1994, Managing Transboundary Resources: Lessons from Ambos Nogales in Environment, Volume 36, Number 4, May 1994, pp. 6 - 38.

Liverman, Diana, Varady, R. G., Chavez, Octavio, and Sanchez, Roberto, 1999, Environmental issues along the U.S.-Mexico border: drivers of change and responses of citizens and institutions, Udall Center Publications, working paper 99-2, Tucson, Arizona, pp. 1-35.

Maidment, David R., 1993, GIS and hydrologic modeling, in Environmental Modeling with GIS, Editors: Michael F. Goodchild, Bradley O. Parks, Louis T. Steyaert. NY, Oxford, Oxford University Press. pp. 147- 167.

Melton, M. A., 1957, An analysis of the relations among elements of climate, surface properties and geomorphology: Office of Naval Research, Geography Branks Project NR 389-042, Technical Report 11, Columbia University.

Naiman, R. J. and others, 1992, Fundamental Elements of Ecologically Healthy Watersheds in the Pacific Northwest Coastal Ecoregion In Naiman (ed.), Watershed Management, Springer- Verlag, New York

Nelson, F.J., 1963, The geology of the Peña Blanca and Walker Canyon areas, Santa Cruz County, Arizona: Tucson, University of Arizona M.S. Thesis, 82 p.

Ortega-Gutierrez, F., Mitre-Salazar, L.M., Roldan-Quintana, J., Aranda- Gomez, J., Moran-Zentendo, D., Alaniz-Alvarez, S., and Nieto-Samaniego, A., 1992, Carta geologica de la Republica Mexicana {Geologic map of Mexico}: Consejo de Recursos Minerales de Mexico and Instituto de Geologia de la UNAM (Universidad Nacional Autonoma de Mexico), scale 1:2,000,000.

Reynolds, S.J., 1988, Geologic Map of Arizona: Arizona Geological Survey Map 26, scale 1:1,000,000.

Sanchez, R., and Lara, F., 1992, Manejo transfronterizo del agua en los dos Nogales: estrategias para Nogales, Sonora. Report to the Ford Foundation, Tijuana: El Colegia de Frontera Norte.

Saxton, K. E., and Shiau, S.Y., 1990, Surface waters of North America; Influence of land and vegetation on streamflow in Wolman, M.G., and Riggs, H.C., eds., The geology of North America Vol. 0-1, Surface water hydrology, The geologic society of America, pp. 55-80.

Simons, F.S., 1974, Geologic map and sections of the Nogales and Lochiel quadrangles, Santa Cruz County, Arizona: U.S. Geological Survey Miscellaneous Investigations Series Map I-762, scale 1:48,000, pamphlet 9 p.

Strahler, A.N, 1957, Quantitative analysis of watershed geomorphology, Transactions of the American Geophysical Union. vol. 8, no. 6, pp. 913- 920.

Strahler, A.N., 1964, Quantitative geomorphology of drainage basins and channel networks, in Chow, V.T., ed., Handbook of Applied Hydrology, section 4- II: New York, McGraw-Hill, p. 4.39-4.76.

USDA-NRCS, 1995, State soil geographic (STATSGO) data base data use information. USDA NRCS, National Soil Survey Center. Miscellaneous publication number 1492.

U.S. Environmental Protection Agency, 1996, US-Mexico border XXI program, framework document, EPA 160-R-96-003, October 1996.

U.S. Environmental Protection Agency, 2000, Upper Santa Cruz Watershed, http://www.epa.gov/surf3/hucs/15050301/

Varady, R.G., and Mack, M.D., 1995, Transboundary water resources and public health in the US-Mexico region in Journal of Environmental Health, 57(8): 8-14.

Varady, R.G., Ingram, Helen, and Milich, Lenard, 1995, The Sonoran Pimería Alta: Shared Environmental Problems and Challenges in Journal of the Southwest Volume 37, Number 1, Spring 1995, pp. 102-122.

Wright, R. and Winckell, A., 1999, Harmonizing framework and resource data across political boundaries: The Tijuana River Watershed GIS in Morain, Stan, Ed., GIS solutions in natural resource management, Onward Press, Santa Fe, New Mexico, p. 71-93.

Wynn, Jeff, and Gettings, Mark, 1997, A Preliminary Interpretation of the 1997 Airborne ElectroMagnetic (EM) Survey over Fort Huachuca, Arizona and Upper San Pedro River Basin, Open File Report 97-XXX



Appendix

ArcView Script for Watershed Delineation developed by Laura M. Brady & Scott Miller (Brady and others, in prep.).
theView = av.GetActiveDoc
theDisplay = theView.GetDisplay
theGridTheme = theView.GetActiveThemes.Get(0)
theGrid = theGridTheme.GetGrid
thePoint = theDisplay.ReturnUserPoint
mPoint = MultiPoint.Make({thePoint})
theSrcGrid = theGrid.ExtractByPoints(mPoint,Prj.MakeNull,FALSE)
theFlowDir = theView.FindTheme("Flowdir").GetGrid
theAccum = theView.FindTheme("Flowacc2").GetGrid
theWater = theFlowDir.Watershed(theSrcGrid.SnapPourPoint(theAccum,240))
' create a theme
theGTheme = GTheme.Make(theWater)
' check if output is ok
if (theWater.HasError) then
return NIL
end
' add theme to the view
theView.AddTheme(theGTheme)




Author Information

Laura Margaret Brady
Geographer/GIS Specialist
U.S. Geological Survey
Geologic Division
Southwest Field Office
520 N. Park Avenue, Suite # 355
Tucson, AZ 85719
ph. (520) 670-5510
fax. (520) 670-5571
lmbrady@swfo.arizona.edu