Evolution of porosity and geochemistry in Marcellus Formation black shale during weathering
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi-year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage-limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.
This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.
Nuclear reactor operators can expand safety margins with more precise information about how materials behave inside operating reactors. INL's new simulation platform makes such studies easier & more informative by letting researchers "plug-n-play" their mathematical models, skipping years of computer code development.
Explore the source record for details and available documents.
Abstract As bedrock weathers to regolith – defined here as weathered rock, saprolite, and soil – porosity grows, guides fluid flow, and liberates nutrients from minerals. Though vital to terrestrial life, the processes that transform bedrock into soil are poorly understood, especially in deep regolith, where direct observations are difficult. A 65-m-deep borehole in the Calhoun Critical Zone Observatory, South Carolina, provides unusual access to a complete weathering profile in an Appalachian granitoid. Co-located geophysical and geochemical datasets in the borehole show a remarkably consistent picture of linked chemical and physical weathering processes, acting over a 38-m-thick regolith divided into three layers: soil; porous, highly weathered saprolite; and weathered, fractured bedrock. The data document that major minerals (plagioclase and biotite) commence to weather at 38 m depth, 20 m below the base of saprolite, in deep, weathered rock where physical, chemical and optical properties abruptly change. The transition from saprolite to weathered bedrock is more gradational, over a depth range of 11–18 m. Chemical weathering increases steadily upward in the weathered bedrock, with intervals of more intense weathering along fractures, documenting the combined influence of time, reactive fluid transport, and the opening of fractures as rock is exhumed and transformed near Earth’s surface.
The US Army Corps of Engineers (USACE) is addressing the Formerly Utilized Sites Remedial Action Program (FUSRAP) Middlesex Sampling Plant Site (MSP) in accordance with the Comprehensive Environmental Response, Compensation, and Liability Act (CERCLA). The MSP has been an active FUSRAP site since the early 1980's and was listed on the National Priority List (NPL) in 1999. Two challenges that the project team are currently addressing are Ground Water (GW) movement through fractured bedrock and Vicinity Properties (VP) assessments. USACE and the US Environmental Protection Agency (USEPA) signed a Record of Decision (ROD) for soils remedial action in 2005 and USACE completed the soils remediation at MSP in 2008, however, VPs were not a part of the soils ROD. Additionally, the State of New Jersey recently identified a potential contaminated VP of the site. VPs are any property in the vicinity of a CERCLA 'facility.' The contaminated VPs of the MSP were addressed prior to the adoption of the CERCLA process to address FUSRAP sites. Since the 1980's many changes to potential Applicable or Relevant and Appropriate Requirements have occurred. As such USACE has begun reassessing the previously addressed VPs for inclusion under the soils ROD as part of the CERCLA process. USACE does expect to have to revisit and potentially remediate up to 12 VPs. This paper will discuss the approach, some challenges, and lessons learned in adequately addressing VPs both from a field work and a post ROD regulatory perspective. USACE is currently in the Feasibility Study phase of the CERCLA process for GW. The project team has encountered and overcome significant challenges with regard to contaminant plume bounding and groundwater movement through the secondary porosity of dipping bedding plane partings and fractured bedrock of the Passaic Aquifer beneath the site. This paper will discuss the challenges, the approaches used, and lessons learned while delineating these contaminants. This paper's discussions are useful for anyone addressing an environmental site with soils contamination on VPs or with delineating contaminants in fractured bedrock. (authors)
Batch sorption and column experiments were conducted to investigate and compare sorption and transport behavior of 90Sr on the assumption of seawater intrusion at nuclear power plants. Batch sorption experiments were carried out on fractured rocks and bedrocks using synthetic groundwater and seawater both spiked with 90Sr. In general, higher 90Sr sorption occurred on fractured rock samples than on bedrocks, because of the presence of weathered secondary minerals (iron oxide and clay) on fractured rock surfaces. However, one particular bedrock sample (WSP-B) which has higher porosity and carbon amount than fractured rock samples also showed the higher 90Sr sorption than its comparable fractured rocks. For all batch sorption studies, 90Sr sorption distribution coefficient, Kd decreased from groundwater to seawater environment due to the higher ionic strength of seawater (6.4×10-1–7.7×10-1 M) compared to groundwater (4.0×10-3–6.0×10-3 M). The three different ionic strength solutions were used in column experiments, and the results showed that transport behavior of Sr through a fractured rock had similar sorption trend to batch sorption results. The highest mobility (or least retardation) for Sr was found for 100% seawater solution compared to the highest retardation (or least mobility) for 100% groundwater solution. These sorption and transport data of Sr on solid materials contacted with various ionic strength solutions corroborate empirically defensible information for assessment of radioactive contamination in groundwater below the NPP sites located nearby shores. In addition, the experimental data will be incorporated to improve transport models of 90Sr in the subsurface environment for severe nuclear accidents.
In October 2018, we conducted an exploratory seismic survey of the East River Watershed Function Science Focus Area, near Crested Butte, CO. Our survey goals were to map seismic bedrock, regolith and sediment properties and relate them to geomorphic and hydrologic processes. We acquired ~12 km of seismic data along roads with a 72-channel, 1.25m spaced streamer, an average shot spacing of 2.5m and an average acquisition rate of 3 km/day. We also obtained ~3km of planted geophone data with geophone and shot spacing between 2.5 and 5 meters and an acquisition rate of 0.5 km/day. Vp results show a sharp transition between regolith and bedrock along hillslopes suggesting a similarly abrupt change in porosity. Vp anisotropy in the Mancos Shale is observable in the lower reaches of the watershed, but is absent in regions near crystalline intrusions where contact metamorphism has likely modified pore structure within the shale. Alluvium thickness varies along the river axis from 0 to > 20 meters reflecting variations in depositional processes as function of stream gradient. We show that when strong discontinuities at the regolith base are present, inhomogeneous waves that travel along subsurface interfaces are generated and can be measured on vertical component geophones. We present a summary manuscript in the report appendix.
We explore the contribution of fractures (joints) in controlling the rate of weathering advance for a low-porosity rock by using methods of homogenization to create averaged weathering equations. The rate of advance of the weathering front can be expressed as the same rate observed in non-fractured media (or in an individual block) divided by the volume fraction of non-fractured blocks in the fractured parent material. In the model, the parent has fractures that are filled with a more porous material that contains only inert or completely weathered material. The low-porosity rock weathers by reaction-transport processes. As observed in field systems, the model shows that the weathering advance rate is greater for the fractured as compared to the analogous non-fractured system because the volume fraction of blocks is <1. The increase in advance rate is attributed both to the increase in weathered material that accompanies higher fracture density, and to the increase in exposure of surface of low-porosity rock to reaction-transport. For constant fracture aperture, the weathering advance rate increases when the fracture spacing decreases. Equations describing weathering advance rate are summarized in the “List of selected equations”. If erosion is imposed at a constant rate, the weathering systems with fracture-bounded bedrock blocks attain a steady state. In the erosional transport-limited regime, bedrock blocks no longer emerge at the air-regolith boundary because they weather away. In the weathering-limited (or kinetic) regime, blocks of various size become exhumed at the surface and the average size of these exposed blocks increases with the erosion rate. For convex hillslopes, the block size exposed at the surface increases downslope. This model can explain observations of exhumed rocks weathering in the Luquillo mountains of Puerto Rico.
The prediction of the geothermal system efficiency is strong linked to the character of the flow system that connects injector and producer wells. If water flow develops channels or “short circuiting” between injection and extraction wells thermal sweep is poor and much of the reservoir is left untapped. The purpose of this project was to understand how channelized flow develops in fracture geothermal reservoirs and how it can be measured in the field. We explored two methods of assessing channelization: hydraulic connectivity tests and tracer tests. These methods were tested at a field site using two verification methods: ground penetrating radar (GPR) images of saline tracer and heat transfer measurements using distributed temperature sensing (DTS). The field site for these studies was the Altona Flat Fractured Rock Research Site located in northeastern New York State. Altona Flat Rock is an experimental site considered a geologic analog for some geothermal reservoirs given its low matrix porosity. Because soil overburden is thin, it provided unique access to saturated bedrock fractures and the ability image using GPR which does not effectively penetrate most soils. Five boreholes were drilled in a “five spot” pattern covering 100 m2 and hydraulically isolated in a single bedding plane fracture. This simple system allowed a complete characterization of the fracture. Nine small diameter boreholes were drilled from the surface to just above the fracture to allow the measurement of heat transfer between the fracture and the rock matrix. The focus of the hydraulic investigation was periodic hydraulic testing. In such tests, rather than pumping or injection in a well at a constant rate, flow is varied to produce an oscillating pressure signal. This pressure signal is sensed in other wells and the attenuation and phase lag between the source and receptor is an indication of hydraulic connection. We found that these tests were much more effective than constant pumping tests in identifying a poorly connected well. As a result, we were able to predict which well pairs would demonstrate channelized flow. The focus of the tracer investigation was multi-ionic tests. In multi-ionic tests several ionic tracers are injected simultaneously and the detected in a nearby pumping well. The time history of concentration, or breakthrough curve, will show a separation of the tracers. Anionic tracers travel with the water but cationic tracer undergo chemical exchange with cations on the surface of the rock. The degree of separation is indicative of the surface area exposed to the tracer. Consequently, flow channelization will tend to decrease the separation in the breakthrough. Estimation of specific surface area (the ration of fracture surface area to formation volume) is performed through matching the breakthrough curve with a transport model. We found that the tracer estimates of surface area were confirmed the prediction of channelized flow between well pairs produced by the periodic hydraulic tests. To confirm that the hydraulic and tracer tests were correctly predicting channelize flow, we imaged the flow field using surface GPR. Saline water was injected between the well pairs which produced a change in the amplitude and phase of the reflected radar signal. A map was produced of the migration of saline tracer from these tests which qualitatively confirmed the flow channelization predicted by the hydraulic and tracer tests. The resolution of the GPR was insufficient to quantitatively estimate swept surface area, however. Surface GPR is not applicable in typical geothermal fields because the penetration depths do not exceed 10’s of meters. Nevertheless, the method of using of phase to measure electrical conductivity and the assessment of antennae polarization represent a significant advancement in the field of surface GPR. The effect of flow character on fracture / rock thermal exchange was evaluated using heated water as a tracer. Water elevated 30 degrees C above the formation water was circulated between two wells pairs. One well pair had been identified in hydraulic and tracer testing as well connected and the other poorly connected. Temperature rise was measured in the adjacent rock matrix using coiled fiber optic cable interrogated for temperature using a DTS. This experimental design produced over 4000 temperature measurements every hour. We found that heat transfer between the fracture and the rock matrix was highly impacted by the character of the flow field. The strongly connected wells which had demonstrated flow channelization produced heat rise in a much more limited area than the more poorly connected wells. In addition, the heat increase followed the natural permeability of the fracture rather than the induced flow field. The primary findings of this work are (1) even in a single relatively planar fracture, the flow field can be highly heterogeneous and exhibit flow channeling, (2) channeling results from a combination of fracture permeability structure and the induced flow field, and (3) flow channeling leads to reduced heat transfer. Multi-ionic tracers effectively estimate relative surface area but an estimate of ion-exchange coefficients are necessary to provide an absolute measure of specific surface area. Periodic hydraulic tests also proved a relative indicator of connectivity but cannot prove an absolute measure of specific surface area.
Significance To understand groundwater in weathered bedrock aquifers, we explore (bio)geochemical reactions using geophysical measurements. The result leads to a conceptual model for the three-dimensional distribution of geochemical alteration in the subsurface—well beyond the typical approach of inferring (bio)geochemical reactions in the subsurface from pits, boreholes, and outcrops. With this geophysical technique we discovered the prevalence of gas beneath the water table at this site. This geophysical viewing of the subsurface highlights patterns of water movement, reactions of clays, and microbial and abiotic gas-producing reactions recorded in the subsurface and thus illuminates storage, flow, and reactions of groundwater in bedrock systems. Protecting the groundwater resource beneath our feet will rely on such enhanced capacities to see into the subsurface.
Models developed to assess the effects of megaregolith on lunar thermal evolution are discussed. It is confirmed that the two sites where lunar heat flow was measured are probably unrepresentative, with heat flows about 25% higher than regional averages, due to focussing of heat flow towards regions with thin megaregolith. Numerous lines of evidence indicate that the megaregolith is generally 2 to 3 km thick under highlands (which cover about 83% of the total lunar surface), and 1 km thick under maria. In most models, megaregolith thickness is assumed to be roughly 6x greater over highlands than over maria. Based on sparse data for porosity among lunar rock types, and the correlation between thermal conductivity and porosity, it is assumed that megaregolith conductivity is roughly 20 kiloerg s(-1)cm(-1)K(-1), and bedrock conductivity is roughly 7x greater. It is also found that insulation by megaregolith exacerbates the problem of reconciling modest temperatures inferred for the (present) matle with a high rate of heat production; an upper limit of 30 ng/g for the bulk-Moon U content can be derived from this constraint alone.
Models developed to assess the effects of megaregolith on lunar thermal evolution are discussed. It is confirmed that the two sites where lunar heat flow was measured are probably unrepresentative, with heat flows about 25 percent higher than regional averages, due to focussing of heat flow towards regions with thin megaregolith. Numerous lines of evidence indicate that the megaregolith is generally 2 to 3 km thick under highlands (which cover about 83 percent of the total lunar surface), and 1 km thick under maria. In most models, megaregolith thickness is assumed to be roughly 6x greater over highlands than over maria. Based on sparse data for porosity among lunar rock types, and the correlation between thermal conductivity and porosity, it is assumed that megaregolith conductivity is roughly 20 kiloerg s(-1)cm(-1)K(-1), and bedrock conductivity is roughly 7x greater. It is also found that insulation by megaregolith exacerbates the problem of reconciling modest temperatures inferred for the (present) mantle with a high rate of heat production; an upper limit of 30 ng/g for the bulk-Moon U content can be derived from this constraint alone.
This dataset includes geoelectrical monitoring data acquired between October 2021 and November 2022, soil moisture and temperature data, groundwater data obtained from borehole SNIB covering the period from June 2021 to September 2022, and hydrological modelling results. The data were acquired to investigate how variations in bedrock type and topography, and vegetation cover control subsurface flow dynamics. To provide insights into the subsurface flow dynamics and their controls, a monitoring transect was installed at the Chicken Bone Meadow, Mt. Snodgrass, Crested Butte, CO, measuring the spatio-temporal variations of soil moisture, soil and snow temperature, subsurface electrical resistivity variations, and groundwater dynamics. Field data are organized in a folder structure, with Electrical Resistivity Tomography (ERT) data being provided as one file per measurement, and data of the soil moisture and temperature sensors being provided as text files covering the entire monitoring period. The ‘Locations.csv’ file contains the location of all sensors, given in NAD83 – UTM Zone 13N. ERT monitoring data has been processed to filter data based on reciprocal errors (data with errors > 30% were removed), a linear error model was fitted to each survey, and to ensure a constant set of measurements for time-lapse inversion, filtered data were interpolated and assigned a 100% measurement error. Soil moisture and temperature data were acquired at 15 min intervals, and averaged to provide 1h data. Weather data and borehole data (groundwater depth, conductivity and temperature) were acquired at 30 min intervals, and are provided as daily measurements; all measurements are averaged, except of precipitation values, which are given as daily accumulation. The hydrological model was set up along the ERT monitoring transect, and net infiltration was used as surface boundary condition and derived from the weather data. Four different results are provided, (1) results for a parameterization using hydraulic permeability and porosity as derived from the ERT data through petrophysical relationships, and (2) three simplified model results, using 1 to 3 geological layers above the bedrock. Modelling was performed using PFLOTRAN, and for each model the PFLOTRAN input files are provided. The result files include weekly hydrological modelling results (e.g., saturation, velocities, pressures), as well as the model parameterization. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type.
Abstract Bedrock weathering regulates nutrient mobilization, water storage, and soil production. Relative to the mobile soil layer, little is known about the relationship between topography and bedrock weathering. Here, we identify a common pattern of weathering and water storage across a sequence of three ridges and valleys in the sedimentary Great Valley Sequence in Northern California that share a tectonic and climate history. Deep drilling, downhole logging, and characterization of chemistry and porosity reveal two weathering fronts. The shallower front is ∼7 m deep at the ridge of all three hillslopes, and marks the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front marks the extent of open fractures and discoloration. This front is 11 m deep under two ridges of similar ridge‐valley spacing, but 17.5 m deep under a ridge with nearly twice the ridge‐valley spacing. Hence, at ridge tops, the fraction of the hillslope relief that is weathered scales with hillslope length. In all three hillslopes, below this second weathering front, closed fractures and unweathered bedrock extend about one‐half the hilltop elevation above the adjacent channels. Neutron probe surveys reveal that seasonally dynamic moisture is stored to approximately the same depth as the shallow weathering front. Under the channels that bound our study hillslopes, the two weathering fronts coincide and occur within centimeters of the ground surface. Our findings provide evidence for feedbacks between erosion and weathering in mountainous landscapes that result in systematic subsurface structuring and water routing.
This report summarizes the research completed during the project entitled "The weathered bedrock vadose zone: A hidden control on water availability in the western United States" (Award: DE-SC0018039). The report covers work completed between 9/15/17 and 9/14/20. The objective of the exploratory project work was to evaluate the overarching hypothesis that the subsurface underlying actively eroding landscapes is systematically structured as a result of a balance between weathering and erosion and that this structure leads to dynamic moisture storage within weathered rock. To accomplish project objectives, we characterized hillslope weathering profile structure and conducted subsurface moisture monitoring across four study sites with comparable bedrock types but differing erosion rate and hydroclimate: The Angelo Coast Range Reserve (ACRR), Sagehorn-Russell Ranch (SRR), Antelope Valley Ranch (AVR), and the LBNL Watershed Function Scientific Focus Area (SFA). New characterization and monitoring infrastructure at the AVR and SFA sites was established through this project and enables long-term monitoring of hydrologic and carbon fluxes at the sites for comparison with long-term monitoring in other sites. Deep drilling, downhole logging, and characterization of chemistry and porosity of rock recovered from drilling reveals two weathering fronts across the sites. A shallower front is marked by the extent of pervasive fracturing and oxidation of pyrite and organic carbon. A deeper weathering front is marked by the appearance of open fractures and discoloration. At the AVR site, multiple ridge tops of different hillslope length were investigated to establish the relationship between the weathering front and topography. Drilling reveals that the fraction of the hillslope relief that is weathered scales with hillslope length, consistent with the theoretical framework that guided the project relating groundwater drainage with landscape evolution. With the establishment of new boreholes, we have documented significant water storage in weathered bedrock using low-field borehole nuclear magnetic resonance (NMR) across all four sites. In the upper 4 m of bedrock, the average water contents at ACRR, AVR, SRR, and SFA are approximately 25%, 16%, 25%, 23% respectively. At the ACRR site, where long-term monitoring via neutron probe surveys and time-domain transmission sensors are available, we have demonstrated that NMR reliably records moisture content and importantly, moisture content dynamics in weathered bedrock. Across all sites, our NMR monitoring has also revealed that a considerable fraction of the water storage occurs in fractures or pores that are significantly larger than the fine-grained bedrock matrix of the parent rock. To identify the proportion of water storage occurring in fractures and large pores, we analyzed the NMR signal by evaluating both the sum of echoes and an inversion of the NMR signal for a distribution of T2 relaxation times. Both methods reveal that at least 20% of the water storage occurs in fractures and larger pores across all sites. To constrain the interpretation of the in-situ field-scale NMR measurements, we compare the NMR response in the field to independent pore-structure information obtained via helium pycnometry. Results support the inference that dynamic seasonal water storage is dominantly restricted to the fracture network, which can reach 7% of the total volume. This dynamic range in water content is comparable to that of some soils, underscoring the need to understand how this region of the Earth System is structured, and functions as a hydrologic reservoir.
Radar backscatter intensity as measured by calibrated synthetic aperture radar (SAR) systems is primarily controlled by three factors: local incidence angle, wavelength-scale roughness, and dielectric permittivity of surface materials. Radar observations may be of limited use for geological investigations of surface composition, unless the relationships between lithology and the above characteristics can be adequately understood. In arid terrains, such as the Southwest U.S., weathering signatures (e.g. soil development, fracturing, debris grain size and shape, and hill slope characteristics) are controlled to some extent by lithologic characteristics of the parent bedrock. These textural features of outcrops and their associated debris will affect radar backscatter to varying degrees, and the multiple-wavelength capability of the JPL Airborne SAR (AIRSAR) system allows sampling of textures at three distinct scales. Diurnal temperature excursions of geologic surfaces are controlled primarily by the thermal inertia of surface materials, which is a measure of the resistance of a material to a change in temperature. Other influences include albedo, surface slopes affecting insolation, local meteorological conditions and surface emissivity at the relevant thermal wavelengths. To first order, thermal inertia variations on arid terrain surfaces result from grain size distribution and porosity differences, at scales ranging from micrometers to tens of meters. Diurnal thermal emission observations, such as those made by the JPL Thermal Infrared Multispectral Scanner (TIMS) airborne instrument, are thus influenced by geometric surface characteristics at scales comparable to those controlling radar backscatter. A preliminary report on a project involving a combination of field, laboratory and remote sensing observations of weathered felsic-to basaltic volcanic rock units exposed in the southern part of the Lunar Crater Volcanic Field, in the Pancake Range of central Nevada is presented. Focus is on the relationship of radar backscatter cross sections at multiple wavelengths, apparent diurnal temperature excursions identified in multi-temporal TIMS images, surface geometries related to weathering style, and parent bedrock lithology.