Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Maps”

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 919 records · Page 51Linked to original sources
USGS2010

Reduced channel conveyance on the Wichita River at Wichita Falls, Texas, 1900-2009

Recent floods on the Wichita River at Wichita Falls, Texas, have reached higher stages compared to historical floods of similar magnitude discharges. The U.S. Geological Survey (USGS) has operated streamflow-gaging station 07312500 Wichita River at Wichita Falls, Tex., since 1938 and flood measurements near the location of the present gage were first made in 1900. Floods recorded in 2007 and 2008 at this gaging station, including the record flood of June 30, 2007, reached higher stages compared to historical floods before 1972 of similar peak discharges. For flood measurements made at stages of more than 18 feet, peak stages were about 1 to 3 feet higher compared to peak stages of similar peak discharges measured before 1972. Flood measurements made at stages of more than 18 feet also indicate a decrease in the measured mean velocity from about 3.5 to about 2.0 feet per second from 1941 to 2008. The increase in stage and decrease in streamflow velocity for similar magnitude floods indicates channel conveyance has decreased over time. A study to investigate the causes of reduced channel conveyance in the Wichita River reach from Loop 11 downstream to River Road in Wichita Falls was done by the USGS in cooperation with the City of Wichita Falls. Historical photographs indicate substantial growth of riparian vegetation downstream from Loop 11 between 1950 and 2009. Aerial photographs taken between 1950 and 2008 also indicate an increase in riparian vegetation. Twenty-five channel cross sections were surveyed by the USGS in this reach in 2009. These cross sections were located at bridge crossings or collocated with channel cross sections previously surveyed in 1986 for use in a floodplain mapping study by the Federal Emergency Management Agency. Four channel cross sections 3,400 to 11,900 feet downstream from Martin Luther King Jr. Boulevard indicate narrowing of the channel. The remaining channel cross sections surveyed in 2009 by the USGS compared favorably with cross sections surveyed in 1986 for the Federal Emergency Management Agency, with no substantial differences noted. Comparison of channel cross sections surveyed in 2009 to those from historic bridge plans indicate no change in cross section has occurred at most of the bridges from Loop 11 downstream to River Road in Wichita Falls, except for obstructions noted at the Scott Avenue bridge and Martin Luther King Jr. bridge. Although obstructions in the channel at these bridges only partially block flow, they could also be contributing to reduced channel conveyance. Step-backwater profiles were used by the USGS to verify channel roughness. The main channel roughness coefficients (Manning's n values) from 2009 surveys were virtually unchanged from those used in a 1991 hydraulic model done for the Federal Emergency Management Agency. The average overbank roughness coefficient (Manning's n value) was 0.15, more than double the value of 0.06 used in the 1991 hydraulic model. Increased overbank vegetation has resulted in higher stages conveying the same amount of discharge, particularly for discharges more than 4,000 cubic feet per second.

Texas↗
USGS2018

Population estimates of the Endangered Hawaiʻi ʻĀkepa (Loxops coccineus) in different habitats on windward Mauna Loa

Endangered Hawai‘i ʻĀkepas ( Loxops coccineus ) are endemic to Hawai‘i island, where they occur in five spatially distinct populations. Data concerning the status and population trends of these unique Hawaiian honeycreepers are crucial for assessing the effectiveness of recovery and management actions. In 2016, we used point‐transect distance sampling to estimate the abundance of Hawai‘i ʻĀkepas in portions of Hawai‘i Volcanoes National Park (HAVO) and the Kaʻū Forest Reserve (KFR) on Mauna Loa volcano. We then compiled the survey data from four other populations to provide a global population estimate. In our HAVO and KFR study area, we mapped habitat classes to determine the population densities in each habitat. Densities were highest (1.03 birds/ha) in open‐canopy montane ʻōhiʻa ( Metrosideros polymorpha ) woodland. In contrast, densities of the largest ʻĀkepa population on Mauna Kea volcano were highest in closed‐canopy ʻōhiʻa and koa ( Acacia koa ) forest where the species is dependent on nest cavities in tall (> 15 m), large (> 50‐cm diameter at breast height) trees. We surveyed potential nesting habitat in HAVO and KFR and found only one cavity in the short‐stature montane ʻōhiʻa woodland and five cavities in the tall‐stature forest. Differences in densities between the Mauna Kea and Mauna Loa populations suggest that Hawai‘i ʻĀkepas may exhibit different foraging and nesting behaviors in the two habitats. The estimated overall population density in the HAVO and KFR study area was 0.52 birds/ha, which equates to 3663 (95% CI 1725–6961) birds in their 11,377‐ha population range. We calculated a global population of 16,428 (95% CI 10,065–25,198) birds, which is similar to an estimate of 13,892 (95% CI 10,315–17,469) birds made in 1986. Our results suggest that populations are stable to increasing in the two largest populations, but the three other populations are smaller (range = 77–1443 birds) and trends for those populations are unknown.

Hawai'i↗
USGS2024

RegionGrow3D: A deterministic analysis for characterizing discrete three-dimensional landslide source areas on a regional scale

Regional-scale characterization of shallow landslide hazards is important for reducing their destructive impact on society. These hazards are commonly characterized by (a) their location and likelihood using susceptibility maps, (b) landslide size and frequency using geomorphic scaling laws, and (c) the magnitude of disturbance required to cause landslides using initiation thresholds. Typically, this is accomplished through the use of inventories documenting the locations and triggering conditions of previous landslides. In the absence of comprehensive landslide inventories, physics-based slope stability models can be used to estimate landslide initiation potential and provide plausible distributions of landslide characteristics for a range of environmental and forcing conditions. However, these models are sometimes limited in their ability to capture key mechanisms tied to discrete three-dimensional (3D) landslide mechanics while possessing the computational efficiency required for broad-scale application. In this study, the RegionGrow3D (RG3D) model is developed to broadly simulate the area, volume, and location of landslides on a regional scale (≥1,000 km 2 ) using 3D, limit-equilibrium (LE)-based slope stability modeling. Furthermore, RG3D is incorporated into a susceptibility framework that quantifies landsliding uncertainty using a distribution of soil shear strengths and their associated probabilities, back-calculated from inventoried landslides using 3D LE-based landslide forensics. This framework is used to evaluate the influence of uncertainty tied to shear strength, rainfall scenarios, and antecedent soil moisture on potential landsliding and rainfall thresholds over a large region of the Oregon Coast Range, USA.

Journal of Geophysical Research: Earth Surface↗
USGS2003

Exploration and discovery in Yellowstone Lake: Results from high-resolution sonar imaging, seismic reflection profiling, and submersible studies

Discoveries from multi-beam sonar mapping and seismic reflection surveys of the northern, central, and West Thumb basins of Yellowstone Lake provide new insight into the extent of post-collapse volcanism and active hydrothermal processes occurring in a large lake environment above a large magma chamber. Yellowstone Lake has an irregular bottom covered with dozens of features directly related to hydrothermal, tectonic, volcanic, and sedimentary processes. Detailed bathymetric, seismic reflection, and magnetic evidence reveals that rhyolitic lava flows underlie much of Yellowstone Lake and exert fundamental control on lake bathymetry and localization of hydrothermal activity. Many previously unknown features have been identified and include over 250 hydrothermal vents, several very large (>500 m diameter) hydrothermal explosion craters, many small hydrothermal vent craters (∼1–200 m diameter), domed lacustrine sediments related to hydrothermal activity, elongate fissures cutting post-glacial sediments, siliceous hydrothermal spire structures, sublacustrine landslide deposits, submerged former shorelines, and a recently active graben. Sampling and observations with a submersible remotely operated vehicle confirm and extend our understanding of the identified features. Faults, fissures, hydrothermally inflated domal structures, hydrothermal explosion craters, and sublacustrine landslides constitute potentially significant geologic hazards. Toxic elements derived from hydrothermal processes also may significantly affect the Yellowstone ecosystem.

Wyoming↗
USGS2000

Clostridium perfringens in Long Island Sound sediments: An urban sedimentary record

Clostridium perfringens is a conservative tracer and an indicator of sewage-derived pollution in the marine environment. The distribution of Clostridium perfringens spores was measured in sediments from Long Island Sound, USA, as part of a regional study designed to: (1) map the distribution of contaminated sediments; (2) determine transport and dispersal paths; (3) identify the locations of sediment and contaminant focusing; and (4) constrain predictive models. In 1996, sediment cores were collected at 58 stations, and surface sediments were collected at 219 locations throughout the Sound. Elevated concentrations of Clostridium perfringens in the sediments indicate that sewage pollution is present throughout Long Island Sound and has persisted for more than a century. Concentrations range from undetectable amounts to 15,000 spores/g dry sediment and are above background levels in the upper 30 cm at nearly all core locations. Sediment focusing strongly impacts the accumulation of Clostridium perfringens spores. Inventories in the cores range from 28 to 70,000 spores/cm2, and elevated concentrations can extend to depths of 50 cm. The steep gradients in Clostridium perfringens profiles in muddier cores contrast with concentrations that are generally constant with depth in sandier cores. Clostridium perfringens concentrations rarely decrease in the uppermost sediment, unlike those reported for metal contaminants. Concentrations in surface sediments are highest in the western end of the Sound, very low in the eastern region, and intermediate in the central part. This pattern reflects winnowing and focusing of Clostridium perfringens spores and fine-grained sediment by the hydrodynamic regime; however, the proximity of sewage sources to the westernmost Sound locally enhances the Clostridium perfringens signals.

Connecticut, Massachusetts, New Hampshire, New Jer↗
USGS2000

On the modified Mercalli intensities and magnitudes of the 1811-1812 New Madrid earthquakes

We reexamine original felt reports from the 1811–1812 New Madrid earthquakes and determine revised isoseismal maps for the three principal mainshocks. In many cases we interpret lower values than those assigned by earlier studies. In some cases the revisions result from an interpretation of original felt reports with an appreciation for site response issues. Additionally, earlier studies had assigned modified Mercalli intensity (MMI) values of V-VII to a substantial number of reports that we conclude do not describe damage commensurate with intensities this high. We investigate several approaches to contouring the MMI values using both analytical and subjective methods. For the first mainshock on 02∶15 LT December 16, 1811, our preferred contouring yields M ω7.2–7.3 using the area-moment regressions of Johnston [1996]. For the 08∶00 LT on January 23, 1812, and 03∶45 LT on February 7, 1812, mainshocks, we obtain M ω7.0 and M ω7.4–7.5, respectively. Our magnitude for the February mainshock is consistent with the established geometry of the Reelfoot fault, which all evidence suggests to have been the causative structure for this event. We note that the inference of lower magnitudes for the New Madrid events implies that site response plays a significant role in controlling seismic hazard at alluvial sites in the central and eastern United States. We also note that our results suggest that thrusting may have been the dominant mechanism of faulting associated with the 1811–1812 sequence.

Arkansas, Illinois, Kentucky, Missouri, Tennessee↗
USGS2012

Tectonics and conductivity structures in the Southern Washington Cascades

The tectonic setting of the southern Washington Cascades has been studied with the aid of magnetotelluric (MT) and other geophysical data. The main feature of interest in the geophysical data is a broad high-conductivity anomaly mapped with MT and geomagnetic variation (GMV) data. This anomaly is located roughly within the triangle formed by the volcanoes Mount Rainier, Mount St. Helens, and Mount Adams but exceeds beyond Mount Rainier to the northwest. We interpret the cause of the anomaly to be conductive rocks with resistivities of 1–4 ohm m and thicknesses possibly greater than 15 km. These conductive rocks are found 2–8 km beneath the overlying less conductive volcanic and sedimentary rocks at the surface. Two aeromagnetic lows follow the trend of the conductivity anomaly, and linear belts of strike-slip seismicity are coincident with both these magnetic lows. One of the aeromagnetic lows is coincident with the western margin of the conductivity anomaly. The geophysical data appear to outline a suture zone of probable Eocene age caused by accretion of a large seamount complex (Siletzia) and that may contain large thicknesses of compressed forearc basin and accretionary prism sedimentary rocks of Cretaceous to Eocene age. Part of the shallower conductive rocks may be associated with carbonaceous continental and transitional marine sedimentary rocks of the Puget Group. The contact between the hypothesized compressed basin and the accreted terranes to the west may localize the release of shear stresses in this region of oblique subduction. Several possible explanations for the conductivity anomaly are considered in addition to the compressed basin hypothesis.

Washington↗
USGS2023

Conserving habitat for migratory ungulates: How wide is a migration corridor?

Conserving migratory ungulates relies on the analysis of GPS collar data and associated maps of migration corridors to inform management and policy actions. Current methods for identifying migratory corridors use complex statistical models designed to account for movement uncertainty rather than estimating the amount of space required by animals to migrate. Furthermore, such methods can complicate conservation efforts by producing highly variable corridor widths and non-contiguous corridors that do not fully connect seasonal ranges. To remedy, we propose an intuitive line buffer approach for delineating individual migration corridors that is simple to implement and focuses on the functional corridor widths needed by migratory ungulates. By buffering a line that connects successive GPS locations, we can delineate individual migration corridors with consistent widths that are robust to variable parameters (GPS fix rate, travel speed, tortuosity) and provide contiguous connection between seasonal ranges. Using a combination of expert knowledge, simulation and 10-min GPS collar data collected from mule deer ( Odocoileus hemionus ) and pronghorn ( Antilocapra americana ), we suggest 400–600 m are reasonable estimates of functional migration corridor widths for individuals of those species. Synthesis and applications . Our line buffer approach is intended to simplify migration corridor delineation, improve transparency and encourage a broader discussion of functional corridor widths. These considerations help advance efforts to conserve habitat within migration corridors and prioritize conservation efforts within a single corridor or across multiple corridors.

Journal of Applied Ecology↗
USGS2020

A two-stage step-wise framework for fast optimization of well placement in coalbed methane reservoirs

Coalbed methane (CBM) has emerged as a clean energy resource in the global energy mix, especially in countries such as Australia, China, India and the USA. The economical and successful development of CBM requires a thorough evaluation and optimization of well placement prior to field-scale exploitation. This paper presents a two-stage, step-wise optimization framework to obtain the optimal placement of wells for large-scale development of CBM reservoirs. In the first stage, an optimal uniform well pattern is obtained by optimizing well pattern description parameters with the particle swarm optimization (PSO) algorithm. Subsequently, the location and status (active/inactive) of each well are perturbed and optimized within the patterns through the integration of the generalized pattern search (GPS) algorithm and a quality map (QM) representing the production potential. This framework was tested in a synthetic anthracite CBM reservoir in the Qinshui basin (with high gas content and low permeability) and a real field high volatile bituminous reservoir in the Illinois basin (with low gas content and high permeability). The results show that: (i) significant variations in the net present value (NPV) exist with respect to different uniform well patterns (even for cases where the total number of wells are identical), the optima of which can be efficiently determined by the PSO within 100 numerical simulation runs; (ii) the optimization of well perturbations by the GPS results in a more noticeable improvement in NPVs for the synthetic (12.3%) than for the real field model (4.6%); (iii) for the low permeable synthetic model with narrow optimal well spacings (320 m × 200 m), the contribution of the optimization of well perturbation to the NPV increment is heavily dependent on the uniform well placement solution; (iv) for the high permeable real field model with large optimal well spacings (1300 m × 1300 m), the initial uniform well placement has a very minor effect on the subsequent well perturbation solutions in terms of NPV; (v) the proposed framework significantly outperforms the conventional well-by-well concatenation procedure in terms of computational efficiency, robustness and optimal criteria set for production potential.

International Journal of Coal Geology↗
USGS2012

Geophysical constraints on Washington convergent margin structure

Gravity and magnetic maps of western Washington reveal the lateral structure and fabric of the Washington Coast Range, Puget Basin, and southern Washington Cascade Range. The magnetic and gravity maps show large amplitude positive anomalies associated with the shallow but largely buried section of Washington Coast Range mafic rocks which are separated by negative anomalies over deep sedimentary basins. The positive anomalies indicate that the Coast Range mafic basement extends farther east than previously thought, at least as far east as the longitude of Seattle. Linear and steep gravity and magnetic gradients indicate many unmapped, often buried faults in the Washington Coast Range Province. Magnetic highs are also associated with mapped batholiths in the Cascade arc. Several magnetic highs observed east of the Coast Range rocks and west of these batholiths may be associated with buried Tertiary plutons or ophiolites. Two-dimensional gravity and magnetic modeling constrained with geological and other geophysical data indicate that the Coast Range Province rocks are about 1 km thick at the coast, thickening to as much as 30 km near their postulated eastern edge. A maximum boundary on the average density of the upper 15–20 km of the rocks that compose the Coast Range Province of 2920 kg/m 3 was established by the modeling, suggesting a composition largely of basalt and gabbro with little interbedded sediments. Under these rocks may be mantle or a subduction complex composed of dense mafic, ultramafic, and sedimentary rocks like that proposed to underlie Vancouver Island. Previous gravity models of the Washington margin include lower densities for the proposed subduction complex than for Vancouver Island, suggesting a lower component of mafic and ultramafic rocks than the rocks underlying Vancouver Island. However, my Washington model requires that the proposed subduction complex be more dense than the trench sediments and, therefore, that material denser than sediments be incorporated within it. The absence of continental mantle and the modeled wedge shape of the Coast Range Province upper crust suggest that erosion of the bottom of the overriding plate by subduction processes may have occurred.

Washington↗
USGS2013

Consideration of vertical uncertainty in elevation-based sea-level rise assessments: Mobile Bay, Alabama case study

The accuracy with which coastal topography has been mapped directly affects the reliability and usefulness of elevationbased sea-level rise vulnerability assessments. Recent research has shown that the qualities of the elevation data must be well understood to properly model potential impacts. The cumulative vertical uncertainty has contributions from elevation data error, water level data uncertainties, and vertical datum and transformation uncertainties. The concepts of minimum sealevel rise increment and minimum planning timeline, important parameters for an elevation-based sea-level rise assessment, are used in recognition of the inherent vertical uncertainty of the underlying data. These concepts were applied to conduct a sea-level rise vulnerability assessment of the Mobile Bay, Alabama, region based on high-quality lidar-derived elevation data. The results that detail the area and associated resources (land cover, population, and infrastructure) vulnerable to a 1.18-m sea-level rise by the year 2100 are reported as a range of values (at the 95% confidence level) to account for the vertical uncertainty in the base data. Examination of the tabulated statistics about land cover, population, and infrastructure in the minimum and maximum vulnerable areas shows that these resources are not uniformly distributed throughout the overall vulnerable zone. The methods demonstrated in the Mobile Bay analysis provide an example of how to consider and properly account for vertical uncertainty in elevation-based sea-level rise vulnerability assessments, and the advantages of doing so.

Alabama↗
USGS2023

Mapping landslide susceptibility over large regions with limited data

Landslide susceptibility maps indicate the spatial distribution of landslide likelihood. Modeling susceptibility over large or diverse terrains remains a challenge due to the sparsity of landslide data (mapped extent of known landslides) and the variability in triggering conditions. Several different data sampling strategies of landslide locations used to train a susceptibility model are used to mitigate this challenge. However, to our knowledge, no study has systematically evaluated how different sampling strategies alter a model's predictor effects (i.e., how a predictor value influences the susceptibility output) critical to explaining differences in model outputs. Here, we introduce a statistical framework that examines the variation in predictor effects and the model accuracy (measured using receiver operator characteristics) to highlight why certain sampling strategies are more effective than others. Specifically, we apply our framework to an array of logistic regression models trained on landslide inventories collected at sub-regional scales over four terrains across the United States. Results show significant variations in predictor effects depending on the inventory used to train the models. The inconsistent predictor effects cause low accuracies when testing models on inventories outside the domain of the training data. Grouping test and training sets according to physiographic and ecological characteristics, which are thought to share similar triggering mechanisms, does not improve model accuracy. We also show that using limited landslide data distributed uniformly over the entire modeling domain is better than using dense but spatially isolated data to train a model for applications over large regions.

California, Kentucky, North Carolina, West Virgini↗
USGS2005

Geology of five small Australian impact craters

Here we present detailed geological maps and cross-sections of Liverpool, Wolfe Creek, Boxhole, Veevers and Dalgaranga craters. Liverpool crater and Wolfe Creek Meteorite Crater are classic bowlshaped, Barringer-type craters, Liverpool was likely formed during the Neoproterozoic and was filled and covered with sediments soon thereafter. In the Cenozoic, this cover was exhumed exposing the crater's brecciated wall rocks. Wolfe Creek Meteorite Crater displays many striking features, including well-bedded ejecta units, crater-floor faults and sinkholes, a ringed aeromagnetic anomaly, rim-skirting dunes, and numerous iron-rich shale balls. Boxhole Meteorite Crater, Veevers Meteorite Crater and Dalgaranga crater are smaller, Odessa-type craters without fully developed, steep, overturned rims. Boxhole and Dalgaranga craters are developed in highly follated Precambrian basement rocks with a veneer of Holocene colluvium. The pre-existing structure at these two sites complicates structural analyses of the craters, and may have influenced target deformation during impact. Veevers Meteorite Crater is formed in Cenozoic laterites, and is one of the best-preserved impact craters on Earth. The craters discussed herein were formed in different target materials, ranging from crystalline rocks to loosely consolidated sediments, containing evidence that the impactors struck at an array of angles and velocities. This facilitates a comparative study of the influence of these factors on the structural and topographic form of small impact craters. ?? Geological Society of Australia.

Boxhole crater, Dalgaranga crater, Liverpool crate↗
USGS1996

Uncertainty of groundwater vulnerability assessments for agricultural regions in Hawaii: Review

There are important challenges associated with assessing potential groundwater vulnerability hazards that may result from regional scale applications of agrochemicals. The increasing availability of Geographic Information System (GIS) software to those involved in assisting with landuse decisions has resulted in the widespread production of multicolored risk management maps for many environmentally sensitive issues. Soil-based GIS's have recently been coupled to various solute-leaching models to make near-surface groundwater vulnerability assessments for guidance in pesticide regulation in several states. In general, these assessments rest on soil, climatic, and chemical data that are extremely sparse and contain considerable uncertainty. It is also important to acknowledge the uncertainty associated with the transport/fate processes that are not accounted for by the modeling approach used to make the assessment. In this paper, we review the results from a series of papers that have focused on characterization of uncertainty in pesticide mobility estimates, using the attenuation and retardation indices (AF and RF), for the Pearl Harbor Basin on the Hawaiian island of Oahu. Relative to data error uncertainties, we discuss the impacts of: (i) soil, climatic, and chemical data base uncertainties, (ii) reductions in data base uncertainties, (iii) extrapolation of soil data base information based on soil taxonomy and soil survey, and (iv) importing information from outside the region of interest. Relative to model error uncertainties, we compare pesticide leaching estimates from the simple AF and RF mobility indices with simulations from the EPA's Pesticide Root Zone Model (PRZM) and field observations. Finally, we outline a Regional Integrated Risk Assessment approach for characterizing regional scale groundwater vulnerability for near-surface nonpoint sources.

Hawaii↗
USGS2008

Mapping and interpretation of Sinlap crater on Titan using Cassini VIMS and RADAR data

Only a few impact craters have been unambiguously detected on Titan by the Cassini-Huygens mission. Among these, Sinlap is the only one that has been observed both by the RADAR and VIMS instruments. This paper describes observations at centimeter and infrared wavelengths which provide complementary information about the composition, topography, and surface roughness. Several units appear in VIMS false color composites of band ratios in the Sinlap area, suggesting compositional heterogeneities. A bright pixel possibly related to a central peak does not show significant spectral variations, indicating either that the impact site was vertically homogeneous, or that this area has been recovered by homogeneous deposits. Both VIMS ratio images and dielectric constant measurements suggest the presence of an area enriched in water ice around the main ejecta blanket. Since the Ku-band SAR may see subsurface structures at the meter scale, the difference between infrared and SAR observations can be explained by the presence of a thin layer transparent to the radar. An analogy with terrestrial craters in Libya supports this interpretation. Finally, a tentative model describes the geological history of this area prior, during, and after the impact. It involves mainly the creation of ballistic ejecta and an expanding plume of vapor triggered by the impact, followed by the redeposition of icy spherules recondensed from this vapor plume blown downwind. Subsequent evolution is then driven by erosional processes and aeolian deposition. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2010

Automated masking of cloud and cloud shadow for forest change analysis using Landsat images

Accurate masking of cloud and cloud shadow is a prerequisite for reliable mapping of land surface attributes. Cloud contamination is particularly a problem for land cover change analysis, because unflagged clouds may be mapped as false changes, and the level of such false changes can be comparable to or many times more than that of actual changes, even for images with small percentages of cloud cover. Here we develop an algorithm for automatically flagging clouds and their shadows in Landsat images. This algorithm uses clear view forest pixels as a reference to define cloud boundaries for separating cloud from clear view surfaces in a spectral-temperature space. Shadow locations are predicted according to cloud height estimates and sun illumination geometry, and actual shadow pixels are identified by searching the darkest pixels surrounding the predicted shadow locations. This algorithm produced omission errors of around 1% for the cloud class, although the errors were higher for an image that had very low cloud cover and one acquired in a semiarid environment. While higher values were reported for other error measures, most of the errors were found around the edges of detected clouds and shadows, and many were due to difficulties in flagging thin clouds and the shadow cast by them, both by the developed algorithm and by the image analyst in deriving the reference data. We concluded that this algorithm is especially suitable for forest change analysis, because the commission and omission errors of the derived masks are not likely to significantly bias change analysis results.

International Journal of Remote Sensing↗
USGS2022

Formation and distribution of coal balls in the Herrin Coal (Pennsylvanian), Franklin County, Illinois Basin, USA

Large areas of concentrated coal balls (permineralized peat) up to 4 m thick obstructed longwall mining in the Herrin Coal at the Old Ben No. 24 mine. The largest coal‐ball area mapped contained >1500 m 3 ; several areas contained >400 m 3 of coal balls. In‐mine mapping established that there were two types of roof (freshwater and marine), and that the coal balls were spatially correlated with the marine roof units. Regional studies and local data revealed that the younger, freshwater Energy Shale (mud) originally covered all of the peat deposited at the mine, but the mud was locally removed during a period of erosion. The great majority of coal balls are found within the coal seam, where they were permineralized in situ by carbonates. Some coal balls are found exposed on the eroded coal surface and others are within the channel fill associated with the erosion, which predates any marine sedimentation. Thus, the mineralization of the coal balls was synchronous with the erosion of the Energy Shale mud. Data from C‐and O‐isotope, geochemical, and mineralogical analyses of coal balls and associated materials were used to refine a depositional model of coal‐ball formation. The concentrated coal‐ball areas were created by the triggered degassing of CO 2 from partially decomposed peat in the presence of cations from fresh waters; 13 C/ 12 C ratios in these coal balls average −23.9‰. 13 C/ 12 C ratios as low as −34‰ occurring in the most concentrated coal balls are consistent with CO 2 produced by anaerobic oxidation of methane. Top‐of‐seam coal balls were formed later and show slight to strong marine influence; their 13 C/ 12 C ratios average −10.8‰.

Illinois↗
USGS2004

Mapping saline groundwater beneath the Sea Galilee and its vicinity using time domain electromagnetic (TDEM) geophysical technique

An extensive time domain electromagnetic (TDEM) survey covering the Sea of Galilee with a dense grid of points has been recently carried out. A total of 269 offshore and 33 supplementary onshore TDEM soundings were performed along six N-S and ten W-E profiles and at selected points both offshore and onshore along the whole coastal line. The interpreted resistivities were calibrated with the direct salinity measurements in the Haon-2 borehole and relatively deep (5 m) cores taken from the lake bottom. It was found that resistivities below 1 ohm-m are solely indicative of groundwater salinity exceeding 10,000 mg Cl/l. Such low resistivities (high salinities) were detected at depths greater than 15 m below almost the entire bottom of the lake. At some parts of the lake, particularly in the south, the saline water was detected at shallower depths, sometimes at a few meters below the bottom. Relatively high resistivity (fresh groundwater) was found along the margins of the lake down to roughly 100 m, the maximum exploration depth of the system. The detected sharp lateral contrasts at the lake margin between high and low resistivities coincide with the faults separating the carbonate and clastic units, respectively. The geometry of the fresh/saline groundwater interface below the central part of the lake is very similar to the shape of the lake bottom, probably due to the diffusive salt transport from the bottom sediments to the lake water. The above geophysical observations suggest differentsalt transport mechanisms from the sediments to the central part of the lake (diffusion) and from regional aquifers to the margins of the lake (advection). ?? 2004 Science From Israel/LPPLtd.

Israel Journal of Earth Sciences↗