Volcanic rifting environments, such as in Iceland, are challenging for conventional seismic reflection methods using active surface seismic sources. This study demonstrates the potential of a novel technique, called Virtual Reflection Seismic Profiling (VRSP) for imaging reflections in geothermal regions, like Hverahlíð, where a dense seismic array recorded a number of local microearthquakes for cross-correlation. Uppsala University, in collaboration with Reykjavik Energy, recorded seismicity in Hverahlíð using both seismometers and geophones. Acquired data were processed using the VRSP method, which applies seismic interferometry only to selected events, in this thesis local microearthquakes. Cross-correlation of the signal from a microearthquake recorded at one of the stations, which acts as a virtual source, with a ghost reflection recorded by the remaining receivers, produces a virtual shot gather. Stacking each station’s result, for all available events, and following a conventional multichannel processing sequence resulted in two stacked seismic images. Potential reflections observed in the obtained sections could be associated with major feed zones identified in the area by the borehole measurements. Eight dynamite explosions were processed with a conventional seismic reflection method, as a complementary source to the microearthquakes. In the produced stacked seismic section two potential reflections could be observed. Results from both passive and active datasets were 3D visualised to verify whether the reflections correspond to each other between sections. Two horizons were traced throughout all three stacked sections. One more interface appeared on two images obtained from processing the passive data. This study shows promising results for using natural sources to image the subsurface in a challenging environment.