Analysis and interpretation of spatial variability of soils is a keystone in site‐specific farming. Soil survey maps may have up to 0.41‐ha inclusions of dissimilar soils within a mapping unit. The objectives of this study were to determine the degree of spatial variability of soil physical properties and variance structure, and to model the sampling interval of alluvial floodplain soils. Soil profiles (n= 209) from 18 parallel transects were sampled with a mean separation distance of 79.4 m. Each profile was classified into surface, subsurface, and deep horizons. Structural analysis of soil bulk density (ρb), sand, clay, saturated hydraulic conductivity (Ks), volumetric water content (θv) at seven pressure potentials (Ψa) (−1, −10, −33, −67, −100, −500, and −1500 kPa) were modeled for the three horizons. Variance of soil physical properties varied from as low as 0.01% (ρb) to as high as 1542% (Ks). The LSD test indicated significant (P< 0.05) differences in sand, clay, ρb,Ks, and θvat various ΨaGeostatistical analyses illustrated that the spatially dependent stochastic component was predominant over the nugget effect. Structured semivariogram functions of each variable were used in generating fine‐scale kriged contour maps. Overall autocorrelation, Moran'sI, indicated a 400‐m sampling range would be adequate for detection of spatial structure of sand, silt, clay, and a 100‐m sampling range for soil hydraulic properties and ρbThe magnitude and spatial patterns soil physical property variability have implications for variable rate applications and design of soil sampling strategies in alluvial floodplain soils.