Systems and methods are disclosed to construct subsurface images from diffuse seismic energy. Various disclosed system embodiments include multiple seismic sensors that each convert received seismic energy into one or more seismic signals. One or more processor combine the seismic signals to determine a subsurface map. As part of determining the map, the processor(s) systematically focus the array of seismic sensors on each bin in the subsurface volume of interest. In this manner each bin becomes a focal point of the array. For each bin, the processor(s) analyze the seismic wave travel time to each seismic sensor and apply a corresponding time shift to align the seismic signals with a uniform travel time. The time-shifted seismic signals are then combined to determine an intensity value for seismic energy radiating from the focal point. A subsurface map can then be derived from the intensity value as a function of position.