A method (600, 630) of estimating a time-varying multipath communication channel, the method comprising: receiving a frequency-domain channel estimate derived from an orthogonal frequency division multiplex, OFDM, data packet, the data packet being transmitted over the multipath communication channel; based on the frequency-domain channel estimate: generating an estimate of tap parameters in a time-domain model of the multipath communication channel, the model characterising the multipath communication channel as a delay line with a plurality of taps (h 1, h 2), each tap having tap parameters for estimation of a delay, a Doppler and a complex gain for the tap, wherein the estimated tap parameters are constant for the data packet; re-estimating the estimate of tap parameters in the time-domain model of the multipath communication channel, said re-estimating comprising, for one or more of a plurality of multipath components: subtracting all of the plurality of multipath components from the estimate of the multipath communication channel except for a respective multipath component; generating a re-estimated estimate of the delay tap parameter and the Doppler tap parameter for the respective multipath component; and performing a least-squares estimation to generate a re-estimated estimate of the gain tap parameter for the respective multipath component; wherein the re-estimating is repeated to generate re-estimated estimates for the plurality of multipath components. receiving a frequency-domain channel estimate derived from an orthogonal frequency division multiplex, OFDM, data packet, the data packet being transmitted over the multipath communication channel; based on the frequency-domain channel estimate: generating an estimate of tap parameters in a time-domain model of the multipath communication channel, the model characterising the multipath communication channel as a delay line with a plurality of taps (h 1, h 2), each tap having tap parameters for estimation of a delay, a Doppler and a complex gain for the tap, wherein the estimated tap parameters are constant for the data packet; re-estimating the estimate of tap parameters in the time-domain model of the multipath communication channel, said re-estimating comprising, for one or more of a plurality of multipath components: subtracting all of the plurality of multipath components from the estimate of the multipath communication channel except for a respective multipath component; generating a re-estimated estimate of the delay tap parameter and the Doppler tap parameter for the respective multipath component; and performing a least-squares estimation to generate a re-estimated estimate of the gain tap parameter for the respective multipath component; generating an estimate of tap parameters in a time-domain model of the multipath communication channel, the model characterising the multipath communication channel as a delay line with a plurality of taps (h 1, h 2), each tap having tap parameters for estimation of a delay, a Doppler and a complex gain for the tap, wherein the estimated tap parameters are constant for the data packet; re-estimating the estimate of tap parameters in the time-domain model of the multipath communication channel, said re-estimating comprising, for one or more of a plurality of multipath components: subtracting all of the plurality of multipath components from the estimate of the multipath communication channel except for a respective multipath component; generating a re-estimated estimate of the delay tap parameter and the Doppler tap parameter for the respective multipath component; and performing a least-squares estimation to generate a re-estimated estimate of the gain tap parameter for the respective multipath component; subtracting all of the plurality of multipath components from the estimate of the multipath communication channel except for a respective multipath component; generating a re-estimated estimate of the delay tap parameter and the Doppler tap parameter for the respective multipath component; and performing a least-squares estimation to generate a re-estimated estimate of the gain tap parameter for the respective multipath component; wherein the re-estimating is repeated to generate re-estimated estimates for the plurality of multipath components.