We have investigated the use of ground-penetrating radar (GPR) as a means of characterizing the heterogeneity of the subsurface. Radar data were collected at several sites in southwestern British Columbia underlain by glaciodeltaic sediments. A cliff face study was conducted in which geostatistical analysis of a digitized photograph of the face and the radar image of the face showed excellent agreement in the maximum correlation direction and the correlation length determined from these two data sets. Other two-dimensional (2-D) sections of radar data were divided into sedimentary architectural elements on the basis of the distinct radar appearance of these sedimentary units. Examples of four sedimentary units were used to obtain semivariograms from the radar data and resulted in maximum correlation lengths between 0.5 and 4.8 m. A 3-D radar survey, collected over a package of gravel and sand foresets, was analyzed to determine the paleoflow direction; a correlation length of 4 m was found in that direction. ¿ 1998 American Geophysical Union