Referring now to FIG. 8, a process flow diagram of a subroutine 800 for sound propagation modeling is shown. In accordance with certain aspects of the present disclosure, subroutine 800 may be implemented or otherwise embodied as a component or subcomponent of a spatial audio processing system; for example, spatial audio processing system 100 as shown and described in FIG. 1. In certain embodiments, subroutine 800 may be a subroutine of routine 600 and/or may comprise one or more sequential or successive steps of routine 600 (as shown and described in FIG. 6). In accordance with an embodiment, subroutine 800 may be initiated by receiving an audio input comprising m-Channels of modeling segment audio 802. The m-Channels are associated with one or more transducers (e.g., microphones) being located within an acoustic space or environment. The one or more transducers may be operably interfaced to comprise an array. In certain specific embodiments, a spatial audio processing system may comprise four or more audio input channels. Subroutine 800 may continue by applying a Fourier Transform to the modeling segment audio, in frames, to convert the modeling segment audio from the time domain to the frequency domain 804. As in routine 600, the Fourier Transform in subroutine 800 may be selected from one or more alternative transform functions, such as Fast Fourier transform, Short Time Fourier transform and/or other window functions or overlap. Subroutine 800 may continue by executed one or more substeps 806, 808, and 810. In certain embodiments, subroutine 800 may proceed by summing (on a per frame basis) the magnitudes of each binary file, or BIN, for each channel of audio 806. The magnitudes of each frame may be sorted in rank order, per BIN 808. Subroutine 800 may apply a magnitude threshold test on the sorted BINs to generate a mask configured to filter silence and stray noise components from the m-Channels of modeling segment audio 810. It is anticipated that alternative techniques to the magnitude threshold test may be employed to generate a temporal and/or spectral mask in substep 810. In certain embodiments, subroutine 800 may continue by applying the mask to the modeling audio segment to obtain only time-frequency BINs containing the source signal 812. Subroutine 800 may continue by calculating the cross power spectral density (CPSD) of the masked modeling audio segment for each BIN, for each of the m-Channels of audio 814. Subroutine 800 may continue by normalizing the CPSD to obtain a frequency domain Green's Function for each BIN 816 to identify an audio propagation model originating from a three-dimensional point source within the audio environment/location. In certain embodiments, the Green's Function data may be continuously updated/refined in response to changing conditions/variables, including tracking a target sound source as it moves to one or more new/different locations within the audio environment/location. Subroutine 800 may conclude by storing/exporting the Green's Function for the point source location within the audio environment 818.