Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geologic model”

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.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Some Opinions on Remote Sensing and Geologic Studies

The principal role of remote sensing data in geologic studies is as a source of geologic information from which meaningful geologic interpretations can be made. Remote sensing data are important in sedimentary basin analysis and other geologic studies as independent and sometimes unique sources of important lithologic and structural information; however, their greatest benefit to exploration-oriented investigations may come when these data are used with other relevant data in a digital database approach to exploration. Modern computer technology facilitates the rapid integration and synthesis of satellite, topographic, gravity, aeromagnetic, geochemical, and other data collected from a given region. Once such data are geometrically registered, they can be digitally processed, within the constraints of defined geologic models, to rapidly identify, and focus further exploration efforts on, target areas that have the greatest potential for success.

Bailey, G. B.↗
NASA NTRS1986

Stress on the surface of Venus: Analysis of data provided by Venera-15 and Venera-16

There are no large-scale sedimentation and erosion processes on Venus, therefore the observed deformation must be formed on the surface or close to it, and in the ellipsoids of stresses one of the main axes must be perpendicular to the surface. According to the deformation visible on photography from Venera-15, 16, the orientations of the ellipsoids of stresses are established in the different sections of the Montes Maxwell and the northern parquet. A geological model is established, according to which the Montes Maxwell is an arched vault and the region to the south a conjugate bend, separated by a shift zone.

Ronka, L. B.↗
NASA NTRS1992

Plate tectonics from VLBI and SLR global data

This study is based on data derived from fifteen years of observations of the SLR (side-looking radar) network and six years of the VLBI (very long baseline interferometry) network. In order to use all available information VLBI and SLR global data sets were combined in a least squares fashion to calculate station horizontal velocities. All significant data pertaining to a single site contribute to the station horizontal motion. The only constraint on the solution is that no vertical motion is allowed. This restriction does not greatly affect the precision of the overall solution given the fact that the expected vertical motion for most stations, even those experiencing post glacial uplift, is well under 1 cm/yr. Since the average baseline is under 4,000 km, only a small fraction of the station vertical velocity is translated into baseline rates so that the error introduced in the solution by restricting up-down station movement is minimal. As a reference, station velocities were then compared to the ones predicted by the NUVEL-1 geological model of DeMets et al. (1990). The focus of the study is on analyzing these discrepancies for global plate tectonics as well as regional tectonic settings. The method used also allows us not only to derive horizontal motion for individual stations but also to calculate Euler vectors for those plates that have enough stations located on the stable interior like North America, Pacific, Eurasia, and Australia.

Harrison, Christopher G. A.↗
NASA NTRS1991

Wave interactions with continuous fractal layers

Many natural structures possess self-similar multiscales which can be characterized by power law spectra. Under appropriate conditions, knowledge of the strength of these scale sizes provides information on the physical processes which formed these objects. In this paper, we investigate wave interactions with continuous fractal layers which model geological and variegated structures. Since fractal characteristics of the layers are embedded in the scattered field, they can be retrieved under appropriate conditions. This inversion can be performed in either the frequency or the time domain as desired.

Kim, Y.↗
NASA NTRS1996

An Integrated Approach To Payload System Simulation

This paper describes a payload system simulation implemented at JPL as part of a comprehensive mission simulation facility. The flight software function includes communication with other process modules, instrument control, and data management. The payload system simulation software consists of: a camera subsystem, a virtual world, and a mission visualization toolset.

dynamics model photometric model geological model ↗
NASA NTRS2013

The Multi-Stage History of Mt. Sharp

The Curiosity rover is exploring Gale crater and Mt. Sharp, Gale's 5‐km high central mound. We are investigating the history of alteration and erosion of Mt. Sharp using orbital imagery, spectroscopy and rover observations. Our results suggest a significant time gap between emplacement of the upper and lower sections of the mound. Crater counts show that the lower mound was formed soon after Gale itself, and that it contains distinct units ranging in altitude from approximately ‐4,500 to ‐1,800 m. Spectral data suggest that many units contain phyllosilicates. We found that these clay‐bearing rocks occur in distinct layers concentrated below ‐2,900 m. Parts of the lower mound exhibit a transition from clays to sulfates with increasing altitude. The lower mound shows evidence of flowing water, including canyons and inverted channels. Wind erosion produced km‐scale yardangs and scalloped cliffs. Our mapping shows that many yardangs in the lower mound are clay‐bearing, with a predominant orientation of around N‐S. Curiosity's ground‐level images show myriad fine‐scale, mainly horizontal layers in the lower mound. The rover has found stream beds and conglomerates, indicating that water once flowed on the crater floor. Drilling near the deepest point in Gale produced abundant clay, providing additional evidence of aqueous alteration. Upper mound units range in altitude from ‐2,100 m to +500 m, and mantle the lower mound above an angular unconformity. Most upper mound units are composed of layers. The formation age of the upper mound is unknown, since few craters are preserved. Clay‐bearing layers are detectable in several locations, mainly at altitudes near ‐2,000 m. There is no evidence of water flow, but wind erosion has scalloped the surfaces and edges of layers, and fine‐scale yardangs are common. Correlations between yardangs and clay spectra are apparent only in the lowermost units of the upper mound. Yardang orientations vary, and include N‐S, NW‐SE, and NE‐SW. Upper mound units resemble the planet‐wide Medusae Fossae formation, dated as Hesperian and argued to be composed of ignimbrites. Medusae Fossae layers are easily eroded by wind, and our mapping demonstrates their resemblance to upper mound fine‐scale yardangs. The history of Mt. Sharp started with deposition and lithification of sediments shortly after crater formation. Some lower mound layers were partially altered to clays and sulfates, and water formed streams and canyons. Wind erosion of the lower mound produced large‐scale yardangs, particularly in clay‐rich layers, oriented generally N‐S. Upper mound units were emplaced following a considerable period of wind erosion. The absence of water flow on the upper mound suggests that these units were emplaced after atmospheric loss rendered water unstable at the surface. The shift in dominant wind direction, as indicated by yardang orientations, also argues for a time gap between erosion of the lower and upper mound. These observations are consistent with upper mound units being related to the Hesperian Medusae Fossae formation. During 2014 Curiosity is expected to reach the foot of Mt. Sharp and ascend through the clay‐rich layers, into the sulfate‐rich layers, and possibly past the interface with the upper mound. This will be a unique opportunity to field check geologic models on the surface of Mars.

Allen, C.↗
DOE OSTI2013

3D characterization of the Astor Pass geothermal system, Nevada

The Astor Pass geothermal system resides in the northwestern part of the Pyramid Lake Paiute Reservation, on the margins of the Basin and Range and Walker Lane tectonic provinces in northwestern Nevada. Seismic reflection interpretation, detailed analysis of well cuttings, stress field analysis, and construction of a 3D geologic model have been used in the characterization of the stratigraphic and structural framework of the geothermal area. The area is primarily comprised of middle Miocene Pyramid sequence volcanic and sedimentary rocks, nonconformably overlying Mesozoic metamorphic and granitic rocks. Wells drilled at Astor Pass show a ~1 km thick section of highly transmissive Miocene volcanic reservoir with temperatures of ~95°C. Seismic reflection interpretation confirms a high fault density in the geothermal area, with many possible fluid pathways penetrating into the relatively impermeable Mesozoic basement. Stress field analysis using borehole breakout data reveals a complex transtensional faulting regime with a regionally consistent west-northwest-trending least principal stress direction. Considering possible strike-slip and normal stress regimes, the stress data were utilized in a slip and dilation tendency analysis of the fault model, which suggests two promising fault areas controlling upwelling geothermal fluids. Both of these fault intersection areas show positive attributes for controlling geothermal fluids, but hydrologic tests show the ~1 km thick volcanic section is highly transmissive. Thus, focused upwellings along discrete fault conduits may be confined to the Mesozoic basement before fluids diffuse into the Miocene volcanic reservoir above. This large diffuse reservoir in the faulted Miocene volcanic rocks is capable of sustaining high pump rates. Understanding this type of system may be helpful in examining large, permeable reservoirs in deep sedimentary basins of the eastern Basin and Range and the highly fractured volcanic geothermal reservoirs in the Snake River Plain.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

EERC Final Topical Report – Findings on Subtask 3.2 – Produced Water Management Through Geologic Homogenization, Conditions, and Reuse

The Energy & Environmental Research Center (EERC) was awarded a contract by the North Dakota Industrial Commission (NDIC) Oil and Gas Research Program (NDIC No. G-051-101) to conduct a study on the recycling of water used in oil and gas operations, also known as produced water, from oil- and gas-producing regions of North Dakota as directed by Section 19 of North Dakota House Bill 1014. This final report provides a compilation of results of the study, which include regulatory, scientific, technological, and feasibility methods and considerations associated with North Dakota produced water management. The report provides a synopsis of this project’s previously submitted produced water assessment report entitled “Produced Water Management and Recycling Options in North Dakota” (Energy & Environmental Research Center, 2020), with updated values provided as appropriate. The report provides the results from the investigation of a novel produced water management strategy, referred to as geologic homogenization, conditioning, and reuse (GHCR), which aims to use a subsurface geologic formation as a natural medium for managing produced water recycling and reuse. Water management is a significant technical and economic challenge for sustainable oil and gas production, and water volumes are intrinsically linked to oil production volumes. North Dakota oil production rose to over 1.5 million barrels (MMbbl)/day in 2019, and despite a downturn in oil price in early 2020, North Dakota oil production has recovered to 1.1 MMbbl/day as of August 2021. Bakken petroleum system development between 2008 and 2020 has resulted in a nearly fourfold increase in produced water volumes to 642 MMbbl/yr in 2020 after peaking at 740 MMbbl/yr in 2019 and a fivefold increase in saltwater disposal (SWD) volumes to 565 MMbbl/yr in 2020 after a peak of 682 MMbbl/yr in 2019. Produced water and SWD volumes are forecasted to double by 2030. SWD is the primary method of produced water management used in North Dakota, with approximately 95% of the SWD volume occurring through subsurface injection into sandstones of the Dakota Group (Dakota). Localized pressurization of the Dakota resulting from SWD and projected increases in produced water volumes could impact the economics of North Dakota oil production. As a result, there is an emerging need to pursue alternative produced water management approaches, including recycling and reuse. While produced water recycling is not yet widespread, commercial operators are making strides in overcoming the technical challenges of using high salinity produced water in completion operations (Marathon Oil, 2020). As water management continues to be a key focal point in companies’ environmental, social, and governance (ESG) initiatives, focus on water management, including recycling, will likely continue to increase. Laboratory column testing, field sample collection, geologic modeling and numerical simulation, and techno-economic analysis all indicate that GHCR could feasibly be implemented as a potential water management option. Laboratory column testing and field sample collection indicate that the Inyan Kara sandstone and native formation fluid are capable of homogenizing with the Bakken produced water to a point where the fluid composition appears to stabilize. Extracting that stabilized fluid could be considered homogeneous and capable of providing individual batches of hydraulic fracturing fluid. Numerical simulation results indicate that extraction of fluids from the Inyan Kara in a GHCR implementation scenario is capable of reducing formation pressure, which would help ease localized pressurization of the Inyan Kara and extend the available capacity for nearby existing SWD wells. Economic analysis indicates that there are scenarios where GHCR implementation can be a competitive or even lower-cost option than a conventional water management approach. Site-specific conditions will dictate the economic potential of GHCR, but potentially attractive sites for GHCR implementation will be those that are located above a pressurized zone of the Inyan Kara, need six or more Bakken infill wells, and face high costs for conventional SWD and/or freshwater. Based on the regulatory review, drilling into the Inyan Kara for SWD and to harness as a source water for industrial use have precedent in North Dakota, and a workable regulatory solution for GHCR seems likely. However, restrictions in the state regarding surface storage and transport of produced fluids may limit some activities, which will affect how GHCR could ultimately be implemented. In summary, this study reveals pursuing GHCR can be a viable approach to water management in North Dakota. The GHCR concept addresses some of the challenges that hinder the more traditional approaches to recycling in the industry. Furthermore, an assessment of the current landscape of water management within the state reveals the ongoing trend of increasing volumes of produced water and SWD. Projections reveal that the volumes of produced water that need to be managed are expected to double over the next decade (Energy & Environmental Research Center, 2020). With the continued development of the Bakken and continuing driving factors related to ESG initiatives, implementing a practice such as GHCR is a feasible approach to adding recycling of produced water to industry within the state. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission Oil and Gas Research Program. References Energy & Environmental Research Center, 2020, Produced water management and recycling options in North Dakota: Final Report for North Dakota Legislative Management Energy Development and Transmission Committee and North Dakota Industrial Commission. Marathon Oil, 2020, Sustainability report: https://cdn.sanity.io/files/ghcnw9z2/website/ 91744eb6ef8fbe59505a911c6b8d2e8dd9a537fa.pdf?dl (accessed November 2021).

02 PETROLEUM↗
DOE OSTI2023

Moment tensor reconstruction

A seismic monitoring system includes a plurality of seismic monitors and a processing device operatively coupled to the plurality of seismic monitors. The processing device receives recordings of waveforms of motion detected at the plurality of seismic detectors in a geographic area. The processing device applies the respective recordings to corresponding positions of the seismic detectors in a three-dimensional geological model that describes its elastic attributes and tests a plurality of moment tensors at a plurality of locations. Based on the testing, the processing device determines a globally convergent source location and moment tensor in the three-dimensional model based on the testing.

Petrov, Petr↗
DOE OSTI2024

Geothermal Reservoir Simulation Analysis in Support of Electricity Co-Production Feasibility Study at the Blackburn Oil Field, Nevada: Preprint

Geothermal electricity co-production is a viable option for oil reservoirs producing large water cuts with elevated wellhead-observed temperatures. Repurposing existing oil wells significantly reduces initial investment costs historically associated with geothermal resource utilization. The National Renewable Energy Laboratory (NREL), partnering with Gradient Geothermal, Inc. (formerly known as Transitional Energy) and Grant Canyon Oil and Gas, has been tasked to evaluate the feasibility of geothermal electricity co-production at the Blackburn Oil Field with Organic Rankine Cycle (ORC) generators. The Devonian steady-state reservoir has historically been producing high water cuts of 240 degrees F (115.6 degrees C) observed at the wellhead without documented pressure drawdown or thermal breakthrough. An estimated initial reservoir temperature of approx. 260degreesF (126.7 degrees C) has been observed in the field and history-matched in a wellbore production analysis and reservoir simulation. Our objective was to develop a conceptual geological model of the subsurface, simulate a natural-state reservoir, model production scenarios, and complete a technical feasibility analysis to accomplish this task. Through extensive modeling and the use of available proprietary and public data, it was possible simulate three scenarios that indicated minimal thermal decline over the duration of a simulated ten-year production and re-injection scheme.

Blackburn Nevada↗
DOE OSTI2025

Geothermal Reservoir Simulation Analysis in Support of Electricity Co-Production Feasibility Study at the Blackburn Oil Field, Nevada

Geothermal electricity co-production is a viable option for oil reservoirs producing large water cuts with elevated wellhead temperatures. Repurposing existing oil wells significantly reduces initial investment costs historically associated with geothermal resource utilization. The National Renewable Energy Laboratory (NREL), partnering with Gradient Geothermal, Inc. (formerly known as Transitional Energy) and Grant Canyon Oil and Gas, has been tasked to evaluate the feasibility of geothermal electricity co-production at the Blackburn Oil Field with Organic Rankine Cycle (ORC) generators. The Devonian steady-state reservoir has historically been producing high water cuts of 240 degrees F (115.6 degrees C) observed at the wellhead without documented pressure drawdown or thermal breakthrough. An estimated initial reservoir temperature of approx. 260 degrees F (126.7 degrees C) has been observed in the field and history-matched in a wellbore production analysis and reservoir simulation. Our objective was to develop a conceptual geological model of the subsurface, simulate a natural-state reservoir, model production scenarios, and complete a technical feasibility analysis to accomplish this task. Through extensive modeling and the use of available proprietary and public data, it was possible simulate three scenarios that indicated minimal thermal decline over the duration of a simulated ten-year production and re-injection scheme.

Blackburn Nevada↗
DOE OSTI2013

Sensitivity of injection costs to input petrophysical parameters in numerical geologic carbon sequestration models

Numerical simulations are widely used in feasibility studies for geologic carbon sequestration. Accurate estimates of petrophysical parameters are needed as inputs for these simulations. However, relatively few experimental values are available for CO2-brine systems. Hence, a sensitivity analysis was performed using the STOMP numerical code for supercritical CO2 injected into a model confined deep saline aquifer. The intrinsic permeability, porosity, pore compressibility, and capillary pressure-saturation/relative permeability parameters (residual liquid saturation, residual gas saturation, and van Genuchten alpha and m values) were varied independently. Their influence on CO2 injection rates and costs were determined and the parameters were ranked based on normalized coefficients of variation. The simulations resulted in differences of up to tens of millions of dollars over the life of the project (i.e., the time taken to inject 10.8 million metric tons of CO2). The two most influential parameters were the intrinsic permeability and the van Genuchten m value. Two other parameters, the residual gas saturation and the residual liquid saturation, ranked above the porosity. These results highlight the need for accurate estimates of capillary pressure-saturation/relative permeability parameters for geologic carbon sequestration simulations in addition to measurements of porosity and intrinsic permeability.

Cost estimate↗
DOE OSTI2018

Final Report: Optimal Model Complexity in Geological Carbon Sequestration: A Response Surface Uncertainty Analysis

The critical component of a risk assessment study in evaluating GCS is an analysis of uncertainty in CO2 modeling. In such analyses, direct numerical simulation of CO2 flow and leakage requires many time-consuming model runs. Alternatively, analytical methods have been developed which allow fast and efficient estimation of CO2 storage and leakage, although restrictive assumptions on formation rock and fluid properties are employed. In this study, an intermediate approach is proposed based on the Design of Experiment and Response Surface methodology, which consists of using a limited number of numerical simulations to estimate a prediction outcome as a combination of the most influential uncertain site properties. The methodology can be implemented within a Monte Carlo framework to efficiently assess parameter and prediction uncertainty while honoring the accuracy of numerical simulations. The choice of the uncertain properties is flexible and can include geologic parameters that influence reservoir heterogeneity, engineering parameters that influence gas trapping and migration, and reactive parameters that influence the extent of fluid/rock reactions. The method was tested and verified on modeling long-term CO2 flow, non-isothermal heat transport, and CO2 dissolution storage by coupling two-phase flow with explicit miscibility calculation using an accurate equation of state that gives rise to convective mixing of formation brine variably saturated with CO2. All simulations were performed using three-dimensional high-resolution models including a target deep saline aquifer, overlying caprock, and a shallow aquifer. To evaluate the uncertainty in representing reservoir permeability, sediment hierarchy of a heterogeneous digital stratigraphy was mapped to create multiple irregularly shape stratigraphic models of decreasing geologic resolutions: heterogeneous (reference), lithofacies, depositional environment, and a (homogeneous) geologic formation. To ensure model equivalency, all the stratigraphic models were successfully upscaled from the reference heterogeneous model for bulk flow and transport predictions (Zhang & Zhang, 2015). GCS simulation was then simulated with all models, yielding insights into the level of parameterization complexity that is needed for the accurate simulation of reservoir pore pressure, CO2 storage, leakage, footprint, and dissolution over both short (i.e., injection) and longer (monitoring) time scales. Important uncertainty parameters that impact these key performance metrics were identified for the stratigraphic models as well as for the heterogeneous model, leading to the development of reduced/simplified models at lower characterization cost that can be used for the reservoir uncertainty analysis. All the CO2 modeling was conducted using PFLOTRAN – a massively parallel, multiphase, multi-component, and reactive transport simulator developed by a multi-laboratory DOE/SciDAC (Scientific Discovery through Advanced Computing) project (Zhang et al., 2017, in review). Within the uncertainty analysis framework, increasing reservoir depth were investigated to explore its effect on the uncertainty outcomes and the potential for developing gravity-stable injection with increased storage security (Dai et al., 20126; Dai et al., 2017, in review). Finally, to accurately model CO2 fluid-rock reactions and resulting long-term storage as secondary carbonate minerals, a modified kinetic rate law for general mineral dissolution and precipitation was proposed and verified that is invariant to a scale transformation of the mineral formula weight. This new formulation will lead to more accurate assessment of mineral storage over geologic time scales (Lichtner, 2016).

58 GEOSCIENCES↗
DOE OSTI2023

Unsaturated alluvium disposal modelling with improved geological realism

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD). The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and system performance assessment (PA) models in several potential host-rock environments, and to simulate and analyze these generic repository concepts and models using the GDSA Framework toolkit and other tools as needed. The specific GDSA goal addressed in this report is reference case development, simulation, and analysis for the unsaturated alluvium (UZ), one of the four potential host-rocks considered by the GDSA. Further, we aim to exercise the simulation tools and methodologies under development by GDSA for PA modelling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Identification of Faults Susceptible to Induced Seismicity (Final Report)

Central to the work documented in this report is the capability of geocellular models to represent the geologic conceptual model updated with fault identification from machine learning and joint inversion modeling of microseismic data measured and recorded as a consequence of CO 2 injection at a field demonstration site: the Illinois Basin - Decatur Project (IBDP). This work required seven unique geocellular models with 100s of simulated variations to gain a very high degree of confidence in the identification of geologic features present that contributed to induced microseismicity at IBDP. All forward modeling: pressure modeling, stress modeling, and seismic modeling used the same geologic conceptual model and representations of that model at different scales. The pressure modeling and poroelastic modeling created “snapshots” of pore pressure and stress field changes at different times during CO 2 injection, in which microseismic events were clustered (in time). These pressure and stress snapshots, within the framework and architecture of the geologic conceptual model via the geocellular model, informed the single fault and fault network models to ascertain the likelihood of fault movement (seismic or aseismic). The outcomes of the pressure, stress, and fault/fault network (seismic) modeling confirmed that the faults in the geologic conceptual model in Task 2 were likely the source of microseismic events measured at IBDP and acted as conduits for pressure to be transmitted from the injection interval into the Precambrian crystalline basement rock. This closely coordinated and integrated unique modeling approach was conducted to prove the viability of our proposed workflow 1) to better resolve crystalline basement faults, 2) detect subseismic faults that could be activated by injection, 3) increase the certainty in fault detection and their susceptibility to release seismic energy, and 4) understand transmission of pressure vertically from the well to the underlying fractured crystalline basement. The proposed methodology was effective in guiding an iterative process of calibrating forward modeling results based on similar geocellular models while honoring the geologic conceptual model (i.e., characterization data and knowledge of regional geology); this led to higher level of certainty in the identification of fault/faults zones to control seismicity and transmission of pressure to the regions of recorded and located injection induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2012

Development of an Advanced Performance Assessment Modeling Capability for Geologic Disposal of Nuclear Waste: Methodology and Requirements

This report describes the planning and initial development of an advanced disposal system PA modeling capability to facilitate the science-based evaluation of disposal system performance for a range of fuel cycle alternatives in a variety of geologic media and generic disposal system concepts. The advanced modeling capability will provide a PA model framework that facilitates PA model development, execution, and evaluation within a formal PA methodology. The PA model framework will provide a formalized structure that enables (a) representation and implementation of a range of generic geologic disposal system options, (b) representation of subsystem processes and couplings at varying levels of complexity in an integrated disposal system model, (c) flexible, modular representation of multi-physics processes, including the use of legacy codes, (d) evaluation of system- and subsystem-level performance, (e) uncertainty and sensitivity analyses to isolate key subsystem processes and components, (f) data and configuration management functions, and (g) implementation in HPC environments.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2018

Constraining the Climate and Ocean PH of the Early Earth with a Geological Carbon Cycle Model

The early Earth's environment is controversial. Climatic estimates range from hot to glacial, and inferred marine pH spans strongly alkaline to acidic. Better understanding of early climate and ocean chemistry would improve our knowledge of the origin of life and its coevolution with the environment. Here, we use a geological carbon cycle model with ocean chemistry to calculate self-consistent histories of climate and ocean pH. Our carbon cycle model includes an empirically justified temperature and pH dependence of seafloor weathering, allowing the relative importance of continental and seafloor weathering to be evaluated. We find that the Archean climate was likely temperate (0-50deg C) due to the combined negative feedbacks of continental and seafloor weathering. Ocean pH evolves monotonically from 6.6(sup + 0.6) (sub -0.4) (2 sigma) at 4.0 Ga to 7.0(sup +0.7) (sub-0.5) (2 sigma) at the Archean-Proterozoic boundary, and to 7.9(sup +0.1) (sub -0.2) (2 sigma) at the Proterozoic-Phanerozoic boundary. This evolution is driven by the secular decline of pCO2, which in turn is a consequence of increasing solar luminosity, but is moderated by carbonate alkalinity delivered from continental and seafloor weathering. Archean seafloor weathering may have been a comparable carbon sink to continental weathering, but is less dominant than previously assumed, and would not have induced global glaciation. We show how these conclusions are robust to a wide range of scenarios for continental growth, internal heat flow evolution and outgassing history, greenhouse gas abundances, and changes in the biotic enhancement of weathering.

Krissansen-Totton, Joshua↗