A methodology based on compressive sampling is developed for incomplete wind time-histories reconstruction and extrapolation in a single spatial dimension, as well as for related stochastic field statistics estimation. This relies on l1-norm minimization in conjunction with an adaptive basis re-weighting scheme. Indicatively, the proposed methodology can be employed for monitoring of wind turbine systems, where the objective relates to either reconstructing incomplete time-histories measured at specific points along the height of a turbine tower, or to extrapolating to other locations in the vertical dimension where sensors and measurement records are not available. Further, the methodology can be used potentially for environmental hazard modeling within the context of performance-based design optimization of structural systems. Unfortunately, a straightforward implementation of the aforementioned approach to account for two spatial dimensions is hindered by significant, even prohibitive in some cases, computational cost. In this regard, to address computational challenges associated with higher-dimensional domains, a methodology based on low rank matrices and nuclear norm minimization is developed next for wind field extrapolation in two spatial dimensions. The efficacy of the proposed methodologies is demonstrated by considering various numerical examples. These refer to reconstruction of wind time-histories with missing data compatible with a joint wavenumber-frequency power spectral density, as well as to extrapolation to various locations in the spatial domain.