SEARCH · Search USGS

Results for “Studies in Maine Geology”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68Linked to original sources

USGS2001

Reconnaissance geologic map of the Dixonville 7.5' quadrangle, Oregon

The Dixonville 7.5 minute quadrangle is situated near the edge of two major geologic and tectonic provinces the northernmost Klamath Mountains and the southeastern part of the Oregon Coast Ranges (Figure 1). Rocks of the Klamath Mountains province that lie within the study area include ultramafic, mafic, intermediate and siliceous igneous types (Diller, 1898, Ramp, 1972, Ryberg, 1984). Similar rock associations that lie to the southwest yield Late Jurassic and earliest Cretaceous radiometric ages (Dott, 1965, Saleeby, et al., 1982, Hotz, 1971, Harper and Wright, 1984). These rocks, which are part of the Western Klamath terrane (Western Jurassic belt of (Irwin, 1964), are considered to have formed within an extensive volcanic arc and rifted arc complex (Harper and Wright, 1984) that lay along western North America during the Late Jurassic (Garcia, 1979, Garcia, 1982, Saleeby, et al., 1982, Ryberg, 1984). Imbricate thrust faulting and collapse of the arc during the Nevadan orogeny, which ranged in age between about 150 to 145 Ma in the Klamath region (Coleman, 1972, Saleeby, et al., 1982, Harper and Wright, 1984) was syntectonic with, or closely followed by deposition of the volcano-lithic clastic rocks of the Myrtle Group. The Myrtle Group consists of Upper Jurassic and Lower to middle Cretaceous turbidity and mass flow deposits considered to be either arc basin and/or post-orogenic flysh basins that were syntectonic with the waning phases of arc collapse (Imlay et al., 1959, Ryberg, 1984, Garcia, 1982, Roure.and Blanchet, 1983). The intermediate and mafic igneous rocks of the Rogue arc and the pre-Nevadan sedimentary cover (the Galice Formation, (Garcia, 1979) are intruded by siliceous and intermediate plutonic rocks principally of quartz diorite and granodiorite composition (Dott, 1965, Saleeby, et al., 1982, Garcia, 1982, Harper and Wright, 1984). The plutonic rocks are locally tectonized into amphibolite, gneiss, banded gneiss and augen gneiss. Similar metamorphic rocks have yielded metamorphic ages of 165 to 150 Ma (Coleman, 1972, Hotz, 1971, Saleeby, et al., 1982, Coleman and Lanphere, 1991). The Jurassic arc rocks and sedimentary cover occur as a tectonic outlier in this region (Figure 2) as they are bound to the northwest and southeast by melange, broken formation and semi-schists of the Dothan Formation and Dothan Formation(?) that are considered part of a late Mesozoic accretion complex (Ramp, 1972, Blake, et al., 1985). The plutonism that accompanied arc formation and tectonic collapse of the arc does not intrude the structurally underlying Dothan Formation, indicating major fault displacements since the Early Cretaceous. Semischistose and schistose rocks of the accretion complex have yielded metamorphic ages of around 125-140 Ma where they have been studied to the southwest (Coleman and Lanphere, 1971, Dott, 1965, Coleman, 1972). These rocks were unroofed and unconformably overlain by marine deposits by late early Eocene time (Baldwin, 1974). The early Tertiary history of this region is controversial. The most recent interpretation is that during the Paleocene and early Eocene the convergent margin was undergoing transtension or forearc extension as suggested by the voluminous extrusion of pillow basalt and related dike complexes (Wells, et al., 1984, Snavely, 1987). This episode was followed shortly by thrust and strike-slip faulting in the late early Eocene (Ryberg, 1984). During the Eocene, the Mesozoic convergent margin association of arc, clastic basin, and accretion complex was partly unroofed and faulted against early Cenozoic rocks of the Oregon Coast Ranges (Ramp, 1972, Baldwin, 1974, Champ, 1969, Ryberg, 1984). Faults that are typical of this period of deformation include high-angle reverse faults with a very strong component of strike-slip displacement characterized by a low-angle rake of striae. Thrust and oblique-slip faults are ubiquitous in early Tertiary rocks to the northwest (Ryberg, 1984, Niem and Niem, 1990). The late Mesozoic and early Cenozoic arc and forearc rocks are unconformably overlain to the east by the late Eocene and younger, mainly continental fluvial deposits and pyroclastic flows of the Cascade arc (Peck, et al., 1964, Baldwin, 1974, Walker and MacLeod, 1991). Minor fossiliferous shallow marine sandstone is locally present. The volcanic sequence consists of a homoclinal section of about 1 to 2 kilometers of andesitic to rhyolitic flows and ash flow tuff. The section is gently east-tilted and is slightly disrupted by NE trending faults with apparent normal separation.

Oregon↗

USGS2013

Estimating irrigation water use in the humid eastern United States

Accurate accounting of irrigation water use is an important part of the U.S. Geological Survey National Water-Use Information Program and the WaterSMART initiative to help maintain sustainable water resources in the Nation. Irrigation water use in the humid eastern United States is not well characterized because of inadequate reporting and wide variability associated with climate, soils, crops, and farming practices. To better understand irrigation water use in the eastern United States, two types of predictive models were developed and compared by using metered irrigation water-use data for corn, cotton, peanut, and soybean crops in Georgia and turf farms in Rhode Island. Reliable metered irrigation data were limited to these areas. The first predictive model that was developed uses logistic regression to predict the occurrence of irrigation on the basis of antecedent climate conditions. Logistic regression equations were developed for corn, cotton, peanut, and soybean crops by using weekly irrigation water-use data from 36 metered sites in Georgia in 2009 and 2010 and turf farms in Rhode Island from 2000 to 2004. For the weeks when irrigation was predicted to take place, the irrigation water-use volume was estimated by multiplying the average metered irrigation application rate by the irrigated acreage for a given crop. The second predictive model that was developed is a crop-water-demand model that uses a daily soil water balance to estimate the water needs of a crop on a given day based on climate, soil, and plant properties. Crop-water-demand models were developed independently of reported irrigation water-use practices and relied on knowledge of plant properties that are available in the literature. Both modeling approaches require accurate accounting of irrigated area and crop type to estimate total irrigation water use. Water-use estimates from both modeling methods were compared to the metered irrigation data from Rhode Island and Georgia that were used to develop the models as well as two independent validation datasets from Georgia and Virginia that were not used in model development. Irrigation water-use estimates from the logistic regression method more closely matched mean reported irrigation rates than estimates from the crop-water-demand model when compared to the irrigation data used to develop the equations. The root mean squared errors (RMSEs) for the logistic regression estimates of mean annual irrigation ranged from 0.3 to 2.0 inches (in.) for the five crop types; RMSEs for the crop-water-demand models ranged from 1.4 to 3.9 in. However, when the models were applied and compared to the independent validation datasets from southwest Georgia from 2010, and from Virginia from 1999 to 2007, the crop-water-demand model estimates were as good as or better at predicting the mean irrigation volume than the logistic regression models for most crop types. RMSEs for logistic regression estimates of mean annual irrigation ranged from 1.0 to 7.0 in. for validation data from Georgia and from 1.8 to 4.9 in. for validation data from Virginia; RMSEs for crop-water-demand model estimates ranged from 2.1 to 5.8 in. for Georgia data and from 2.0 to 3.9 in. for Virginia data. In general, regression-based models performed better in areas that had quality daily or weekly irrigation data from which the regression equations were developed; however, the regression models were less reliable than the crop-water-demand models when applied outside the area for which they were developed. In most eastern coastal states that do not have quality irrigation data, the crop-water-demand model can be used more reliably. The development of predictive models of irrigation water use in this study was hindered by a lack of quality irrigation data. Many mid-Atlantic and New England states do not require irrigation water use to be reported. A survey of irrigation data from 14 eastern coastal states from Maine to Georgia indicated that, with the exception of the data in Georgia, irrigation data in the states that do require reporting commonly did not contain requisite ancillary information such as irrigated area or crop type, lacked precision, or were at an aggregated temporal scale making them unsuitable for use in the development of predictive models. Confidence in the reliability of either modeling method is affected by uncertainty in the reported data from which the models were developed or validated. Only through additional collection of quality data and further study can the accuracy and uncertainty of irrigation water-use estimates be improved in the humid eastern United States.

Eastern United States↗

USGS1934

Zoning of ore deposits in and adjoining the Leadville District, Colorado

The ore deposits in the heart of the Leadville district , Colorado , are divisible, as shown by Emmons, Irving, and Loughlin,2 into two main classes on the basis of mineral composition: "contact-metamorphic" silicates accompanied by iron oxides, which adjoin an obscure igneous stock, and mixed sulphides accompanied by manganosiderite and quartz gangue, which cut the "contact-metamorphic" deposits and form a broad zone extending for 2 1/2 miles to the west and for shorter distances in other directions. A detailed consideration of the mineralogy of the mixed sulphide ore deposits clearly shows that they in turn are divisible into groups, referred to as hotter and intermediate mesothermal deposits . These groups, however, because of complications in structural control, are not ideally arranged from the standpoint of zoning .Emmons, Irving, and Loughlin noted that barite was a conspicuous gangue mineral in certain outlying deposits on the west slope of the Mosquito Range. The junior author during 1928-1931 had an opportunity to study the outlying deposits over an area of twenty square miles, which includes mines near the heads of Iowa Gulch and Empire Gulch, and the Weston Pass district (Fig. 1), and found that these baritic deposits form a somewhat distinct outer zone in which the primary ore consists of galena and light-colored zinc blende, with small amounts of silver and practically no gold, in a gangue of barite, jasperoid, and crystalline quartz. There is also recognized a still more remote though poorly defined zone in which the same sulphides are present but jasperoid is less conspicuous and barite is practically absent, whereas dolomite is the characteristic gangue mineral. The deposits in these two outer zones are fissure-fillings and comparatively small replacement bodies. They contrast strongly with the sulphide deposits in the heart of the Leadville district and are assigned to the outer mesothermal and epithermal zones. As they have not been described before, general features of these deposits are given in some detail, followed by their relations to the deposits within the Leadville district . Analogies are also drawn with similar deposits in other districts. In the preparation of this paper the authors have been impressed with the confusion that may result from too simple a conception of zoning . For a proper conception of Leadville , the geologist must have sufficient regard for structural conditions, must make due allowance for successive stages of mineral deposition within a small block of ground or even within a single ore shoot, and must also allow for the reactions between minerals formed at an early stage and solutions introduced during a late stage. The effects of such complications in the Leadville district might give the impression that certain minerals found in one of the inner zones were formed at relatively high temperatures, but close study of paragenesis indicates late deposition at relatively low temperatures. Again, as already stated, there is sufficient inconsistency in the distribution of certain varieties of deposits to give a confused impression of zonal arrangement unless the influence of structure on the circulation of ore -forming solutions is appreciated. The field facts give a far more complex picture than might be anticipated from a simple statement of the theory of zoning .

Colorado↗

USGS2017

Runoff and water-quality characteristics of three Discovery Farms in North Dakota, 2008–16

Agricultural producers in North Dakota are aware of concerns about degrading water quality, and many of the producers are interested in implementing conservation practices to reduce the export of nutrients from their farms. Producers often implement conservation practices without knowledge of the water quality of the runoff from their farm or if conservation practices they may implement have any effect on water quality. In response to this lack of information, the U.S. Geological Survey, in cooperation with North Dakota State University Extension Service and in coordination with an advisory group consisting of State agencies, agricultural producers, and commodity groups, implemented a monitoring study as part of a Discovery Farms program in North Dakota in 2007. Three data-collection sites were established at each of three farms near Underwood, Embden, and Dazey, North Dakota. The purpose of this report is to describe runoff and water-quality characteristics using data collected at the three Discovery Farms during 2008–16. Runoff and water-quality data were used to help describe the implications of agricultural conservation practices on runoff and water-quality patterns. Runoff characteristics of monitoring sites at the three farms were determined by measuring flow volume and precipitation. Runoff at the Underwood farm monitoring sites generally was controlled by precipitation in the area, antecedent soil moisture conditions, and, after 2012, possibly by the diversion ditch constructed by the producer. Most of the annual runoff was in March and April each year during spring snowmelt. Runoff characteristics at the Embden farm are complex because of the mix of surface runoff and flow through two separate drainage tile systems. Annual flow volumes for the drainage tiles sites (sites E2 and E3) were several orders of magnitude greater than measured at the surface water site E1. Site E1 generally only had runoff briefly in March and April during spring snowmelt and during only a few large rain events throughout 2009–16. Flow was somewhat continuous at sites E2 and E3 throughout the year during years of increased precipitation, such as in 2010 and 2011. At Dazey farm, annual flow volumes at the most downstream site D3 for 2010–15 ranged from 88 acre-feet (2012) to 12,060 acre-feet (2010). The largest monthly runoff volumes at D1 (most upstream site; combination of data from site D1a [original site] and site D1b [relocated site]) and D3 were in March and April during spring snowmelt runoff and rain events. At Underwood farm, total ammonia and total phosphorus had the highest concentrations at the most upstream site (U1) and decreased sequentially at sites U2 and U3 downstream. Total ammonia and total phosphorus concentrations at the sites for Underwood farm also generally were higher than measured at sites for the Dazey and Embden farms. At Embden farm, nitrate plus nitrite concentrations were lowest at site E1 (surface-water site) and highest at sites E2 and E3 (drainage tile sites). Nitrate plus nitrite concentrations at sites E2 and E3 also were the highest among all the sites at all three farms. Median total nitrate plus nitrite concentrations for sites E1, E2, and E3 were 0.22, 13, and 10 milligrams per liter as nitrogen, respectively. Nutrient concentrations generally were greater at site D1 (most upstream site) compared to site D3 (most downstream site) at Dazey farm. Higher concentrations at site D1, which is farther upstream and closer to potential sources of nutrients, compared to lower concentrations at site D3, which is farther downstream and receives more runoff, indicates that dilution may be the reason concentrations decrease downstream. Annual loads for chloride at all three Underwood sites were the greatest in 2011 and the least in 2012, which coincided with years of the greatest and least annual flow volume, respectively. Total ammonia had a similar pattern at the three sites. Nitrate plus nitrite loads displayed a different pattern than chloride and total ammonia, indicating possible different sources. Chloride, total ammonia, total phosphorus, and suspended sediment were transported past site U1 mostly in March and the least from July through October. Monthly nitrate plus nitrite loads had a different pattern than the other constituents, indicating other possible sources such as fertilizer application in the surrounding cropland. Annual loads for Embden farm were considerably greater at sites E2 and E3 compared to site E1. Annual yields for all constituents also were substantially greater at sites E2 and E3 compared to site E1, mainly because of a combination of higher flow volumes and small contributing drainage areas at sites E2 and E3 compared to site E1. The greatest annual loads at Dazey farm site D3 for chloride, nitrate plus nitrite, and suspended sediment were in 2010 and 2011, and zero loads were estimated for 2012 because no flow was measured at the site. Mean monthly loads generally were greatest for most constituents in March and April at sites D1 and D3 except for suspended sediment that had the greatest monthly loads in May. To mitigate runoff and water-quality effects of their operations, the producers implemented various agricultural conservation practices before and during the Discovery Farms monitoring. Even though it was difficult to quantify the effects of the agricultural conservation practices implemented at the farms, the data collected from the Discovery Farms program provided a better understanding of some of the variables that affect runoff and water quality.

North Dakota↗

USGS2005

Reconnaissance study of late quaternary faulting along Cerro Goden fault zone, western Puerto Rico

The Cerro Goden fault zone is associated with a curvilinear, continuous, and prominent topographic lineament in western Puerto Rico. The fault varies in strike from northwest to west. In its westernmost section, the fault is ∼500 m south of an abrupt, curvilinear mountain front separating the 270- to 361-m-high La Cadena de San Francisco range from the Rio Añasco alluvial valley. The Quaternary fault of the Añasco Valley is in alignment with the bedrock fault mapped by D. McIntyre (1971) in the Central La Plata quadrangle sheet east of Añasco Valley. Previous workers have postulated that the Cerro Goden fault zone continues southeast from the Añasco Valley and merges with the Great Southern Puerto Rico fault zone of south-central Puerto Rico. West of the Añasco Valley, the fault continues offshore into the Mona Passage (Caribbean Sea) where it is characterized by offsets of seafloor sediments estimated to be of late Quaternary age. Using both 1:18,500 scale air photographs taken in 1936 and 1:40,000 scale photographs taken by the U.S. Department of Agriculture in 1986, we identified geomorphic features suggestive of Quaternary fault movement in the Añasco Valley, including aligned and deflected drainages, apparently offset terrace risers, and mountain-facing scarps. Many of these features suggest right-lateral displacement. Mapping of Paleogene bedrock units in the uplifted La Cadena range adjacent to the Cerro Goden fault zone reveals the main tectonic events that have culminated in late Quaternary normal-oblique displacement across the Cerro Goden fault. Cretaceous to Eocene rocks of the La Cadena range exhibit large folds with wavelengths of several kms. The orientation of folds and analysis of fault striations within the folds indicate that the folds formed by northeast-southwest shortening in present-day geographic coordinates. The age of deformation is well constrained as late Eocene–early Oligocene by an angular unconformity separating folded, deep-marine middle Eocene rocks from transgressive, shallow-marine rocks of middle-upper Oligocene age. Rocks of middle Oligocene–early Pliocene age above unconformity are gently folded about the roughly east-west–trending Puerto Rico–Virgin Islands arch, which is well expressed in the geomorphology of western Puerto Rico. Arching appears ongoing because onshore and offshore late Quaternary oblique-slip faults closely parallel the complexly deformed crest of the arch and appear to be related to extensional strains focused in the crest of the arch. We estimate ∼4 km of vertical throw on the Cerro Goden fault based on the position of the carbonate cap north of the fault in the La Cadena de San Francisco and its position south of the fault inferred from seismic reflection data in Mayaguez Bay. Based on these observations, our interpretation of the kinematics and history of the Cerro Goden fault zone includes two major phases of motion: (1) Eocene northeast-southwest shortening possibly accompanied by left-lateral shearing as determined by previous workers on the Great Southern Puerto Rico fault zone; and (2) post–early Pliocene regional arching of Puerto Rico accompanied by normal offset and right-lateral shear along faults flanking the crest of the arch. The second phase of deformation accompanied east-west opening of the Mona rift and is inferred to continue to the present day.

Puerto Rico↗

USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗

USGS1965

Hydrology of the Little Plover River basin, Portage County, Wisconsin, and the effects of water resource development

The Little Plover River basin is in the sand-plain area of central Wisconsin. The basin and the surrounding sand-plain area provide a good fish and wildlife habitat and is a popular locale for sport fishing. Good yields may be obtained in the area from irrigated crops, and the irrigated acreage has been increasing rapidly in recent years. Sportsmen and conservationists are concerned about the effects of increased development of the water resources on the streams as trout habitat. In the past, many political and legal conflicts among water users have arisen from erroneous opinions as to the behavior of water. Many of these conflicts would be diminished or eliminated if the participants were cognizant of fundamental hydrologic principles. This study was made to demonstrate the extent and nature of the interrelation of ground water and surface water and the fundamental hydrologic principles governing water movement. The study was also made to determine the hydrologic changes that might occur following development, to provide information that might be used as a basis for planning water development, and for drafting legislation that recognizes the relation between ground water and surface water. Water has been developed in the Little Plover River basin for industry, for domestic and stock supplies, and for irrigation. Irrigated acreage is increasing in the area, and the use of water for irrigation may alter the hydrology of the basin somewhat. About 4,000-4,500 acres of land within the basin, or 50-60 percent of the basin area, is suitable for irrigated farming, but probably no more than 2,500 acres will be under irrigation in any one year, unless present crop-rotation practices are changed. Most of the Little Plover River basin is underlain by from 40 to 100 feet of glacial outwash consisting of highly permeable sand and gravel. The glacial outwash is the main aquifer in the area and is capable of yielding large quantities of water to wells. An aquifer test in the area indicated that the coefficient of transmissibility of the glacial outwash is about 140,000 gallons per day per foot. The specific yield of the outwash is about 20 percent, as determined from water-level and streamfiow data. Morainal deposits occur locally with the glacial outwash. These deposits transmit water readily and do not form barriers to ground water in the outwash. Relatively impermeable crystalline rocks underlie the glacial deposits, and a sandstone ridge of low permeability impedes the movement of ground water from the basin by underflow. The glacial outwash and morainal deposits are recharged by infiltration of 9-10 inches of the 31 inches of precipitation that falls on the area in an average year. If it is not withdrawn by wells for consumptive use or by phreatophytes, water that infiltrates the sand and gravel discharges later into the Little Plover River. This ground-water discharge constitutes 90-95 percent of the total flow of the Little Plover River. Annual evapotranspiration varies considerably, but generally ranges from 2 to 8 inches less than the potential evapotranspiration of 24 inches. Consumptive use of irrigation water averages about 4 inches per year. Most of the water pumped from wells otherwise would be discharged to the stream, and consumptive use of irrigation water will deplete streamflow by the amount of evapotranspiration. Pumping wells have little effect on the water level in the highly permeable sand and gravel. Significant interference between wells would occur only if large capacity wells were within a few tens of feet of each other. Ground water and surface water are closely interrelated in the sand-plain area and ground-water withdrawals near the Little Plover River may cause a measurable streamflow depletion. In a test, a well that was pumping about 1,120 gpm (gallons per minute) and that was 300 feet from the stream derived about 30 percent of its flow from the stream after 3 days of pumping. For this study, the effects of increased ground-water development were evaluated from a hypothetical development schedule, for which it was assumed that 500 acres were irrigated the first year and that an additional 50 acres were irrigated in each succeeding year for 10 years. It also was assumed that the average annual consumptive-use requirement for irrigation water would be one- third of an acre-foot per acre. Calculations indicate that the maximum monthly rate of depletion due to the consumptive use of 4 inches of ground water per year on 500 acres would be about 0.4 cfs (cubic feet per second) the first year and 0.5 cfs after 10 years of pumping. Other computations indicate that the maximum monthly rate of depletion due to irrigating 500 acres the first year and 50 additional acres each year for 10 years would be about 0.8 cfs. Maximum depletion would occur during the summer months, concurrent with the irrigation withdrawals. Because of the close interrelation between ground and surface water, surface- water withdrawals will cause an increased inflow of ground water to the stream and a decline in ground-water levels near the stream. These effects were demonstrated by pumping from the stream. After 29 hours of pumping, a depletion of 1,120 gpm at a site 7,000 feet downstream was about 200 gpm less than the diversion at the pump. Most of the 200 gpm was supplied from the stream-banks, and ground-water levels near the stream declined as much as 0.3 foot. Computations indicated that ground-water inflow, following a streamflow diver- sion that lowered the stage 0.15 foot, would be 0.14 cfs after 5 days and 0.06 cfs after 30 days. The demonstration of the quantitative relation between ground water and surface water, as given by this study, should provide a sound basis for planning water development to minimize conflicts of interest. The demonstrations also should provide a basis for drafting legislation that recognizes the interrelation of ground water and surface water. Because the geology and the hydrology are relatively uniform throughout the sand plains, many of the methods and hydrologic values determined for this detailed study of the Little Plover River basin may be applied to other basins in the sand-plain area.

Wisconsin↗

USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗

USGS2018

High-water marks from Hurricane Sandy for coastal areas of Connecticut, Rhode Island, and Massachusetts, October 2012

Because coastal areas in Connecticut, Rhode Island, and Massachusetts were heavily affected by Hurricane Sandy in October 2012, the U.S. Geological Survey (USGS), under a mission agreement with the Federal Emergency Management Agency, collected storm tide high-water marks in those coastal areas. This effort was undertaken to better understand the areal extent and impact of storm tides resulting from strong storms. On October 27–29, 2012, Hurricane Sandy moved up the eastern coast of the United States after passing over the Bahamas. On October 29 at about 7:30 p.m. eastern daylight time, Hurricane Sandy made landfall its final time near Brigantine, New Jersey, with recorded wind speeds of about 80 miles per hour. The damages from Hurricane Sandy exceeded $50 billion in total, making it the second most costly Atlantic hurricane at that time, second only to Hurricane Katrina in 2005. Hurricane Sandy also resulted in 147 deaths, and about 650,000 homes and many businesses being damaged along the eastern coast of the United States. The severity of Hurricane Sandy’s effects resulted in presidential disaster declarations being declared in 10 States from Virginia to Massachusetts and the District of Columbia in the months following Hurricane Sandy; the list of States affected included Connecticut, Rhode Island, and Massachusetts. In response to the approach of Hurricane Sandy, the USGS deployed 60 temporary storm tide sensors and 2 temporary real-time rapid deployment gages to collect tide elevation data during the storm along the coastal areas of Connecticut, Rhode Island, and Massachusetts. This activity was done from Virginia to Maine before the storm. Following Hurricane Sandy, in October and November 2012, 371 storm tide high-water marks were identified and flagged in the coastal areas of Connecticut, Rhode Island, and Massachusetts. High-water marks near USGS temporary storm tide sensors, real-time rapid deployment gages, and streamgages affected by the tides as well as high-water marks on Block Island, R.I., and Martha’s Vineyard and Nantucket, Mass., were surveyed at the same time the high-water marks were identified and flagged in October and November 2012. The remaining high-water marks flagged during October and November 2012 were surveyed from December 2013 through June 2014 and in December 2016. Elevations of all high-water marks were referenced to the North American Vertical Datum of 1988 and horizontal coordinates to the North American Datum of 1983 using the Global Navigation Satellite System, survey-grade Digital Global Positioning System receivers, and total station surveying equipment. Of the 371 storm tide high-water marks flagged following Hurricane Sandy, only 364 high-water marks were surveyed; the remaining 7 could not be found or had been destroyed when locations were revisited to conduct surveys. The 157 high-water marks surveyed in Connecticut had elevations that ranged from 2.5 to 12.2 feet (ft) with an average elevation of 8.1 ft and a median elevation of 8.3 ft. The 76 high-water marks in Rhode Island had elevations that ranged from 3.6 to 16.2 ft and averaged 7.1 ft with a median of 6.6 ft. The 131 high-water marks in Massachusetts had elevations that ranged from 2.8 to 22.7 ft and averaged 7.3 ft with a median of 6.6 ft. Individual information on the location, type, accuracy, and elevation of the 371 high-water marks can be found in an accompanying USGS data release and at the USGS Flood Event Viewer website for Hurricane Sandy ( https://stn.wim.usgs.gov/fev/#Sandy ). The high-water marks along the coast line of Connecticut and eastern Massachusetts, including Nantucket, generally had higher storm tide elevations than the coast line of Rhode Island including Block Island and southern Massachusetts, including Martha’s Vineyard. The high-water mark elevations compare well with recorded peak-storm tide data at USGS temporary storm tide sensors and real-time rapid deployment gages deployed for Hurricane Sandy in Connecticut, Rhode Island, and Massachusetts. High-water mark data collected following Hurricane Sandy will be used by Federal, State, and local government agencies, nongovernmental organizations, universities, and the public for better understanding the areal extent and impact of the storm tides. Additionally, these data can be used for such activities as land-use planning, flood risk studies, flood resiliency studies, and coastal modeling. These data from this historic storm can be compared with other regional hurricanes and tropical storms for planning into the future.

Connecticut, Massachusetts, Rhode Island↗

USGS2013

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗

USGS1990

The Ag-Mn-Pb-Zn vein, replacement, and skarn deposits of Uchucchacua, Peru; studies of structure, mineralogy, metal zoning, Sr isotopes, and fluid inclusions

Uchucchacua is an Ag-Mn-Pb-Zn vein, replacement, and skarn mineral district in the central Andes of Peru. Host rocks are massive Jumasha Formation shelf limestones of Turonian age that have been folded into an asymmetric northeast-verging anticline of Andean trend. Strata near the fold crest are cut by minor dacitic intrusions and have been displaced by a conjugate set of steep wrench faults that strike northwest-southeast and northeast-southwest. Most ore occurrences are restricted to host rocks that lie below marly limestone at the top of the middle Jumasha Formation. Vein ores located along the fracture system have formed by fissure infill and by replacement of limestone wall rocks. Larger sheetlike replacement orebodies are parallel and adjacent to a large fault. Replacement is in zones of brecciation adjacent to fault bends which were dilatent during sinistral slip of the master fracture. Such fracture belts may have been subject to paleokarst solution processes before mineralization.Four paragenetic stages have been identified. Fe, Mn, and Si were introduced at the exoskarn stage (I) as the anhydrous silicates ferroan tephroite, johannsenite, rhodonite, and bustamite. During the early main stage (II) ferroan tephroite was replaced by friedelite and magnetite under oxidizing conditions at a relatively low pH. Pb, Zn, Fe, Cu, and B were introduced; principal sulfides are pyrrhotite, Fe-rich sphalerite, Mn-rich wurtzite, alabandite, galena, chalcopyrite, and tetrahedrite. Pyrrhotite was replaced by other sulfides during later stage II. Main gangue minerals were calcite, kutnohorite, rhodochrosite, and quartz. Ag, As, and Sb were introduced during the late stage (III) in the form of sulfosalts, principally pyrargyrite. Redistribution of metals introduced at stage II resulted in the growth of Fe-poor sphalerite and alabandite accompanied by calcite gangue. Decreasing Fe contents of alabandite and sphalerite during late stage II and stage III, together with the appearance of pyrite, indicate an increase in sulfur fugacity and/or decrease in temperature over this period. The supergene stage (IV) affects the upper 30 to 150 m of most veins and involves the growth of Mn hydroxides, goethite, orpiment, marcasite, cerussite, and siderite.Distribution patterns of metal ratios and high metal values define ore bands, with a succession of antiforms and synforms. Ore-band locations are determined by vein width, itself a function of wall-rock reactivity, and fracture permeability. Metal ratio distributions in the Luz vein show symmetric dispersion of metals at right angles to the ore-band axis. Ag shows the least dispersion with progressively greater dispersion of Pb and Zn.Geologic considerations indicate mineralization at a minimum depth of 1,600 m. Most primary fluid inclusions in calcite of probable stage II origin were trapped over a temperature range of about 200 degrees to 337 degrees C (hydrostatic pressure) or about 225 degrees to 362 degrees C (lithostatic). Fluid inclusions from late stage II quartz were trapped at about 165 degrees to 205 degrees C (hydrostatic) or 190 degrees to 230 degrees C (lithostatic). Salinity ranged from 0.5 to 29.7 equiv wt percent NaCl-CaCl 2 with little relation to temperature. Calcite 87 Sr/ 86 Sr values of 0.70743 to 0.71122 extend well outside the range of values both determined for magmatic rocks of the region and estimated for the host limestones but are probably compatible with older sedimentary rocks of the mine district. Salinity levels and NaCl/CaCl 2 ratios are similar to basinal brines associated with Mississippi Valley-type deposits. It is concluded that basinal brine, expelled from the deforming Mesozoic sequence, was a significant component of the ore fluid. This fluid was probably heated by dacitic magmas or cooling dacitic intrusions; a polygenetic high-salinity magmatic-basinal fluid then rose along existing strike-slip fractures toward the mineralization site where it mixed with heated low-salinity ground water.

Uchucchacua↗

USGS2013

Detections, concentrations, and distributional patterns of compounds of emerging concern in the San Antonio River Basin, Texas, 2011-12

During 2011–12, the U.S. Geological Survey, in cooperation with the San Antonio River Authority, evaluated detections, concentrations, and distributional patterns of selected compounds of emerging concern (hereinafter referred to as “CECs”) from water-quality samples (hereinafter referred to as “samples”) collected at a total of 20 sampling sites distributed throughout the San Antonio River Basin, Texas. Of the 54 wastewater compounds analyzed, 32 were detected in at least one sample collected from the San Antonio River Basin, and 22 of those compounds were not detected in any samples. The flame retardants tris (2-chloroethyl) phosphate and tris (dichloroisopropyl) phosphate, both possible endocrine disruptors, were the most frequently detected wastewater compounds with 28 of the 33 samples analyzed for wastewater compounds having measureable concentrations of those compounds. Of the 13 analyzed pharmaceuticals, 4 compounds were detected in a least one sample. Carbamazepine, an anticonvulsant, was the most frequently detected prescription pharmaceutical with 24 detections in 34 samples analyzed for pharmaceuticals. Of the 17 steroidal hormones, 4 were detected in at least one sample from the San Antonio River Basin. Estrone was detected in 9 of 34 samples analyzed for steroidal hormones, making it the most frequently detected steroidal hormone. Of the 4 sterols, all 4 were detected in at least one sample from the San Antonio River Basin. Cholesterol, detected in 19 of 34 samples analyzed for sterols, was the most frequently detected sterol. Three synoptic sampling events were completed as part of this study. The first and second synoptic sampling events included samples collected at the same 12 sampling sites. During the first and second synoptic sampling events, the lowest number of detections (2 and 0, respectively) and the lowest total concentrations of all measured compounds (0.62 and not measureable, respectively) occurred in samples collected at the Macdona site (Medina River near Macdona, Tex.). The highest number of detections (21 and 23, respectively) and highest total concentrations of all measured compounds (7.75 and 3.97 micrograms per liter [µg/L], respectively) occurred in samples collected at the SAR Elmendorf site (San Antonio River near Elmendorf, Tex.). The third synoptic sampling event included samples collected at seven sites that were added to the study after the first two synoptic sampling events were completed. During the third synoptic sampling event, the lowest number of detections (two) and the lowest total concentration (0.14 µg/L) of compounds were measured in samples collected at the North Prong site (North Prong Medina River above confluence Wallace Creek near Medina, Tex.). The highest number of detections (21) occurred at the SAR Mitchell site (San Antonio River at Mitchell Street, San Antonio, Tex.). The Dos Rios site (the Dos Rios wastewater treatment plant outfall at San Antonio, Tex.) had the highest total concentration of all measured compounds (4.37 µg/L) in the third synoptic sampling event. Because Ecleto Creek flows only intermittently at the Ecleto site (Ecleto Creek near Runge, Tex.), samples from the Ecleto site were collected at different times than were samples from the other sites and were not included in a synoptic sampling event. The presence of wastewater compounds at the Ecleto site indicates that at least some wastewater compounds can be introduced into surface waters in rural parts of the San Antonio River Basin during runoff or because of onsite wastewater system seepage. The steroidal hormone and sterols detected at the Ecleto site, including estrone, cholesterol, beta -sitosterol, and beta -stigmastanol, likely were derived from cattle waste rather than from wastewater effluent. The distributional patterns of detections and concentrations of individual compounds and compound classes show the influence of wastewater-treatment plant (WWTP) outfalls on the quality of water in the San Antonio River Basin. In the Medina River Subbasin, the minimal influence of wastewater is evident as far downstream as the Macdona site. Downstream from the Macdona site, the Medina River receives treated municipal wastewater from both the Medio Creek Water Recycling Center site from an unnamed tributary at the plant and the Leon Creek Water Recycling Center site from Comanche Creek at the plant, and corresponding increases in both the number of detections and the total concentrations of all measured compounds at all downstream sampling sites were evident. Similarly, the San Antonio River receives treated municipal wastewater as far upstream as the SAR Witte site (San Antonio River at Witte Museum, San Antonio, Tex.) and additional WWTP outfalls along the Medina River upstream from the confluence of the Medina and San Antonio Rivers. Consequently, all samples collected along the main stem of the San Antonio River had higher concentrations of CECs in comparison to sites without upstream WWTPs. Sites in urbanized areas without upstream WWTPs include the Leon 35 site (Leon Creek at Interstate Highway 35, San Antonio, Tex.), the Alazan site (Alazan Creek at Tampico Street, San Antonio, Tex.), and the San Pedro site (San Pedro Creek at Probandt Street, at San Antonio, Tex.). The large number of detections at sites with no upstream wastewater source demonstrated that CECs can be detected in streams flowing through urbanized areas without a large upstream source of treated municipal wastewater. A general lack of detection of pharmaceuticals in streams without upstream outfalls of treated wastewater appears to be typical for streams throughout the San Antonio River Basin and may be a useful indicator of point-source versus nonpoint-source contributions of these compounds in urban streams. Observations of lower concentrations of compounds at the furthest downstream sampling sites in the basin indicate some natural attenuation of these compounds during transport; however, a more focused assessment is needed to make this determination.

Texas↗

USGS2015

Climate and streamflow characteristics for selected streamgages in eastern South Dakota, water years 1945–2013

Upward trends in precipitation and streamflow have been observed in the northeastern Missouri River Basin during the past century, including the area of eastern South Dakota. Some of the identified upward trends were anomalously large relative to surrounding parts of the northern Great Plains. Forcing factors for streamflow trends in eastern South Dakota are not well understood, and it is not known whether streamflow trends are driven primarily by climatic changes or various land-use changes. Understanding the effects that climate (specifically precipitation and temperature) has on streamflow characteristics within a region will help to better understand additional factors such as land-use alterations that may affect the hydrology of the region. To aid in this understanding, a study was completed by the U.S. Geological Survey, in cooperation with the East Dakota Water Development District and James River Water Development District, to assess trends in climate and streamflow characteristics at 10 selected streamgages in eastern South Dakota for water years (WYs) 1945–2013 (69 years) and WYs 1980–2013 (34 years). A WY is the 12-month period, October 1 through September 30, and is designated by the calendar year in which it ends. One streamgage is on the Whetstone River, a tributary to the Minnesota River, and the other streamgages are in the James, Big Sioux, and Vermillion River Basins. The watersheds for two of the James River streamgages extend into North Dakota, and parts of the watersheds for two of the Big Sioux River streamgages extend into Minnesota and Iowa. The objectives of this study were to document trends in streamflow and precipitation in these watersheds, and characterize the residual streamflow variability that might be attributed to factors other than precipitation. Residuals were computed as the departure from a locally-weighted scatterplot smoothing (LOWESS) model. Significance of trends was based on the Mann-Kendall nonparametric test at a 0.10 significance level. Of the 10 streamgages selected, only the Elm River at Westport (in the upper part of James River Basin) did not have a significant upward trend in annual mean streamflow for WYs 1945–2013, whereas only one-half of the streamgages had significant upward trends in annual mean streamflow for WYs 1980–2013. Mean and 7-day minimum streamflows also had upward trends for the spring runoff period (March–May) for most of the streamgages during WYs 1945–2013 and for one streamgage during WYs 1980–2013. Magnitudes of increases in streamflow were as great as 30 cubic feet per second per year for the streamgage on the James River near Scotland during WYs 1980–2013. Precipitation trends for WYs 1945–2013 were not necessarily significant for the watersheds of streamgages with a significant streamflow trend. Annual total precipitation had a significant upward trend for the watersheds of 4 of the 10 streamgages during WYs 1945–2013 and no significant trends for WYs 1980–2013. The most widespread precipitation increase was for September–November, with significant upward trends for the watersheds of 8 of the 10 streamgages during WYs 1945–2013; however, no trends in September– November precipitation were significant for WYs 1980–2013. The greatest magnitude of increase in precipitation was for the December–May season during WYs 1980–2013, which had a mean increase of 0.106 inch per year in the watersheds of streamgages with significant trends. The correlation between streamflow and precipitation metrics was low as indicated by the mean coefficient of determination (R 2 ) of 0.18 for all pairs considered. The highest locally-weighed scatterplot smoothing (LOWESS) correlation was between annual precipitation (by water year) and annual mean streamflow (by water year), which had a mean R 2 of 0.47 for all streamgages and was as high as 0.72 for one streamgage. The correlation between annual precipitation and March–May mean streamflow had a mean R 2 of 0.33 for all streamgages and was as high as 0.52 for one streamgage. Other metrics had R 2 values for LOWESS correlations that were less than 0.3 and were not further considered for analyses of residuals. For annual precipitation as a predictor of annual mean flow, precipitation-removed streamflow had significant upward trends during WYs 1945–2013 for one-half of the streamgages. Upward trends in residual annual mean streamflow were indicated for the Whetstone River and lower part of the Big Sioux River Basin, indicating that other factors are contributors to streamflow variability during WYs 1945–2013. In contrast, most of the streamgages in the James and Vermillion River Basins had no trends in residual annual mean streamflow, indicating that streamflow trends can be explained primarily by precipitation. Precipitation-removed streamflow had fewer trends during the more recent analysis period of WYs 1980–2013 than WYs 1945–2013 for all streamgages in eastern South Dakota. Upward trends in residuals for March– May mean streamflow were indicated for Skunk Creek at Sioux Falls and the Big Sioux River at Akron, but trends in residuals were not significant at the remaining streamgages. For the streamgages with significant trends in residual streamflow (such as the streamgage on the Whetstone River and streamgages in the Big Sioux River Basin), land-use changes likely are minor factors, with the main factors probably being changes in the timing and frequency of large precipitation events and persistently wetter antecedent conditions. Changes in the relation between precipitation and streamflow since 1945 were evident when considering the runoff efficiency of the watershed. For example, the streamflow response to annual precipitation of 25 inches for the James River near Scotland increased from approximately 1,000 cubic feet per second for WYs 1945–1990 to about 2,500 cubic feet per second for WYs 1991–2013. The importance of antecedent conditions on annual mean streamflow also was indicated by the significance of the multiple linear regression coefficients of annual mean streamflow and precipitation from preceding water years for all but one streamgage. In addition, rising groundwater levels are present in wells in eastern South Dakota, particularly since the 1980s.

South Dakota↗

USGS1998

Detailed study of selenium and selected constituents in water, bottom sediment, soil, and biota associated with irrigation drainage in the San Juan River area, New Mexico, 1991-95

In response to increasing concern about the quality of irrigation drainage and its potential effects on fish, wildlife, and human health, the U.S. Department of the Interior began the National Irrigation Water Quality Program (NIWQP) to investigate these concerns at irrigation projects sponsored by the Department. The San Juan River area in northwestern New Mexico was one of the areas designated for study. Study teams composed of scientists from the U.S. Geological Survey, the U.S. Fish and Wildlife Service, the Bureau of Reclamation, and the Bureau of Indian Affairs collected water, bottom-sediment, soil, and biological samples at 61 sites in the San Juan River area during 1993-94. Supplemental data collection conducted during 1991-95 by the Bureau of Indian Affairs and its contractor extended the time period and sampling sites available for analysis. Analytical chemistry performed on samples indicated that most potentially toxic elements other than selenium generally were not high enough to be of concern to fish, wildlife, and human health. Element concentrations in some water, bottom-sediment, soil, and biological samples exceeded applicable standards and criteria suggested by researchers in current literature. Selenium concentrations in water samples from 28 sites in the study area exceeded the 2-microgramper-liter (lg/L) wildlife-habitat standard. Vanadium concentrations in water exceeded the 100-Kg/L standard for livestock-drinking water at one site. In biota, selenium and aluminum concentrations regularly equaled or exceeded avian dietary threshold concentrations. In bottom sediment and soil, element concentrations above the upper limit of the baseline range for western soils were: selenium, 24 exceedances; lead, 2 exceedances; molybdenum, 2 exceedances;strontium, 4 exceedances; and zinc, 4 exceedances. Concentrations of total selenium in bottom-sediment and soil samples were significantly greater for Cretaceous than for non-Cretaceous soil types in the study area and were generally similar for habitats within and outside irrigation-affected areas. Mean and median total-selenium concentrations in samples from areas with Cretaceous soil types were 4.6 and 2.2 micrograms per gram (ps/g), respectively. Mean and median total-selenium concentrations in samples from areas with non-Cretaceous soil types were 0.6 and 0.15 pg/g, respectively. Samples from the study area had low concentrations of organic constituents. Organochlorine pesticides and polychlorinated biphenyls were detected in a few biological samples at low concentrations. Polycyclic aromatic hydrocarbon (PAH) compounds were not detected in whole-water samples collected using conventional water-sampling techniques. In tests involving the use of semipermeable-membrane devices to supplement conventional water assays for PAH's, low concentrations of PAH's were found at several locations in the Hammond Irrigation Supply Canal, but were not detected in the Hammond ponds at the downstream reach of the Hammond irrigation service area. PAH compounds do not appear to reach the San Juan River through the Hammond Canal. Data indicate that water samples from irrigation-drainage-affected habitats had increased mean selenium concentrations compared with samples from irrigation-delivery habitat. The mean selenium concentration in water was greatest at seeps and tributaries draining irrigated land (17 μg/L); less in irrigation drains and in ponds on irrigated land (61.tg/L); and least in backwater, the San Juan River, and irrigation-supply water (0.5 - 0.6 μg/L). Statistical tests imply that irrigation significantly increases selenium concentrations in water samples when a Department of the Interior irrigation project is developed on selenium-rich sediments. Water samples from sites with Cretaceous soils had significantly greater selenium concentrations than water samples from sites with non-Cretaceous soils. Water samples from Department of the Interior project irrigation-drainage sites developed on Cretaceous soils contained a mean selenium concentration about 10 times greater than those in samples from Department of the Interior project sites developed on non-Cretaceous soils. Selenium was much less concentrated in water than in bottom sediment, soil, or biota in the study area. The range in concentrations of dissolved selenium in water was less than 1 ptg/L to 37 1.1g/L (less than 1 to 37 parts per billion). The range in concentrations of total selenium in bottom sediment and soil was less than 0.1 to 23lig/g (less than 100 to 23,000 parts per billion). The range in concentration of selenium in biota was less than 0.1 to 24.0 fig/g (less than 100 to 24,000 parts per billion). Data indicated that bioaccumulation and leaching from soil were the important processes at the study area that lead to elevated levels of selenium. Other processes examined included: (1) evapoconcentration of selenium; (2) atmospheric deposition of aerosols containing selenium; and (3) contamination of surface water by point-source or non-point-source discharges. Selenium concentrations in biological samples were evaluated by a number of variables including: (1) media sampled (emergent and submergent plants, nektonic and benthic invertebrates, omnivore/herbivore and carnivore fish, and terrestrial and aquatic amphibians); (2) habitat (San Juan River main-stem reaches, backwaters, tributary reaches, irrigation delivery or drainage canals, and ponds); (3) irrigation project area and reference sites; and (4) soil type (non-Cretaceous or Cretaceous soils). Graphical techniques and nonparametric statistical tests were applied to determine the influence of selected physiographic variables on selenium concentrations in biological samples collected in the San Juan River area. Species of sucker and of smaller fish contained significantly higher selenium concentrations in the upstream portion of the river where a productive community of plants and animals is found that is associated with warming, nutrient-rich waters discharged from an upstream reservoir. Selenium concentrations in algae, odonates, and mosquitofish collected from both irrigation-drain and pond habitats underlain by Cretaceous soils were significantly greater than in those collected from similar habitats underlain by non-Cretaceous soils. Investigators conclude that the major factor affecting the variability of selenium accumulation in biota at aquatic habitats was the presence of underlying Cretaceous soils. Median selenium concentrations were less than 2 lAg/g for plant samples, less than 7 μg/g for invertebrate samples, and less than 6 lAg/g for whole-fish samples collected from aquatic habitats underlain by non-Cretaceous soils. Similar samples collected from aquatic habitats underlain by Cretaceous soils contained median selenium concentrations two to five times greater. Leaching of selenium from Cretaceous soils in the San Juan River area increases the accumulation of selenium concentrations in the biota and thereby increases the exposure and potential health risks associated with selenium to migratory birds, fish, and other wildlife that use these aquatic habitats extensively. Aquatic habitats presenting the greatest average exposure to excess selenium concentrations in the diets of resident wildlife are from consumption of plants, invertebrates, and fish at irrigation-drain habitats underlain by Cretaceous soils. Of the irrigation projects evaluated in the San Juan River area, the highest median selenium concentrations in algae, cattail leaves, odonate nymphs, mosquitofish, and leopard frog samples from the study area were collected from the east hogback irrigation drain.

New Mexico↗

USGS2019

Tectono-magmatic evolution of porphyry belts in the central Tethys region of Turkey, the Caucasus, Iran, western Pakistan, and southern Afghanistan

Exploration in the central Tethys region of Turkey, Armenia, Azerbaijan, Georgia, Iran, and weste rn Pakistan has led to the identification of the giant Reko Diq (24 Mt Cu and 1300 t Au), Sar Cheshmeh (8.9 Mt Cu and 0.46 Mt Mo), Sungun (5.1 Mt Cu and 0.20 Mt Mo), and Kadjaran (4.6 Mt Cu, 0.94 Mt Mo, and 1100 t Au), and 10 other large (1–2 Mt Cu) porphyry deposits including Saindak, Cevizlidere, Teghout, Meiduk, and Halilağa. Continued exploration efforts have also resulted in the development of porphyry-related gold deposits such as Kişladağ (9.6 Moz Au), Çöpler (3.7 Moz Au), Aği Daği (1.7 Moz Au), and Sary Gunay (3.0 Moz Au), and in the generation of several other promising exploration projects. The distribution in space and time of porphyry deposits in the central Tethys region was shaped by complex pre- to post-mineral tectonic, igneous, collisional, uplift and burial events. These events are represented by a partially-overlapping and variably exhumed and covered collage of twenty-six Early Jurassic to Holocene magmatic belts permissive for the occurrence of porphyry deposits (porphyry tracts and sub-tracts). Twelve tracts or sub-tracts are characterized by compressional continental arcs that formed on drifting terranes or continental margins, 10 developed in compressional to extensional intra-oceanic arc and backarc-rift settings, and 4 formed in extensional post-collisional environments over amalgamated terranes. Eight of these belts were variably affected by coeval and younger metamorphic, fold-and-thrust, and extensional faulting events. Fifty-four porphyry Au-(Cu), Cu-Au, Cu-Mo, Mo-Cu deposits, 15 porphyry-related Au, Au-(Mo) and W-(Mo-Au) deposits, 239 porphyry prospects, and 68 other porphyry-related mineral sites were identified in the study region. Of the 376 porphyry and porphyry-related sites, about 11% formed in island arc, 42% in continental arc, 20% in backarc, and 27% in post-collisional settings. Of the 69 porphyry and porphyry-related deposits, 7% developed in intra-oceanic arc, 41% in continental arc, 27% in backarc, and 25% in post-collisional settings. The largest occur in either compressional continental arc (18 deposits including the Reko Diq and Sar Cheshmeh giants) or post-collisional (13 deposits including the Kadjaran and Sungun giants) environments. Ninety percent of the largest porphyry or porphyry-related deposits occur in only 9 of the 26 permissive porphyry tracts or sub-tracts. Moreover, 88, 90, and 77% of the identified Cu, Mo, and Au resources are contained in porphyry deposits that occur in only 4 of these 9 tracts. Of these 4 tracts, 3 outline arc settings, and one delimits a post-collisional environment. The compositional diversity of porphyry intrusions in these tectono-magmatic environments generally varies from island arc settings with the most restricted range (partly alkaline but mainly calc-alkaline dioritic to granodioritic-tonalitic), to continental arc (calc-alkaline dioritic-quartz dioritic, granodioritic, quartz monzonitic-granitic, and less commonly mildly alkaline), to backarc (mildly alkaline and calc-alkaline dioritic to granitic), to post-collisional settings with the most expansive range (alkaline and calc-alkaline mafic to felsic, and weakly peraluminous). Metal associations also vary broadly as a function of porphyry intrusion composition from weakly peraluminous to metaluminous felsic (Mo[±W ± Cu]; <2% of porphyry-related systems [i.e., Tyrnyauz]), to metaluminous felsic and intermediate (Cu-Mo[±Au]; 85% [i.e., Cevizlidere, Haft Cheshmeh, Kahang, Sar Cheshmeh, Sungun, Teghout, Reko Diq, Saindak]), to mildly alkaline felsic and intermediate (Cu-Au[±Mo] [i.e., Agarak, Kadjaran, Kale Kafi]) and mafic (Au-Cu; 12% [i.e., Çöpler]), and to alkaline felsic (Au-Mo; 1% of porphyry-related systems [i.e., Kişladağ). Tectonic changes were critical in triggering the formation of large porphyry deposits in the region. Large porphyry deposits were preferentially emplaced in continental arc settings shortly before major collisional events (Dar Alu, Kahang, Meiduk, Now Chun, and the giant Sar Cheshmeh and Reko Diq deposits), or in post-subduction environments shortly after collision (Bakirçay, Güzelyayla, Haft Cheshmeh, Masjed Daghi, and the giant Kadjaran and Sungun deposits) or during periods of prominent extension (Aği Daği, Halilağa, Kişladağ, Sari Gunay, and Zarshuran porphyry-related deposits). Collision-induced uplift, erosion, and removal of coeval volcanic rocks favorably exposed the hypabyssal level of subduction-related porphyry deposits. Extensional structures that developed parallel and orthogonal to the compressional principal stress component along transtensional or transpressional strike-slip faults or in pull-apart basins commonly controlled porphyry-related deposits in post-collisional settings. The latter deposits typically exhibit shallow epithermal levels of emplacement because of preservation by burial. Seventeen porphyry deposits and one porphyry-related deposit in the study region are reported to contain significant supergene resources. Relatively mature levels of secondary copper enrichment in dominantly granodioritic to granitic porphyry deposits occur in areas where large pyrite-rich quartz-sericite alteration zones have been preserved and exposed to surface oxidation (Güzelyayla and Ulutaş in northeastern Turkey; Agarak, Ankavan, Dastakert, Kadjaran, and Teghout in Armenia; Ali Javad in northern Iran; Kale Kafi in central Iran; Darreh Zar, Meiduk, Now Chun, and Sar Cheshmeh in southeastern Iran; and Tanjeel in southwestern Pakistan). Chalcocite blankets also developed over porphyry deposits in regions where significant post-mineral faulting has occurred (Muratdere and Sarıçayıryayla in western Turkey). Normal faulting also enhanced secondary enrichment of gold in the Halilağa porphyry and Sary Gunay porphyry-related deposits located respectively in western Turkey and northern Iran. Evaluation of provincial as well as local controls strongly suggests that continued exploration in the region will lead to the identification of additional porphyry and porphyry-related deposits. These deposits will likely be found under younger cover formations in porphyry belts that are already known, and in association with superjacent high- and intermediate-sulfidation epithermal deposits, or increasingly peripheral skarn, carbonate-replacement, and sediment-hosted deposits. Application of suitable exploration techniques to detect concealed and/or deformed deposits in porphyry belts that remain under-explored may also prove productive.

Tethys region↗

USGS2006

Evaluation of some software measuring displacements using GPS in real-time

For the past decade, the USGS has been monitoring deformation at various locations in the western United States using continuous GPS. The main focus of these measurements are estimates of displacement averaged over one day. Essentially, these consist of recording at 30 seconds intervals the carrier-frequency phase-data (equivalent to travel-time) between a GPS receiver and the GPS satellite network. In turn, these observations, which are converted to pseudo—ranges, are processed using one of the “research grade” programs (GIPSY, Zumberge et al., or GAMIT, wwwgpsg.mit.edu/~simon/gtgk) to estimate the position of the GPS receiver averaged over 24 hours. However, it is possible and desirable to estimate the position of the receiver (actually the antenna) more frequently and to do this within a few seconds of the time actual measurement (known as real-time). A recent example, the 2004 Magnitude 6, Parkfield, California earthquake, demonstrated that having GPS estimates of position more frequently than simply a daily average is required if one requires discrimination between co-seismic and post-seismic deformation (Langbein et al., 2006). The high-rate estimates of position obtained at Parkfield show that post-seismic deformation started less than one-hour after the mainshock and that this deformation was roughly the same magnitude as the co-seismic deformation. The high-rate solutions for Parkfield were done by others including Yehuda Bock at UCSD and Kristine Larson at U. of Colorado, but not the USGS. The Parkfield experience points out the need for an in-house capability by the USGS to be able to accurately measure co-seismic displacements and other rapid, deformation signals using GPS. This applies to both the Earthquake and Volcano Hazard programs. Although at many locations where we monitor deformation, we have strainmeters and tiltmeters in addition to GPS which, in principle, are far more sensitive to rapid deformation over periods of less than a day (Langbein and Bock, 2004). But, not all locales include strain and tiltmeters. Thus, having the capability to extract signals with periods of less than a day is desirable since the distribution of GPS is more extensive than strain and tilt. At both Parkfield and Long Valley, the USGS has been using other software packages to process the GPS data at sub-daily intervals and in real-time. The underlying goal of these types of measurements is to detect any deformation event as it evolves; the 24 hour processing might not provide timely results if such a deformation event is precursory to a geologic hazard (an earthquake for Parkfield and either a volcanic event or an earthquake for Long Valley). In Long Valley, We use the software package called 3DTracker (http://www.3dtracker.com, http://www.condorearth.com) to estimate the changes of in position of a remote site relative to a “fixed” site. The 3DTracker software uses double difference GPS code measurements and receiversatellite-time triple differences from one epoch to the next of the GPS phase data (a proxy for travel-time measurements) and employs a Kalman filter to obtain stability in the estimate of position. That is, the estimate of the current position depends upon the estimate of the prior position. Hence, a time series of position looks fairly smooth depending upon the coefficient selected for the Kalman filter. With triple differences, the sometimes troublesome initial integer cycle ambiguity terms cancel (number of wavelengths between the receiver and each satellite), but only the incremental change in position is calculated. This triple difference Kalman filter solution is slow to converge and less accurate than a double difference (e.g., RTD, Track) solution, but it is robust and computationally efficient (Remondi and Brown, 2000). 3D-Tracker allows use of various single-frequency and dual-frequency GPS phase and code observables including the ionospheric-free combinations (known as LC or L3 and P(L3)) formed from an linear combination of the L1 and L2 carrier phase and code data. The lowest noise observable is the L1 carrier, but it is biased by ionospheric refraction that has amplitudes of about 1 to 10 ppm. This results in a systematic scale error in the relative positions. The L3 phase noise is about 3 times greater than the L1 phase noise, but it is generally used to solve for all but the shortest baselines (< 5 km). In addition, the software does output the position changes is a standard format that can be used for other analysis. At Parkfield, we use the software package called RTD (http://www.geodetics.com). The RTD software has been described in the literature (Bock et al., 2000) but basically, it estimates the position without the constraint of a Kalman filter. It uses double differences (in our studies the LC or ionospheric free observable is used) and the integer ambiguities are resolved independently for each 1-second measurement; Most GPS software that use double-differences require several epochs of measurements to resolve the integer ambiguities. The data files use a proprietary format and can not be read by me or others; rather, Yehuda Bock at UCSD (and author of RTD) translates these files into a standard format that can be read by me. Recently, Tom Herring of MIT has modified the GAMIT software to process kinematically GPS data (www-gpsg.mit.edu/~simon/gtgk/tutorial/Lecture_13.pdf). At this time, the software, known as TRACK, does not process the observations in real-time. Consequently, the latency between the time of the observation and the time when a position estimate is available depends upon the frequency that the data are downloaded and the speed of actually processing the observations; there could be a delay of an hour or two before the a position estimates are available. Unlike RTD and 3DTracker, TRACK comes with GAMIT (which is distributed freely) and is currently operating in a test mode at the USGS office in Pasadena. The LC or ionosphere free observable is used in our TRACK solutions. JPL has a version of their GIPSY software called “Real-time GIPSY (RTG)” (gipsy.jpl.nasa.gov/orms/rtg), which, like TRACK, can process the pseudo-range data “off—line”. However, this software is not freely distributed. Instead, at least one company, NAVCOM, has teamed with JPL to integrate RTG with GPS receivers and telemetry that yields positions in realtime. Kristine Larson of University of Colorado has modified the original GIPSY to estimate positions kinematically. Again, like TRACK, the positions are estimated off—line. Much of her research is described in Larson et al. (2003), and Choi et al. (2004). For Long Valley, out of the 17 GPS sites, we monitor 5 baselines within the caldera at 5 second intervals relative to the Bald Mountain site at the edge of the caldera using 3DTracker. The baseline measurement using 3DTracker consists of determination of the 3 dimensional positions of the 5 remote points (GPS receivers) relative to a GPS site at Bald. A second, independent system collects and downloads once a day the 30-second data used for the 24-hour solutions for the 12 sites not monitored with 3DTracker. For the sites monitored with 3DTracker, the pseudo—range data are decimated to 30 seconds and converted to a form used for the 24-hour solutions. Both sets of telemetry employ 900 MHz spread spectrum radios which require line of site between all of the links. The telemetry for the 3DTracker sites require a dedicated radios at each end and intermediate repeaters as needed, while the telemetry required for the other sites use a single master radio, repeaters as needed, and a radio at each remote site. (The 5 sites being monitored with 3DTracker require 13 radios.) At Parkfield, RTD is used to measure the position changes all 12 baselines at 1 second intervals relative to a site, Pomm, adjacent to the San Andreas Fault. The complete RTD package (hardware and software) collects all of the data and determines the position of each site relative to Pomm. In addition, the system stores both the 1-second and 30-second pseudo-range data for later downloading which are ultimately used in the 24-hour solutions. To do this, each site has a 2.4 GHz radio and a telemetry buffer. The telemetry buffer holds 24-hours of data (in the event that the telemetry link is broken) and converts the RS232 data stream from the GPS receiver into a form compatible with an IP (Internet protocol) network connection. In contrast with the Long Valley system, the telemetry link for GPS at Parkfield consists of a single radio at each remote sites and a single radio at the central site. Although position estimates are produced within 1-second of the observations, these results are not immediately available because there is no high speed Internet connection to Parkfield. Instead, the data are stored on a removable disk and sent to UCSD once per month. Below, I describe the results of a simple experiment to examine the response of some of these systems to simulated deformation that could be an analogue of a tectonic or volcanic event. In many engineering applications, the system response is tested by inputting a step to the system and measuring the output of the system. Essentially, this is what I've done. The experiment described below moves the GPS antenna from its original position to a new position within 1 second; the software tracks the translation. These measurements were conducted in August 2004 with the RTD software at Parkfield, and twice in Long Valley. The first Long Valley test was conducted in September 2004 using 3DTracker on a single baseline. The test was repeated in September 2005 using 3DTracker on two baselines and, importantly, saving the RINEX files of the data so that the data could be replayed through 3DTracker using other options in the program and, using other software packages including TRACK. In addition, we observed a short-term event at the Three Sisters volcano in Oregon. This event was snow melt at a remote GPS site which gave an apparent 15 cm displacement in vertical in less than one-day. 3DTracker is used to monitor this site, and the event was captured with this software. In addition, with the assistance of others, I got additional estimates of position using other software packages; those results are presented. Finally, the precision of both 3DTracker and RTD are compared using a power spectrum. Those results would suggest that 3DTracker using appropriate Kalman filter coefficients would have better precision than RTD; instead, the lower noise level from 3DTracker is a result of smoothing from the Kalman filter. Given the results described in this report, high-rate GPS is certainly capable of accurately measuring displacements of 1 centimeter with a high degree of statistical confidence. Plotting these results show that the time of the displacement can be visually determined to that of the sampling interval of the data. However, especially with small amplitude signals, any of the software packages can yield erroneous deformation “signals” that are either due excess travel-time of the GPS carrier frequency from multipath or a limitation in the software. Thus, the time series of displacements must be viewed with caution and knowledge of external circumstances that might cause a change in position. The casual reader should continue with the next section describing the methods then jump to the last two sections for the discussion and conclusions. I have made some recommendations there.

Open-File Report↗

USGS2019

India National Gas Hydrate Program Expedition-02: Operational and technical summary

The India National Gas Hydrate Program is being steered by the government of India's Ministry of Petroleum and Natural Gas (MoPNG) with participation of Directorate General of Hydrocarbons (DGH), Oil and Natural Gas Corporation Limited (ONGC), and the National Oil Companies and Research Institutes of India. The India National Gas Hydrate Program Expedition 01 (NGHP-01) established the presence of gas hydrate in the Krishna Godavari (KG) and Mahanadi Basins and in the offshore area of the Andaman Sea Basin. However, the gas hydrates discovered during NGHP-01 were mainly distributed as fracture-filling material in fine-grained clay-rich sediments. The India National Gas Hydrate Program Expedition 02 (NGHP-02) was carried out with an objective to discover gas hydrate in sand-rich sediment along the eastern offshore margin of India. ONGC planned and executed NGHP-02 on the behalf of the MoPNG. NGHP-02 started on March 3, 2015 and was completed on July 28, 2015 (total 147 days) using the Japanese scientific Drilling Vessel Chikyu (D/V Chikyu). During NGHP-02, 42 holes at 25 sites were drilled, cored, and/or surveyed with downhole logging tools. These sites were located in four areas along the eastern margin of India and formally named Area A (Mahanadi Basin, three sites), Area B (northern part of the KG-Basin, twelve sites), Area C (central part of the KG-Basin, six sites), and Area E (southern part to the KG-Basin, four sites). All 25 sites established during NGHP-02 were first drilled and logged with logging-while-drilling (LWD) tools and an additional 17 holes were then drilled and/or cored with conventional coring tools (HPCS/ESCS) or pressure coring tools (PCTB). Wireline logging was conducted in 10 holes and formation tests using a dual packer Modular Formation Dynamics Tester (MDT) tool were carried out in two holes. The onboard science team used the laboratory facilities on the D/V Chikyu to examine and analyse the physical properties, geochemistry, and sedimentology of all the cores collected during the expedition. Core samples were also analysed in additional post-expedition shore-based studies conducted in numerous domestic and international gas hydrate research laboratories. The NGHP-02 sediment cores were archived at the National Gas Hydrate Core Repository in Mumbai (India), which is associated with the ONGC Gas Hydrate Research and Technology Centre (GHRTC). The necessary data for characterizing the occurrence of gas hydrate, such as interstitial water chlorinities, core-derived gas chemistry, physical and sedimentological properties, thermal images of the recovered cores, pressure core and downhole measured logging data (LWD and/or conventional wireline log data), were obtained from most of the drill sites established during NGHP-02. Almost all the drill sites yielded evidence for the occurrence of gas hydrate; however, the inferred in situ concentration of gas hydrate varied substantially from site to site. For the most part, the interpretation of downhole logging data, core thermal images, interstitial water analyses, and pressure core images from the sites established during NGHP-02 indicate that the occurrence of concentrated gas hydrate is mostly associated with coarser grained (sand-rich) sediments. This paper presents the operational and technical summary of NGHP-02. NGHP-02 started on March 3, 2015 and was completed on July 28, 2015 (total 147 days) using the Japanese scientific Drilling Vessel Chikyu (D/V Chikyu). During NGHP-02, 42 holes at 25 sites were drilled, cored, and/or surveyed with downhole logging tools. These sites were located in four areas along the eastern margin of India and formally named Area A (Mahanadi Basin, three sites), Area B (northern part of the KG-Basin, twelve sites), Area C (central part of the KG-Basin, six sites), and Area E (southern part to the KG-Basin, four sites). All 25 sites established during NGHP-02 were first drilled and logged with logging-while-drilling (LWD) to

Journal of Marine and Petroleum Geology↗

USGS2024

Growth, survival, and cohort formation of juvenile Lost River (Deltistes luxatus) and shortnose suckers (Chasmistes brevirostris) in Upper Klamath Lake, Oregon, and Clear Lake Reservoir, California—2021–22 monitoring report

Executive Summary The work reported in this publication provides updated data and interpretation for sampling years 2015 and 2022 of the juvenile monitoring project. The study objectives, background, study area, species description, and methods remained the same or similar throughout the years, while the executive summary, results, and discussion were updated each year. Therefore much of this paper was originally presented in previous reports (Bart and others 2020a, b; Bart and others, 2021; Burdick and others, 2016; Burdick and others, 2018; Martin and others, 2022) and is repeated here for the reader’s convenience. Populations of federally endangered Lost River ( Deltistes luxatus ) and shortnose suckers ( Chasmistes brevirostris ) in Upper Klamath Lake, Oregon, and Clear Lake Reservoir (hereinafter, Clear Lake), California, are experiencing long-term decreases in abundance. Upper Klamath Lake populations are decreasing not only because of adult mortality, which is relatively low, but also because they are not being balanced by recruitment of young adult suckers into adult spawning aggregations. Long-term monitoring of juvenile sucker populations is conducted to (1) determine if there are annual and species-specific differences in production, survival, and growth; (2) better understand when juvenile sucker mortality is greatest; and (3) identify potential causes of high juvenile sucker mortality particularly in Upper Klamath Lake. The U.S. Geological Survey (USGS) monitoring program, begun in 2015, tracks cohorts through summer months and among years in Upper Klamath and Clear Lakes. Data on juvenile suckers captured in trap nets are used to provide information on annual variability in age-0 sucker production, juvenile sucker apparent survival, growth, species composition, and health. Upper Klamath Lake indices of year-class strength suggest that the 2022 age-0 cohort is the lowest since standardized monitoring began. The 2021 cohort, like most cohorts, had moderately low catch rates their first year of life, with a steep drop off during the second year. Although the 2020 cohort persisted through the September 2022 sampling, this cohort was sparsely represented after the first year with no representatives from this cohort captured from July 2021 through July 2022. Despite apparently low fall through spring apparent survival, the relatively large 2019 cohort persisted in our 2020–21 samples, but has not been detected since June 2021. Klamath largescale ( Catostomus snyderi ) and shortnose suckers were only differentiated from each other starting in 2020. Shortnose suckers dominated the age-1 catch in 2020 and 2022, whereas age-1 Klamath largescale suckers were slightly more prevalent in 2021. Although there were occasionally age-2 and older suckers captured, none of these fish were Lost River suckers. Except for 2015, 2017, and 2021, there were more age-0 Lost River suckers than presumed shortnose suckers in Upper Klamath Lake. However, in all years sampled, there were more age-1 presumed shortnose suckers than Lost River suckers. Age distribution of suckers captured in Clear Lake indicates greater juvenile survival than in Upper Klamath Lake. Most juvenile suckers captured throughout the years were from the 2016 and 2017 cohorts; however, by 2022 most of these fish were no longer susceptible to standard trap nets and were not as prevalent in 2022 juvenile catches, and these suckers presumedly recruited to the adult population. As the 2016 and 2017 cohorts catches declined, so did the catch in overall numbers of suckers. Excluding age-0 catches, the 2016 cohort catches peaked at age-3 and the 2017 catches peaked at age-2. In 2022, the majority of the catch was composed of age-3 to age-5 suckers. The majority of suckers captured in Clear Lake during this multiyear project were classified as the combination of Klamath largescale suckers and shortnose suckers from the Lost River Basin, from the 2016 and 2017 cohorts. The few suckers identified as Lost River or definitive shortnose suckers were from the 2016 and 2017 cohorts. A lack of age-0 suckers captured in Clear Lake during years with low spawning tributary inflow or lake levels suggested that low water prevented spawning and year class formation. However, recent data indicate that some cohorts with Klamath largescale and shortnose sucker genetics that were not captured as age-0 suckers were detected in later years at age-1 or age-2. This finding indicates that juvenile suckers in Clear Lake may spend one or more years in the tributaries and that these cohorts may primarily be represented by Klamath largescale suckers. The first 7 years of this monitoring program indicated different patterns in recruitment and survival of juvenile suckers between Upper Klamath and Clear Lakes. Since the monitoring program began in 2015, age-0 sucker catch rates, interpreted as indices of year-class strength, were greatest in Upper Klamath Lake in 2016 and 2019. In those years, Lost River suckers made up the majority of age-0 sucker catches. However, in 2017 and 2020, the age-1 sucker catches from these cohorts were mainly composed of shortnose suckers or suckers with genetic markers of both Klamath largescale and shortnose suckers, indicating a low first year survival for Lost River suckers even when age-0 catches were high. Age-0 suckers do not fully recruit to our sampling gear in Upper Klamath Lake until August, experience high mortality by September, and are almost undetectable in subsequent years. In Clear Lake, suckers are often not captured until age-1 or age-2 and juvenile annual survival appears much greater; however, there does appear to be a drop-off in catch rates as the suckers age and become less susceptible to the fishing gear.

California, Oregon↗