Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Inverse Problems”

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 163 records · Page 9Linked to original sources
USGS2003

Time‐lapse imaging of saline‐tracer transport in fractured rock using difference‐attenuation radar tomography

Accurate characterization of fractured‐rock aquifer heterogeneity remains one of the most challenging and important problems in groundwater hydrology. We demonstrate a promising strategy to identify preferential flow paths in fractured rock using a combination of geophysical monitoring and conventional hydrogeologic tests. Cross‐well difference‐attenuation ground‐penetrating radar was used to monitor saline‐tracer migration in an experiment at the U.S. Geological Survey Fractured Rock Hydrology Research Site in Grafton County, New Hampshire. Radar data sets were collected every 10 min in three adjoining planes for 5 hours during each of 12 tracer tests. An innovative inversion method accounts for data acquisition times and temporal changes in attenuation during data collection. The inverse algorithm minimizes a combination of two functions. The first is the sum of weighted squared data residuals. Second is a measure of solution complexity based on an a priori space‐time covariance function, subject to constraints that limit radar‐attenuation changes to regions of the tomograms traversed by high difference‐attenuation ray paths. The time series of tomograms indicate relative tracer concentrations and tracer arrival times in the image planes; from these we infer the presence and location of a preferential flow path within a previously identified zone of transmissive fractures. These results provide new insights into solute channeling and the nature of aquifer heterogeneity at the site.

New Hampshire↗
USGS1977

SIPT: a seismic refraction inverse modeling program for timeshare terminal computer systems

SIPB is an interactive Fortran computer program that was developed for use with a timeshare computer system with program control information submitted from a remote terminal, and output data displayed on the terminal or printed on a line printer. The program is an upgraded version of FSIPI (Scott, Tibbetts, and Burdick, 1972) with several major improvements in addition to .its adaptation to timeshare operation. The most significant improvement was made in the procedure for handling data from in-line offset shotpoints beyond the end shotpoints of the geophone spread. The changes and improvements are described, user's instructions are outlined, examples of input and output data for a test problem are presented, and the Fortran program is listed in this report. An upgraded batch-mode program, SIPB, is available for users who do not have a timeshare computer system available (Scott, 1977).

Open-File Report↗
USGS2007

Integrating laboratory creep compaction data with numerical fault models: A Bayesian framework

[1] We developed a robust Bayesian inversion scheme to plan and analyze laboratory creep compaction experiments. We chose a simple creep law that features the main parameters of interest when trying to identify rate-controlling mechanisms from experimental data. By integrating the chosen creep law or an approximation thereof, one can use all the data, either simultaneously or in overlapping subsets, thus making more complete use of the experiment data and propagating statistical variations in the data through to the final rate constants. Despite the nonlinearity of the problem, with this technique one can retrieve accurate estimates of both the stress exponent and the activation energy, even when the porosity time series data are noisy. Whereas adding observation points and/or experiments reduces the uncertainty on all parameters, enlarging the range of temperature or effective stress significantly reduces the covariance between stress exponent and activation energy. We apply this methodology to hydrothermal creep compaction data on quartz to obtain a quantitative, semiempirical law for fault zone compaction in the interseismic period. Incorporating this law into a simple direct rupture model, we find marginal distributions of the time to failure that are robust with respect to errors in the initial fault zone porosity.

Journal of Geophysical Research B: Solid Earth↗
USGS1988

Cobalt in ferromanganese crusts as a monitor of hydrothermal discharge on the Pacific sea floor

Ferromanganese oxide crusts, which accumulate on unsedimented surfaces in the open ocean 1–6 , derive most of their metal content from dissolved and particulate matter in ambient bottom water 7,8 , in proportions modified by the variable scavenging efficiency of the oxide phase for susceptible ions 9 . They differ in this respect from abyssal nodules, much of whose metals are remobilized from host sediments. Here we present maps of cobalt concentration and inferred accumulation rate of ferromanganese crusts from the Pacific Ocean. We propose that depletion of cobalt in Pacific crusts measures the location and intensity of submarine hydrothermal discharge. Use of the 'cobalt chronometer', an algorithm inversely relating cobalt content and crust growth rate, permits mapping of the accumulation rate of ferromanganese crusts with only indirect recourse to radioactivity-based dating methods. These maps show that crusts in hydrothermal areas grow from two to more than four orders of magnitude faster than in the Central Pacific Ocean. Cobalt-enriched crusts are found where water masses are most isolated from continental-coastal and hydrothermal sources of metals, now and in the past. This relationship can resolve the problem of cobalt enrichment in crusts without recourse to hypotheses invoking special cobalt sources or enrichment mechanisms.

Nature↗
USGS2024

Physics-based satellite-derived bathymetry (SDB) using Landsat OLI images

The estimation of depth in optically shallow waters using satellite imagery can be efficient and cost-effective. Active sensors measure the distance traveled by an emitted laser pulse propagating through the water with high precision and accuracy if the bottom peak intensity of the waveform is greater than the noise level. However, passive optical imaging of optically shallow water involves measuring the radiance after the sunlight undergoes downward attenuation on the way to the sea floor, and the reflected light is then attenuated while moving back upward to the water surface. The difficulty of satellite-derived bathymetry (SDB) arises from the fact that the measured radiance is a result of a complex association of physical elements, mainly the optical properties of the water, bottom reflectance, and depth. In this research, we attempt to apply physics-based algorithms to solve this complex problem as accurately as possible to overcome the limitation of having only a few known values from a multispectral sensor. Major analysis components are atmospheric correction, the estimation of water optical properties from optically deep water, and the optimization of bottom reflectance as well as the water depth. Specular reflection of the sky radiance from the water surface is modeled in addition to the typical atmospheric correction. The physical modeling of optically dominant components such as dissolved organic matter, phytoplankton, and suspended particulates allows the inversion of water attenuation coefficients from optically deep pixels. The atmospheric correction and water attenuation results are used in the ocean optical reflectance equation to solve for the bottom reflectance and water depth. At each stage of the solution, physics-based models and a physically valid, constrained Levenberg–Marquardt numerical optimization technique are used. The physics-based algorithm is applied to Landsat Operational Land Imager (OLI) imagery over the shallow coastal zone of Guam, Key West, and Puerto Rico. The SDB depths are compared to airborne lidar depths, and the root mean squared error (RMSE) is mostly less than 2 m over water as deep as 30 m. As the initial choice of bottom reflectance is critical, along with the bottom reflectance library, we describe a pure bottom unmixing method based on eigenvector analysis to estimate unknown site-specific bottom reflectance.

Remote Sensing↗
USGS2012

Approaches in highly parameterized inversion - GENIE, a general model-independent TCP/IP run manager

GENIE is a model-independent suite of programs that can be used to generally distribute, manage, and execute multiple model runs via the TCP/IP infrastructure. The suite consists of a file distribution interface, a run manage, a run executer, and a routine that can be compiled as part of a program and used to exchange model runs with the run manager. Because communication is via a standard protocol (TCP/IP), any computer connected to the Internet can serve in any of the capacities offered by this suite. Model independence is consistent with the existing template and instruction file protocols of the widely used PEST parameter estimation program. This report describes (1) the problem addressed; (2) the approach used by GENIE to queue, distribute, and retrieve model runs; and (3) user instructions, classes, and functions developed. It also includes (4) an example to illustrate the linking of GENIE with Parallel PEST using the interface routine.

Techniques and Methods↗
USGS1997

Neural network-based nonlinear model predictive control vs. linear quadratic gaussian control

One problem with the application of neural networks to the multivariable control of mineral and extractive processes is determining whether and how to use them. The objective of this investigation was to compare neural network control to more conventional strategies and to determine if there are any advantages in using neural network control in terms of set-point tracking, rise time, settling time, disturbance rejection and other criteria. The procedure involved developing neural network controllers using both historical plant data and simulation models. Various control patterns were tried, including both inverse and direct neural network plant models. These were compared to state space controllers that are, by nature, linear. For grinding and leaching circuits, a nonlinear neural network-based model predictive control strategy was superior to a state space-based linear quadratic gaussian controller. The investigation pointed out the importance of incorporating state space into neural networks by making them recurrent, i.e., feeding certain output state variables into input nodes in the neural network. It was concluded that neural network controllers can have better disturbance rejection, set-point tracking, rise time, settling time and lower set-point overshoot, and it was also concluded that neural network controllers can be more reliable and easy to implement in complex, multivariable plants.

Mining, Metallurgy & Exploration (MME)↗
USGS2003

Lateral line pore diameters correlate with the development of gas bubble trauma signs in several Columbia River fishes

Gas bubble trauma (GBT) caused by gas supersaturation of river water continues to be a problem in the Columbia River Basin. A common indicator of GBT is the percent of the lateral line occluded with gas bubbles; however, this effect has never been examined in relation to lateral line morphology. The effects of 115, 125 and 130% total dissolved gas levels were evaluated on five fish species common to the upper Columbia River. Trunk lateral line pore diameters differed significantly ( P <0.0001) among species (longnose sucker>largescale sucker>northern pikeminnow≥chinook salmon≥redside shiner). At all supersaturation levels evaluated, percent of lateral line occlusion exhibited an inverse correlation to pore size but was not generally related to total dissolved gas level or time of exposure. This study suggests that the differences in lateral line pore diameters between species should be considered when using lateral line occlusion as an indicator of gas bubble trauma.

Coumbia river basin↗
USGS2001

Effects of pool formation and flash flooding on relative abundance of young-of-year flannelmouth suckers in the Paria River, Arizona

Flannelmouth sucker, Catostomus latipinnis , a fish endemic to the Colorado River basin in the western United States, appears to experience poor recruitment to adult size in the Colorado River, downstream of Glen Canyon Dam. Lack or impermanence of rearing areas for young-of-year (YOY) fish is hypothesized to be the problem. Knowing the importance of tributary mouths as rearing areas in other river systems, we studied use of the mouth of the Paria River, a tributary of the Colorado River, by YOY flannelmouth suckers, and the availability of rearing area in the mouth at different flow levels in the Colorado River in 1996 and 1997. We also examined the relationship between flash floods in the Paria River and catch-per-unit-effort (CPUE) of YOY in the Paria River between 1991 and 1996. Maximum mean daily discharge in the Paria River was inversely correlated with CPUE of YOY flannelmouth suckers (Spearman Rho=−0.9856, p =0.0003) during their critical rearing period (15 March–30 June). Thus, it appears that YOY flannelmouth suckers rear longer in the Paria River in years when flash flooding is minimal. Recruitment of YOY flannelmouth suckers at the Paria River may also be improved by enhancing pool formation during spring and summer rearing seasons. YOY flannelmouth sucker was captured in a pool created by high Colorado River flows (≥336 m 3 /s) that inundated the mouth of the Paria River during spring and summer, 1996. In 1997, high flows (about 550–750 m 3 /s) in the Colorado River during winter and spring initially inundated the Paria River and formed a pool in the mouth. However, these high flows eventually caused 0.5–1.0 m of suspended sediment from the incoming Paria River to deposit in the mouth. Thus, despite higher flows than 1996, the slackwater area formed only occasionally in 1997. Differences in pool formation between 1996 and 1997 demonstrate that pool formation cannot be inferred solely from Colorado River flows.

Arizon↗
USGS1999

The influence of sampling interval on the accuracy of trail impact assessment

Trail impact assessment and monitoring (IA&M) programs have been growing in importance and application in recreation resource management at protected areas. Census-based and sampling-based approaches have been developed in such programs, with systematic point sampling being the most common survey design. This paper examines the influence of sampling interval on the accuracy of estimates for selected trail impact problems. A complete census of four impact types on 70 trails in Great Smoky Mountains National Park was utilized as the base data set for the analyses. The census data were resampled at increasing intervals to create a series of simulated point data sets. Estimates of frequency of occurrence and lineal extent for the four impact types were compared with the census data set. The responses of accuracy loss on lineal extent estimates to increasing sampling intervals varied across different impact types, while the responses on frequency of occurrence estimates were consistent, approximating an inverse asymptotic curve. These findings suggest that systematic point sampling may be an appropriate method for estimating the lineal extent but not the frequency of trail impacts. Sample intervals of less than 100 m appear to yield an excellent level of accuracy for the four impact types evaluated. Multiple regression analysis results suggest that appropriate sampling intervals are more likely to be determined by the type of impact in question rather than the length of trail. The census-based trail survey and the resampling-simulation method developed in this study can be a valuable first step in establishing long-term trail IA&M programs, in which an optimal sampling interval range with acceptable accuracy is determined before investing efforts in data collection.

Landscape and Urban Planning↗
USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗
USGS1983

Water resources of the Tulalip Indian Reservation, Washington

Water will play a significant role in the future development of the Tulalip Indian Reservation. Ground-water resources are sufficient to supply several times the 1978 population. Potential problems associated with increased ground-water development are saltwater encroachment in the coastal areas and septic-tank contamination of shallow aquifers. There are sufficient good-quality surface-water resources to allow for significant expansion of the tribe)s fisheries activities. The tribal well field is the only place where the ground-water system has been stressed) resulting in declining water levels (1,5 feet per year), The well field has a useful life of at least 1.5-20 years, This can be increased by drilling additional wells to expand the present well field, Inflow of water to the reservation is in the form of precipitation (103 cubic feet per second) ft3/s)) surface-water inflow (13 ft3/s)) and ground-water inflow (4 ft3/s), Outflow is as evapotranspiration (62 ft3/s)) surface-water outflow (40 ft3/s)) and ground-water outflow (18 ft3/s), Total inflow and outflow are equal (120 ft3/s). Ground water is generally suitable for domestic use without treatment) but a serious quality problem is the presence of coliform bacteria in some shallow wells, High values of turbidity and color and large concentrations of iron and manganese are common problems regarding the esthetic quality of the water, In a few places, large concentrations of chloride and dissolved solids indicate the possibility of saltwater encroachment, but no ongoing trend has been identified, Surface waters have been observed to contain undesirably high concentrations of total phosphorus and total and fecal-coliform bacteria) and to have temperatures too high for fish-rearing. The concentration of nutrients appears to be related to flow conditions. Nitrate and total nitrogen are greater in wet-season runoff than during low-flow periods) and total phosphorus shows an inverse relationship. Total phosphorus and ammonia concentrations are greatest in dry-season storm runoff. Generally) surface-water quality is adequate for fish-rearing and (with treatment) for public supply,

Open-File Report↗
USGS2016

Geology, selected geophysics, and hydrogeology of the White River and parts of the Great Salt Lake Desert regional groundwater flow systems, Utah and Nevada

The east-central Great Basin near the Utah-Nevada border contains two great groundwater flow systems. The first, the White River regional groundwater flow system, consists of a string of hydraulically connected hydrographic basins in Nevada spanning about 270 miles from north to south. The northernmost basin is Long Valley and the southernmost basin is the Black Mountain area, a valley bordering the Colorado River. The general regional groundwater flow direction is north to south. The second flow system, the Great Salt Lake Desert regional groundwater flow system, consists of hydrographic basins that straddle the Utah-Nevada border, with a length of about 150 miles from north to south. The general regional groundwater flow direction is from south to north towards the Great Salt Lake Desert. For 15 years with support from the Southern Nevada Water Authority (SNWA), hydrologists, geologists, and geophysicists studied the basin connections and the groundwater resources in these and adjacent flow systems over an area of about 25,000 square miles. A major first part of the SNWA study was constructing a 3-dimensional digital hydrogeologic framework based on geologic maps and cross sections at 1:250,000 scale. This framework documents the presence of three major aquifers: (1) Paleozoic carbonate rocks, (2) Eocene to Miocene volcanic rocks, and (3) Miocene to Holocene basin-fill sediments, as well as confining units that constrain flow. We interpret that movement of most groundwater through and across basins is by fracture-dominated flow along faults/fractures, yet in most places flow is prevented or retarded across faults, so mapping structures gives a first approximation to conduits and barriers to flow. The most important structures by far are high-angle normal faults of the basin-range episode of east-west extensional deformation. This event began at about 20 Ma, although most deformation and the formation of the present topography took place between 10 Ma and present. This topography consists of north-trending basins (mostly grabens) that alternate with north- trending ranges (mostly horsts); erosion of the ranges filled the basins with clastic alluvial basin-fill deposits. Geophysics provides data on the third dimension (cross sections) of the hydrogeologic framework. Audiomagnetotelluric profiles and gravity inversion located faults and enabled us to estimate thicknesses of basin-fill deposits. To this framework, hydrologic studies addressed precipitation, surface water, and springs, as well as groundwater levels, volumes, geochemistry, water budgets, and monitoring. At nearly the same time as our study, the Utah Geological Survey (UGS) and U.S. Geological Survey (USGS) addressed the same issues in many of the same areas, and publication of the efforts by all three agencies reveals a surprising similarity of conclusions, with some critical exceptions, which therefore demonstrates the great value of many scientists independently studying the same complex scientific problem. The differences in conclusions include directions and volumes of some ground- water flow paths, such as one proposed by the USGS of unlikely groundwater flow from Steptoe Valley to southern Snake Valley, and another proposed by the UGS of unlikely significant groundwater recharge flow from the Snake Range to the Fish Springs complex.

Utah↗
USGS2003

Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California

We use land-subsidence observations from repeatedly surveyed benchmarks and interferometric synthetic aperture radar (InSAR) in Antelope Valley, California, to estimate spatially varying compaction time constants, ??, and inelastic specific skeletal storage coefficients, Skv*, in a previously calibrated regional groundwater flow and subsidence model. The observed subsidence patterns reflect both the spatial distribution of head declines and the spatially variable inelastic skeletal storage coefficient. Using the nonlinear parameter estimation program UCODE we estimate compaction time constants between 3.8 and 285 years. The Skv* values are estimated by linear estimation and range from 0 to almost 0.09. We find that subsidence observations over long time periods are necessary to constrain estimates of the large compaction time constants in Antelope Valley. The InSAR data used in this study cover only a three-year period, limiting their usefulness in constraining these time constants. This problem will be alleviated as more SAR data become available in the future or where time constants are small. By incorporating the resulting parameter estimates in the previously calibrated regional model of groundwater flow and land subsidence we can significantly improve the agreement between simulated and observed land subsidence both in terms of magnitude and spatial extent. The sum of weighted squared subsidence residuals, a common measure of model fit, was reduced by 73% with respect to the original model. However, the ability of the model to adequately reproduce the subsidence observed over only a few years is impaired by the fact that the simulated hydraulic heads over small time periods are often not representative of the actual aquifer hydraulic heads. Errors in the simulated hydraulic aquifer heads constitute the primary limitation of the approach presented here.

Water Resources Research↗
USGS2012

Use of focal mechanisms to determine stress: A control study.

To allow focal mechanisms to be inverted for the stress field requires a different inversion algorithm than for slickenside data because focal mechanisms do not represent fault slip data unless one can decide which nodal plane is the fault plane. If one can decide which nodal plane is the fault plane, then the focal mechanisms can be inverted with the slickenside inversion algorithm. This decision cannot always be made, so algorithms for inverting focal mechanisms for the stress field are studied. These algorithms either use both of the possible fault planes or attempt to choose the correct fault plane while determining the stress tensor. Simulated focal mechanisms are made from slickenside data and used to provide a control study for the focal mechanism inversion algorithms. The results of this control study show that focal mechanisms can be inverted to find the best stress tensor, but the resolution is decreased unless the fault planes can be picked a priori. The resolution can also be increased by including constraints on the magnitude of the tangential traction on the fault plane. Therefore, using focal mechanisms to study small variations in the stress field requires that other data (e.g., studies of the hypocenters, surface faulting, or structural information concerning the region) be introduced to pick which of the nodal planes is the fault plane. This study also introduces the method of bootstrap resampling to the statistics of this problem. The non-Gaussian nature of the data makes the nonparametric formulation of the bootstrap approach ideal for this problem.

Journal of Geophysical Research Solid Earth↗
USGS1996

Calculation of earthquake rupture histories using a hybrid global search algorithm: Application to the 1992 Landers, California, earthquake

A method is presented for the simultaneous calculation of slip amplitudes and rupture times for a finite fault using a hybrid global search algorithm. The method we use combines simulated annealing with the downhill simplex method to produce a more efficient search algorithm then either of the two constituent parts. This formulation has advantages over traditional iterative or linearized approaches to the problem because it is able to escape local minima in its search through model space for the global optimum. We apply this global search method to the calculation of the rupture history for the Landers, California, earthquake. The rupture is modeled using three separate finite-fault planes to represent the three main fault segments that failed during this earthquake. Both the slip amplitude and the time of slip are calculated for a grid work of subfaults. The data used consist of digital, teleseismic P and SH body waves. Long-period, broadband, and short-period records are utilized to obtain a wideband characterization of the source. The results of the global search inversion are compared with a more traditional linear-least-squares inversion for only slip amplitudes. We use a multi-time-window linear analysis to relax the constraints on rupture time and rise time in the least-squares inversion. Both inversions produce similar slip distributions, although the linear-least-squares solution has a 10% larger moment (7.3 ?? 1026 dyne-cm compared with 6.6 ?? 1026 dyne-cm). Both inversions fit the data equally well and point out the importance of (1) using a parameterization with sufficient spatial and temporal flexibility to encompass likely complexities in the rupture process, (2) including suitable physically based constraints on the inversion to reduce instabilities in the solution, and (3) focusing on those robust rupture characteristics that rise above the details of the parameterization and data set.

Physics of the Earth and Planetary Interiors↗
USGS2014

Simulation of water-table aquifers using specified saturated thickness

Simulating groundwater flow in a water-table (unconfined) aquifer can be difficult because the saturated thickness available for flow depends on model-calculated hydraulic heads. It is often possible to realize substantial time savings and still obtain accurate head and flow solutions by specifying an approximate saturated thickness a priori, thus linearizing this aspect of the model. This specified-thickness approximation often relies on the use of the “confined” option in numerical models, which has led to confusion and criticism of the method. This article reviews the theoretical basis for the specified-thickness approximation, derives an error analysis for relatively ideal problems, and illustrates the utility of the approximation with a complex test problem. In the transient version of our complex test problem, the specified-thickness approximation produced maximum errors in computed drawdown of about 4% of initial aquifer saturated thickness even when maximum drawdowns were nearly 20% of initial saturated thickness. In the final steady-state version, the approximation produced maximum errors in computed drawdown of about 20% of initial aquifer saturated thickness (mean errors of about 5%) when maximum drawdowns were about 35% of initial saturated thickness. In early phases of model development, such as during initial model calibration efforts, the specified-thickness approximation can be a very effective tool to facilitate convergence. The reduced execution time and increased stability obtained through the approximation can be especially useful when many model runs are required, such as during inverse model calibration, sensitivity and uncertainty analyses, multimodel analysis, and development of optimal resource management scenarios.

Ground Water↗
USGS2026

Unraveling protracted modification of Archean and Paleoproterozoic crust in central Laurentia, Penokean orogen, with garnet and accessory mineral geochronology and microstructural analysis

Proterozoic metamorphism and deformation of the southern margin of the Superior craton in the Lake Superior region is attributed to the Penokean orogeny (1890−1830 Ma). This model includes a period of crustal inversion in which Archean basement blocks were exhumed through overlying Paleoproterozoic strata, producing the corridor of gneiss domes that parallels the trend of the Penokean orogen across the northern Midcontinent, USA. However, recent geologic mapping and 40 Ar/ 39 Ar geochronology challenge this interpretation, suggesting instead that the gneiss dome structures reflect younger episodes of tectonic activity along the southern margin of Laurentia. In absence of integrated pressure-temperature-time-deformation constraints for these rocks, interpretations are largely limited to their final cooling history, making it difficult to both identify the tectonic forces that shaped the architecture of the Penokean orogenic belt and assess the extent to which later Proterozoic tectonism modified the southern Superior craton. We address this problem with an approach joining thermodynamic modeling, garnet and accessory mineral geochronology, and microstructural analysis for several metamorphic rocks across the gneiss dome corridor. The U-Pb ages of titanite reveal that the Proterozoic geometries of exhumed basement gneiss domes are governed by preexisting Archean structures. Garnet Lu-Hf geochronology constrains the timing of prograde-to-peak metamorphism in the Penokean orogenic belt. Granulite facies metamorphism is related to the final stages of the Penokean orogeny at 1837 Ma and localized in a belt of high-grade rocks near a major Penokean suture. Garnet Lu-Hf ages of samples adjacent to gneiss domes reflect regional metamorphism following the accretionary phase of the Penokean orogeny, between 1825 Ma and 1782 Ma, which we suggest reflects continued crustal thickening related to convergence farther south during this time interval. Combination of garnet microstructures and Sm-Nd ages reflects later exhumation of gneiss domes and buried metasedimentary rocks by ca. 1750 Ma, consistent with previously published 40 Ar/ 39 Ar cooling ages across the region. Reset Lu-Hf and Sm-Nd garnet ages and U-Pb ages of syn-kinematic titanite reflect reactivation of primary Penokean structures during this period of basement uplift. These data document significant modification of the Penokean orogen and the Archean crust of the southern Superior province between 1800 Ma and 1700 Ma. Tectonic activity during this interval coincides with collisional events recognized in western Laurentia, suggesting that the period immediately following the Penokean orogeny may be a broadly important time for crustal growth and modification in proto-North America.

Michigan, Wisconsin↗