1. A method comprising the steps of: determining a grid of G tentative source locations r g to localize more than one acoustic sources, wherein a position of one of the acoustic sources at time t is denoted by {r k (t)ε d } k=1 K , where d=[2, 3]; using an acoustic model to compute replicas that correspond to simulated acoustic sources at locations in , wherein the replicas are model-predicted Fourier coefficient vectors at F frequencies {ω f } f=1 F corresponding to the acoustic pressure field as sampled by an acoustic sensor array having N sensors, where G>>max{KF, N} and K is a presumed number of acoustic sources; collecting, using the acoustic sensor array, time series data of actual acoustic measurements at each sensor of the acoustic sensor array caused by the acoustic sources; using a discrete-time Fourier transform (DFT) on the collected time series data to compute Fourier coefficients at frequencies {ω f } f=1 F for all N sensors; constructing DFT coefficient vectors y f (t):=[y f 1 (t),. .. , y f N (t)]′ε N using the Fourier coefficients previously obtained via the DFT, where y f n (t) denotes the Fourier-coefficient estimate corresponding to ω f at time t obtained from data gathered by the n-th sensor; modeling DFT coefficient vectors at ω f as y f (t)=Σ g=1 G s g,f (t)p g,f +ε f (t), f=1,. .. , F, where s g,f (t) denotes the Fourier coefficient at ω f of the spectrum corresponding to the acoustic signature of the acoustic source located at r g ε at time t, p g,f ε N is the replica for ω f corresponding to a source located at r g ε normalized so that ∥p g,f ∥ 2 =1, and ε f (t)ε N denotes a zero-mean additive noise component; obtaining, at time t, an estimate Ŝ(t) according to the equation arg ⁢ ⁢ min S ∈ ℂ G × F ⁢ 1 2 ⁢ ∑ f = 1 F ⁢ ⁢  y f ⁡ ( t ) - P f ⁢ s f  2 2 + λ 2 ⁢ ∑ f = 1 F ⁢ ⁢  s f ⁡ ( t ) - s ^ f ⁡ ( t - 1 )  2 2 + μ ⁢ ∑ g = 1 G ⁢ ⁢  ϛ g  2 , where S:=[s 1 ,. .. , s F ], g ′ is the g-th row of S, P f :=[p 1,f ,. .. , p G,f ]ε N×G is the matrix of replicas for ω f , and μ, λ>0 are tuning parameters; and constructing one or more source location maps (SLMs) over per frequency ω f using Ŝ(t), wherein each location on a particular SLM is associated with its corresponding acoustic gain estimate ŝ g,f (t), wherein estimates of the actual locations of the K presumed acoustic sources, {{circumflex over (r)} k (t)} kεκ , correspond to the locations of the K-largest coefficients |ŝ g,f (t)| depicted in the SLM.