Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Basic Data Report”

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.

At least 757 records · Page 42Linked to original sources
USGS2023

Southern (California) sea otter population status and trends at San Nicolas Island, 2020–2023

The population of southern sea otters ( Enhydra lutris nereis ) at San Nicolas Island, California, has been monitored annually since the translocation of 140 southern sea otters to the island was completed in 1990. Monitoring efforts have varied in frequency and type across years. In 2017, the U.S. Navy and the U.S. Fish and Wildlife Service initiated a southern sea otter monitoring and research plan to determine the effects of military readiness activities on the growth or decline of the southern sea otter population at San Nicolas Island. The southern sea otter is the only subspecies of sea otter in California (hereafter, “sea otter"). The monitoring program, at its basic level, includes seasonal surveys of population abundance, distribution, and foraging activity. From 2020 to 2023, we measured a 10-percent per annum increase in population abundance (95-percent confidence interval =0–20 percent), with 146 total individuals as of April 2023. Coinciding with the recent population growth, the sea otter distribution, which previously tended to concentrate on the island’s west end during 2003–2006 before shifting toward more use in the north and south sides during 2017–2019, appears to have shifted again during 2020–2023 to concentrate at the island’s east end. Forage data were collected between February 2020 and April 2023. There was a total of 773 forage dives in 60 forage bouts, with most of the identified prey on successful dives (n=401) recorded as sea urchins (66 percent), followed by bivalves (15 percent), snails (12 percent), and crabs (5.2 percent). Two lobsters and three abalone also were identified among the sea otter prey. Estimates of energy intake rates averaged 14.0 kilocalories per minute (95-percent confidence interval =10.8–17.2 kilocalories per minute). Monitoring data from the past two decades indicate that sea otters at San Nicolas Island have maintained a steady pattern of energy intake and population growth characteristic of a robust population, including a sixfold growth between 2000 and 2023. There was no conclusive evidence of density-dependent effects based on these patterns; however, estimates of energy intake rates for 2020–2023 were slightly lower than previous estimates from 2017 to 2019. Additionally, subtidal monitoring results at four sites around San Nicolas Island indicated that counts of purple sea urchins ( Strongylocentrotus purpuratus ) have increased between 2003 and 2023, whereas sea otter foraging surveys completed during the same period revealed that some sea otters have shifted toward higher consumption of purple sea urchins and bivalves compared to red sea urchins ( S. fransicanus ), which generally are the preferred larger prey of sea otters. These results contribute to the understanding of population dynamics and to the conservation and planning of future monitoring and research of sea otters at San Nicolas Island.

California↗
USGS1947

Electrochemical stabilization as a means of preventing ground failure in railroads

Laboratory and field data on electrochemical stabilization of clays, by three Russian authors, are here presented in translation. Abstracts of the Russian papers were published in May 1947 issue of the Engineering News Record (pp. 100-101). There exists also a small body of literature, in German and English, dealing with the electrochemical stabilization and related subjects. Elements of the electrochemical process were patented by Casagrande in Germany, shortly before the last war. Results of the Russians and of others, including the German patent, appear to be sound and interesting accordingly. Mechanism of the electrochemical stabilization, however, appears to be surmised rather than established. Unless the mechanism of such stabilization is understood in detail, little progress may be expected in field applications of the electrochemical method. Electroosmosis, a poorly reversible coagulation of the soil colloids, and introduction of exchangeable aluminum into the clay complex have been given credit for the ground-stabilizing effects of direct electrical current. Much remains to be done, as the reader may see, in developing further the theory of the method. A critical study is indicated, in this connection, by agencies or individuals qualified and equipped for basic research in soil physics. Optimum schedules for field treatments need be ascertained with particular care, to suit any given kind of material and environment. A wide range of variation in such schedules, is most certainly to be encountered in dealing with materials as diverse in their composition and properties as are clays. Any generalization on relationships between soil, electrolytes, moisture, and current could be premature if based on the Russian work alone. Stabilization of ground is a major engineering geologic problem of national interest. Needless to say, perhaps, that failures are to be expected, in laboratory and in the field, in this as well as in any other kind of research. To minimize probabilities of such failures, it may be recommended that investigators develop the electrochemical stabilization problem not merely against the relatively narrow background of soil mechanics, but with a certain feeling for geology, mineralogy, pedology, soil physics, and soil chemistry.

Open-File Report↗
USGS2014

Integrated synoptic surveys of the hydrodynamics and water-quality distributions in two Lake Michigan rivermouth mixing zones using an autonomous underwater vehicle and a manned boat

The U.S. Geological Survey (USGS), in cooperation with the National Monitoring Network for U.S. Coastal Waters and Tributaries, launched a pilot project in 2010 to determine the value of integrated synoptic surveys of rivermouths using autonomous underwater vehicle technology in response to a call for rivermouth research, which includes study domains that envelop both the fluvial and lacustrine boundaries of the rivermouth mixing zone. The pilot project was implemented at two Lake Michigan rivermouths with largely different scales, hydrodynamics, and settings, but employing primarily the same survey techniques and methods. The Milwaukee River Estuary Area of Concern (AOC) survey included measurements in the lower 2 to 3 miles of the Milwaukee, Menomonee, and Kinnickinnic Rivers and inner and outer Milwaukee Harbor. This estuary is situated in downtown Milwaukee, Wisconsin, and is the most populated basin that flows directly into Lake Michigan. In contrast, the Manitowoc rivermouth has a relatively small harbor separating the rivermouth from Lake Michigan, and the Manitowoc River Watershed is primarily agricultural. Both the Milwaukee and Manitowoc rivermouths are unregulated and allow free exchange of water with Lake Michigan. This pilot study of the Milwaukee River Estuary and Manitowoc rivermouth using an autonomous underwater vehicle (AUV) paired with a manned survey boat resulted in high spatial and temporal resolution datasets of basic water-quality parameter distributions and hydrodynamics. The AUV performed well in these environments and was found primarily well-suited for harbor and nearshore surveys of three-dimensional water-quality distributions. Both case studies revealed that the use of a manned boat equipped with an acoustic Doppler current profiler (ADCP) and multiparameter sonde (and an optional flow-through water-quality sampling system) was the best option for riverine surveys. To ensure that the most accurate and highest resolution velocity data were collected concurrently with the AUV surveys, the pilot study used a manned boat equipped with an ADCP. Combining the AUV and manned boat datasets resulted in datasets that are essentially continuous from the fluvial through the lacustrine zones of a rivermouth. Whereas the pilot studies were completed during low flows on the tributaries, completion of surveys at higher flows using the same techniques is possible, but the use of the AUV would be limited to areas with relatively low velocities (less than 2 feet per second) such as the harbors and nearshore zones of Lake Michigan. Overall, this pilot study aimed at evaluation of AUV technology for integrated synoptic surveys of rivermouth mixing zones was successful, and the techniques and methods employed in this pilot study should be transferrable to other sites with similar success. The use of the AUV provided significant time savings compared to traditional sampling techniques. For example, the survey of outer Milwaukee Harbor using the AUV required less than 7 hours for approximately 600 profiles compared to the 150 hours it would have taken using traditional methods in a manned boat (a 95 percent reduction in man-hours). The integrated datasets resulting from the AUV and manned survey boat are of high value and present a picture of the mixing and hydrodynamics of these highly dynamic, highly variable rivermouth mixing zones from the relatively well-mixed fluvial environment through the rivermouth to the stratified lacustrine receiving body of Lake Michigan. Such datasets not only allow researchers to understand more about the physical processes occurring in these rivermouths, but they provide high spatial resolution data required for interpretation of relations between disparate point samples and calibration and validation of numerical models.

Wisconsin↗
USGS2013

Effects of Canada goose herbivory on the tidal freshwater wetlands in Anacostia Park, 2009-2011

Herbivory has played a major role in dictating vegetation abundance and species composition at Kingman Marsh in Anacostia Park, Washington, D.C., since restoration of this tidal freshwater wetland was initiated in 2000. The diverse and robust vegetative cover that developed in the first year post-reconstruction experienced significant decimation in the second year, after the protective fencing was removed, and remained suppressed throughout the five-year study period. In June 2009 a herbivory study was initiated to document the impacts of herbivory by resident and nonmigratory Canada geese (Branta canadensis) to vegetation at Kingman Marsh. Sixteen modules consisting of paired fenced plots and unfenced control plots were constructed. Eight of the modules were installed in vegetated portions of the restoration site that had been protected over time by pre-existing fencing, while the remaining eight modules were placed in portions of the site that had not been protected over time and were basically unvegetated at the start of the experiment. Exclosure fencing was sufficiently elevated from the substrate level to allow access to other herbivores such as fish and turtles, while hopefully excluding mature Canada geese. The study was designed with an initial exclosure elevation of 20 cm. This elevation was chosen based on the literature, as adequate to exclude mature Canada geese, while maximizing access to other herbivores such as fish and turtles. Repeated measures analysis of variance (ANOVA) was used to analyze the differences between paired fenced and unfenced control plots for a number of variables including total vegetative cover. Differences in total vegetative cover were not statistically significant for the baseline data collected in June 2009. By contrast, two months after the old protective fencing was removed from the initially-vegetated areas to allow Canada geese access to the unfenced control plots, total vegetative cover had declined dramatically in the initially-vegetated unfenced control plots, and differences between paired fenced and unfenced control plots were statistically significant. These differences have remained steady and significant throughout the remainder of these first three years of the study. Total vegetative cover has followed a somewhat different path in the initially-unvegetated modules, where cover in the fenced plots did not significantly exceed cover in the unfenced control plots until the August 2010 sampling event. In spite of the slow start in the initially-unvegetated modules, differences between paired fenced plots and unfenced control plots have remained significant and even increased significantly over time. This indicates that total vegetative cover in the initially-unvegetated fenced plots and unfenced control plots is continuing to diverge over time as vegetation increases in the protected plots compared to the basically unvegetated unfenced control plots. Total vegetative cover has been composed almost entirely of native species during the first three years of the study, with cover by exotics averaging less than 1% during each sampling event. Species richness did not differ significantly between fenced plots and unfenced control plots during 2009, the first year of the study. Since August 2010, species richness has remained significantly greater in the fenced plots than in the unfenced control plots. These differences have remained relatively steady over time for both the initially-vegetated and initially unvegetated modules. During the study it became apparent that our elevated fence plots were more accessible to mature geese than we had expected. Even after lowering the exclosure fencing to 15 cm in 2010 and 10 cm in 2011, we documented geese inside exclosures in both years. Nonetheless the data indicate that even at 10 cm, we have limited the numbers of mature geese entering the fenced plots, rather than totally preventing their access through low spots in the uneven substrate surface. At an exclosure elevation of 10 cm and with a soft, mucky substrate, we are assuming that non-goose herbivores such as fish and turtles still have free access to the fenced plots. Annual wildrice (Zizania aquatica), known from previous studies to be especially palatable to Canada geese, has seen the greatest impact from partial access to the fenced plots by mature geese, moving from an overwhelming dominant in the initially-vegetated plots to a minor presence there by August 2011. Interestingly, pickerelweed (Pontederia cordata), also known to be highly palatable to Canada geese, has so far shown only minor herbivory in the fenced plots. By August 2011, pickerelweed had actually increased to significantly greater cover levels in the fenced plots compared to the unfenced control plots. In conclusion, the first three years of data document that vegetation exposed to full herbivory by resident and nonmigratory Canada geese for three years in the unfenced control plots showed significantly lower total vegetative cover and species richness compared to the vegetation in the fenced plots, which experienced reduced herbivory by resident and nonmigratory Canada geese. These effects were documented for modules located in both initially-vegetated and initially-unvegetated habitats.

Maryland↗
USGS1983

Evaluating environmental and economic consequences of alternative pest management strategies: results of modeling workshops

The U.S. Environmental Protection Agency (EPA) needs a comprehensive method to evaluate the human health and environmental effects of alternative agricultural pest management strategies. This project explored the utility of Adaptive Environmental Assessment (AEA) techniques for meeting this need. The project objectives were to produce models for environmental impact analysis, improve communications, identify research needs and data requirements, and demonstrate a process for resolving conflicts. The project was structured around the construction (in an initial 2 1/2-day workshop) and examination (in a second 2 1/2-day workshop) of a simulation model of a corn agroecosystem. The model conceptualized at the first workshop simulates the effect of corn agrecosystem decisions on crop production, economic returns, and environmental indicators. The model is composed of five interacting submodels: 1) a Production Strategies submodel which makes decisions concerning tillage, planting, fertilizer and pesticide applications, and harvest; 2) a Hydrology/Chemical Transport submodel which represents soil hydrology, erosion, and concentrations of fertilizers and pesticides in the soil, runoff, surface waters, and percolation; 3) a Vegetation submodel which simulates growth of agricultural crops (corns and soybeans) and weeds; 4) a Pests submodel which calculates pest population levels and resulting crop damage; and 5) an Environmental Effects submodel which calculates indicators of potential fish kills, human health effects, and wildlife habitat. The most persistent data gaps encountered in quantifying the model were coefficients to relate environmental consequences to alternative pest management strategies. While the model developed in the project is not yet accurate enough to be used for real-world decisions about the use of pesticides on corn, it does contain the basic structure upon which such a model could be built. More importantly at this stage of development, the project has shown that very complex systems can be modeled in short periods of time and that the process of building such models increases understanding among disciplinary specialists and between diverse institutional interests. This process can be useful to EPA as the agency cooperates with other institutions to meet its responsibilities in less costly ways. Activities at the second 2 1/2-day workshop included a review of the model, incorporation of necessary corrections, simulation of policy scenarios, and examination of techniques to address remaining institutional conflicts. Participants were divided into three groups representing environmental, production or industry, and regulatory interests. Each group developed scenarios that would be most appealing to their particular interest and the scenarios were simulated by the agroecosystem computer model. Negotiators from each of the interest groups decided whether a hypothetical herbicide should be relabeled and if certain restrictions should be imposed on its use. Other participants functioned as experts and consultants on caucus teams. A solution to the hypothetical problem was successfully negotiated. Workshop participants and project staff agreed that the model and processes developed during the project should be used in training students, extension specialists, farmers, researchers, and chemical producers in collaborative problem solving methods. More productive research can be planned, and more realistic models of complex systems can be built in this way. More importantly, greater trust of decisionmakers in computer models, better understanding by technical experts about disciplines other than their own, and improved cooperation between institutional interests can be achieved. This trust, understanding, and cooperation are critical ingredients in solving problems that are too complex to be resolved by independent disciplinary activity and unilateral decision authority.

Report↗
USGS1982

Radioactivity and geochemistry of selected mineral-spring waters in the Western United States; basic data and multivariate statistical analysis

Multivariate statistical analyses were performed on data from 156 mineral-spring sites in nine Western States to analyze relationships among the various parameters measured in the spring waters. Correlation analysis and R-mode factor analysis indicate that three major factors affect water composition in the spring systems studied: (1) duration of water circulation, (2) depth of water circulation, and (3) partial pressure of carbon dioxide. An examination of factor scores indicates that several types of hydrogeologic systems were sampled. Most of the samples are (1) older water from deeper circulating systems having relatively high salinity, high temperature, and low Eh or (2) younger water from shallower circulating systems having relatively low salinity, low temperature, and high Eh. The rest of the samples are from more complex systems. Any of the systems can have a relatively high or low content of dissolved carbonate species, resulting in a low or high pH, respectively. Uranium concentrations are commonly higher in waters of relatively low temperature and high Eh, and radium concentrations are commonly higher in waters having a relatively high carbonate content (low pH) and, secondarily, relatively high salinity. Water samples were collected and (or) measurements were taken at 156 of the 171 mineral-spring sites visited. Various samples were analyzed for radium, uranium, radon, helium, and radium-228 as well as major ions and numerous trace elements. On-site measurements for physical properties including temperature, specific conductance, pH, Eh, and dissolved oxygen were made. All constituents and properties show a wide range of values. Radium concentrations range from less than 0.01 to 300 picocuries per liter; they average 1.48 picocuries per liter and have an anomaly threshold value of 171 picocuries per liter for the samples studied. Uranium concentrations range from less than 0.01 to 120 micrograms per liter and average 0.26 micrograms per liter; they have an anomaly threshold value of 48.1 micrograms per liter. Radon content ranges from less than 10 to 110,000 picocuries per liter, averages 549 picocuries per liter and has an anomaly threshold of 20,400 picocuries per liter. Helium content ranges from -1,300 to +13,000 parts per billion relative to atmospheric helium; it averages +725 parts per billion and has an anomaly threshold of 10,000 parts per billion. Radium-228 concentrations range from less than 2.0 to 33 picocuries per liter; no anomaly threshold was determined owing to the small number of samples. All of the anomaly thresholds may be somewhat high because the sampling was biased toward springs likely to be radioactive. The statistical variance in radium and uranium concentrations unaccounted for by the identified factors testifies to the complexity of some hydrogeologic systems. Unidentified factors related to geologic setting and the presence of uranium-rich rocks in the systems also affect the observed concentrations of the radioactive elements in the water. The association of anomalous radioactivity in several springs with nearby known uranium occurrences indicates that other springs having anomalous radioactivity may also be associated with uranium occurrences as yet undiscovered.

Open-File Report↗
USGS2022

Water-quality data and trends in the Rapid Creek Basin, South Dakota, 1970–2020

Surface-water-quality data in the Rapid Creek Basin in South Dakota were compiled to assess basic trends in the water quality of Rapid Creek. Spatial and temporal patterns in water quality were described for major ions, sediment, total suspended solids, nutrients, field measurements, bacteria, and select metals for the period of 1970–2020, and a water-quality trend analysis was completed for sites with enough data for selected constituents. Major ions and total suspended solids had higher median concentrations in the lower basin (downstream from the city of Rapid City) relative to the upper and middle basins. Nutrient concentrations were generally low, and increased concentrations were only detected at the sites downstream from the City of Rapid City Water Reclamation Facility. Fecal indicator bacteria ( Escherichia coli and fecal coliform) concentrations were highest downstream from the main urbanized area of Rapid City. Water-quality trends were analyzed for total dissolved solids, specific conductance, calcium, magnesium, total suspended solids, total phosphorus, dissolved phosphorus, and total Kjeldahl nitrogen for the period of 1979–2019. Concentrations for major ions and total dissolved solids typically changed by less than 15 percent. Total dissolved solids concentrations upstream from Rapid City were generally decreasing, whereas concentrations downstream were generally increasing. The flow-averaged geometric mean concentration of total dissolved solids at three sites upstream from Rapid City decreased overall by 3–5 percent, and concentrations at two sites downstream from Rapid City increased by at least 7 percent between 1979 and 2019. Trends in specific conductance in the Rapid Creek Basin were mixed with alternating increasing and decreasing trends at many of the sites between 1979 and 2014. Total suspended solids concentrations were observed to be decreasing at two sites analyzed for trends. Concentrations in total phosphorus were observed to be decreasing at every site analyzed for trends between 1989 and 2014. Significant downward trends in total Kjeldahl nitrogen were observed at two sites in the lower Rapid Creek Basin for the trend period of 1999–2019. The decreases in total suspended solids and nutrient concentrations in the Rapid Creek Basin could be related to several processes such as the implementation of a stormwater management plan in Rapid City, improvements to the water reclamation facility downstream from Rapid City, and residual climatic effects.

South Dakota↗
USGS1993

A detailed taxonomy of Upper Cretaceous and lower Tertiary Crassatellidae in the eastern United States: An example of the nature of extinction at the boundary

Current theories on the causes of extinction at the Cretaceous-Tertiary boundary have been based on previously published data; however, few workers have stopped to ask the question, 'How good is the basic data set?' To test the accuracy of the published record, a quantitative and qualitative analysis of the Crassatellidae (Mollusca, Bivalvia) of the Gulf and Mid-Atlantic Coastal Plains of the United States for the Upper Cretaceous and lower Tertiary was conducted. Thirty-eight species names and four generic names are used in publications for the Crassatellidae within the geographic and stratigraphic constraints of this analysis. Fourteen of the 38 species names are represented by statistically valid numbers of specimens and were tested by using canonical discriminant analysis. All 38 names, with the exception of 1 invalid name and 4 names for which no representative specimen could be located, were evaluated qualitatively. The results show that the published fossil record is highly inaccurate. Only 8 valid, recognizable species exist in the Crassatellidae within the limits of this study, 14 names are synonymized, and 11 names are represented by indeterminate molds or poorly preserved specimens. Three of the four genera are well founded; the fourth is based on the juvenile of another genus and therefore synonymized. This detailed taxonomic analysis of the Crassatellidae illustrates that the published fossil record is not reliable. Calculations of evolutionary and paleobiologic significance based on poorly defined, overly split fossil groups, such as the Crassatellidae, are biased in the following ways: Rates of evolution and extinction are higher, Faunal turnover at mass extinctions appears more catastrophic, Species diversity is high, Average species durations are shortened, and Geographic ranges are restricted. The data on the taxonomically standardized Crassatellidae show evolutionary rates one-quarter to one-half that of the published fossil record; faunal change at the Cretaceous-Tertiary boundary that was not catastrophic; a constant number of species on each side of the Cretaceous-Tertiary boundary; a decrease in abundance in the Tertiary; and lower species diversity, longer average species durations, and expanded geographic ranges. Similar detailed taxonomic studies need to be conducted on other groups of organisms to test the patterns illustrated for the Crassatellidae and to determine the extent and direction of the bias in the published fossil record. Answers to our questions about evolutionary change cannot be found in the literature but rather with the fossils themselves. Evolution and extinction occur within small populations of species groups, and it is only through detailed analysis of these groups that we can achieve an understanding of the causes and effects of evolution and extinction.

Professional Paper↗
USGS2013

Digital database of the Holocene tephras of the Mono-Inyo Craters, California

This digital product comprises a collection of age and isopach data for the Holocene tephras of the Mono-Inyo Craters, California. Data on the most recent eruptions from this volcanic chain are relatively comprehensive, getting less so the further back in time. For the most recent eruptions to about 1,500 years ago, tephra beds within separate eruptive sequences have been studied and isopached. Before this, from about 2,000 years ago to about 5,000 years ago, there are insufficient data for isopaching. However, one isolated tephra of about 9,000 years ago was studied and isopached in detail. Regarding ages, there are many tens of radiocarbon ages that have been obtained on the Holocene Mono-Inyo volcanic products. The vast majority of these radiocarbon dates are associated with tephras at locales that can be considered distal (basically where the primary tephra is less than several centimeters (cm) thick). These dates represent carbon that was sequestered perhaps within several hundred years of the eruption but do not represent the ages of separate eruptive pulses. There are two reasons for this. In some cases, it is clear that the dated material is not associated with the eruption products. This is the case in some lake strata where carbon is either not physically close to a given tephra layer or where an age for a tephra layer was obtained by interpolation assuming a sedimentation rate. In other cases, it is not clear that a given tephra layer represents a primary tephra; in such cases the layer could instead be redeposited. At most distal localities (beyond about 5 kilometers (km) from the chain), there was no record made of whether tephra was primary or redeposited, and at these distances where tephra is thin, it is generally redeposited during later events such as fires or thunderstorms. These age data are not appropriate for use in dating the eruptive history of the volcanic chain, and are therefore not included in the present contribution. The carbon age data in the present contribution were obtained by careful consideration of the material being collected. In the best instances, carbon was collected from new growth on plants that were probably killed by an eruption event through burning and burial. Slightly poorer data were collected from burned and buried forest duff that is renewed frequently. Finally, some dates for older Holocene tephra layers at Black Lake, Nevada, downwind of the Mono-Inyo Craters, appear to allow correlation of the layers to proximal occurrences. In cases where these poorer data were collected but yielded ages statistically indistinguishable from better data, the poorer data were included in the analysis. In the most difficult cases, usually the furthest back in time, poorer data that were nevertheless statistically indistinguishable were weighted together to generate the age estimate. There are some known Holocene eruptions from the Mono-Inyo Craters that are not included in this tabulation, as so far a tephra has not been associated with the eruptions. A good example of this is the Java blocks. The Java block eruption, from a vent underlying the northwestern corner of Negit Island in Mono Lake, expelled numerous blocks that were rafted within the lake and that are mostly deposited on the southwestern and northern lakeshore. No tephra that can be correlated to this deposit has been found, and therefore the eruption is not included in this tabulation.

California↗
USGS2014

Water quality and algal community dynamics of three deepwater lakes in Minnesota utilizing CE-QUAL-W2 models

Water quality, habitat, and fish in Minnesota lakes will potentially be facing substantial levels of stress in the coming decades primarily because of two stressors: (1) land-use change (urban and agricultural) and (2) climate change. Several regional and statewide lake modeling studies have identified the potential linkages between land-use and climate change on reductions in the volume of suitable lake habitat for coldwater fish populations. In recent years, water-resource scientists have been making the case for focused assessments and monitoring of sentinel systems to address how these stress agents change lakes over the long term. Currently in Minnesota, a large-scale effort called “Sustaining Lakes in a Changing Environment” is underway that includes a focus on monitoring basic watershed, water quality, habitat, and fish indicators of 24 Minnesota sentinel lakes across a gradient of ecoregions, depths, and nutrient levels. As part of this effort, the U.S. Geological Survey, in cooperation with the Minnesota Department of Natural Resources, developed predictive water quality models to assess water quality and habitat dynamics of three select deepwater lakes in Minnesota. The three lakes (Lake Carlos in Douglas County, Elk Lake in Clearwater County, and Trout Lake in Cook County) were assessed under recent (2010–11) meteorological conditions. The three selected lakes contain deep, coldwater habitats that remain viable during the summer months for coldwater fish species. Hydrodynamics and water-quality characteristics for each of the three lakes were simulated using the CE-QUAL-W2 model, which is a carbon-based, laterally averaged, two-dimensional water-quality model. The CE-QUAL-W2 models address the interaction between nutrient cycling, primary production, and trophic dynamics to predict responses in the distribution of temperature and oxygen in lakes. The CE-QUAL-W2 models for all three lakes successfully predicted water temperature, on the basis of the two metrics of absolute mean error and root mean square error, using measured inputs of water temperature and nutrients. One of the main calibration tools for CE-QUAL-W2 model development was the vertical profile temperature data, available for all three lakes. For all three lakes, the absolute mean error and root mean square error were less than 1.0 degree Celsius and 1.2 degrees Celsius, respectively, for the different depth ranges used for vertical profile comparisons. In Lake Carlos, simulated water temperatures compared better to measured water temperatures in the epilimnion than in the hypolimnion. The reverse was true for the other two lakes, Elk Lake and Trout Lake, where the simulated results were slightly better for the hypolimnion than the epilimnion. The model also was used to approximate the location of the thermocline throughout the simulation periods, approximately April to November, in all three lake models. Deviations between the simulated and measured water temperatures in the vertical lake profile commonly were because of an offset in the timing of thermocline shifts rather than the simulated results missing thermocline shifts altogether.

Minnesota↗
USGS2006

Kendall-Theil Robust Line (KTRLine--version 1.0)-A Visual Basic Program for Calculating and Graphing Robust Nonparametric Estimates of Linear-Regression Coefficients Between Two Continuous Variables

The Kendall-Theil Robust Line software (KTRLine-version 1.0) is a Visual Basic program that may be used with the Microsoft Windows operating system to calculate parameters for robust, nonparametric estimates of linear-regression coefficients between two continuous variables. The KTRLine software was developed by the U.S. Geological Survey, in cooperation with the Federal Highway Administration, for use in stochastic data modeling with local, regional, and national hydrologic data sets to develop planning-level estimates of potential effects of highway runoff on the quality of receiving waters. The Kendall-Theil robust line was selected because this robust nonparametric method is resistant to the effects of outliers and nonnormality in residuals that commonly characterize hydrologic data sets. The slope of the line is calculated as the median of all possible pairwise slopes between points. The intercept is calculated so that the line will run through the median of input data. A single-line model or a multisegment model may be specified. The program was developed to provide regression equations with an error component for stochastic data generation because nonparametric multisegment regression tools are not available with the software that is commonly used to develop regression models. The Kendall-Theil robust line is a median line and, therefore, may underestimate total mass, volume, or loads unless the error component or a bias correction factor is incorporated into the estimate. Regression statistics such as the median error, the median absolute deviation, the prediction error sum of squares, the root mean square error, the confidence interval for the slope, and the bias correction factor for median estimates are calculated by use of nonparametric methods. These statistics, however, may be used to formulate estimates of mass, volume, or total loads. The program is used to read a two- or three-column tab-delimited input file with variable names in the first row and data in subsequent rows. The user may choose the columns that contain the independent (X) and dependent (Y) variable. A third column, if present, may contain metadata such as the sample-collection location and date. The program screens the input files and plots the data. The KTRLine software is a graphical tool that facilitates development of regression models by use of graphs of the regression line with data, the regression residuals (with X or Y), and percentile plots of the cumulative frequency of the X variable, Y variable, and the regression residuals. The user may individually transform the independent and dependent variables to reduce heteroscedasticity and to linearize data. The program plots the data and the regression line. The program also prints model specifications and regression statistics to the screen. The user may save and print the regression results. The program can accept data sets that contain up to about 15,000 XY data points, but because the program must sort the array of all pairwise slopes, the program may be perceptibly slow with data sets that contain more than about 1,000 points.

Techniques and Methods↗
USGS1974

Seismic hazards and land-use planning

Basic earth-science data are necessary for a realistic assessment of seismic hazards and as a basis for limiting corrective land-use controls only to those areas of greatest hazard. For example, the location, character, and amount of likely displacement and activity of surface faulting can be predicted if detailed geologic maps and seismic data are available and are augmented by field studies at critical localities. Because few structures can withstand displacement of their foundations, they should be located off active fault traces, the distance varying with the character of faulting, the certainty with which fault traces are known, and the importance of the structure. Recreational activities and other nonoccupancy land uses should be considered for fault zone areas where land is under pressure for development; elsewhere, such areas should remain as open space. Two methods of predicting ground shaking effects have applications to land-use decisions: (1) Relative earthquake effects can be related to firmness of the ground and can be used in a gross way to allocate population density in the absence of more sophisticated analyses; and (2) intensity maps, based on, (a) damage from former earthquakes, or (b) a qualitative analyses of geologic units added to a design earthquake, can be helpful both for general and specific plans. Theoretical models are used with caution to predict ground motion for critical structures to be located at specific sites with unique foundation conditions. Fully adequate methods of assessing possible shaking remain to be developed. Where land-use decisions do not reflect likely ground shaking effects, stringent building codes are needed, particularly for important structures. Ground failure (landsliding, ground cracking and lurching, differential settlement, sand boils, and subsidence) commonly results from liquefaction, loss of soil strength, or compaction. Areas suspected of being most likely to fail should not be developed unless detailed site studies can demonstrate the hazard does not exist or can be overcome. Various methods can be used to reduce the high, long-term public costs that follow development of unstable ground. However, areas subject to tectonic deformation generally cannot be predicted nor can effects of such deformation be minimized. Large water waves, such as produced by tsunamis, seiches, and dam failure or overtopping, can be anticipated in many places. Their effects can be lessened by land-use regulations similar to flood-plain zoning, restrictions on location of critical structures, and appropriate warning systems. Many local, state, and federal government agencies, universities, and private consultants may be able to assist planners by advising them of pertinent data and where those data can be obtained. Interpretation of the data for an evaluation of seismic risk commonly requires a team of planners, geologists, and soil and structural engineers.

Circular↗
USGS1947

Major winter and nonwinter floods in selected basins in New York and Pennsylvania

The scientific design of flood-control works is based on an evaluation of the hydrologic factors basic to flood events, particularly how rainfall and snow runoff, soil conditions, and channel influences can combine to produce greater or lesser floods. For this purpose an analysis of the pertinent hydrologic data is needed. The methods of analysis adopted should conform as closely as possible to those already in use and must be adapted to the quality of the available information. Maximum floods in 8 basins in New York and Pennsylvania during the winter and nonwinter months were studied, a total of 21 floods. The most outstanding winter flood of record in the North Atlantic region was that of March 1936. Rainfall plus snow melt in the basins studied ranged between 3.04 and 6.87 inches, and associated volumes of direct runoff from 1.88 to 5.63 inches. Winter floods have a common characteristic in their relation to freezing temperature. The antecedent periods, representing a period of snow accumulation and frost penetration, are below freezing, and the flood itself is contemporaneous with a period of above-freezing temperatures, usually associated with rain, during which the previously accumulated snow is melted. A second common characteristic of major winter floods is their tendency to be associated with widespread causal meteorologic conditions. There was a more complete conversion of rainfall and snow melt into runoff during the winter storms studied than during the wettest nonwinter flood. Snow melt during winter floods ranged from 0.04 to 0.07 inch per degree-day above 32° F. The depth of mean areal rainfall produced by the nonwinter storms studied ranged from 3.05 to 4.96 inches. The maximum 24-hour quantity at single stations was 14 inches, which was measured during the storm of July 1935 in New York. The volume of direct runoff ranged between 1.39 and 3.41 inches. The portion of rainfall that was converted into runoff varied in accordance with the rate of antecedent base flow, expressed in second-feet per square mile, and emphasized the influence of antecedent conditions. The average volume of direct runoff during winter floods was 4.24 inches, and the average during nonwinter floods was 2.44 inches. The latter, however, were more concentrated as to time, tending to compensate for large volume of runoff in winter, so that the crest rates of direct runoff averaged 0.056 inches per hour during the winter and 0.051 inches during the nonwinter period.

New York, Pennsylvania↗
USGS1994

High-resolution single-channel seismic reflection surveys of Orange Lake and other selected sites of north central Florida

The potential fluid exchange between lakes of north central Florida and the Floridan aquifer and the process by which exchange occurs is of critical concern to the St. Johns Water Management District. High-resolution seismic tools with relatively new digital technology were utilized in collecting geophysical data from Orange, Kingsley, Lowry and Magnolia Lakes, and the Drayton Island area of St. Johns River. The data collected shows the application of these techniques in understanding the formation of individual lakes, thus aiding in the management of these natural resources by identifying breaches or areas where the confining units are thin or absent between the water bodies and the Floridan aquifer. Orange Lake, the primary focus of the study, is a shallow flooded plain that was formed essentially as an erosional depression in the clayey Hawthorn formation. The primary karstic features identified in the lake were cover subsidence, cover collapse and buried sinkholes structures in various sizes and stages of development. Orange Lake was divided into three areas southeast, southwest, and north-central. Karst features within the southeast area of Orange Lake are mostly cover subsidence sinkholes and associated features. Many of the subsidence features found are grouped together to form larger composite sinkholes, some greater than 400 m in diameter. The size of these composite sinkholes and the number of buried subsidence sinkholes distinguish the southeast area from the others. The potential of lake waters leaking to the aquifer in the southeast area is probably controlled by the permeability of the cover sediments or by fractures that penetrate the lake floor. The lake bottom and subsurface of the north-central areas are relatively subsidence sinkholes that have no cover sediments overlying them, implying that the sinks have been actively subsiding with some seepage into the aquifer from the lake in this area due to the possible presence of the active subsidence and faulting. The largest and most important features in the lake are the collapse sinkholes found along the southwestern shore that provide conduits for exchange between the lake and subsurface aquifer. There are two basic differences between the southwest and other areas of the lake: (1) the features found towards the central part of the lake are smaller in scale (1to 10 m across) and tend to be singular structures compare to the southwest area where features combined to form larger sinkholes (>400 m), and; (2) the southwest area is the only site where collapse dolines were identified. These dolines are located along the southwestern shoreline adjacent to Heagy-Burry Park. The comparison of seismic profiles from the several other selected lake and river sites to the Orange Lake profiles showed that other study areas were constructed of one or two large subsidences or a combination of sinkholes to form one large sinkhole. Aside from the difference in scale the basic characteristics of the subsidence sinkholes were similar.

Open-File Report↗
USGS2007

Estimated water use and availability in the East Narragansett Bay study area, Rhode Island, 1995-99

Water availability became a concern in Rhode Island during a drought in 1999, and further investigation was needed to assess the current demands on the hydrologic system from withdrawals during periods of little to no precipitation. The low ground-water levels and streamflows measured in Rhode Island prompted initiation of a series of studies on water use and availability in each major drainage area in Rhode Island for the period 1995–99. The investigation of the East Narragansett Bay area is the last of these studies. The East Narragansett Bay study area (130.9 square miles) includes small sections of the Ten Mile and Westport River Basins in Rhode Island. The area was divided into three regions (islands and contiguous land areas separated by the bay) within each of which the freshwater water use and availability were assessed. During the study period from 1995 through 1999, three major public water suppliers in the study area withdrew 7.601 million gallons per day (Mgal/d) from ground-water and surface-water reservoirs. The estimated water withdrawals by minor public water suppliers during the study period were 0.063 Mgal/d. Total self-supply domestic, industrial, commercial, and agricultural withdrawals from the study area averaged 1.891 Mgal/d. Total water use in the study area averaged 16.48 Mgal/d, of which about 8.750 Mgal/d was imported from other basins. The average return flow to freshwater within the basin was 2.591 Mgal/d, which included effluent from permitted facilities and septic systems. The average return flow to saltwater (Narragansett Bay) outside of the basin was about 45.21 Mgal/d and included discharges by permitted facilities (wastewater-treatment plants and Rhode Island Pollutant Discharge Elimination Systems). The PART program, a computerized hydrographseparation application, was used for the data collected at two selected index stream-gaging stations in the East Narragansett Bay study area to determine water availability on the basis of the 75th, 50th, and 25th percentiles of the total base flow; the base flow for the 7-day, 10-year low-flow scenario; and the base flow for the Aquatic Base Flow scenario for both stations. Base flows in the study area were lowest in September for the 75th, 50th, and 25th percentiles. The safe yields determined for the surface-water reservoirs (14.10 Mgal/d) were added to the estimated available ground water (gross yield) in the Southeastern Narragansett and East Narragansett Islands regions to give the total available water. The water availability in the study area at the 50th percentile ranged from 33.18 Mgal/d in September to 94.62 Mgal/d in June, water availability for the 7-day, 10-year low-flow scenario at the 50th percentile ranged from 21.87 Mgal/d in September to 83.03 Mgal/d in June, and water availability for the Aquatic Base Flow scenario at the 50th percentile ranged from 14.10 Mgal/d in August and September to 65.48 Mgal/d in June. Because water withdrawals and use are greater during the summer than at other times of the year, water availability in June, July, August, and September was compared to water withdrawals in the three regions. For the study period, the withdrawals in July were higher than in the other summer months. For the 50th percentile, the ratios of water withdrawn to water available were close to one in August for the estimated basic and Aquatic Base Flow scenarios and in September for the estimated 7-day, 10-year low-flow scenario. For the 25th percentile, the ratios were close to one in August for the estimated basic and for the 7-day, 10-year low-flow scenario, and were close to one in July for the estimated Aquatic Base Flow scenario. A long-term water budget was calculated for the East Narragansett Bay study area to identify and assess inflows and outflows by region. The water withdrawals and return flows used in the budget were from 1995 through 1999. Total inflow and outflow were calculated separately for each region. Inflow was assumed to equal outflow; the total water budget was 292.1 Mgal/d for the study area. Precipitation and return flow were 99 and less than 1 percent of the total estimated inflow to the study area, respectively. Evapotranspiration, streamflow, and water withdrawals were 47, 49, and 3 percent of the total outflow from the study area, respectively.

Rhode Island↗
USGS1995

U.S. Geological Survey Subsidence Interest Group conference, Edwards Air Force Base, Antelope Valley, California, November 18-19, 1992; Abstracts and summary

Land subsidence, the loss of surface elevation as a result of the removal of subsurface support, affects every state in the United States. More than 17,000 mi 2 of land in the United States has been lowered by the various processes that produce land subsidence with annual costs from resulting flooding and structural damage that exceed $125 million. It is estimated that an additional $400 million is spent nationwide in attempts to control subsidence. Common causes of land subsidence include the removal of oil, gas, and water from underground reservoirs; dissolution of limestone aquifers (sinkholes); underground mining activities; drainage of organic soils; and hydrocompaction (the initial wetting of dry soils). Overdrafting of aquifers is the major cause of areally extensive land subsidence, and as ground-water pumping increases, land subsidence also will increase. Land subsidence and its effects on engineering structures have been recognized for centuries, but it was not until this century that the processes that produce land subsidence were identified and understood. In 1928, while working with field data from a test of the Dakota Sandstone aquifer, O.E. Meinzer of the U.S. Geological Survey recognized the compressibility of aquifers. Around the same time, Karl Terzaghi, a soil scientist working at Harvard University, developed the one-dimensional consolidation theory that provided a quantitative means of predicting soil compaction resulting from the drainage of compressible soils. Thus, with the recognition of the compressibility of aquifers (Meinzer), and the development of a quantitative means of predicting soil compaction as a consequence of the reduction of intergranular pore pressure (Terzaghi), the theory of aquifer-system compaction was formed. With the widespread availability of electric power in rural areas, and the advent of the deep turbine pump, ground-water withdrawals increased dramatically throughout the country in the 1940's and 1950's. Along with this unprecedented increase in pumpage, substantial amounts of land subsidence were observed in several areas of the United States, most notably in Arizona, California, and Texas. Beginning in 1955, under the direction of Joseph Poland, the Geological Survey began the "Mechanics of Aquifers Project," which focused largely on the processes that resulted in land subsidence due to the withdrawal of ground water. This research team gained international renown as they advanced the scientific understanding of aquifer mechanics and land-subsidence theory. The results of field studies by members of this research group not only verified the validity of the application of Terzaghi's consolidation theory to compressible aquifers, but they also provided definitions, methods of quantification, and confirmation of the interrelation among hydraulic head declines, aquifer-system compaction, and land subsidence. In addition to conducting pioneering research, this group also formed a "center of expertise," providing a focal point within the Geological Survey for the dissemination of technology and scientific understanding in aquifer mechanics. However, when the "Mechanics of Aquifers Project" was phased out in 1984, the focal point for technology transfer no longer existed. Interest among various state and local agencies in land subsidence has persisted, and the Geological Survey has continued to participate in a broad spectrum of cooperative and Federally funded projects in aquifer mechanics and land subsidence. These projects are designed to identify and monitor areas with the potential for land subsidence, to conduct basic research in the processes that control land subsidence and the development of earth fissures, as well as to develop new quantitative tools to predict aquifer-system deformation. In 1989 an ad hoc "Aquifer Mechanics and Subsidence Interest Group" (referred to herein as the "Subsidence Interest Group") was formed

California↗
USGS2002

Crustal structure of the coastal and marine San Francisco Bay region, California

As of the time of this writing, the San Francisco Bay region is home to about 6.8 million people, ranking fifth among population centers in the United States. Most of these people live on the coastal lands along San Francisco Bay, the Sacramento River delta, and the Pacific coast. The region straddles the tectonic boundary between the Pacific and North American Plates and is crossed by several strands of the San Andreas Fault system. These faults, which are stressed by about 4 cm of relative plate motion each year, pose an obvious seismic hazard. We have many ways to study earthquake faults. Where faults break the land surface, we may learn valuable information needed for hazard assessment, such as cumulative offset, slip rate, and earthquake history. However, many of the major faults in the region are partly submerged beneath San Francisco and Monterey Bays. Although this situation poses problems in gathering observational data for hazard assessment, bay-region waterways provide an opportunity to study faultzone structure by using marine subsurface-imaging techniques, which are easier and cheaper than equivalent studies on land. In 1993, the U.S. Geological Survey (USGS) launched a 5-year project aimed at unearthing the basic science of the submerged San Andreas strike-slip fault system in the San Francisco Bay region with its many interacting strands. Primary project goals were structural, such as to discover how the San Andreas and Hayward Faults are connected or related at depth, to learn how the complex of faults in the San Andreas stepover zone on the Golden Gate platform functions, and to locate previously unknown faults. This volume thus contains mostly structural information about the San Francisco Bay region, much of it gathered through exploratory geophysical experiments. The volume is organized “top down,” from studies in the shallowest crust to the base of the crust. The first three chapters are linked through their use of novel geophysical techniques to study earthquake effects, coseismic slip, and shallow stratigraphy. Kayen and others examine crustal structure at very high resolution and demonstrate the use of ground-penetrating-radar tomography to measure the liquefaction potential of coastal sedimentary deposits. McGann and others use microfossils from drill cores along the San Francisco-Oakland Bay Bridge to determine a more detailed late Pleistocene stratigraphy of San Francisco Bay than was previously available. Geist and Zoback use the historical record of a small local tsunami generated by the great 1906 San Francisco earthquake to model the rupture process of that earthquake. The last four chapters are dedicated to studies of fault related structure of the seismogenic crust in and around the San Andreas Fault system in the San Francisco Bay region. Jachens and others compile an aeromagnetic anomaly map from new high-resolution flights across the bay region. Some of these anomalies mark the positions of offshore faults, and others are offset by faults, providing constraints on cumulative slip. Hart and others concisely summarize the marine seismic data recorded in and around San Francisco Bay, map the coverage, and provide archival information for those interested in acquiring data. The last two chapters present the results of the seismic data that have been analyzed. Bruns and others present their analysis of high-quality intermediate-resolution (~5-km penetration) seismic-reflection data gathered over the complex San Andreas-San Gregorio Fault junction. This junction, which is thought to be where the 1906 San Francisco earthquake originated (see Geist and Zoback, this volume), contains an apparent extensional right stepover in the San Andreas Fault. Finally, Parsons and others review and summarize the results of deep-crustal seismic-reflection experiments and local-earthquake tomographic studies, including previously unpublished data, and provide additional support and discussion for already-published studies. In summary, these studies were carried out in an environment where background information on faults in the San Francisco Bay region was sought. Much of the structural information presented here comes from experiments of a style unlikely to be conducted by the USGS in the near future. Together, the chapters in this volume provide a structural framework for a major part of a complex strike-slip fault system.

California↗
USGS2011

Hydrogeology and water quality of the Floridan aquifer system and effects of Lower Floridan aquifer pumping on the Upper Floridan aquifer at Fort Stewart, Georgia

Test drilling, field investigations, and digital modeling were completed at Fort Stewart, GA, during 2009?2010, to assess the geologic, hydraulic, and water-quality characteristics of the Floridan aquifer system and evaluate the effect of Lower Floridan aquifer (LFA) pumping on the Upper Floridan aquifer (UFA). This work was performed pursuant to the Georgia Environmental Protection Division interim permitting strategy for new wells completed in the LFA that requires simulation to (1) quantify pumping-induced aquifer leakage from the UFA to LFA, and (2) identify the equivalent rate of UFA pumping that would produce the same maximum drawdown in the UFA that anticipated pumping from LFA well would induce. Field investigation activities included (1) constructing a 1,300-foot (ft) test boring and well completed in the LFA (well 33P028), (2) constructing an observation well in the UFA (well 33P029), (3) collecting drill cuttings and borehole geophysical logs, (4) collecting core samples for analysis of vertical hydraulic conductivity and porosity, (5) conducting flowmeter and packer tests in the open borehole within the UFA and LFA, (6) collecting depth-integrated water samples to assess basic ionic chemistry of various water-bearing zones, and (7) conducting aquifer tests in new LFA and UFA wells to determine hydraulic properties and assess interaquifer leakage. Using data collected at the site and in nearby areas, model simulation was used to assess the effects of LFA pumping on the UFA. Borehole-geophysical and flowmeter data indicate the LFA at Fort Stewart consists of limestone and dolomitic limestone between depths of 912 and 1,250 ft. Flowmeter data indicate the presence of three permeable zones at depth intervals of 912-947, 1,090-1,139, and 1,211?1,250 ft. LFA well 33P028 received 50 percent of the pumped volume from the uppermost permeable zone, and about 18 and 32 percent of the pumped volume from the middle and lowest permeable zones, respectively. Chemical constituent concentrations increased with depth, and water from all permeable zones contained sulfate at concentrations that exceeded the U.S. Environmental Protection Agency secondary maximum contaminant level of 250 milligrams per liter. A 72-hour aquifer test pumped LFA well 33P028 at 740 gallons per minute (gal/min), producing about 39 ft of drawdown in the pumped well and about 0.4 foot in nearby UFA well 33P029. Simulation using the U.S. Geological Survey finite-difference code MODFLOW was used to determine long-term, steady-state flow in the Floridan aquifer system, assuming the LFA well was pumped continuously at a rate of 740 gal/min. Simulated steady-state drawdown in the LFA was identical to that observed in pumped LFA well 33P028 at the end of the 72-hour test, with values larger than 1 ft extending 4.4 square miles symmetrically around the pumped well. Simulated steady-state drawdown in the UFA resulting from pumping in LFA well 33P028 exceeded 1 ft within a 1.4-square-mile circular area, and maximum drawdown in the UFA was 1.1 ft. Leakage from the UFA through the Lower Floridan confining unit contributed about 98 percent of the water to the well; lateral flow from specified-head model boundaries contributed about 2 percent. About 80 percent of the water supplied to LFA well 33P028 originated from within 1 mile of the well, and 49 percent was derived from within 0.5 mile of the well. Vertical hydraulic gradients and vertical leakage are progressively higher near the LFA pumped well which results in a correspondingly higher contribution of water from the UFA to the pumped well at distances closer to the pumped well. Simulated pumping-induced interaquifer leakage from the UFA to the LFA totaled 725 gal/min (1.04 million gallons per day), whereas simulated pumping at 205 gal/min (0.3 million gallons per day) from UFA well 33P029 produced the equivalent maximum drawdown as pumping LFA well 33P028 at 740 gal/min during the aquifer test. This equivalent pumpin

Georgia↗