“…These include: (i) how best to address, in a numerically efficient and physically realistic manner, the handling of layer debonding and fluid invasion along layer interfaces with associated stunting of fracture height growth in shallower wells-relatively limited progress has been made in this area [94][95][96]; (ii) how to appropriately adjust current (linear elastic) simulators to enable modeling of the propagation of hydraulic fractures in highly cleated coal bed seams (for the extraction of methane) [97]; (iii) how to appropriately adjust current (linear elastic) simulators to enable modeling of the propagation of hydraulic fractures in weakly consolidated and unconsolidated ''soft'' sandstones, such as are found in the Gulf of Mexico-limited progress has been made in this area [98,99]; (iv) laboratory and field observations demonstrate that mode III fracture growth does occur [100], and this needs to be further researched; (v) related to (iii), the effect of the invaded zone ahead of the fracture tip needs to be further researched-a criterion to switch from a fluid lag based approach to an invaded zone based approach in a numerical model is required; (vi) suitable models for the propagation of hydraulic fractures in naturally fractured reservoirs that result in complex (non-planar) geometric configurations requires development [101]; and (vii) how to efficiently model 3D or ''out of plane'' effects, such as fracture re-alignment (when the fracture initiates following an orientation that is not perpendicular to the minimum in situ stress and then tries to re-align itself), which could be a cause of near-wellbore tortuosity or even ''pinching,'' a factor that usually determines the success or failure of hydraulic fracturing treatments [102].…”