A method of determining characteristics (310, 512, 570, 572, 574) of a subterranean body in reservoir characterization, comprising: performing a base seismic survey operation (302, 550); after performing the base seismic survey operation (302, 550), performing pressure testing (304, 552) in a well, wherein the pressure testing (304, 552) comprises drawing down pressure in the well and measuring pressure data in the well during the pressure testing (304, 552); performing a second seismic survey operation (306, 556) coincidentally with the pressure testing (304, 552), wherein seismic data obtained by the second seismic survey operation (306, 556) are affected by pressure changes in the subterranean body due to the pressure testing (304, 552); comparing seismic data of the base seismic survey operation (302, 550) with the seismic data of the second seismic survey operation (308, 556); providing the pressure data and the result of comparing of the seismic data to a processing system for processing to determine the characteristics (310, 512, 570, 572, 574) of the subterranean body; wherein the seismic data obtained from the base and second seismic survey operations make up time-lapse seismic data; characterized in that the method further comprises: processing the time-lapse seismic data to detect pressure changes; and determining a four-Dimensional spatio-temporal evolution of pressure in the subterranean body based on the measured pressure data and the processing of the time-lapse seismic data to capture movement of pressure fronts as a function of both time and space. performing a base seismic survey operation (302, 550); after performing the base seismic survey operation (302, 550), performing pressure testing (304, 552) in a well, wherein the pressure testing (304, 552) comprises drawing down pressure in the well and measuring pressure data in the well during the pressure testing (304, 552); performing a second seismic survey operation (306, 556) coincidentally with the pressure testing (304, 552), wherein seismic data obtained by the second seismic survey operation (306, 556) are affected by pressure changes in the subterranean body due to the pressure testing (304, 552); comparing seismic data of the base seismic survey operation (302, 550) with the seismic data of the second seismic survey operation (308, 556); providing the pressure data and the result of comparing of the seismic data to a processing system for processing to determine the characteristics (310, 512, 570, 572, 574) of the subterranean body; wherein the seismic data obtained from the base and second seismic survey operations make up time-lapse seismic data; characterized in that the method further comprises: processing the time-lapse seismic data to detect pressure changes; and determining a four-Dimensional spatio-temporal evolution of pressure in the subterranean body based on the measured pressure data and the processing of the time-lapse seismic data to capture movement of pressure fronts as a function of both time and space. processing the time-lapse seismic data to detect pressure changes; and determining a four-Dimensional spatio-temporal evolution of pressure in the subterranean body based on the measured pressure data and the processing of the time-lapse seismic data to capture movement of pressure fronts as a function of both time and space.