SEARCH · Engineering Papers
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.
The project established a Tight Oil Field Laboratory to address technical challenges in developing stacked unconventional reservoirs in the Powder River Basin. Key activities included data compilation, subsurface mapping, drilling, logging, coring, deployment of fiber optics and microseismic, completion and stimulation optimization, and well performance evaluation.
Abstract The Congo rainforest located in central equatorial Africa is an important, yet understudied part of the globe surrounded by complex orographic features. A primitive understanding of precipitation processes such as mesoscale convective dynamics magnifies uncertainties in the future climate projections of the hydrological cycle over the Congo. Furthermore, the effects of orography, which is an important forcing for convection and precipitation, are poorly resolved by climate models, and ill‐conceptualized over the Congo. To address this knowledge gap, perturbed orographic forcing experiments are conducted using the Weather Research and Forecasting (WRF) mesoscale numerical model in a high‐resolution convection‐permitting model setup. The model simulated selected dates in November 2014. The thunderstorms and rainfall simulated in the control run for the case study analyzed in this article compared reasonably well to satellite‐derived brightness temperature and rainfall data. The results from this case study show that the dynamical impact of increasing the height of the East African Highlands is the blocking of the lower‐tropospheric tropical easterlies. This weakening of the lower‐tropospheric zonal winds increases the windshear over the Congo Basin resulting in slower propagating, more intense mesoscale convective systems with enhanced rainfall.
Abstract Hatchery production of Chinook Salmon Oncorhynchus tshawytscha in the Columbia River basin comprises most of the anadromous salmonid production in this region. Hatchery facilities and programs serve to mitigate for impacts to salmonids due to the construction and operation of hydropower dams and habitat impacts from development in addition to the conservation and restoration of natural populations. A genetic method referred to as parentage-based tagging (PBT) enables highly reliable detection of hatchery-origin fish and inference of multigeneration pedigrees. This study compiles 11 years of PBT data from nearly 125,000 interior stream-type Chinook Salmon from 24 spawning hatcheries located in tributaries of the mid- and upper Columbia River as well as the Salmon, Clearwater, and Grande Ronde River subbasins. Multigenerational pedigrees allowed investigation of the proportion of natural- and hatchery-origin broodstock (pNOB and pHOB, respectively) for each hatchery and enumeration of the scale of production between segregated and integrated programs. We then compared how pHOB and the scale of production influenced the number of stray fish observed, genetic diversity, relatedness, and age-class compositions within broodstocks. Over 91.0% of hatchery broodstock could be assigned back to their parents, and there was an overall rate of less than 1.0% of broodstock that were unintentionally incorporated into nonnatal hatchery programs. We evaluated 14 segregated programs with 0.0–10.0% pNOB, 7 intermediate programs (10.1–50.0% pNOB), and 3 integrated programs (>50.0% pNOB). There was no correlation between the scale of production or pNOB with the level of genetic diversity, but as production size increased, so did the effective number of breeders. This study demonstrates the utility of PBT as a monitoring tool for hatchery broodstocks, and results suggest that segregated and integrated programs have tradeoffs that generally align with their intended broodstock management purpose of providing fish for harvest and/or fish for supplementation or reintroduction.
Abstract Hatchery production of Chinook Salmon Oncorhynchus tshawytscha in the Columbia River basin comprises most of the anadromous salmonid production in this region. Hatchery facilities and programs serve to mitigate for impacts to salmonids due to the construction and operation of hydropower dams and habitat impacts from development in addition to the conservation and restoration of natural populations. A genetic method referred to as parentage-based tagging (PBT) enables highly reliable detection of hatchery-origin fish and inference of multigeneration pedigrees. This study compiles 11 years of PBT data from nearly 125,000 interior stream-type Chinook Salmon from 24 spawning hatcheries located on tributaries of the mid- and upper Columbia River and in the Salmon, Clearwater, and Grande Ronde River subbasins. Multigenerational pedigrees allowed for investigation of the proportions of natural- and hatchery-origin broodstock (pNOB and pHOB, respectively) for each hatchery and enumeration of the scale of production between segregated and integrated programs. We then compared how pHOB and the scale of production influenced the number of stray fish observed, genetic diversity, relatedness, and age-class compositions within broodstocks. Over 91.0% of hatchery broodstock could be assigned back to their parents, and overall less than 1.0% of broodstock consisted of fish that were unintentionally incorporated into nonnatal hatchery programs. We evaluated 11 segregated programs with 0.0–10.0% pNOB, 9 intermediate programs (10.1–50.0% pNOB), and 3 integrated programs (>50.0% pNOB). There was no correlation between the scale of production or pNOB with the level of genetic diversity, but as production size increased, so did the effective number of breeders. This study demonstrates the utility of PBT as a monitoring tool for hatchery broodstocks, and results suggest that segregated and integrated programs have tradeoffs that generally align with their intended broodstock management purpose of providing fish for harvest and/or fish for supplementation or reintroduction.
Riparian enhancement is a common restoration technique in the Columbia River Basin (CRB) and the Pacific Northwest. However, relatively few studies have evaluated its success and even fewer studies include long-term monitoring. Forty-one riparian planting projects located in the CRB, each with a paired treatment and control reach, were evaluated in the summer of 2018 and 2019 using an extensive post-treatment (EPT) design. At each reach, we quantified woody plant abundance, richness, diversity, and vegetation and canopy cover as riparian response variables, and measured terrace height and documented individuals with bud browse, deceased, and with predator protection as potential explanatory variables. Species richness, woody plant abundance (all height classes combined and within shrub height class), and the proportion of woody plants with bud browse or deceased were all higher in treatment than control reaches. However, no significant improvements were observed for any other riparian response or explanatory variables. Implementation techniques, including invasive species removal efforts, follow-up planting, watering, and the type of restoration project (floodplain restoration and planting, in-stream restoration and planting, or solely planting) were also considered as potential explanatory variables and were found to have significant impacts on restoration response for abundance and cover variables. Further, our results suggest that planting implementation methods, site level physical factors, and time since planting all influence the success of riparian planting projects and if not addressed in restoration design and implementation can contribute to a lack of detectable response.
Underground Hydrogen Storage (UHS) in depleted oil and gas reservoirs could provide a cost-effective solution to balance seasonal fluctuations in renewable energy generation. However, data and knowledge on UHS at subsurface conditions are limited so it is difficult to estimate how effective this type of storage could be. In this study, we perform high pressure experiment to measure the effectiveness of cyclic hydrogen (H 2 ) storage in a specimen of Temblor sandstone retrieved from the San Joaquin Basin of California. Our experiment mimics reservoir pressure conditions to measure H 2 -brine relative permeability and fluid-rock interactions over the course of ten charging and discharging cycles. Initial gas breakthrough occurred at 15 % to 25 % H2 saturation in the specimen with 3 % NaCl brine as the resident fluid. Continuing injecting to 4 pore volumes (PV) of H 2 yielded an asymptotic H 2 saturation of 38 % to 41 %, a level often referred to as the irreducible gas saturation based on two-phase flow. The boundary condition in this study mimics the near wellbore region, which experiences bi-directional H 2 flow. This bi-directional flow led to evaporative drying of the specimen resulting in 94 % H 2 saturation at the end of 10th cycle. This indicates that cyclic flow and evaporative drying can lead to more efficient reservoir storage where a larger fraction of the reservoir porosity is usable to store H 2 . The produced gas stream consisted of H 2 mixed with 8 % to 22 % H 2 O, indicating formation dry-out by evaporation. Meanwhile, produced water chemistry indicated calcite and silicate dissolution, with calcite sourced from fossil fragments. This led to a loss of cementation and weakened the rock sample. Combined, our results indicate dry-out, compaction, increased H 2 saturation, rock weakening, and permeability loss during cyclic UHS. Overall, we anticipate that the combined effects should lead to higher than anticipated UHS storage efficiency per volume of sandstone reservoir rock.
Active reservoir management (ARM) through brine extraction can reduce pressure buildup during large-scale implementation of carbon capture and storage (CCS) projects. This study used machine learning (ML)-assisted approaches to analyze bottomhole pressure (BHP) responses to various brine injection and extraction scenarios. Field monitoring data were collected over a 2-year operation period at two injection wells and one extraction well (about 400 meters away) as part of a Brine Extraction and Storage Test (BEST) in the North Dakota portion of the Williston Basin. Injection activities increased the BHPs at the injection wells by around 0.70 MPa (~100 psi) during the operation period. Extraction activities demonstrated the capability to decrease the BHPs at the injection wells by approximately 0.21–0.34 MPa (30–50 psi) depending on the ratio of the extraction and injection well flow rates (the “extraction ratio” – a normalization procedure used in the analysis). The pressure reduction provided by the extraction well equated to 30%–50% of the pressure buildup at the injection well. Furthermore, this work shows how ML analytics can play a key role in estimating reservoir pressure responses for complex injection and extraction activities during pressure management practices of CCS projects.
Oil production in the Permian Basin gives rise to large volumes of produced water contaminated by silt, emulsified oil, and additives used for enhanced oil recovery. There is intense interest in the design of membrane modules as sustainable alternatives for produced water treatment to enable the reuse of produced water for agricultural applications, injection into aquifers, and redeployment in oil recovery. Here, we report a hierarchically textured cement-based membrane exhibiting orthogonal wettability, specifically, superhydrophilic and underwater superoleophobic characteristics. The in situ formation of ettringite needles accompanied by embedding of glass spheres imbues multiscale texturation to stainless-steel mesh membranes, enabling the separation of silt and oil from produced water at high flux rates (1600 L h -1 ∙m -2 , at ca. 2.7 bar). Oil concentration is reduced as low as 1 ppb with an overall separation efficiency of 99.7% in singlepass filtration. The membranes show outstanding mechanical resilience and retention of performance across multiple cycles.
Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.
Lake Erie is threatened by eutrophication and harmful algal blooms due to excess nutrient loading from agricultural sources. Agricultural conservation practices (ACPs) have been developed and implemented to reduce nutrient losses but estimating ACP effectiveness is challenging. The Soil and Water Assessment Tool (SWAT) has been used to investigate ACP effectiveness for water quality improvement. Many SWAT applications have been developed by different investigators to evaluate ACP effectiveness for reducing nutrient, particularly phosphorus (P), loading in the agriculturally-dominated Western Lake Erie Basin (WLEB). Our objective is to establish what has been achieved by past modeling research and make suggestions for future applications and improvements. We synthesized the findings of 28 SWAT modeling studies within the WLEB. Models generally performed satisfactorily against accepted criteria for streamflow and sediment, but performance for P loads, like soluble reactive P, was mostly “unsatisfactory”. The “unsatisfactory” performance maybe due to imperfections and idealizations in model formulations and/or parameterization. Thus, simulations of P transport and transformation processes need improvement. In addition, model parameter selection is the key part of model set-up. Most SWAT modeling studies used default values during initial set-up, then performed calibration and validation. It was found that the calibrated P related parameter values varied widely across different studies, even within the same watershed with some values unrealistic for the study areas. The phenomena of different combinations of model parameters producing similar outputs indicates equifinality. Equifinality in the baseline model may impact results when ACPs are incorporated. Furthermore, the unrealistic values used in ACP assessment undermine the credibility of ACP effectiveness. Future model applications should try to re-examine the calibrated P parameters and make sure they are realistic for the study area as well as reduce equifinality by constraining the model with characterization of watershed conditions, better understanding of hydrologic processes, and parameter values based on real-world observations. As a result, future model applications should focus on improving P transport and transformation processes, using measured watershed characteristics for parameterization, and improving reflections of climate change, which could result in more accurate assessments of ACP effectiveness to meet targeted goals.
Accurate estimates of nutrient loads are necessary to identify critical source areas and quantify the impact of management practices on pollutant export. Previous studies have investigated nutrient load estimate uncertainty, but they often focus on nutrient loads estimated using an interpolation method for large-scale watersheds with short-term datasets. The study objective was to quantify uncertainty in soluble reactive phosphorus (SRP), total phosphorus (TP), and suspended solids (SS) load estimates from two small (<10 3 km 2 ) agricultural watersheds in the western Lake Erie Basin resulting from different sampling frequencies. Each watershed had high temporal resolution datasets of discharge (15 min) and nutrient concentration (1 to 3 samples per day) collected over a 30-year period (1990–2020). Firstly, SRP, TP, and SS loads were calculated using the high temporal resolution datasets, which was assumed as “true loads”. Secondly, the high temporal concentration data were decomposed to semiweekly, weekly, biweekly, and monthly sampling and annual loads were estimated using four common load estimation methods to assess the effect of sampling frequency and load estimation method on load estimate error. Across the four different methods, the composite method had the lowest relative root mean square and absolute bias, but the rectangular interpolation method was the most precise. Furthermore, even with semiweekly sampling, the composite method resulted in an unacceptable level of precision (average imprecision = 39 %), while the interpolation method resulted in an unacceptable bias (average absolute bias = 16 %). Because neither method could provide acceptable accuracy and precision at the lowest decrease in sampling (e.t. semiweekly sampling), continued daily sampling is recommended in these watersheds.
Plant transpiration plays a critical role in global water and energy cycles, requiring better process understanding as climate change intensifies drought stress and alters plant responses. Most hydrological models such as the widely-used SWAT lack representation of plant hydraulics, the mechanistic processes controlling plant water regulation and transpiration. Here, this study developed SWAT-PHS by integrating a plant hydraulics scheme (PHS) into SWAT hydrological model, enabling explicit simulation of root water uptake, sap flow, storage and transpiration at 30-minute timescales for watershed-scale modeling. In the Hanjiang River Basin, SWAT-PHS mitigated overestimation of runoff during the rainy season and underestimation during the dry season, reducing the overall simulation error by 29% across the entire simulation period. The model can simulate reasonable plant water dynamics, including diurnal transpiration patterns and drought responses showing declining transpiration flux, hydraulic buffering through stem water storage, and depth-dependent root water uptake strategies. Sensitivity analysis shows that SWAT-PHS captured mechanistic relationships between plant hydraulic traits and transpiration, with root distribution and stem capacitance positively affecting annual transpiration while vulnerability parameters showed negative effects. This work provides a pathway for improving hydrologic modeling and water resource management by better representing plant water regulation under climate change and expected intensifying water stress conditions.
The complex topography of submarine canyons may result in different composition of benthic foraminifera assemblages. To understand how trophic, hydrological and sedimentological conditions in submarine canyons can influence the distribution of benthic foraminifera, and to use this information to corroborate paleoenvironmental interpretations for the Holocene, we investigated recent benthic foraminiferal assemblages (total fauna >63 μm) and sedimentological data in two canyons (CANWN and CAND) in the Espírito Santo Basin (ESB) between 18°20' and 21°20' S. Surface sediment samples (0–2 cm) were collected inside the canyons (150 to 1300 m water depth) and in shelf-slope adjacent transects (50 to 1300 m water depth). The density, taxonomic diversity, and assemblage composition of benthic foraminifera change with depth and location. The distinct ecological preferences of the most abundant taxa allowed us to recognize five benthic foraminiferal groups. Three groups (V, III, and I) are present in different bathymetric sectors; Group V: outer shelf (50 m), Group III: upper, and middle – lower slope (150 to 400 m), and Group I: middle – lower slope (1000 to 1300 m). Groups II and IV show no characteristic bathymetric distribution and are present only in CAND and in CANWN, respectively. Group V consists of Hanzawaia boueana, Peneroplis planatus, and Quinqueloculina lamarckiana; Group III is dominated by Globocassidulina rossensis and Trifarina spp.; Group I consists of Globocassidulina crassa, Bolivina lowmani, Gavelinopsis versiformis, Alabaminella weddellensis, and Epistominella exigua. The main species in Group II (CAND, 150, 1000 to 1300 m) are Trifarina angulosa, Globocassidulina subglobosa, and Discorbis vilardeboanus. Group IV (middle – lower CANWN, 1000 to 1300 m), consists mainly of agglutinated species Glomospira charoides, Rhabdammina abyssorum, and Psammosphaera fusca. Further, our data suggest that the quantity (and quality) of food supply, hydrodynamic conditions and sediment properties are the main drivers controlling the bathymetric distribution of benthic foraminiferal assemblages in both canyons. The middle – lower CANNW revealed unstable trophic conditions, related to terrigenous sediment input due to turbidity currents. In CAND, the foraminiferal assemblages ecology indicated sufficient organic matter supply that favors species establishment and diversity, indicating a more productive and less unstable environment than in CANWN.
This study reports the application of Soil andWater Assessment Tool (SWAT) within the Hydrologic andWater Quality System(HAWQS) on-line platform, for the UpperMississippi River Basin (UMRB). TheUMRB is an important ecosystem located in the north central U.S. that is experiencing a range of ecological stresses. Specifically, testing of SWAT was performed for: (1) Hargreaves (HG) and Penman-Monteith (PM) PET methods, and (2) Livneh, National Climatic Data Center (NCDC) and Parameter-elevation Regressions on Independent Slopes Model (PRISM) climate datasets. The Livneh-PM combination resulted in the highest average annual water yield of 380.6 mm versus the lowest estimated water yield of 193.9 mm for the Livneh-HG combination, in response to 23-year uncalibrated simulations. Higher annual ET and PET values were predicted with HG method versus the PM method for all three weather datasets in response to the uncalibrated simulations, due primarily to higher HG-based estimates during the growing season. Based on these results, it was found that the HG method is the preferred PET option for the UMRB. Initial calibration of SWAT was performed using the Livneh data and HGmethod for threeMississippi Rivermain stem gauge sites, which was followed by spatial validation at 10 other gauge sites located within the UMRB streamnetwork. Overall satisfactory results were found for the calibration and validation gauge sites, with the majority of R2 values ranging between 0.61 and 0.82, Nash- Sutcliffe modeling efficiency (NSE) values ranging between 0.50 and 0.79, and Kling-Gupta efficiency (KGE) values ranging between 0.61 and 0.84. The results of an additional experimental suite of six scenarios,which represented different combinations of climate data sets and calibrated parameters, revealed that suggested statistical criteria were again satisfied by the different scenario combinations. Overall, the PRISM data exhibited the strongest reliability for the UMRB.
We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.
Despite the clear ecological significance of the microbiomes inhabiting groundwater and connected ecosystems, our current understanding of their habitats, functionality, and the ecological processes controlling their assembly have been limited. In this study, an efficient pipeline combining geochemistry, high-throughput Fluidigm TM functional gene amplification and sequencing was developed to analyze the suspended and attached microbial communities inhabiting five groundwater monitoring wells in the Illinois Basin, USA. The dominant taxa in the suspended and the attached microbial communities exhibited significantly different spatial and temporal changes in both alpha- and beta-diversity. Further analyses of representative functional genes affiliated with N 2 fixation (nifH), methane oxidation (pmoA), and sulfate reduction (dsrB, and aprA), suggested functional redundancy within the shallow aquifer microbiomes. While more diversified functional gene taxa were observed for the suspended microbial communities than the attached ones except for pmoA, different levels of changes over time and space were observed between these functional genes. Notably, deterministic and stochastic ecological processes shaped the assembly of microbial communities and functional gene reservoirs differently. While homogenous selection was the prevailing process controlling assembly of microbial communities, the neutral processes (e.g., dispersal limitation, drift and others) were more important for the functional genes. The results suggest complex and changing shallow aquifer microbiomes, whose functionality and assembly vary even between the spatially proximate habitats and fractions. As a result, this research underscored the importance to include all the interface components for a more holistic understanding of the biogeochemical processes in aquifer ecosystems, which is also instructive for practical applications.
We report aircraft observations of extreme levels of HCl and the dihalogens Cl 2 , Br 2 , and BrCl in an industrial plume near the Great Salt Lake, Utah. Complete depletion of O 3 was observed concurrently with halogen enhancements as a direct result of photochemically produced halogen radicals. Observed fluxes for Cl 2 , HCl, and NO x agreed with facility-reported emissions inventories. Bromine emissions are not required to be reported in the inventory, but are estimated as 173 Mg year –1 Br 2 and 949 Mg year –1 BrCl, representing a major uncounted oxidant source. A zero-dimensional photochemical box model reproduced the observed O 3 depletions and demonstrated that bromine radical cycling was principally responsible for the rapid O 3 depletion. Inclusion of observed halogen emissions in both the box model and a 3D chemical model showed significant increases in oxidants and particulate matter (PM 2.5 ) in the populated regions of the Great Salt Lake Basin, where winter PM 2.5 is among the most severe air quality issues in the U.S. The model shows regional PM 2.5 increases of 10%–25% attributable to this single industrial halogen source, demonstrating the impact of underreported industrial bromine emissions on oxidation sources and air quality within a major urban area of the western U.S.
Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.