Even Castles made of sand, fall into the sea, eventually. - Jimi Hendrix
As Chek Jawa is approximately 1 km by 660 metres, sampling plan and sampling frame has to be thought out carefully so that the interpolation of a logistically possible amount of samples would be as close to the actual scenario as possible.
The distribution of the grain sizes along the shore is affected by many physical variables. Wind direction, the longshore drift and strength of incoming waves affect the distribution of grain sizes on the coast. However, two general trends had been observed about the distribution of grain sizes along the coast. Firstly, largest sediment particles are generally located in the zone of most intense wave breaking and decreases in grain size toward deeper waters and shoreward across the surf and swash zones. (Komar, 1977) Secondly, finer grain sizes tend to be further downcoast in the direction of the longshore drift. This implies that the grain sizes are not equally distributed along the coast. Hence the probability of obtaining a sample of a particular grain size will vary across the coast.
Simple random sampling requires that each possible sample that is to be obtained from a population has to have the same probability of being selected (Petersen and Calvin, 1965). However, based on the two general trends mentioned above, it can be inferred that the probability of selecting any sample with a particular grain size will not hold true. There are many possible variations which may not be detected due to possible clustering of sampling points. Therefore we avoided simple random sampling.
Systematic sampling is another popular method used in the field of soil sampling. This method requires the field to be divided equally such that each unit in the field is regularly spaced out from each other (Cochran, 1953). It requires the sampler to obtain samples from all the strata at a close-enough distance to capture variations, which implies investment of a certain amount of time and effort. Judging from low tide time window for us to collect our samples, we decided that based on these limitations systematic sampling is not feasible due to time and logistical constraints. Another possible disadvantage is a probability of a zero estimate if the sampling interval between systematic samples is larger than the possible period of occurence of significant variations (Christman, 2000; cited in Loh, 2008).
Within these constraints, we used stratified random sampling for our data collection (white squares in "toogle analysis"), based on transect lines which were used in earlier studies, so that it is possible for our raster layer to be superimposed with past data with minimal error from interpolation. Three random sampling points were generated for each 100 metres. This method is used considering time and logistical limitations for data collection and sieving; importantly, it is a compromise between random and systematic sampling, to mitigate extreme clustering and systematic bias.
The properties of the soil are very much determined by the composition of the sizes of its particles. By separating a soil sample, we can obtain its composition. Today many methods exists to segregate the soil particles, to name a few, sieving, sedimentation and pipette sampling methods. The sieving method is amongst the simplest and most commonly used.
Sieving method involve vigorously shaking of sets of differently sized sieves to separate the sample into various grain sizes. It is a convenient method for segregating particles coarser than 0.05mm. (Day,1965) We adopted the recognized Udden-Wentworth grain size classification, which takes into account the characteristics of terrigenous sediments (Wentworth, 1922) (see figure below).

Based on the Udden-Wentworth grain-size classification, we expected that most of our sample grain size to fall within 500 μm to 65 μm, and selected sieves with the following gap sizes: 2mm, 1mm, 500m, 320 μm, 250 μm, 125 μm, 100 μm, 65 μm, 45 μm and sieving pan (to collect grain sizes less than 45 μm ) for segregating the various grain sizes of our samples.
For particle-size analyses, the distribution of a soil is expressed in portions of the various sizes of particles which it contains. The proportions are commonly represented by relative weights of such classes (Day, 1965).
Samples were dried at 400 degrees Celcius for approximately 3hrs to remove excessive moisture content that is held within the sample. This will ensure that all the samples will have approximately the same amount of moisture when it is sieved. A mechanical sieve shaker to ensure that each sample will be subjected to the same length of time and amplitude of shaking. Sieve brushes were used to clean the sieves before segregating the next sample. This is to ensure that earlier segregated samples do not contaminate the later samples as well as keep the sieves dry so as not to collate the sediments and interfere with the results of the segregation.
The sieving method has several limitations. The probability of a particle passing a given sieve in a given time of shaking depends upon the nature of the particle and properties of the sieve (Day, 1965). The dimensions of the particle might hinder the particle from passing through a given sieve opening, even though the particle has the correct mean size to pass though. Thus it will require extensive shaking before the particle is orientated into the “correct” orientation to pass through the sieve opening. Sieve openings are generally unequal in size, requiring extensive shaking before all particles have the opportunity of approaching the largest openings. A "complete" sieving can rarely be met in practical times of shaking, thus we decided to designate our own shaking time (15 minutes). We obtained the 15 minutes sieving time by performing a few trials sieving to see what is the minimum time required to achieve sufficient segregation of the samples without compromising the quality of the data.
As our data collection is discrete by nature, there is a need to interpolate the points to "fill in the gaps" as grain distribution is continuous by nature. Application of interpolation in the field of Geography can often be seen during the estimation of rainfall, temperature and vegetation species. There are several well known techniques around, including inverse distance weighted (IDW) interpolation and kriging.
IDW uses the idea of weighting each data by the inverse of the distance to the estimation location. This implies that the points further away from the estimation location will have an effect of only a power p of the inverse distance, thus not affecting as much as points that are nearer to the estimated point.
IDW assumes a homogeneous weight of influence by all the data points. Kriging, on the other hand, is preceded by a modeling of the spatial structure of the data; unlike IDW, kriging assumes that each data point has a different weight with respect to the estimated point. In addition, points that are behind another sample point will be "screened off" by lowering of weight, thus reducing the effect of this "blocked" point on the estimated point. This effect is known as the Screen effect (Wackernagel, 1995). Kriging estimates an underlying spatial variogram model to represent spatial variability. Furthermore, in kriging, as a statistical method, information such as the estimation of error can be known. This is obtained from the kriging standard deviation. With these considerations, we used kriging as our interpolation technique of choice. Our krigged raster layers of grain sizes were generated in ESRI’s ArcGIS.