Hydraulic diffusivity is an essential parameter in understanding and predicting an aquifer's response to perturbations. Common techniques of aquifer hydraulic properties estimation, such as slug tests and pumping tests, require an artificial flow and pressure perturbation. This can be expensive and extremely invasive, and therefore sometimes impossible to apply due to environmental risks or resource limitations. However, many hydrogeologic systems have time-periodic behavior due to natural and anthropogenic influences (Depner & Rasmussen, 2016). A common example are aquifers which are connected to naturally oscillating water bodies that cause tide pulses propagating deep into the reservoir as diffusive waves. This is the case for the Meghna River, which has strong ocean-driven tides with semi-diurnal amplitudes of more than 0.5 m even 200 km inland of the Bay of Bengal. The Meghna River receives groundwater which carries high levels of naturally occurring, dissolved arsenic. Tides may drive mixing of oxic river water with anoxic groundwater, forcing the precipitation of iron oxides and the accumulation (sorption) of dissolved arsenic in those oxides, which pro-Abstract This note introduces the estimation of aquifer diffusivity (D) through simultaneous inversion of the attenuation and lag of multiple head fluctuation frequencies due to a dynamic source or boundary, that is, a river. Spectral analysis, with optimized moving time window length and step size, was used to extract the dominant constituents and their attenuation through space; the cross-power spectral density method was used to determine time lags. The Jacob-Ferris analytical model was then used for inverting for D. Unlike most similar applications to date, here we propose using all frequencies with robust signal-to-noise ratios (five total in our test cases) and both the amplitude attenuation and time lag in the inversion. The method was implemented using observations from wells in the banks of the fluctuating Meghna River in Bangladesh that is connected with a semi-confined sandy alluvial aquifer. The estimated D using the technique provides estimates that are very similar to those from pumping tests. The estimates are more accurate compared to previous implementation of the Jacob-Ferris model on the same data that used only a dominant frequency's amplitude attenuation or time lag. The workflow and codes for the analysis are provided for straightforward implementation of the robust and cost-effective method.Plain Language Summary There is a rich tradition of inferring aquifer diffusivity (D) (transmissivity divided by storativity) by measuring the speed and decay of periodic pressure pulses. The analytical equations that describe these pulses, although simple to apply, are suspected of being inaccurate since they do not explicitly account for aquifer complexity. Our goal in this study is to evaluate if a simple model can make accurate predictions when modern signal processing tools are used. We utilize five natural solar and lunar tide signals to prod...