A computer-implemented method, comprising: receiving faults geometry data of a subsurface earth volume; receiving at least one data set associated with the subsurface earth volume; simulating a plurality of models for a subsurface earth volume based on the received faults geometry data and the received dataset, the plurality of models comprising: one or more linearly independent far field stress models, one or more discontinuity pressure models including the pressure inside a 3D discontinuity and one or more discontinuity pressure shift models including a salt-pressure-shift parameter (Sps) corresponding to an intersection of a gradient pressure linear function as a function of depth from the Earth's surface, wherein pressure inside discontinuity at depth z is given by P Z = Sps + P r | z |; using the solutions of six simulation results, to invert for pressure gradient inside 3D discontinuities, the gradient regional stress regime and the salt-pressure-shift parameter, the inversion comprising a transformation matrix: I 6 α = σ xx σ xy σ yy σ zz Sps Pr = σ R where I 6 represents the identity matrix of dimension six, Sps represents a salt-pressure-shift, Pr represents the pressure inside discontinuities, and α represents six linear coefficients representing parameters of the far field stress; predicting the stress attribute of the subsurface earth volume, based on the stress value, the discontinuity pressure shift at the surface of the Earth, the strain value, the displacement value, or the combination thereof computed through inversion of the models; and generating at least one parameter to configure a control signal based upon the predicted stress attribute, wherein the control signal is configured to cause equipment to perform a physical operation, and wherein the equipment comprises drilling equipment, exploration equipment, a well control injector, a valve, or a combination thereof. receiving faults geometry data of a subsurface earth volume; receiving at least one data set associated with the subsurface earth volume; simulating a plurality of models for a subsurface earth volume based on the received faults geometry data and the received dataset, the plurality of models comprising: one or more linearly independent far field stress models, one or more discontinuity pressure models including the pressure inside a 3D discontinuity and one or more discontinuity pressure shift models including a salt-pressure-shift parameter (Sps) corresponding to an intersection of a gradient pressure linear function as a function of depth from the Earth's surface, wherein pressure inside discontinuity at depth z is given by P Z = Sps + P r | z |; using the solutions of six simulation results, to invert for pressure gradient inside 3D discontinuities, the gradient regional stress regime and the salt-pressure-shift parameter, the inversion comprising a transformation matrix: I 6 α = σ xx σ xy σ yy σ zz Sps Pr = σ R where I 6 represents the identity matrix of dimension six, Sps represents a salt-pressure-shift, Pr represents the pressure inside discontinuities, and α represents six linear coefficients representing parameters of the far field stress; predicting the stress attribute of the subsurface earth volume, based on the stress value, the discontinuity pressure shift at the surface of the Earth, the strain value, the displacement value, or the combination thereof computed through inversion of the models; and generating at least one parameter to configure a control signal based upon the predicted stress attribute, wherein the control signal is configured to cause equipment to perform a physical operation, and wherein the equipment comprises drilling equipment, exploration equipment, a well control injector, a valve, or a combination thereof. one or more linearly independent far field stress models, one or more discontinuity pressure models including the pressure inside a 3D discontinuity and one or more discontinuity pressure shift models including a salt-pressure-shift parameter (Sps) corresponding to an intersection of a gradient pressure linear function as a function of depth from the Earth's surface, wherein pressure inside discontinuity at depth z is given by P Z = Sps + P r | z |; using the solutions of six simulation results, to invert for pressure gradient inside 3D discontinuities, the gradient regional stress regime and the salt-pressure-shift parameter, the inversion comprising a transformation matrix: I 6 α = σ xx σ xy σ yy σ zz Sps Pr = σ R where I 6 represents the identity matrix of dimension six, Sps represents a salt-pressure-shift, Pr represents the pressure inside discontinuities, and α represents six linear coefficients representing parameters of the far field stress; predicting the stress attribute of the subsurface earth volume, based on the stress value, the discontinuity pressure shift at the surface of the Earth, the strain value, the displacement value, or the combination thereof computed through inversion of the models; and generating at least one parameter to configure a control signal based upon the predicted stress attribute, wherein the control signal is configured to cause equipment to perform a physical operation, and wherein the equipment comprises drilling equipment, exploration equipment, a well control injector, a valve, or a combination thereof.