Kira Shonkwiler

Dr. Kira Shonkwiler


Research Focus & Metrics

Academic Career Length

2Years (2025-Present)

(Calculated once Scilit publication data is synced: current year − year of first publication + 1)

Representative Publications

View All Publications
2 July 2025 in Atmosphere
Verified
Journal Article
Open access
PDF file
AAAbhinav Anand
Abhinav Anand
SRStuart N. Riddick
Stuart N. Riddick
KSKira Shonkwiler
Kira Shonkwiler
AUAashish Upreti
Aashish Upreti
MMMichael Moy
Michael Moy
DZDaniel Zimmerle
Daniel Zimmerle
Recent studies estimate that emissions from oil and gas production facilities contribute between 20 and 50% of the total methane (CH4) emitted in the US; therefore, quantifying and reducing these emissions are crucial for achieving climate goals. Methane quantification depends on both measuring methane concentrations and converting them to emissions through a modeling framework. Currently, simple atmospheric dispersion models are primarily used to quantify emissions and concentrations, but these estimates are highly uncertain when quantifying emissions from complex aerodynamic sources, such as oil and gas facilities. This investigation used a CFD modeling approach, which can account for aerodynamic complexity but has hitherto not been used to model methane concentrations downwind of a methane release of a known rate, and compared it against in situ measurements. High-time-resolution (1 Hz) methane concentration and meteorological data were measured during experiments conducted at the METEC on 21 March and 11 July 2024. The METEC site configuration, measured wind data, and controlled emission rates were used as input for the CONVERGE CFD model to model downwind CH4 concentration. The modeling was carried out between 20 and 70 m, from two different points of release in two separate controlled-release experiments, one from a separator and another from a wellhead. In these experiments, we found that the CFD model could predict the CH4 concentrations downwind of the release to a good degree. The model was evaluated on multiple metrics to assess its performance in estimating methane concentrations at typical fence line distances (∼30 m). These results help us to understand external flows and the ability of CFD models to predict downwind concentrations in aerodynamically complex environments.

29 January 2026 in Preprints.org
Verified
Preprint
Open access
PDF file
JMJacob Mdigo
Jacob Mdigo
ASArthur Santos
Arthur Santos
GDGerald Duggan
Gerald Duggan
PVPrajay Vora
Prajay Vora
KSKira Shonkwiler
Kira Shonkwiler
DZDaniel Zimmerle
Daniel Zimmerle
This work presents a mechanistic modeling approach for simulating methane emissions from tri-ethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (< 0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

20 April 2026 in Atmosphere
Verified
Journal Article
Open access
PDF file
AUAashish Upreti
Aashish Upreti
KSKira Shonkwiler
Kira Shonkwiler
SRStuart N. Riddick
Stuart N. Riddick
DZDaniel Zimmerle
Daniel Zimmerle
Methane (CH4) is a major component of natural gas and a potent greenhouse gas. Increasing atmospheric methane concentrations are attributed to emissive anthropogenic activities by an average of 13 ppb per yr since 2020 and are linked to a changing global climate. Mitigating CH4 emissions from oil and gas production sites has recently become a target to reduce overall greenhouse gas emissions; however, monitoring the efficacy of mitigation strategies depends on accurate quantification of CH4 emissions at the facility-level. Near-field quantification of methane (CH4) emissions from oil and gas (O&G) facilities remains challenging due to the effects of atmospheric variability and sensor configuration on atmospheric dispersion models. This study evaluates the performance of two atmospheric dispersion models, the Gaussian plume (GP) and backward Lagrangian stochastic (bLS), by comparing calculated CH4 emissions to controlled single-point emissions between 0.4 and 5.2 kg CH4 h−1. Emissions were calculated by both models using 121 individual sets of measurements comprising five-minute averaged downwind methane mixing ratios and matching meteorological data. The comparison shows that the bLS approach achieved a higher proportion of emission estimates within a factor of two (FAC2) of the known emission rates compared to the GP approach. The emissions calculated by the bLS model also had a lower multiplicative error and reduced bias relative to GP. Other error-based metrics further confirmed the bLS model performed better, as it yielded lower RMSE and MAE than GP. Statistical analysis of the emission data shows that the lateral and vertical alignment of the source and the sensor plays a critical role in emission estimations, as measurements made closer to the plume centerline and at a distance between 40 and 80 m downwind yielded the best FAC2 agreement. High wind meander degraded the ability of both approaches to generate representative emissions, particularly with the GP approach, as it violates the modeling approach’s assumption of steady-state emissions. Data suggest emissions calculated by the bLS model are comprehensively in better agreement, but the computational demands of the modeling approach and integration into fenceline systems limit real-time applicability. While these results provide insight into model performance under controlled near-field conditions, their applicability to more complex or heterogeneous oil and gas production environments (e.g., the regions Marcellus or Unita Basins) remains limited and uncertain.

7 April 2026 in Fuels
Verified
Journal Article
Open access
PDF file
JMJacob Mdigo
Jacob Mdigo
ASArthur Santos
Arthur Santos
GDGerald Duggan
Gerald Duggan
PVPrajay Vora
Prajay Vora
KSKira Shonkwiler
Kira Shonkwiler
DZDaniel Zimmerle
Daniel Zimmerle
This work presents a mechanistic modeling approach for simulating methane emissions from triethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (<0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

4 March 2026 in Preprints.org
Verified
Preprint
Open access
PDF file
AUAashish Upreti
Aashish Upreti
KSKira Shonkwiler
Kira Shonkwiler
SRStuart N. Riddick
Stuart N. Riddick
DZDaniel Zimmerle
Daniel Zimmerle
Methane (CH4) is a major component of natural gas and a potent greenhouse gas. Increasing atmospheric methane concentrations are attributed to emissive anthropogenic activities by an average of 13 ppb per yr since 2020 and are linked to a changing global climate. Mitigating CH4 emissions from oil and gas production sites has recently become a target to reduce overall greenhouse gas emissions, however, monitoring the efficacy of mitigation strategies depends on accurate quantification of CH4 emissions at the facility-level. Near-field quantification of methane (CH4) emissions from oil and gas (O&G) facilities remains challenging due to the effects of atmospheric variability and sensor configuration on atmospheric dispersion models. This study evaluates the performance of two atmospheric dispersion models, the Gaussian Plume (GP) and backward Lagrangian Stochastic (bLS), by comparing calculated CH4 emissions to controlled single-point emissions of between 0.4 and 5.2 kg CH4 h-1. Emissions were calculated by both models using 121 individual sets of measurements comprising five-minute averaged downwind methane mixing ratios and matching meteorological data. Comparison shows the bLS approach showed better predictive performance with twice as many emission estimates were within a factor of two (FAC2) of the known emission rates compared to those calculated using the GP approach. The emissions calculated by the bLS model also had a lower multiplicative error and reduced bias relative to GP. Other error-based metrics further confirmed the bLS model performed better, as it yielded lower RMSE and MAE than GP. Statistical analysis of the emission data shows the lateral and vertical alignment of source and sensor plays a critical role in emission estimations as measurements made closer to the plume centerline and at a distance between 40 to 80 m downwind yielded the best FAC2 agreement. High wind meander degraded ability of both approaches to generate representative emissions particularly with the GP approach as it violates the modelling approach’s assumption of steady-state emissions. Data suggest emissions calculated by the bLS model are comprehensively in better agreement but the computational demands of the modeling approach and integration into fenceline systems limit real-time applicability. While it is likely that the results presented here are suitable for informing leak detection technology in relatively flat unvegetated environments, it is currently unknown if these findings will be applicable in more vertiginous or heavily vegetated oil and gas producing regions of the Marcellus or Uinta Basins.