From 2009 to 2017, parts of Central America experienced marked increase in the number of small to moderate-sized earthquakes. For example, three significant earthquakes (∼MW5) occurred near Prague, Oklahoma, the United States in 2011. On 6 Nov 2011, an Mw 5.7 earthquake occurred in Prague, central Oklahoma with a sequence of aftershocks. The seismic activity has been attributed to slip on the Wilzetta fault system. In this study, we provide a 3D fully coupled poroelastic analysis (using FLAC3D) of the Wilzetta fault system response to saltwater injection in the underpressured subsurface layers, especially Arbuckle group and the basement, to evaluate the conditions that might have led to the increased seismicity. In view of the data-limited nature of the problem, we have considered multiple plausible scenarios, and use the available data to evaluate the hydromechanical response of the faults of interest in the study area. Numerical simulations show that the injection of large volumes of fluid into Arbuckle group tends to bring the part of the Wilzetta faults in Arbuckle group and basement into near critical conditions.