Large-eddy simulations were employed over half-ice and half-water surfaces, with varying surface temperatures, wind speeds, directions, as to test if the atmospheric interaction with the heterogeneous surface can be predicted via a heterogeneity Richardson number. This dataset was used to determine that surface heat fluxes over ice, water, and the aggregate surface seem to be captured reasonably well by the wind direction and the heterogeneity Richardson number, but the mean wind and turbulent kinetic energy (TKE) profiles were not, suggesting that not only the difference in stability between the two surface, but also the individual stabilities over each surface influence the dynamics.
Large-eddy simulations were employed over five different sea ice patterns, with a constant ice fraction, to test if the overlying atmospheric boundary layer (ABL) dynamics and thermodynamics differs. The results of these simulations were used to determine that there were differences in vertical heat flux, momentum flux, and horizontal wind speed, and that more surface information is needed to predict the ABL over the sea ice surface. To see what other surface information is needed, twenty-two landscape metrics were calculated over forty-four different maps at differing resolutions, using the FRAGSTATs program. The results of that analysis are available in a .csv file in this dataset.
Guo, Xuehui; Pan, Da; Daly, Ryan; Chen, Xi; Walker, John; Tao, Lei; McSpiritt, James; Zondlo, Mark
Gas-phase ammonia (NH3), emitted primarily from agriculture, contributes significantly to reactive nitrogen (Nr) deposition. Excess deposition of Nr to the environment causes acidification, eutrophication, and loss of biodiversity. The exchange of NH3 between land and atmosphere is bidirectional and can be highly heterogenous when underlying vegetation and soil characteristics differ. Direct measurements that assess the spatial heterogeneity of NH3 fluxes are lacking. To this end, we developed and deployed two fast-response, quantum cascade laser-based open-path NH3 sensors to quantify NH3 fluxes at a deciduous forest and an adjacent grassland separated by 700 m in North Carolina, United States from August to November, 2017. The sensors achieved 10 Hz precisions of 0.17 ppbv and 0.23 ppbv in the field, respectively. Eddy covariance calculations showed net deposition of NH3 (-7.3 ng NH3-N m−2 s−1) to the forest canopy and emission (3.2 ng NH3-N m−2 s−1) from the grassland. NH3 fluxes at both locations displayed diurnal patterns with absolute magnitudes largest midday and with smaller peaks in the afternoons. Concurrent biogeochemistry data showed over an order of magnitude higher NH3 emission potentials from green vegetation at the grassland compared to the forest, suggesting a possible explanation for the observed flux differences. Back trajectories originating from the site identified the upwind urban area as the main source region of NH3. Our work highlights the fact that adjacent natural ecosystems sharing the same airshed but different vegetation and biogeochemical conditions may differ remarkably in NH3 exchange. Such heterogeneities should be considered when upscaling point measurements, downscaling modeled fluxes, and evaluating Nr deposition for different natural land use types in the same landscape. Additional in-situ flux measurements accompanied by comprehensive biogeochemical and micrometeorological records over longer periods are needed to fully characterize the temporal variabilities and trends of NH3 fluxes and identify the underlying driving factors.
This dataset contains supplementary materials for Chapter 4 and Chapter 5 of Yiheng Tao's PhD dissertation (2022). The dissertation’s abstract is provided here:
Carbon capture, utilization, and storage (CCUS) mitigates climate change by capturing carbon dioxide (CO2) emissions from large point sources, or CO2 from the ambient air, and subsequently reusing the captured CO2 or injecting it into deep geological formations for long-term and secure storage. Almost all current decarbonization pathways include large-scale CCUS, on the order of a billion tonnes (Gt) of CO2 captured and stored each year globally starting in 2030, yet the actual deployment has lagged far behind (around 0.04 Gt CO2 was captured in 2021). In this dissertation, I contribute to several aspects of largescale deployment of CCUS by (1) developing and applying efficient numerical models to simulate geological CO2 storage and (2) identifying key policies to address the bottlenecks of overall CCUS deployment. This dissertation concerns the United States, China, and the Belt and Road Initiative (BRI) region through research projects that are consistent with each location’s current development stage of CCUS.
Chapters 2 and 3 contain computational modeling studies. In Chapter 2, I develop a new series of vertical-equilibrium (VE) models in the dual-continuum modeling framework to simulate CO2 injection and migration in fractured geological formations. Those models are shown to be effective and efficient when properties of the formation allow for the VE assumption. In Chapter 3, I apply a VE model to simulate basin-scale CO2 injection in the Junggar Basin of Northwestern China. The results show that current regional emissions of more than 100 million tonnes of CO2 per year can be stored effectively, thereby confirming the great potential of the Junggar Basin for early CCUS deployment.
Chapters 4 and 5 contain policy analyses. In Chapter 4, I propose a dynamic system consisting of new CO2 pipelines and novel Allam-cycle power plants in the Central United States, and examine how government policies, including an extended Section 45Q tax credit, may improve the economic feasibility of this system. Lastly, in Chapter 5, I investigate and quantify CO2 emissions implications of power plant projects associated with the BRI. I also propose a “greenness ratio” to measure the level of environmental sustainability of BRI in the power sector.
Pan, Da; Gelfand, Ilya; Tao, Lei; Abraha, Michael; Sun, Kang; Guo, Xuehui; Chen, Jiquan; Robertson, G. Philip; Zondlo, Mark A.
This dataset contains spectroscopic simulations, experimental results for the 2202 cm-1 N2O absorption line, and N2O flux measurements shown in "A New Open-path Eddy Covariance Method for N2O and Other Trace Gases that Minimizes Temperature Corrections" by Da Pan, Ilya Gelfand, Lei Tao, Michael Abraha, Kang Sun, Xuehui Guo, Jiquan Chen, G. Philip Robertson, and Mark A. Zondlo. The HITRAN Application Programming Interface (HAPI) with HITRAN 2016 was used for spectroscopic simulations. Experiments were conducted to quantify H2O-broadened half-width at half maximum and validate spectroscopic simulations. N2O flux was measured with both eddy covariance and static chamber methods.
Elevated reactive nitrogen (Nr) deposition is a concern for alpine ecosystems, and dry NH3 deposition is a key contributor. Understanding how emission hotspots impact downwind ecosystems through dry NH3 deposition provides opportunities for effective mitigation. However, direct NH3 flux measurements with sufficient temporal resolution to quantify such events are rare. Here, we measured NH3 fluxes at Rocky Mountain National Park (RMNP) during two summers and analyzed transport events from upwind agricultural and urban sources in northeastern Colorado. We deployed open-path NH3 sensors on a mobile laboratory and an eddy covariance tower to measure NH3 concentrations and fluxes. Our spatial sampling illustrated an upslope event that transported NH3 emissions from the hotspot to RMNP. Observed NH3 deposition was significantly higher when backtrajectories passed through only the agricultural region (7.9 ng m-2 s-1) versus only the urban area (1.0 ng m-2 s-1) and both urban and agricultural areas (2.7 ng m-2 s-1). Cumulative NH3 fluxes were calculated using observed, bidirectional modeled, and gap-filled fluxes. More than 40% of the total dry NH3 deposition occurred when air masses were traced back to agricultural source regions. More generally, we identified that 10 (25) more national parks in the U.S. are within 100 (200) km of an NH3 hotspot, and more observations are needed to quantify the impacts of these hotspots on dry NH3 depositions in these regions.
Wang, Rui; Guo, Xuehui; Pan, Da; Kelly, James; Bash, Jesse; Sun, Kang; Paulot, Fabien; Clarisse, Lieven; Van Damme, Martin; Whitburn, Simon; Coheur, Pierre-François; Clerbaux, Cathy; Zondlo, Mark
Monthly, high resolution (~2 km) ammonia (NH3) column maps from the Infrared Atmospheric Sounding Interferometer (IASI) were developed across the contiguous United States and adjacent areas. Ammonia hotspots (95th percentile of the column distribution) were highly localized with a characteristic length scale of 12 km and median area of 152 km2. Five seasonality classes were identified with k-means++ clustering. The Midwest and eastern United States had a broad, spring maximum of NH3 (67% of hotspots in this cluster). The western United States, in contrast, showed a narrower mid-summer peak (32% of hotspots). IASI spatiotemporal clustering was consistent with those from the Ammonia Monitoring Network. CMAQ and GFDL-AM3 modeled NH3 columns have some success replicating the seasonal patterns but did not capture the regional differences. The high spatial-resolution monthly NH3 maps serve as a constraint for model simulations and as a guide for the placement of future, ground-based network sites.
Natural gas vehicles (NGVs) have been promoted in China to mitigate air pollution, yet our measurements and analyses show that NGV growth in China may have significant negative impacts on climate change. We conducted real-world vehicle emission measurements in China and found high methane emissions from heavy-duty NGVs (90% higher than current emission limits). These emissions have been ignored in previous emission estimates, leading to biased results. Applying our observations to life-cycle analyses, we found that switching to NGVs from conventional vehicles in China has led to a net increase in greenhouse gas (GHG) emissions since 2000. With scenario analyses, we also show that the next decade will be critical for China to reverse the trend with the upcoming China VI standard for heavy-duty vehicles. Implementing and enforcing the China VI standard is challenging, and the method demonstrated here can provide critical information regarding the fleet-level CH4 emissions from NGVs.
Yang, Yuan; Pan, Ming; Beck, Hylke; Fisher, Colby; Beighley, R. Edward; Kao, Shih-Chieh; Hong, Yang; Wood, Eric
Conventional basin-by-basin approaches to calibrate hydrologic models are limited to gauged basins and typically result in spatially discontinuous parameter fields. Moreover, the consequent low calibration density in space falls seriously behind the need from present-day applications like high resolution river hydrodynamic modeling. In this study we calibrated three key parameters of the Variable Infiltration Capacity (VIC) model at every 1/8° grid-cell using machine learning-based maps of four streamflow characteristics for the conterminous United States (CONUS), with a total of 52,663 grid-cells. This new calibration approach, as an alternative to parameter regionalization, applied to ungauged regions too. A key difference made here is that we tried to regionalize physical variables (streamflow characteristics) instead of model parameters whose behavior may often be less well understood. The resulting parameter fields no longer presented any spatial discontinuities and the patterns corresponded well with climate characteristics, such as aridity and runoff ratio. The calibrated parameters were evaluated against observed streamflow from 704/648 (calibration/validation period) small-to-medium-sized catchments used to derive the streamflow characteristics, 3941/3809 (calibration/validation period) small-to-medium-sized catchments not used to derive the streamflow characteristics) as well as five large basins. Comparisons indicated marked improvements in bias and Nash-Sutcliffe efficiency. Model performance was still poor in arid and semiarid regions, which is mostly due to both model structural and forcing deficiencies. Although the performance gain was limited by the relative small number of parameters to calibrate, the study and results here served as a proof-of-concept for a new promising approach for fine-scale hydrologic model calibrations.
Since 1850 the concentration of atmospheric methane (CH4), a potent greenhouse gas, has more than doubled. Recent studies suggest that emission inventories may be missing sources and underestimating emissions. To investigate whether offshore oil and gas platforms leak CH4 during normal operation, we measured CH4 mole fractions around eight oil and gas production platforms in the North Sea which were neither flaring gas nor off-loading oil. We use the measurements from summer 2017, along with meteorological data, in a Gaussian plume model to estimate CH4 emissions from each platform. We find CH4 mole fractions of between 11 and 370 ppb above background concentrations downwind of the platforms measured, corresponding to a median CH4 emission of 6.8 g CH4 s-1 for each platform, with a range of 2.9 to 22.3 g CH4 s-1. When matched to production records, during our measurements individual platforms lost between 0.04% and 1.4% of gas produced with a median loss of 0.23%. When the measured platforms are considered collectively, (i.e. the sum of platforms’ emission fluxes weighted by the sum of the platforms’ production), we estimate the CH4 loss to be 0.19% of gas production. These estimates are substantially higher than the emissions most recently reported to the National Atmospheric Emission Inventory (NAEI) for total CH4 loss from United Kingdom platforms in the North Sea. The NAEI reports CH4 losses from the offshore oil and gas platforms we measured to be 0.13% of gas production, with most of their emissions coming from gas flaring and offshore oil loading, neither of which were taking place at the time of our measurements. All oil and gas platforms we observed were found to leak CH4 during normal operation and much of this leakage has not been included in UK emission inventories. Further research is required to accurately determine total CH4 leakage from all offshore oil and gas operations and to properly include the leakage in national and international emission inventories.
Geochemical and geomechanical perturbations of the subsurface caused by the injection of fluids present the risk of leakage and seismicity. This study investigated how flow of acidic fluids affects hydraulic and frictional properties of fractures using experiments with 3.8 cm-long specimens of Eagle Ford shale, a laminated shale with carbonate-rich strata. In low-pressure flow cells, one set of samples was exposed to an acidic brine and another set was exposed to a neutral brine. X-ray computed tomography and x-ray fluorescence analysis revealed that samples exposed to the acidic brine were calcite-depleted and had developed a porous altered layer, while the other set showed little evidence of alteration. After reaction, samples were compacted and sheared in a triaxial cell that supplied normal stress and differential pore pressure at prescribed sliding velocities, independently measuring friction and permeability. During the initial compaction, the porous altered layer collapsed into fine particles that filled the fracture aperture. This effectively impeded flow and sealed the fracture, resulting in a decrease in fracture permeability by 1 to 2 orders of magnitude relative to the compressed unaltered fractures. During shear, the collapsed layer of fine-grained particles prevented the formation of interlocking micro-asperities resulting in lower frictional strength. With regard to subsurface risks, this study showcases how coupled geochemical and geomechanical processes could favorably seal fractures to inhibit leakage, but also could increase the likelihood of induced seismicity. These findings have important implications for geological carbon sequestration, pressurized fluid energy storage, geothermal energy, and other subsurface technologies.
Fractures in geological formations may enable migration of environmentally relevant fluids, as in leakage of CO2 through caprocks in geologic carbon sequestration. We investigated geochemically induced alterations of fracture geometry in Indiana Limestone specimens. Experiments were the first of their kind, with periodic high-resolution imaging using X-ray computed tomography (xCT) scanning while maintaining high pore pressure (100 bar). We studied two CO2-acidified brines having the same pH (3.3) and comparable thermodynamic
disequilibrium but different equilibrated pressures of CO2 (PCO2 values of 12 and 77 bar). High-PCO2 brine has a faster calcite dissolution kinetic rate because of the accelerating effect of carbonic acid. Contrary to expectations, dissolution extents were comparable in the two experiments. However, progressive xCT
images revealed extensive channelization for high PCO2, explained by strong positive feedback between ongoing flow and reaction. The pronounced channel increasingly directed flow to a small region of the fracture, which explains why the overall dissolution was lower than expected. Despite this, flow simulations revealed large increases in permeability in the high-PCO2 experiment. This study shows that the permeability evolution of dissolving fractures will be larger for faster-reacting fluids. The overall mechanism is not because more rock dissolves, as would be commonly assumed, but because of accelerated fracture channelization.
This three-year project, performed by Princeton University in partnership with the University of Minnesota and Brookhaven National Laboratory, examined geologic carbon sequestration in regard to CO2 leakage and potential subsurface liabilities. The research resulted in basin-scale analyses of CO2 and brine leakage in light of uncertainties in the characteristics of leakage processes, and generated frameworks to monetize the risks of leakage interference with competing subsurface resources. The geographic focus was the Michigan sedimentary basin, for which a 3D topographical model was constructed to represent the hydrostratigraphy. Specifically for Ottawa County, a statistical analysis of the hydraulic properties of underlying sedimentary formations was conducted. For plausible scenarios of injection into the Mt. Simon sandstone, leakage rates were estimated and fluxes into shallow drinking-water aquifers were found to be less than natural analogs of CO2 fluxes. We developed the Leakage Impact Valuation (LIV) model in which we identified stakeholders and estimated costs associated with leakage events. It was found that costs could be incurred even in the absence of legal action or other subsurface interference because there are substantial costs of finding and fixing the leak and from injection interruption. We developed a model framework called RISCS, which can be used to predict monetized risk of interference with subsurface resources by combining basin-scale leakage predictions with the LIV method. The project has also developed a cost calculator called the Economic and Policy Drivers Module (EPDM), which comprehensively calculates the costs of carbon sequestration and leakage, and can be used to examine major drivers for subsurface leakage liabilities in relation to specific injection scenarios and leakage events. Finally, we examined the competitiveness of CCS in the energy market. This analysis, though qualitative, shows that financial incentives, such as a carbon tax, are needed for coal combustion with CCS to gain market share. In another part of the project we studied the role of geochemical reactions in affecting the probability of CO2 leakage. A basin-scale simulation tool was modified to account for changes in leakage rates due to permeability alterations, based on simplified mathematical rules for the important geochemical reactions between acidified brines and caprock minerals. In studies of reactive flows in fractured caprocks, we examined the potential for permeability increases, and the extent to which existing reactive transport models would or would not be able to predict it. Using caprock specimens from the Eau Claire and Amherstburg, we found that substantial increases in permeability are possible for caprocks that have significant carbonate content, but minimal alteration is expected otherwise. We also found that while the permeability increase may be substantial, it is much less than what would be predicted from hydrodynamic models based on mechanical aperture alone because the roughness that is generated tends to inhibit flow.
Complete dataset of pore water chemical parameters measured at the Marsh Resource Meadowlands Mitigation Bank, a tidal marsh within the New Jersey Meadowlands, from March 2011 to April 2012. Analytes measured include dissolved methane, sulfate, dissolved organic carbon, temperature, salinity, and pH. Measurements were conducted using porewater dialysis samplers, and water was sampled from the surface to a depth of 60 cm.