Search USGSSearch

SEARCH · Search USGS

Results for “Space Technology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

225 records · Page 13Linked to original sources

Chemical loading into surface water along a hydrological, biogeochemical, and land use gradient: A holistic watershed approach

Identifying the sources and impacts of organic and inorganic contaminants at the watershed scale is a complex challenge because of the multitude of processes occurring in time and space. Investigation of geochemical transformations requires a systematic evaluation of hydrologic, landscape, and anthropogenic factors. The 1160 km 2 Boulder Creek Watershed in the Colorado Front Range encompasses a gradient of geology, ecotypes, climate, and urbanization. Streamflow originates primarily as snowmelt and shows substantial annual variation. Water samples were collected along a 70-km transect during spring-runoff and base-flow conditions, and analyzed for major elements, trace elements, bulk organics, organic wastewater contaminants (OWCs), and pesticides. Major-element and trace-element concentrations were low in the headwaters, increased through the urban corridor, and had a step increase downstream from the first major wastewater treatment plant (WWTP). Boron, gadolinium, and lithium were useful inorganic tracers of anthropogenic inputs. Effluent from the WWTP accounted for as much as 75% of the flow in Boulder Creek and was the largest chemical input. Under both hydrological conditions, OWCs and pesticides were detected in Boulder Creek downstream from the WWTP outfall as well as in the headwater region, and loads of anthropogenic-derived contaminants increased as basin population density increased. This report documents a suite of potential endocrine-disrupting chemicals in a reach of stream with native fish populations showing indication of endocrine disruption.

Environmental Science & Technology

Three-dimensional geologic mapping to assess geothermal potential: Examples from Nevada and Oregon

Geologic structure plays an important role in controlling fluid flow in geothermal systems. In particular, very complex structural settings, consisting of many closely spaced and intersecting faults, host many geothermal systems. To elucidate the key geologic factors that affect fault-controlled geothermal circulation, it is critical to precisely characterize the structural and stratigraphic geometries in these complex settings. Here, we present a methodology and the results of 3D geologic analyses of two geothermal systems in the Basin and Range, USA. This methodology is a quantitative and geologically focused technique that can be used to precisely characterize geothermal areas, in a time when future geothermal growth demands increased exploration precision and efficiency. Surficial and subsurface geologic and geophysical data are synthesized in the construction of detailed 3D geologic maps of geothermal areas. Based on these 3D geologic maps, we examine several geologic attributes that control permeability development and geothermal fluid flow along faults. We use the stress state of faults and the distribution of structural discontinuities (i.e., fault intersections and fault terminations) to identify locations of upflow along faults in these geothermal systems. These results and the methodology presented herein are directly applicable to structurally controlled geothermal fields in the Basin and Range and worldwide. As development focus shifts toward blind geothermal resources, integration of precisely characterized subsurface structural information into exploration methods will be increasingly critical to continued growth in geothermal exploration and development.

Nevada, Oregon

Design and performance of a horizontal mooring for upper-ocean research

This paper describes the design and performance of a two-dimensional moored array for sampling horizontal variability in the upper ocean. The mooring was deployed in Massachusetts Bay in a water depth of 84 m for the purpose of measuring the horizontal structure of internal waves. The mooring was instrumented with three acoustic current meters (ACMs) spaced along a 170-m horizontal cable that was stretched between two subsurface buoys 20 m below the sea surface. Five 25-m-long vertical instrument strings were suspended from the horizontal cable. A bottom-mounted acoustic Doppler current profiler (ADCP) was deployed nearby to measure the current velocity throughout the water column. Pressure sensors mounted on the subsurface buoys and the vertical instrument strings were used to measure the vertical displacements of the array in response to the currents. Measurements from the ACMs and the ADCP were used to construct time-dependent, two-dimensional current fields. The current fields were used as input to a numerical model that calculated the deformation of the array with respect to the nominal zero-current configuration. Comparison of the calculated vertical offsets of the downstream subsurface buoy and downstream vertical instrument string with the pressure measurements were used to verify the numerical code. These results were then used to estimate total deformation of the array due to the passage of the internal waves. Based on the analysis of the three internal wave events with the highest measured vertical offsets, it is concluded that the geometry of the main structure (horizontal cable and anchor legs) was kept to within ±2.0 m, and the geometry of the vertical instrument strings was kept to within ±4.0 m except for one instance when the current velocity reached 0.88 m s −1 .

Massachusetts

Radio-tracking manatees from land and space: tag design, implementation, and lessons learned from long-term study

West Indian manatees (Trichechus manatus) were tracked along the Atlantic coast of Florida and Georgia (N = 83 manatees, n = 439 tag deployments, 1986-1996) and in eastern Puerto Rico (N = 8, n = 43, 1992-1996) using conventional and satellite-based radio-telemetry systems. A floating radio-tag, attached by a flexible tether to a padded belt around the base of the tail, enabled us to track manatees in saltwater environments. The tag incorporated VHF (very high frequency) and ultrasonic transmitters for field tracking and tag recovery, and an Argos satellite-monitored transmitter for remote tracking. We located each animal in the field about twice per week, received more than 60 000 good-quality Argos locations, and recovered tags in over 90% of deployments. The tag was designed to detach from the belt when entangled to prevent injury or drowning, and this often led to premature termination of tracking bouts. We had considerable success, however, in retagging belted manatees without recapture (97% of 392 retagging events). Most individuals were radio-tagged more than once (median = 3.0, maximum = 43) for a median total duration of 7.5 months (maximum = 6.8 yr). Data obtained through Argos have been valuable in addressing questions relating to long-distance movements, site fidelity, and identification of high-use areas. Fine-scale analyses of manatee habitat use and movements may require restricting the data set to the highest location quality or developing new analytical techniques to incorporate locational error. Field tracking provided useful ancillary data on life-history parameters, but sample sizes were small and survival estimates imprecise. Modification of the existing tag design to include Global Positioning System (GPS) functionality, with its finer spatial and temporal resolution, will offer new opportunities to address critical research and management problems facing this endangered species.

Marine Technology Society Journal

Three-dimensional imaging, change detection, and stability assessment during the centerline trench levee seepage experiment using terrestrial light detection and ranging technology, Twitchell Island, California, 2012

A full scale field seepage test was conducted on a north-south trending levee segment of a now bypassed old meander belt on Twitchell Island, California, to understand the effects of live and decaying root systems on levee seepage and slope stability. The field test in May 2012 was centered on a north-south trench with two segments: a shorter control segment and a longer seepage test segment. The complete length of the trench area measured 40.4 meters (m) near the levee centerline with mature trees located on the waterside and landside of the levee flanks. The levee was instrumented with piezometers and tensiometers to measure positive and negative porewater pressures across the levee after the trench was flooded with water and held at a constant hydraulic head during the seepage test—the results from this component of the experiment are not discussed in this report. We collected more than one billion three-dimensional light detection and ranging (lidar) data points before, during, and after the centerline seepage test to assess centimeter-scale stability of the two trees and the levee crown. During the seepage test, the waterside tree toppled (rotated 20.7 degrees) into the water. The landside tree rotated away from the levee by 5 centimeters (cm) at a height of 2 m on the tree. The paved surface of the levee crown had three regions that showed subsidence on the waterside of the trench—discussed as the northern, central, and southern features. The northern feature is an elongate region that subsided 2.1 cm over an area with an average width of 1.35 m that extends 15.8 m parallel to the trench from the northern end of the trench to just north of the trench midpoint, and is associated with a crack 1 cm in height that formed during the seepage test on the trench wall. The central subsidence feature is a semicircular region on the waterside of the trench that subsided by as much as 6.2 cm over an area 3.4 m wide and 11.2 m long. The southern feature is an elongate region that has a maximum subsidence of 3.5 cm over an area 0.75 m wide and 8.1 m long and is associated with a number of small fractures in the pavement that are predominately north-south-trending and parallel to the trench. We determined that there was no significant motion of the levee flank during the last week of the seepage test. We also determined biomorphic parameters for the landside tree, such as the 3D positioning on the levee, tree height, levee parallel/perpendicular cross sectional area, and canopy centroid. These biomorphic parameters were requested to support a University of California Berkeley team studying seepage and stability on the levee. A gridded, 2-cm bare-earth digital elevation model of the levee crown and the landside levee flank from the final terrestrial lidar (T-Lidar) survey provided detailed topographic data for future assessment. Because the T-Lidar was not integrated into the project design, other than an initial courtesy dataset to help characterize the levee surface, our ability to contribute to the overall science goals of the seepage test was limited. Therefore, our analysis focused on developing data collection and processing methodology necessary to align ultra high-resolution T-Lidar data (with an average spot spacing 2–3 millimeters on the levee crown) from several instrument setup locations to detect, measure, and characterize dynamic centimeter-scale deformation and surface changes during the seepage test.

California

Methods of practice and guidelines for using survey-grade global navigation satellite systems (GNSS) to establish vertical datum in the United States Geological Survey

Geodetic surveys have evolved through the years to the use of survey-grade (centimeter level) global positioning to perpetuate and post-process vertical datum. The U.S. Geological Survey (USGS) uses Global Navigation Satellite Systems (GNSS) technology to monitor natural hazards, ensure geospatial control for climate and land use change, and gather data necessary for investigative studies related to water, the environment, energy, and ecosystems. Vertical datum is fundamental to a variety of these integrated earth sciences. Essentially GNSS surveys provide a three-dimensional position x, y, and z as a function of the North American Datum of 1983 ellipsoid and the most current hybrid geoid model. A GNSS survey may be approached with post-processed positioning for static observations related to a single point or network, or involve real-time corrections to provide positioning "on-the-fly." Field equipment required to facilitate GNSS surveys range from a single receiver, with a power source for static positioning, to an additional receiver or network communicated by radio or cellular for real-time positioning. A real-time approach in its most common form may be described as a roving receiver augmented by a single-base station receiver, known as a single-base real-time (RT) survey. More efficient real-time methods involving a Real-Time Network (RTN) permit the use of only one roving receiver that is augmented to a network of fixed receivers commonly known as Continually Operating Reference Stations (CORS). A post-processed approach in its most common form involves static data collection at a single point. Data are most commonly post-processed through a universally accepted utility maintained by the National Geodetic Survey (NGS), known as the Online Position User Service (OPUS). More complex post-processed methods involve static observations among a network of additional receivers collecting static data at known benchmarks. Both classifications provide users flexibility regarding efficiency and quality of data collection. Quality assurance of survey-grade global positioning is often overlooked or not understood and perceived uncertainties can be misleading. GNSS users can benefit from a blueprint of data collection standards used to ensure consistency among USGS mission areas. A classification of GNSS survey qualities provide the user with the ability to choose from the highest quality survey used to establish objective points with low uncertainties, identified as a Level I, to a GNSS survey for general topographic control without quality assurance, identified as a Level IV. A Level I survey is strictly limited to post-processed methods, whereas Level II, Level III, and Level IV surveys integrate variations of a RT approach. Among these classifications, techniques involving blunder checks and redundancy are important, and planning that involves the assessment of the overall satellite configuration, as well as terrestrial and space weather, are necessary to ensure an efficient and quality campaign. Although quality indicators and uncertainties are identified in post-processed methods using CORS, the accuracy of a GNSS survey is most effectively expressed as a comparison to a local benchmark that has a high degree of confidence. Real-time and post-processed methods should incorporate these "trusted" benchmarks as a check during any campaign. Global positioning surveys are expected to change rapidly in the future. The expansion of continuously operating reference stations, combined with newly available satellite signals, and enhancements to the conterminous geoid, are all sufficient indicators for substantial growth in real-time positioning and quality thereof.

Techniques and Methods

Modified GIC estimation using 3-D Earth conductivity

Geomagnetically induced currents (GICs) are quasi-direct current (DC) electric currents that flow in technological conductors during geomagnetic storms. Extreme GICs are hazardous to man-made infrastructure. GICs enter and exit the technological systems, such as the electric power grid, at grounding points, and their magnitudes depend on the currents that flow underground. They are, therefore, a function of the Earth's electrical conductivity, represented at ground level as Earth impedances, as well as the resistance parameters of the power network. Traditional GIC estimation practices are based on Earth impedances obtained from laterally homogeneous or piecewise layered-Earth models. We refer to these methods, collectively, as the 1-D approximation. However, GIC hazard mitigation can be improved with more accurate GIC modeling that takes the spatially heterogeneous Earth's conductivity into account. Here, we propose a modified approximation for GIC estimation that is very similar to the 1-D approximation but is instead derived from empirical 3-D Earth impedances. Our formulation sets up the computation of static, frequency-dependent power line telluric response functions, which, once computed, may be considered part of the power grid system model. These response functions may then be used for historical scenario analysis of GIC hazards and for simplified real-time, albeit approximate, GIC estimation in a power grid. This modest modification to the simpler local field formulation approach avoids real-time integration of geoelectric fields along power lines while taking the realistic 3-D Earth into account in a rigorous manner. Once implemented, the method provides a power grid operator with the benefits of convenience and computational speed for a first look real-time operational GIC hazard assessment. We estimate that the proposed modified 3-D GIC modeling approach produces GIC values that are well within 50% of those obtained with the full-scale power line integration of spatially variable geoelectric fields, for storms comparable in scale to the 2003 Halloween storm, all geological structures, and power lines located in the contiguous United States and other low- to middle-latitude regions.

Space Weather

Special topic—Unoccupied aircraft systems

Introduction Unoccupied aircraft systems (UAS) increasingly support volcano monitoring and eruption response activities in the United States and abroad (James and others, 2020). Advances in UAS platforms and miniaturization of sensors over the past decade have expanded the use of this technology for a wide range of applications within volcanology (Jordan, 2019; James and others, 2020). UAS can greatly enhance existing ground-, aerial-, and satellite-based observation and in situ monitoring networks at volcanoes by providing new avenues for data collection in terms of access, resolution, and timing. UAS can collect data in difficult and hazardous environments, reducing risk to occupied aircraft and (or) ground crews; support the generation of dense time series of data through frequent, low-cost, high-resolution surveys; and provide real-time, on-demand measurements at volcanic systems for indicators such as gas, thermal output, and topographic change without the need to wait for contracted aerial flight services or satellite orbit intervals. During the 2018 response to the Kīlauea eruption on the Island of Hawaiʻi, UAS were used extensively and successfully to monitor, track, investigate, and (or) warn of ongoing volcanic activity (fig. L1; Neal and others, 2019). Throughout the eruption, the UAS team was able to provide data products rapidly to emergency managers for situational awareness and to scientists for quantitative hazard assessment (Diefenbach and others, 2018). Over the course of 4 months, more than 1,200 UAS missions were flown and yielded critical data that included (1) live video to emergency operations centers in Hilo and Honolulu for situational awareness; (2) gas emission rates, compositions, and concentrations; (3) repeat nadir videos over sections of the lava channel to support measurements of lava effusion rate; (4) oblique videos for hazards assessment and outreach; and (5) photogrammetry surveys to create very high-resolution topographic models and orthophoto mosaics (Diefenbach and others, 2018). In coming years, the U.S. Geological Survey (USGS) Volcano Hazards Program (VHP) plans to expand its fleet of UAS, associated sensors, and remote pilots to enhance volcano monitoring and response capabilities. Currently (2023), USGS operational capabilities are restricted to small class UAS (sUAS; less than [<] 55 pounds) that are limited in range, payload capacity, and flight duration. Additionally, USGS-piloted platforms are restricted to the U.S. Department of the Interior Office of Aviation Services approved fleet, which includes a limited number of small and medium multi-rotor aircraft and vertical take-off and landing fixed-wing aircraft ( https://www.doi.gov/aviation/uas/fleet ). Each type of platform has advantages and disadvantages. Small rotor-wing quadcopters are fast to deploy, can be carried in a backpack, and are highly maneuverable, but are typically only equipped with a small camera and have a minimal flight range. Medium rotor-wing hexacopters can carry larger payloads (< 20 kilograms [kg]) and varied sensors, but, with the drawback of minimal flight time (<30 minutes), they typically have similar range capabilities to their smaller counterparts and are not as easily deployable. Fixed-wing platforms provide relatively long endurance (<60 minutes) and range and, with the vertical take-off and landing capabilities, can launch and land in relatively small spaces; however, they have less maneuverability and hovering capability than the rotor-wing platforms. Although the 2018 Kīlauea response showed the benefit of the current UAS fleet, all platforms have limited range [<10 kilometers (km)], such that operators must be stationed relatively close to the region of interest. To expand UAS monitoring capabilities, VHP staff have been working closely with industry partners and the National Aeronautics and Space Administration to develop a next-generation UAS for volcano monitoring (Kern and others, 2020). This ruggedized, mid-range (>20 km), multiparametric (gas and photogrammetry) UAS has been developed to meet volcano monitoring needs, particularly at less accessible, more dangerous stratovolcanoes. It is expected in the coming years that additional UAS platforms with new and smaller sensors will expand our capabilities to meet the Nation’s volcano monitoring objectives.

Scientific Investigations Report

Swatara Creek basin of southeastern Pennsylvania: An evaluation of its hydrologic system

Local concentrations of population in the Swatara Creek basin of Pennsylvania find it necessary to store, transport, and treat water because local supplies are either deficient or have been contaminated by disposal of wastes in upstream areas. Water in the basin is available for the deficient areas and for dilution of the coal-mine drainage in the northern parts and the sewage wastes in the southern parts. Swatara Creek drains 576 square miles just east of Harrisburg, Pa., and is the largest tributary to the Susquehanna River from the north side below Harrisburg. It rises in the southern Pocono Mountains and flows southwestward across the Lebanon Plateau. On an average day Swatara Creek discharges more than 630 million gallons into the Susquehanna River at Middletown, Pa. In a year this amounts to about 23 inches of water over the entire basin and is the residual from an average annual precipitation of 45.5 inches. During an average year the flow in Swatara Creek from the upper third of the basin above Harper Tavern is always greater than 1,300 mgd (million gallons per day) for at least 15 days and is always greater than 25 mgd for at least 350 days. The daily streamflow from the basin averages 1.1 mgd per sq mi, but yields from different areas range from 0.97 to 1.22 mgd per sq mi. These variations are caused chiefly by differences in precipitation and land cover. The area of lowest yield is in the valleys west of Tremont, and the highest yields are in the Upper and Lower Little Swatara Creek subbasins. At high and medium stages the chemical character of the water in the streams is suitable for public and private supplies. At lower stages, defending on the areas and the amounts of contamination by coal-mine drainage and sewage pollution, the natural flow may require some treatment. At low stages the chemical characteristics of the natural flow not affected by man is almost identical with that of the ground water in the area drained by the stream. In general, the total dissolved solids range from about 25 to 400 parts per million and the hardness is as much as about 300 parts per million. The ground-water increment to the base flow of Swatara Creek averages about 240 mgd, or about 8.8 inches annually, for the basin. Generally, ground-water supplies in amounts of less than 0.5 mgd can be developed south of Blue Mountain. Supplies of several million gallons per day have been developed for industrial use from the permeable limestones in the south-central part of the basin. More intensive investigation in other parts of the basin would indicate areas where supplies of more than 0.5 mgd could be developed from properly spaced wells. The chemical character of water from wells depends largely on the host rock. In highly soluble rocks water contains large amount of dissolved solids; in more resistant rocks concentrations are lower. The chemical character of unpolluted ground water generally reflects the composition of the more readily soluble minerals in the local geologic environment. Areas contaminated by septic- tank effluent may have above normal amounts of nitrate and detergent products. Except where polluted, most ground water is suitable for public and industrial uses without extensive treatment. Sites for storage of surface water exist in the part of the basin lying in the valley and ridge area. As much as 30 to 40 percent of the annual flow could be impounded for release as low-flow augmentation for dilution of mine drainage and other wastes in the basin. Low sediment yields of supplying drainage areas would ensure a long life expectancy of reservoirs at these sites. Overbank flooding of the main stem of the Swatara Creek and its tributaries has occurred many times in the past. However, it has not been a hazard because urban development has not encroached on the flood plain. An inundation map of the August 1933 flood provides a basis that urban planners may use to avoid future damage. As water in the Swatara Creek moves downstream to the Susquehanna River, the flow is influenced consecutively by a large annual rainfall on the northern valley and ridge area, the wastes of surface and subsurface coal-mining activities, and less annual rainfall on the part of the basin lying in the Lebanon Plateau area; the flow is supplemented and further influenced by many tributaries and by the industrial and domestic wastes that are carried by these secondary streams. The annual precipitation ranges from 52 inches at the east edge and 49 inches at the west edge of the mountainous part of the basin to about 41 inches at the southwestern part at Middletown. The rainfall generally is adequate during the growing season to mature the crops. The mean annual temperature at Lebanon is about 52&deg;F, and the growing season is about 180 days. In this report the basin has been divided into eight hydrologic zones, leased on runoff, natural use of water, and chemical character of water. Four zones lie in the valley and ridge area, three lie in the Lebanon Plateau area, and one lies in the highland along the southeastern basin boundary. In each of the zones the hydrologic characteristics are virtually the same, but they may be completely different from those in adjacent zones. The boundaries of the zones generally coincide with boundaries between geologic formations, and the areas in each zone include rocks of similar influence on water. Streams in zone 4 at the northeast edge of the plateau have the highest average surface runoff from 1.2 to 1.1 mgd per sq mi whereas those in zone 2 at the northwest edge of the valley and ridge area have the lowest, about 1.0 mgd. Streams in zone 8, along the southeast edge of the basin, have the largest sustained low-flow yield, about 0.26 to 0.19 mgd per sq mi; those in zone 5 overlying the Martinsburg Shale east of Harrisburg have the smallest sustained low-flow yields, 0.03 to 0.01 mgd. Streams in the limestone area of zone 7 have the greatest range in low-flow yields in any one zone from 0.60 to 0 mgd per sq mi. Low-flow yields in zones 1 through 4 range from 0.13 to 0.03 mgd per sq mi. Surface flows from zones 1 and 2 are generally acidic and contain high concentrations of sulfate, iron, and total dissolved solids especially where contaminated with mine wastes. Surface flows from zones 3 and 4 are dilute, slightly alkaline, and suitable for public water supplies. Surface flows from zones 5, 6, and 7 are alkaline and contain moderate concentrations of dissolved solids with waters of highest hardness occurring in zone 7. Surface flows from zone 8 are dilute to moderately mineralized and are relatively high in silica concentration. Nitrate concentrations are high in surf Fee flows below sewage outfalls and in ground water contaminated by septic tank effluent and industrial wastes. Average annual sediment yields of 550 to 650 tons per square mile are characteristic of zones 1 and 2 where strip mining has destroyed the forest cover and coal culm is carried into the streams. From agricultural lands on the Martinsburg Shale in zones 5 and 6, annual sediment yields range from 300 to 350 tons per square mile; but from agricultural lands on the siliceous rocks in zone 8 and zones 3 and 4 in the valley and ridge area, the sediment yield ranges from 200 to 250 tons annually per square mile. Lowest annual sediment yields in the basin are in the forested areas of siliceous rocks in zones 2, 3, 4, and 5, and in the sinkhole topography of the limestones in zone 7 where the yield ranges from 30 to 35 tons and 50 to 60 tons per square mile, respectively. The amount of ground water that can be developed in the basin is dependent on the ability of the underlying rocks to yield water to wells. More than 300 gpm (gallons per minute) can be obtained from wells in alluvial materials in the valley bottoms and in some of the limestones where large solution channels and fractures are penetrated by the wells. From 50 to 300 gpm can be obtained from wells in loosely cemented sandstones and in fractured limestones. From 10 to 50 gpm can be developed from wells in the shales and harder sandstones. The most dense rocks will yield from 1 to 10 gpm from fractures and crevices. Most wells yield water from the upper 350 feet of the formation, for this part contains the most fractures or solution channels. Studies show that the velocity at which a contaminant will move downstream in the basin is related to the discharge of the stream at the time. At a stream discharge of about 400 mgd at Pine Grove, a contaminant in Swatara Creek would require about 40 hours to move from Pine Grove to Middletown. As a result of dispersion and dilution, the maximum concentration of the contaminant at Middletown would be less than 20 percent the concentration at Pine Grove under these conditions. An evaluation of the availability of water in the basin indicates that about I,239 mgd enters as precipitation, 630 mgd leaves as streamflow, 580 mgd is evaporated and transpired, and 56 mgd is diverted for use by man. Not all the diversions for man's use are lost to the basin, as about 27 mgd is returned as sewage for reuse. About one-fourth of the waste water is returned to the ground and the remainder to stream drainageways. Of that diverted by man, 11.6 mgd is used for public supply and 44.4 mgd for industrial and private supplies. Diversions of streamflow furnish 86 percent of the public supply and 27 percent of the industrial supply, and ground-water sources yield the remainder. Municipal and private sewage treatment plants are upgrading the waste water in many places, but no provisions are being made for treatment other than natural dilution and assimilation for the 15 mgd of coal-nine drainage in the northern part of the basin. Technology for economic treatment of mine water is not available at this time, although research in this field is being done. Urbanization eastward from Harrisburg and around Lebanon has increased the population density of the basin. Densities of 500 people per square mile and water use exceeding 2.0 mgd per sq mi can be expected in the future. By the year 2000 the population of the basin may increase 60 percent; and if the per capita rate of use increases 0.5 percent per year the domestic requirements for water will be about two times the present use, or 23 mgd. Similarly, if the present 1:4 ratio of domestic use to industrial use of water continues, at least 89 mgd will be needed for industry in the future. Although an increase to twice the present use of water can be foreseen, or 112 mgd, water for the dilution and assimilation of wastes from treatment systems are not included. Providing water for dilution of wastes from treatment plants has not been a problem, but in the future the amounts needed for this purpose will be greater as the population increases. As water becomes more valuable, treatment of sewage wastes to reduce the biochemical-oxygen-demand load by at least 80 to 90 percent will be necessary to conserve water for more productive uses. As much as 100 mgd may be needed for waste dilution in the basin by year 2000. The present trends in suburban and light industrial development will probably persist in the basin. Problems arising through changes in economic value of water, conflicts in use, and alternatives in development are typical of those confronting the manager of a water-resource system.

Pennsylvania