Geology Reports⌕ Search

SEARCH · Geology Reports

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.

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

At least 1,045 records · Page 58Linked to original sources

Paleoseismic potential of sublacustrine landslide records in a high-seismicity setting (south-central Alaska)

Sublacustrine landslide stratigraphy is considered useful for quantitative paleoseismology in low-seismicity settings. However, as the recharging of underwater slopes with sediments is one of the factors that governs the recurrence of slope failures, it is not clear if landslide deposits can provide continuous paleoseismic records in settings of frequent strong shaking. To test this, we selected three lakes in south-central Alaska that experienced a strong historical megathrust earthquake (the 1964 M w 9.2 Great Alaska Earthquake) and exhibit high sedimentation rates in their main basins (0.2 cm yr − 1 –1.0 cm yr − 1 ). We present high-resolution reflection seismic data (3.5 kHz) and radionuclide data from sediment cores in order to investigate factors that control the establishment of a reliable landslide record. Seismic stratigraphy analysis reveals the presence of several landslide deposits in the lacustrine sedimentary infill. Most of these landslide deposits can be attributed to specific landslide events, as multiple landslide deposits sourced from different lacustrine slopes occur on a single stratigraphic horizon. We identify numerous events in the lakes: Eklutna Lake proximal basin (14 events), Eklutna Lake distal basin (8 events), Skilak Lake (7 events) and Kenai Lake (7 events). The most recent event in each basin corresponds to the historic 1964 megathrust earthquake. All events are characterized by multiple landslide deposits, which hints at a regional trigger mechanism, such as an earthquake (the synchronicity criterion). This means that the landslide record in each basin represents a record of past seismic events. Based on extrapolation of sedimentation rates derived from radionuclide dating, we roughly estimate a mean recurrence interval in the Eklutna Lake proximal basin, Eklutna Lake distal basin, Skilak Lake and Kenai Lake, at ~ 250 yrs., ~ 450 yrs., ~ 900 yrs. and ~ 450 yrs., respectively. This distinct difference in recording can be explained by variations in preconditioning factors like slope angle, slope recharging (sedimentation rate) and the sediment source area: faster slope recharging and a predominance of delta and alluvial fan failures, increase the sensitivity and lower the intensity threshold for slope instability. Also, the seismotectonic setting of the lakes has to be taken into account. This study demonstrates that sublacustrine landslides in several Alaskan lakes can be used as reliable recorders of strong earthquake shaking, when a multi-lake approach is used, and can enhance the temporal and spatial resolution of the paleoseismic record of south-central Alaska.

Alaska↗

Regression equations for estimation of annual peak-streamflow frequency for undeveloped watersheds in Texas using an L-moment-based, PRESS-minimized, residual-adjusted approach

Annual peak-streamflow frequency estimates are needed for flood-plain management; for objective assessment of flood risk; for cost-effective design of dams, levees, and other flood-control structures; and for design of roads, bridges, and culverts. Annual peak-streamflow frequency represents the peak streamflow for nine recurrence intervals of 2, 5, 10, 25, 50, 100, 200, 250, and 500 years. Common methods for estimation of peak-streamflow frequency for ungaged or unmonitored watersheds are regression equations for each recurrence interval developed for one or more regions; such regional equations are the subject of this report. The method is based on analysis of annual peak-streamflow data from U.S. Geological Survey streamflow-gaging stations (stations). Beginning in 2007, the U.S. Geological Survey, in cooperation with the Texas Department of Transportation and in partnership with Texas Tech University, began a 3-year investigation concerning the development of regional equations to estimate annual peak-streamflow frequency for undeveloped watersheds in Texas. The investigation focuses primarily on 638 stations with 8 or more years of data from undeveloped watersheds and other criteria. The general approach is explicitly limited to the use of L-moment statistics, which are used in conjunction with a technique of multi-linear regression referred to as PRESS minimization. The approach used to develop the regional equations, which was refined during the investigation, is referred to as the 'L-moment-based, PRESS-minimized, residual-adjusted approach'. For the approach, seven unique distributions are fit to the sample L-moments of the data for each of 638 stations and trimmed means of the seven results of the distributions for each recurrence interval are used to define the station specific, peak-streamflow frequency. As a first iteration of regression, nine weighted-least-squares, PRESS-minimized, multi-linear regression equations are computed using the watershed characteristics of drainage area, dimensionless main-channel slope, and mean annual precipitation. The residuals of the nine equations are spatially mapped, and residuals for the 10-year recurrence interval are selected for generalization to 1-degree latitude and longitude quadrangles. The generalized residual is referred to as the OmegaEM parameter and represents a generalized terrain and climate index that expresses peak-streamflow potential not otherwise represented in the three watershed characteristics. The OmegaEM parameter was assigned to each station, and using OmegaEM, nine additional regression equations are computed. Because of favorable diagnostics, the OmegaEM equations are expected to be generally reliable estimators of peak-streamflow frequency for undeveloped and ungaged stream locations in Texas. The mean residual standard error, adjusted R-squared, and percentage reduction of PRESS by use of OmegaEM are 0.30log 10 , 0.86, and -21 percent, respectively. Inclusion of the OmegaEM parameter provides a substantial reduction in the PRESS statistic of the regression equations and removes considerable spatial dependency in regression residuals. Although the OmegaEM parameter requires interpretation on the part of analysts and the potential exists that different analysts could estimate different values for a given watershed, the authors suggest that typical uncertainty in the OmegaEM estimate might be about +or-0.10 10 . Finally, given the two ensembles of equations reported herein and those in previous reports, hydrologic design engineers and other analysts have several different methods, which represent different analytical tracks, to make comparisons of peak-streamflow frequency estimates for ungaged watersheds in the study area.

Texas↗

A conceptual site model of contaminant transport pathways from the Bremerton Naval Complex to Sinclair Inlet, Washington, 2011–21

Historical activities on the Bremerton Naval Complex (BNC) in Puget Sound, Washington, have resulted in Sinclair Inlet sediments with elevated concentrations of contaminants, including organic contaminants such as polychlorinated biphenyls and trace elements including mercury. Six U.S. Geological Survey–U.S. Navy datasets have been collected since the last major assessment, in 2013, of soil and groundwater contaminant transport pathways and mercury loading estimates from the BNC to Sinclair Inlet. These include: mercury isotope analysis to support sourcing of mercury in Sinclair Inlet; mercury sampling within the dry dock systems; nearshore thermal surveys to identify potential groundwater discharge locations to Sinclair Inlet; time-series monitoring in nearshore wells to understand the inland extent and dynamics of the tidal mixing zone; tidal studies of mercury in nearshore monitoring wells in an area of contaminated fill material called Site 1; and a spatial survey of trace elements and other parameters in nearshore monitoring wells, pore water, seeps, surface water, and sediment along unwalled shorelines in the western part of the BNC. The results were incorporated into an updated Conceptual Site Model and used to update contaminant load estimates from the terrestrial BNC to Sinclair Inlet. The results from these studies provide data to the U.S. Navy to support prioritization of on-going remediation actions to manage contamination on the BNC that reduce potential impacts to Sinclair Inlet sediment, surface water, and fish and shellfish tissue. Mercury isotope analysis of surface sediments and particulate material indicated that a similar industrial mercury profile is present throughout Puget Sound, including terrestrial and marine BNC samples and in other Sinclair Inlet sediments and persists across regions with low and elevated mercury concentrations. Two sources of mercury at the BNC are Sites 1 and 2 subsurface soils/fill material, with total mercury concentrations in particulates collected from the bottom of monitoring wells drilled in these materials ranging from 18,000 to 44,000 nanograms per gram (as compared to the Washington State Marine Sediment Cleanup Screening Level of 590 nanograms per gram). Contaminants are transported from the terrestrial BNC to Sinclair Inlet via three primary pathways, (1) stormwater outfalls, (2) dry dock discharges, and (3) direct discharge along unwalled shorelines. Previous loading estimates (based on filtered total mercury) ranked stormwater outfalls, particularly outfall PSNS015 in Site 2 soils, as the largest soil and groundwater contaminant transport pathway from the terrestrial BNC to Sinclair Inlet. Updated loading estimates in this report suggest that the dry dock systems may be a larger pathway of mercury from the terrestrial BNC to Sinclair Inlet than previously thought, within the same order of magnitude as the PSNS015 storm-drain system. Trace-element loads via direct shoreline discharge are difficult to estimate due to the large and dynamic tidal mixing zone of groundwater and seawater in the nearshore along unwalled shorelines. However, current best estimated ranges suggest that direct shoreline discharge is one of the three main pathways and may contribute smaller mercury loads than the stormwater and the dry dock systems. Along unwalled shorelines, direct groundwater discharge of terrestrial contaminants may be less important than recirculating seawater in the nearshore mixing zone that can extract contaminants from nearshore subsurface material. Total estimated mercury loads from the terrestrial BNC to Sinclair Inlet range from approximately 40 to 200 grams of filtered total mercury per year and a minimum of 70–350 grams of particulate total mercury per year, for a minimum total of 110–525 grams of whole (filtered plus particulate) total mercury per year. Data gaps are identified that, if filled, would further refine the Conceptual Site Model and contaminant loading estimates from the terrestrial BNC to Sinclair Inlet.

Washington↗

Forest Types in the Lower Suwannee River Floodplain, Florida?-A Report and Interactive Map

A map of forest types in the lower Suwannee River floodplain, Florida, was created during a study conducted from 1996 to 2000 by the U.S. Geological Survey in cooperation with the Suwannee River Water Management District. The map is presented with this report on a compact disc with interactive viewing software. The forest map can be used by scientists for ecological studies in the floodplain based on land cover types and by landowners and management personnel making land use decisions. The study area is the 10-year floodplain of the lower Suwannee River from its confluence with the Santa Fe River to the lower limit of forests near the Gulf of Mexico. The floodplain is divided into three reaches: riverine (non-tidal), upper tidal, and lower tidal, due to changes in hydrology, vegetation, and soils with proximity to the coast. The 10-year floodplain covers about 21,170 hectares; nearly 88 percent of this area (18,580 hectares) is mapped as 14 major forest types. Approximately 29 percent (5,319 hectares) of these forests have been altered by agriculture or development. About 75 percent of the area of major forest types (13,994 hectares) is wetland forests and about 25 percent (4,586 hectares) is upland forests. Tidal wetland forests (8,955 hectares) cover a much greater area than riverine wetland forests (5,039 hectares). Oak/pine upland forests are present in the riverine and upper tidal reaches of the floodplain on elevations that are inundated only briefly during the highest floods. High bottomland hardwoods are present on the higher levees, ridges, and flats of the riverine reach where soils are usually sandy. Low bottomland hardwood forests are present in the riverine reach on swamp margins and low levees and flats that are flooded continuously for several weeks or longer every 1 to 3 years. Riverine swamps are present in the lowest and wettest areas of the non-tidal floodplain that are either inundated or saturated most of the time. Upper tidal bottomland hardwood forests are present on sandy soils on high flats and in transitional areas between upland forests and swamps. Upper tidal mixed forests are found on low levees or between swamps and higher forest types. Upper tidal swamps are present at elevations below median monthly high stage and usually have surface soils that are permanently saturated mucks. Lower tidal hammocks are found on higher elevations that do not receive regular tidal inundation but have a high water table and are briefly inundated by storm surges several times a decade. Lower tidal mixed forests include swamps with numerous small hummocks or less common larger hummocks. Lower tidal swamps are found on deep muck soils that are below the elevation of the median daily or monthly high stage. Seven additional land cover types (2,590 hectares) are mapped. Water in the main channel of the lower Suwannee River (1,767 hectares) was mapped separately from open water in the floodplain (239 hectares). Other land cover types are: seepage slopes (70 hectares), isolated forested wetlands (19 hectares), marshes upstream of the tree line (505 hectares), beds of emergent aquatic vegetation (21 hectares), and floodplain glades (46 hectares)

Water-Resources Investigations Report↗

Structure, metamorphism, and plutonism in the south-central Klamath Mountains, California

In the south-central Klamath Mountains 50 miles of the the north-trending central metamorphic belt and adjacent parts of the eastern Paleozoic and western Paleozoic and Triassic belts have been mapped and studied in detail. Within the central metamorphic belt a sequence of three lithologically distinctive metamorphic units has been recognized (from bottom to top): (1) siliceous metasedimentary rocks and greenstones of the Stuart Fork Formation; (2) the Salmon Hornblende Schist; and (3) siliceous, calcareous, and amphibolitic rocks, predominantly metasedimentary, of the Grouse Ridge Formation. The age of these metamorphic rocks is uncertain; they are known only to predate intrusion of Late Jurassic (Nevadan) granitic rocks. Ultramafic rocks, mainly alpine-type peridotites, were emplaced before the granitic rocks and occur primarily in a single large sheetlike body which separates the central metamorphic belt from the eastern Paleozoic belt. Granitic plutons, including quartz diorites, trondhjemites, granodiorites, diorites, and gabbros, in decreasing order of abundance, range in size from less than 1 to about 80 square miles in area. Two orogenic phases in the central metamorphic belt have been distinguished by structural and textural features. A late deformation uniformly affected the metamorphic terrane and the ultramafic rocks but predated granitic rocks. It was accompanied by some metamorphism in the lower to middle greenchist facies and produced upright folds that trend south and plunge gently. An earlier phase affected the various rock units differentially; it produced widespread recumbent folding and upper greenschist- to amphibolite-facies metamorphism in Salmon and Grouse Ridge rocks, but involved the underlying Stuart Fork Formation less severely, producing at least local recumbent folding and lower greenschist-facies metamorphism. The preferred interpretation of this upward increase in structural complexity and metamorphic grade is that the Salmon-Grouse Ridge sequence is a thrust sheet which overrode the Stuart Fork rocks concurrently with emplacement of the ultramafic rocks during the culmination of early recumbent folding and metamorphism. Thrusting was then followed by upright folding during the waning stages of regional metamorphism. The first deformational phase, and possibly the second, occurred during late Paleozoic time as indicated by recent isotopic ages of Salmon Hornblende Schist.

California↗

Assessment of arsenic concentrations in domestic well water, by town, in Maine 2005-09

Prior studies have established that approximately 10 percent of domestic wells in Maine have arsenic levels greater than the U.S. Environmental Protection Agency maximum contaminant limit (10 micrograms per liter (ug/L)). Of even greater concern are multiple discoveries of wells with very high arsenic levels (> 500 ug/L) in several areas of the State. A study was initiated to assist the Maine Center for Disease Control and Prevention (ME-CDC) in developing a better understanding of the statewide spatial occurrence of wells with elevated arsenic levels at the individual town level, identify areas of the State that should be targeted for increased efforts to promote well-water testing, and generate data for potential use in predicting areas of the State likely to have very high levels of arsenic. The State's Health and Environmental and Testing Laboratory (HETL) annually analyzes samples from thousands of domestic wells for arsenic. Results of arsenic analyses of domestic well water submitted to the HETL from 2005 to 2009 were screened and organized, by town, in order to summarize the results for all towns with samples submitted to the HETL. In order to preserve the privacy of well owners, the screening and organization of samples was conducted in the offices of the ME-CDC, following applicable Maine and United States laws, rules, and privacy policies. After screening, the database contained samples from 531 towns in Maine and from 11,111 individual wells. Of those towns, 385 had samples from 5 or more individual wells, 174 towns had samples from 20 or more individual wells, and 49 towns had samples from 60 or more wells. These samples, because they were submitted by homeowners and were not part of a random sample, may not be representative of all wells in a given area. The minimum, maximum, and median arsenic values for the towns with five or more samples were calculated, and the maximum and median values were mapped for the State. The percentages of samples exceeding 10, 50, 100, and 500 ug/L were calculated for the 174 towns with 20 or more sampled wells, and statewide maps were prepared for each of these categories. More than 25 percent of the sampled wells in 44 towns exceeded 10 ug/L. Many fewer towns had wells with samples that exceeded the 50, 100, or 500 ug/L categories. For 19 towns, more than 10 percent of the sampled wells had arsenic concentrations that exceeded 50 ug/L, and in 45 towns, 1 percent or more exceeded 100 ug/L. Of these, Surry in Hancock County had 120 wells tested, and 23 percent of those wells had arsenic concentrations that exceeded 100 ug/L, which is a much higher rate than for other towns. In only four towns (Danforth in Washington County, Surry and Blue Hill in Hancock County, and Woolwich in Sagadahoc County), 1 percent or more of the sampled wells had arsenic concentrations greater than 500 ug/L during 2005-09. The distribution of high arsenic concentrations in wells follows some geographic patterns, which are generally geologically controlled. There are clusters or belts of towns with high arsenic concentrations (> 50 ug/L), such as in southern coastal areas, the Kennebec County area, and towns along the central coastal part of Maine. In contrast, there are areas of the State with low arsenic concentrations, such as the northernmost towns, as well as towns in the western and west-central areas. There appear to be three distinct large-scale areas of high concentrations of arsenic in groundwater-one in southern coastal areas, one in central Kennebec County, and one in the town of Ellsworth (Hancock County) and the surrounding areas. In addition, several smaller clusters of isolated high concentrations of arsenic in groundwater exist. Earlier testing has identified other clusters of very high arsenic concentrations in groundwater in the towns of Northport, Buxton/Hollis, and Waldoboro, but those samples were collected before 2005 and did not factor in this analysis.

Scientific Investigations Report↗

Hydrologic models and analysis of water availability in Cuyama Valley, California

Changes in population, agricultural development practices (including shifts to more water-intensive crops), and climate variability are placing increasingly larger demands on available water resources, particularly groundwater, in the Cuyama Valley, one of the most productive agricultural regions in Santa Barbara County. The goal of this study was to produce a model capable of being accurate at scales relevant to water management decisions that could be considered in the evaluation of the sustainable water supply. The Cuyama Valley Hydrologic Model (CUVHM) was designed to simulate the most important natural and human components of the hydrologic system, including components dependent on variations in climate, thereby providing a reliable assessment of groundwater conditions and processes that can inform water users and help to improve planning for future conditions. Model development included a revision of the conceptual model of the flow system, construction of a precipitation-runoff model using the Basin Characterization Model (BCM), and construction of an integrated hydrologic flow model with MODFLOW-One-Water Hydrologic Flow Model (MF-OWHM). The hydrologic models were calibrated to historical conditions of water and land use and, then, used to assess the use and movement of water throughout the Valley. These tools provide a means to understand the evolution of water use in the Valley, its availability, and the limits of sustainability. The conceptual model identified inflows and outflows that include the movement and use of water in both natural and anthropogenic systems. The groundwater flow system is characterized by a layered geologic sedimentary sequence that—in combination with the effects of groundwater pumping, natural recharge, and the application of irrigation water at the land surface—displays vertical hydraulic-head gradients. Overall, most of the agricultural demand for water in the Cuyama Valley in the initial part of the growing season is supplied by groundwater, which is augmented by precipitation during wet winter and spring seasons. In addition, the amount of groundwater used for irrigation varies from year to year in response to climate variation and can increase dramatically in dry years. Model simulation results, however, also indicated that irrigation may have been less efficient during wet years. Agricultural pumpage is a major component to simulated outflow that is often poorly recorded. Therefore, an integrated, coupled farm-process model is used to estimate historical pumpage for water-balance subregions that evolved with the development of groundwater in the Valley from 1949 through 2010. The integrated hydrologic model includes these water-balance subregions and delineates natural, municipal, and agricultural land use; streamflow networks; and groundwater flow systems. The redefinition of the geohydrologic framework (including the internal architecture of the sedimentary units) and incorporation of these units into the simulation of the regional groundwater flow system indicated that faults have compartmentalized the alluvial deposits into subregions, which have responded differently to regional groundwater flow, locations of recharge, and the effects of development. The Cuyama Valley comprises nine subregions grouped into three regional zones, the Main, Ventucopa Uplands, and Sierra Madre Foothills, which are fault bounded, represent different proportions of the three alluvial aquifers, and have different water quality. The CUVHM uses MF-OWHM to simulate and assess the use and movement of water, including the evolution of land use and related water-balance regions. The model is capable of being accurate at annual to interannual time frames and at subregional to valley-wide spatial scales, which allows for analysis of the groundwater hydrologic budget for the water years 1950–2010, as well as potential assessment of the sustainable use of groundwater. Simulated changes in storage over time showed that significant withdrawals from storage generally occurred not only during drought years (1976–77 and 1988–92) but also during the early stages of industrial agriculture, which was initially dominated by alfalfa production. Since the 1990s, agriculture has shifted to more water-intensive crops. Measured and simulated groundwater levels indicated substantial declines in selected subregions, mining of groundwater that is thousands to tens of thousands of years old, increased groundwater storage depletion, and land subsidence. Most of the recharge occurs in the upland regions of Ventucopa and Sierra Madre Foothills, and the largest fractions of pumpage and storage depletion occur in the Main subregion. The long-term imbalance between inflows and outflows resulted in simulated overdraft (groundwater withdrawals in excess of natural recharge) of the groundwater basin over the 61-year period of 1949–2010. Changes in storage varied considerably from year to year, depending on land use, pumpage, and climate conditions. Climatically driven factors can greatly affect inflows, outflows, and water use by more than a factor of two between wet and dry years. Although precipitation during inter-decadal wet years previously replenished the basin, the water use and storage depletion have lessened the effects of these major recharge events. Simulated and measured water-level altitudes indicated the presence of large areas where depressed water levels have resulted in large desaturated zones in the younger and Older Alluvium layers in the Main-zone subregions. The results of modeled projection of the base-case scenario 61 years into the future indicated that current supply-and-demand are unsustainable and will result in additional groundwater-level declines and related storage depletion and land subsidence. The reduced-supply and reduced-demand projections reduced groundwater storage depletion but may not allow for sustainable agriculture under current demands, agricultural practices, and land use.

California↗

Hydrology of the San Luis Valley, south-central Colorado

An investigation of the water resources of the Colorado part of the San Luis Valley was begun in 1966 by the U.S. Geological Survey, in cooperation with the Colorado Water Conservation Board. (See index map, fig. 1). The purpose of the investigation is to provide information for planning and implementing improved water-development and management practices. The major water problems in the San Luis Valley include (1) waterlogging, (2) waste of water by nonbeneficial evapotranspiration, (3) deterioration of ground-water chemical quality, and (4) failure of Colorado to deliver water to New Mexico and Texas in accordance with the Rio Grande Compact. This report describes the hydrologic environment, extent of water-resource development, and some of the problems related to that development. Information presented is based on data collected from 1966 to 1968 and on previous studies. Subsequent reports are planned as the investigation progresses. The San Luis Valley extends about 100 miles from Poncha Pass near the northeast corner of Saguache County, Colo., to a point about 16 miles south of the Colorado-New Mexico State line. The total area is 3,125 square miles, of which about 3,000 are in Colorado. The valley is nearly flat except for the San Luis Hills and a few other small areas. The Colorado part of the San Luis Valley, which is described in this report, has an average altitude of about 7,700 feet. Bounding the valley on the west are the San Juan Mountains and on the east the Sangre de Cristo Mountains. Most of the valley floor is bordered by alluvial fans deposited by streams originating in the mountains, the most extensive being the Rio Grande fan (see block diagram, fig. 2 in pocket). Most of the streamflow is derived from snowmelt from 4,700 square miles of watershed in the surrounding mountains. The northern half of the San Luis Valley is internally drained and is referred to as the closed basin. The lowest part of this area is known locally as the "sump." The remainder of the valley is drained by the Rio Grande and its tributaries. The climate of the San Luis Valley is arid, and a successful agricultural economy would not be possible without irrigation. It is characterized by cold winters, moderate summers, and much sunshine. The average annual precipitation on the valley floor ranges from 7 to 10 inches. More than half the precipitation occurs from July to September. Moisture deficiency in the valley is shown by the graph comparing pan evaporation and precipitation {fig. 3}. For the years 1961-67 average pan evaporation for the period April through September was 52.25 inches, but average precipitation for the period was only 5.02 inches. Average annual precipitation was 7.8 inches. Owing to the short growing season (90-120 days), crops a.re restricted mainly to barley, oats, potatoes, and other vegetables.

Colorado↗

Total Phosphorus Loads for Selected Tributaries to Sebago Lake, Maine

The streamflow and water-quality datacollection networks of the Portland Water District (PWD) and the U.S. Geological Survey (USGS) as of February 2000 were analyzed in terms of their applicability for estimating total phosphorus loads for selected tributaries to Sebago Lake in southern Maine. The long-term unit-area mean annual flows for the Songo River and for small, ungaged tributaries are similar to the long-term unit-area mean annual flows for the Crooked River and other gaged tributaries to Sebago Lake, based on a regression equation that estimates mean annual streamflows in Maine. Unit-area peak streamflows of Sebago Lake tributaries can be quite different, based on a regression equation that estimates peak streamflows for Maine. Crooked River had a statistically significant positive relation (Kendall's Tau test, p=0.0004) between streamflow and total phosphorus concentration. Panther Run had a statistically significant negative relation (p=0.0015). Significant positive relations may indicate contributions from nonpoint sources or sediment resuspension, whereas significant negative relations may indicate dilution of point sources. Total phosphorus concentrations were significantly larger in the Crooked River than in the Songo River (Wilcoxon rank-sum test, p<0.0001). Evidence was insufficient, however, to indicate that phosphorus concentrations from medium-sized drainage basins, at a significance level of 0.05, were different from each other or that concentrations in small-sized drainage basins were different from each other (Kruskal-Wallis test, p= 0.0980, 0.1265). All large- and medium-sized drainage basins were sampled for total phosphorus approximately monthly. Although not all small drainage basins were sampled, they may be well represented by the small drainage basins that were sampled. If the tributaries gaged by PWD had adequate streamflow data, the current PWD tributary monitoring program would probably produce total phosphorus loading data that would represent all gaged and ungaged tributaries to Sebago Lake. Outside the PWD tributary-monitoring program, the largest ungaged tributary to Sebago Lake contains 1.5 percent of the area draining to the lake. In the absence of unique point or nonpoint sources of phosphorus, ungaged tributaries are unlikely to have total phosphorus concentrations that differ significantly from those in the small tributaries that have concentration data. The regression method, also known as the rating-curve method, was used to estimate the annual total phosphorus load for Crooked River, Northwest River, and Rich Mill Pond Outlet for water years 1996-98. The MOVE.1 method was used to estimate daily streamflows for the regression method at Northwest River and Rich Mill Pond Outlet, where streamflows were not continuously monitored. An averaging method also was used to compute annual loads at the three sites. The difference between the regression estimate and the averaging estimate for each of the three tributaries was consistent with what was expected from previous studies.

Water-Resources Investigations Report↗

Numerical simulation of streamflow distribution, sediment transport, and sediment deposition along Long Beach Creek in Northeast Missouri

This report presents the results of a study conducted by the U.S. Geological Survey in cooperation with the Missouri Department of Conservation to describe the hydrology, sediment transport, and sediment deposition along a selected reach of Long Branch Creek in Macon County, Missouri. The study was designed to investigate spatial and temporal characteristics of sediment deposition in a remnant forested riparian area and compare these factors by magnitude of discharge events both within and outside the measured range of flood magnitudes. The two-dimensional finite-element numerical models RMA2-WES and SED2D-WES were used in conjunction with measured data to simulate streamflow and sediment transport/deposition characteristics during 2-, 5-, 10-, and 25-year recurrence interval floods. Spatial analysis of simulated sediment deposition results indicated that mean deposition in oxbows and secondary channels exceeded that of the remaining floodplain areas during the 2-, 5-, 10-, and 25-year recurrence interval floods. The simulated mass deposition per area for oxbows and secondary channels was 1.1 to 1.4 centimeters per square meter compared with 0.1 to 0.60 centimeters per square meter for the remaining floodplain. The temporal variability of total incremental floodplain deposition during a flood was found to be strongly tied to sediment inflow concentrations. Most floodplain deposition, therefore, occurred at the beginning of the streamflow events and corresponded to peaks in sediment discharge. Simulated total sediment deposition in oxbows and secondary channels increased in the 2-year through 10-year floods and decreased in the 25- year flood while remaining floodplain deposition was highest for the 25-year flood. Despite increases in sediment inflows from the 2-year through 25-year floods, the retention ratio of sediments (the ratio of floodplain deposition to inflow load) was greatest for the 5-year flood and least for the 25-year flood. The decrease in retention ratio at greater flows is likely the result of higher velocities on the floodplain, resulting in higher bed shear stress, greater suspension time of deposited material, and greater sediment transport through the system. Simulated sediment deposition was most sensitive to sediment inflow concentrations and modification of floodplain roughness—factors that can be controlled through management practices. The increase in floodplain sediment deposition resulting from a simulated increase in vegetation density (increase in floodplain roughness from a Manning's n of 0.11 to 0.12) was 142,000 kilograms, or 6.5 percent for a 10-year recurrence interval flood. This increase was comparable to total oxbow and secondary channel deposition mass in the simulations, but would result in a mean increase in floodplain deposition thickness of only 0.025 centimeter. The hydrodynamic model results show the importance of the secondary channels and meander cutoff channels in this system because these areas quickly bring floodwaters and sediment to areas not close to the main channel. The meander cutoff channels in the simulation also effectively decrease flow and velocities in some main channel sections thereby affecting sediment deposition in the vicinity of these features.

Missouri↗

Streamflow Simulations and Percolation Estimates Using the Soil and Water Assessment Tool for Selected Basins in North-Central Nebraska, 1940-2005

The U.S. Geological Survey, in cooperation with the Upper Elkhorn, Lower Elkhorn, Upper Loup, Lower Loup, Middle Niobrara, Lower Niobrara, Lewis and Clark, and Lower Platte North Natural Resources Districts, used the Soil and Water Assessment Tool to simulate streamflow and estimate percolation in north-central Nebraska to aid development of long-term strategies for management of hydrologically connected ground and surface water. Although groundwater models adequately simulate subsurface hydrologic processes, they often are not designed to simulate the hydrologically complex processes occurring at or near the land surface. The use of watershed models such as the Soil and Water Assessment Tool, which are designed specifically to simulate surface and near-subsurface processes, can provide helpful insight into the effects of surface-water hydrology on the groundwater system. The Soil and Water Assessment Tool was calibrated for five stream basins in the Elkhorn-Loup Groundwater Model study area in north-central Nebraska to obtain spatially variable estimates of percolation. Six watershed models were calibrated to recorded streamflow in each subbasin by modifying the adjustment parameters. The calibrated parameter sets were then used to simulate a validation period; the validation period was half of the total streamflow period of record with a minimum requirement of 10 years. If the statistical and water-balance results for the validation period were similar to those for the calibration period, a model was considered satisfactory. Statistical measures of each watershed model's performance were variable. These objective measures included the Nash-Sutcliffe measure of efficiency, the ratio of the root-mean-square error to the standard deviation of the measured data, and an estimate of bias. The model met performance criteria for the bias statistic, but failed to meet statistical adequacy criteria for the other two performance measures when evaluated at a monthly time step. A primary cause of the poor model validation results was the inability of the model to reproduce the sustained base flow and streamflow response to precipitation that was observed in the Sand Hills region. The watershed models also were evaluated based on how well they conformed to the annual mass balance (precipitation equals the sum of evapotranspiration, streamflow/runoff, and deep percolation). The model was able to adequately simulate annual values of evapotranspiration, runoff, and precipitation in comparison to reported values, which indicates the model may provide reasonable estimates of annual percolation. Mean annual percolation estimated by the model as basin averages varied within the study area from a maximum of 12.9 inches in the Loup River Basin to a minimum of 1.5 inches in the Shell Creek Basin. Percolation also varied within the studied basins; basin headwaters tended to have greater percolation rates than downstream areas. This variance in percolation rates was mainly was because of the predominance of sandy, highly permeable soils in the upstream areas of the modeled basins.

Scientific Investigations Report↗

Origin of concretionary Mn-Fe-oxides in stream sediments of Maine, U.S.A.

Studies of stream and sediment-pore waters largely explain the genesis of concretionary Mn-Fe-oxides in Maine. Waters of two small streams near Jackman, Maine, were studied in terms of pH, Eh, dissolved oxygen, dissolved organic carbon, dissolved Mn, total dissolved Fe, and ferrous and ferric Fe. Pyrite Creek has profuse concretions and coatings of Mn-Fe-oxides, whereas West Pyrite Creek has only sparse Mn-Fe-oxide stains. Pyrite Creek drains boggy terrain and West Pyrite Creek drains well-drained terrain. In West Pyrite Creek, stream and subjacent pore waters have chemical characteristics that do not differ greatly. However, dissolved Mn, ferrous Fe, dissolved oxygen, and in situ Eh measurements show that a steep Eh gradient exists between stream and subjacent pore waters of Pyrite Creek. The steep Eh gradient is manifested by the common zonation of coatings and stains on rocks in stream sediment. The bottom zone has no deposition of oxides, the middle zone is red and consists mostly of Fe-oxides, and the upper zone is black or dark-brown and consists of Mn-oxides with varying amounts of Fe-oxides. The zonation agrees with theoretical predictions of oxide stability as one moves from a reducing to an oxidizing environment. At locations where concretionary Mn-Fe-oxides form, pore waters are depleted of oxygen because of abundant decaying organic material in the stream sediment. The pore waters are charged with dissolved Mn and Fe because mechanically deposited Mn-Fe-oxides are remobilized due to the low-Eh conditions. Groundwaters also contribute dissolved Mn and Fe. Stream waters, on the other hand, are oxygenated and the high-Eh conditions result in low concentrations of dissolved Mn and Fe in stream waters because of the insolubility of Mn-Fe-oxides in high-Eh environments. Therefore, concretionary Mn-Fe-oxides form at the interface between pore and stream waters because Mn- and Fe-rich pore waters, which are undersaturated with respect to Mn-Fe-oxides, mix with oxygen-rich stream waters, which are saturated with respect to Mn-Fe-oxides. ?? 1983.

Chemical Geology↗

Lead isotope compositions as guides to early gold mineralization: The North Amethyst vein system, Creede district, Colorado

The North Amethyst vein system, which is hosted by approximately 27 Ma Carpenter Ridge Tuff and approximately 26 Ma Nelson Mountain Tuff, has two mineral associations separated by brecciation and sedimentation in the veins. The early association consists of quartz, rhodonite, hematite, magnetite, electrum (Au (sub 0.3-0.5) Ag (sub 0.7-0.5)) , and Mn carbonate, Au-Ag sulfide, Ag sulfosalt, and base metal sulfide minerals. The later mineral association cuts the Mn- and Au-bearing assemblages and consists of quartz, calcite, sericite, chlorite, hematite, adularia, fluorite, base metal sulfides, and Ag-bearing tetrahedrite.Our detailed studies show that the Pb isotope compositions of paragenetically early galenas associated with Au-rich mineralization in the North Amethyst vein system are relatively unradiogenic ( 206 Pb/ 204 Pb: 18.826-18.881, 207 Pb/ 204 Pb: 15.588-15.602, and 208 Pb/ 204 Pb: 37.790-37.926) compared to Pb isotope compositions of galenas formed later at about 25 Ma during Ag and base metal mineralization ( 206 Pb/ 204 Pb: 19.041-19.115, 207 Pb/ 204 Pb: 15.627- 15.672, and 208 Pb/ 204 Pb: 37.829-38.057). New Pb isotope data for the central and southern parts of the Creede district, which are located 5 to 7 km south of the North Amethyst area, agree with the results of a regional study by Doe et al. (1979) that included five samples from the main part of the Creede district. Galenas and adularia from the central and southern Creede district and galenas from the Bondholder district, 7 km north of the North Amethyst area, are similar to the Pb isotope compositions of galenas formed later in the North Amethyst area during Ag and base metal mineralization. Galenas from the Alpha-Corsair vein, which was mined for Ag prior to 1910, are isotopically similar to galenas associated with North Amethyst Au-stage mineralization. This isotopic similarity suggests that unexplored segments of the Alpha-Corsair structure may have the same mineralogy as the North Amethyst Au stage; thus, the Alpha-Corsair structure has the potential for high gold contents.Pb isotope compositions from the late stage of the North Amethyst vein system and from the Bondholder and central and southern Creede mining districts are more radiogenic than the host volcanic rocks of the central cluster of the San Juan volcanic field. Our Pb isotope results indicate that early Au mineralization of the North Amethyst area may represent the product of an older and relatively local hydrothermal system distinct from that of the younger base metal and Ag mineralization found throughout the region. Fluids that deposited Au minerals may have derived their Pb isotope composition by a greater degree of interaction with shallow, relatively less radiogenic volcanic wall rocks. The younger, base metal and Ag-rich mineralization that overprints the Au mineralization in the North Amethyst area clearly has a more radiogenic isotopic signature, which implies that the later mineralization derived a greater component of its Pb from Proterozoic source rocks, or sediments derived from them.Paragenetically early sulfide-rich vein assemblages have the least radiogenic galenas and generally also have the highest Au contents. Thus, identification of paragenetically early vein assemblages with relatively unradiogenic Pb isotope compositions similar to those of the North Amethyst area provides an additional exploration tool for Au in the central San Juan Mountains area.

Colorado↗

Availability and distribution of low flow in Anahola Stream, Kauaʻi, Hawaiʻi

Anahola Stream is a perennial stream in northeast Kauaʻi, Hawaiʻi, that supports agricultural, domestic, and cultural uses within its drainage basin. Beginning in the late 19th century, Anahola streamflow was diverted by Makee Sugar Company at altitudes of 840 feet (upper intake) and 280 feet (lower intake) for irrigating sugarcane in the Keālia area. When sugarcane cultivation in the Keālia area ceased in 1988, part of the Makee Sugar Company’s surface-water collection system (Makee diversion system) in the Anahola drainage basin was abandoned. In an effort to better manage available surface-water resources, the State of Hawaiʻi Department of Hawaiian Home Lands is considering using the existing ditches in the Anahola Stream drainage basin to provide irrigation water for Native Hawaiian farmers in the area. To provide information needed for successful management of the surface-water resources, the U.S. Geological Survey investigated the availability and distribution of natural low flow in Anahola Stream and also collected low-flow data in Goldfish Stream, a stream that discharges into Kaneha Reservoir, which served as a major collection point for the Makee diversion system. Biological surveys of Anahola Stream were conducted as part of a study to determine the distribution of native and nonnative aquatic stream fauna. Results of the biological surveys indicated the presence of the following native aquatic species in Anahola Stream: ʻoʻopu ʻakupa (Sandwich Island sleeper) and ʻoʻopu naniha (Tear-drop goby) in the lower stream reaches surveyed; and ʻoʻopu nākea (Pacific river goby), ʻoʻopu nōpili (Stimpson’s goby), and ʻōpae kalaʻole (Mountain shrimp) in the middle and upper stream reaches surveyed. Nonnative aquatic species were found in all of the surveyed stream reaches along Anahola Stream. The availability and distribution of natural low flow were determined using a combination of discharge measurements made from February 2011 to May 2012 at low-flow partial-record and seepage-run stations established at locations of interest along study-area streams. Upstream of the upper intake, the estimated natural (undiverted) median flow in Anahola Stream is 2.7 million gallons per day, and the flow is expected to be greater than or equal to 0.97 million gallons per day 95 percent of the time. About 0.7 mile upstream of the lower intake and downstream from the confluence with Keaʻoʻopu Stream, the estimated natural (undiverted) median flow in Anahola Stream is 6.3 million gallons per day, and the flow is expected to be greater than or equal to 2.7 million gallons per day 95 percent of the time. In Goldfish Stream, about 0.4 mile upstream from the point of discharge into Kaneha Reservoir, the estimated natural median flow is 0.54 million gallons per day, and the flow is expected to be greater than or equal to 0.23 million gallons per day 95 percent of the time. The discharge estimates are representative of low-flow conditions in the study-area streams, and they are applicable to the base period (water years 1961–2011) over which they have been computed. The distribution of natural low flow in Anahola Stream was characterized through data collected during wet- and dry-season seepage runs. Seepage-run results show that Anahola Stream was generally a gaining stream under natural low-flow conditions. During the wet-season seepage run, Anahola Stream at the station located upstream of tributary Kaʻalula Stream had more than five times the flow that was measured upstream from the upper intake. The estimated total gain (including tributary inflow) in the 6.1-mile seepage-run reach was 6.97 million gallons per day; about 42 percent of that gain was groundwater discharge to the main channel of Anahola Stream. During the dry-season seepage run, about 34 percent of the estimated total gain of 3.93 million gallons per day in the same seepage-run reach was groundwater discharge to the main channel of Anahola Stream. A 15-percent seepage loss was estimated in a 0.3-mile reach downstream from the confluence of Anahola and Keaʻoʻopu Streams. The report summarizes scenarios that describe (1) surface-water availability under regulated conditions of Anahola Stream if the upper and lower intakes are restored in the future; and (2) amount of flow available for agricultural use at the upper intake under a variety of potential instream-flow standards that may be established by the State of Hawaiʻi for the protection of instream uses.

Hawai'i↗

Streamflow, groundwater hydrology, and water quality in the upper Coleto Creek watershed in southeast Texas, 2009&ndash;10

The U.S. Geological Survey (USGS), in cooperation with the Goliad County Groundwater Conservation District, Victoria County Groundwater Conservation District, Pecan Valley Groundwater Conservation District, Guadalupe-Blanco River Authority, and San Antonio River Authority, did a study to examine the hydrology and stream-aquifer interactions in the upper Coleto Creek watershed. Findings of the study will enhance the scientific understanding of the study-area hydrology and be used to support water-management decisions to help ensure protection of the Evangeline aquifer and surface-water resources in the study area. This report describes the results of streamflow measurements, groundwater-level measurements, and water quality (from both surface-water and groundwater sites) collected from three sampling events (July&ndash;August 2009, January 2010, and June 2010) designed to characterize groundwater (from the Evangeline aquifer) and surface water, and the interaction between them, in the upper Coleto Creek watershed upstream from Coleto Creek Reservoir in southeast Texas. This report also provides a baseline level of water quality for the upper Coleto Creek watershed. Three surface-water gain-loss surveys&mdash;July 29&ndash;30, 2009, January 11&ndash;13, 2010, and June 21&ndash;22, 2010&mdash;were done under differing hydrologic conditions to determine the locations and amounts of streamflow recharging or discharging from the Evangeline aquifer. During periods when flow in the reaches of the upper Coleto Creek watershed was common (such as June 2010, when 12 of 25 reaches were flowing) or probable (such as January 2010, when 22 of 25 reaches were flowing), most of the reaches appeared to be gaining (86 percent in January 2010 and 92 percent in June 2010); however, during drought conditions (July 2009), streamflow was negligible in the entire upper Coleto Creek watershed; streamflow was observed in only two reaches during this period, one that receives inflow directly from Audilet Spring and another reach immediately downstream from Audilet Spring. Water levels in the aquifer at this time declined to the point that the aquifer could no longer provide sufficient water to the streams to sustain flow. Groundwater-level altitudes were measured at as many as 33 different wells in the upper Coleto Creek watershed during three different survey events: August 4&ndash;7 and 12, 2009; January 12&ndash;14 and 22, 2010; and June 21&ndash;24, 2010. These data were used in conjunction with groundwater-level altitudes from three continuously monitored wells to generate potentiometric surface maps for each of the three sampling events to help characterize the groundwater hydrology of the Evangeline aquifer. The altitudes of potentiometric surface contours from all three sampling events are highest in the northeast part of the study area and lowest in the southwest part of the study area. Groundwater flow direction shifts from southeast to east across the watershed, roughly coinciding with the general flow direction of the main stem of Coleto Creek. Groundwater-level altitudes increased an average of 2.35 inches between the first and third sampling events as drought conditions in summer 2009 were followed by consistent rains the subsequent fall and winter, an indication that the aquifer responds relatively quickly to both the absence and relative abundance of precipitation. A total of 44 water-quality samples were collected at 21 different sites over the course of the three sampling events (August 4&ndash;7, 2009, January 12&ndash;14, 2010, and June 21&ndash;24, 2010). In most cases, samples from each site were analyzed for the following constituents: dissolved solids, major ions, alkalinity, nutrients, trace elements, and stable isotopes (hydrogen, oxygen, and strontium). Major-ion compositions were relatively consistent among most of the samples from the upper Coleto Creek watershed (generally calcium bicarbonate waters, with chloride often making a major contribution). Of the 23 trace elements that were analyzed in water samples as part of this study, only arsenic (in two samples) and manganese (in seven samples) had concentrations that exceeded public drinking-water standards or guidelines. At 3 of the 19 sites sampled&mdash;State wells 79-06-411, 79-14-204, and Audilet Spring&mdash;nitrate concentrations exceeded the threshold (2.0 milligrams per liter) associated with anthropogenic contributions. The majority of the water samples (36 out of 44) that were analyzed for stable isotopes of hydrogen and oxygen during the three sampling events plotted in a relatively tight cluster centered near the global meteoric water line. The eight remaining samples, which include the four surface-water samples collected in June 2010, the sample collected from Coleto Creek Reservoir in January 2010, and all three samples collected at State well 79-15-904, deviate from the global meteoric water line in a way that indicates evaporative losses. The isotopic signatures of the three samples collected at State well 79-15-904, when taken in conjunction with its proximity to Coleto Creek Reservoir, indicate that there is likely a hydraulic connection between the two. When all of the sites are examined as a whole, there is a general pattern in strontium concentrations across the entire watershed that indicates that both the surface-water and groundwater samples derive from a single source (the Evangeline aquifer) with relatively uniform water-rock interactions.

Texas↗

Variability in eddy sandbar dynamics during two decades of controlled flooding of the Colorado River in the Grand Canyon

Sandbars are iconic features of the Colorado River in the Grand Canyon, Arizona, U.S.A. Following completion of Glen Canyon Dam in 1963, sediment deficit conditions caused erosion of eddy sandbars throughout much of the 360 km study reach downstream from the dam. Controlled floods in 1996, 2004, and 2008 demonstrated that sand on the channel bed could be redistributed to higher elevations, and that floods timed to follow tributary sediment inputs would increase suspended sand concentrations during floods. Since 2012, a new management protocol has resulted in four controlled floods timed to follow large inputs of sand from a major tributary. Monitoring of 44 downstream eddy sandbars, initiated in 1990, shows that each controlled flood deposited significant amounts of sand and increased the size of subaerial sandbars. However, the magnitude of sandbar deposition varied from eddy to eddy, even over relatively short distances where main-stem suspended sediment concentrations were similar. Here, we characterize spatial and temporal trends in sandbar volume and site-scale (i.e., individual eddy) sediment storage as a function of flow, channel, and vegetation characteristics that reflect the reach-scale (i.e., kilometer-scale) hydraulic environment. We grouped the long-term monitoring sites based on geomorphic setting and used a principal component analysis (PCA) to correlate differences in sandbar behavior to changes in reach-scale geomorphic metrics. Sites in narrow reaches are less-vegetated, stage changes markedly with discharge, sandbars tend to remain dynamic, and sand storage change dominantly occurs in the eddy compared to the main channel. In wider reaches, where stage-change during floods may be half that of narrow sites, sandbars are more likely to be stabilized by vegetation, and floods tend to aggrade the vegetated sandbar surfaces. In these locations, deposition during controlled floods is more akin to floodplain sedimentation, and the elevation of sandbar surfaces increases with successive floods. Because many sandbars are intermediate to the end members described above, high-elevation bar surfaces stabilized by vegetation often have a more dynamic unvegetated sandbar on the channel-ward margin that aggrades and erodes in response to controlled flood cycles. Ultimately, controlled floods have been effective at increasing averaged sandbar volumes, and, while bar deposition during floods decreases through time where vegetation has stabilized sandbars, future controlled floods are likely to continue to result in deposition in a majority of the river corridor.

Colorado River, Grand Canyon↗

Principal reference section for part of the Eocene Ghazij Formation, Moghal Mine area, Mach coal field, Balochistan, Pakistan

The information presented on this sheet was collected as part of a joint U.S. Geological Survey-Geological Survey of Pakistan program sponsored by the U.S. Agency for International Development. As a project within this program, the coal-bearing Ghazij Formation (Eocene) was investigated in the northeastern part of Balochistan cast and south of the provincial capital of Quetta. Strata exposed in this area range in age from Permian to Holocene and crop out as a belt of folded and thrusted rocks that form a southeast-facing orocline. In this region of Pakistan, the Ghazij can usually be divided into three parts. The lower part is the thickest (probably more than 1,000 m) and consists of gray-weathering calcareous mudrock (shale, mudstone, and impure claystone) and a few tabular bodies of fine-to medium-grained calcareous sandstone. The middle part (27-300 m) consists of gray-weathering calcareous mudrock and tabular to lenticular bodies of fine-to medium-grained calcareous sandstone; beds of carbonaceous shale and coal are common. The upper part (as thick as 533 m) contains reddish-weathering calcareous mudrock that contains scattered lenticular bodies of fine- to medium-grained calcareous sandstone. Fossil plant debris is common in mudrock of the lower and middle parts of the Ghazij, and bivalves and gastropods are common in the middle part; the upper part of the Ghazij is usually unfossiliferous. This three-fold division of the Ghazij is less distinct in the Johan area. Here, the upper part of the formation is clearly identifiable, but rocks below it are poorly exposed and assigning a stratigraphic level that separates the middle and lower parts of the formation is problematic. Below the upper part of the formation is a thick sequence of greenish-gray calcareous mudrock that contains locally abundant plant debris and isolated bodies of brown-weathering sandstone. Rare carbonaceous shale and even rarer coal are present in the upper part of this sequence, and this interval of the formation might correspond to the middle part of the Ghazji exposed in areas to the north. We propose that, in the Johan area, those rocks below the upper part of the formation be referred to as the main body of the Ghazij (for example, main-body Ghazij). Underlying the Ghazij are the carbonate rocks of the Paleocene Dungan Formation (or its equivalent), and overlying the Ghazij are the mostly carbonate rocks of the Eocene Kirthar Formation (or its equivalent). Both contacts can be conformable or unconformable. All of the pre-Neogene rocks in Balochistan are greatly deformed by the collision of India and Asia. The Ghazij is especially susceptible to regional compressional tectonics because it contains a large amount of shale and is sandwiched between two thick carbonate units. As a result, bedding-plane faults and isoclinal folds are common. As part of our study of the Ghazij Formation, five stratigraphic sections were measured: one near Pir Ismail Ziarat, one in the Sor Range, two in the vicinity of Mach, and one near Johan. Each area's section is published separately.

Balochistan↗

Principal reference section for part of the Eocene Ghazij Formation, Sarawan River area, Johan coal field, Balochistan, Pakistan

The information presented on this sheet was collected as part of a joint U.S. Geological Survey-Geological Survey of Pakistan program sponsored by the U.S. Agency for International Development. As a project within this program, the coal-bearing Ghazij Formation (Eocene) was investigated in the northeastern part of Balochistan cast and south of the provincial capital of Quetta. Strata exposed in this area range in age from Permian to Holocene and crop out as a belt of folded and thrusted rocks that form a southeast-facing orocline. In this region of Pakistan, the Ghazij can usually be divided into three parts. The lower part is the thickest (probably more than 1,000 m) and consists of gray-weathering calcareous mudrock (shale, mudstone, and impure claystone) and a few tabular bodies of fine-to medium-grained calcareous sandstone. The middle part (27-300 m) consists of gray-weathering calcareous mudrock and tabular to lenticular bodies of fine-to medium-grained calcareous sandstone; beds of carbonaceous shale and coal are common. The upper part (as thick as 533 m) contains reddish-weathering calcareous mudrock that contains scattered lenticular bodies of fine- to medium-grained calcareous sandstone. Fossil plant debris is common in mudrock of the lower and middle parts of the Ghazij, and bivalves and gastropods are common in the middle part; the upper part of the Ghazij is usually unfossiliferous. This three-fold division of the Ghazij is less distinct in the Johan area. Here, the upper part of the formation is clearly identifiable, but rocks below it are poorly exposed and assigning a stratigraphic level that separates the middle and lower parts of the formation is problematic. Below the upper part of the formation is a thick sequence of greenish-gray calcareous mudrock that contains locally abundant plant debris and isolated bodies of brown-weathering sandstone. Rare carbonaceous shale and even rarer coal are present in the upper part of this sequence, and this interval of the formation might correspond to the middle part of the Ghazji exposed in areas to the north. We propose that, in the Johan area, those rocks below the upper part of the formation be referred to as the main body of the Ghazij (for example, main-body Ghazij). Underlying the Ghazij are the carbonate rocks of the Paleocene Dungan Formation (or its equivalent), and overlying the Ghazij are the mostly carbonate rocks of the Eocene Kirthar Formation (or its equivalent). Both contacts can be conformable or unconformable. All of the pre-Neogene rocks in Balochistan are greatly deformed by the collision of India and Asia. The Ghazij is especially susceptible to regional compressional tectonics because it contains a large amount of shale and is sandwiched between two thick carbonate units. As a result, bedding-plane faults and isoclinal folds are common. As part of our study of the Ghazij Formation, five stratigraphic sections were measured: one near Pir Ismail Ziarat, one in the Sor Range, two in the vicinity of Mach, and one near Johan. Each area's section is published separately.

Balochistan↗