The Pierre Auger Observatory is being used to study cosmic rays with energies larger than 10 EeV. An essential quantity that must be deduced from data is the lateral distribution function (LDF). Knowledge of the LDF is important for the reconstruction of the shower core and the shower direction. Here we describe how the LDF is measured using the large sample of events recorded with the surface detector (SD) array and with a small sample observed with the fluorescence detectors (FD). For hybrid events, in which SD and FD measurements of the same shower are available, the core position is much better constrained than for SD-only events, thus providing an important cross-check on the LDF determined from SD measurements alone.

Evidence of azimuthal asymmetries in the time structure and signal size have been found in non-vertical showers at the Pierre Auger Observatory. It has been previously shown that the asymmetry in time distributions offers a new possibility for the determination of the mass composition. New studies have demonstrated that the dependence of the asymmetry parameter in the rise-time and fall-time distributions with sec(theta) shows a clear peak. Both, the position of the peak, X_asymax, and the size of the asymmetry at X_asymax are sensitive to primary mass composition and have a small dependence on energy. In this paper a study of the discriminating power of the new observables to separate primary species is presented.

Linsley used the Volcano Ranch array to collect data on the lateral distribution of showers produced by cosmic rays at energies above 10^17 eV. Very precise measurements of the steepness of the lateral distribution function were made on 366 events. The current availability of sophisticated hadronic interaction models has prompted an interpretation of the measurements. In this analysis we use the AIRES Monte Carlo code to generate showers, together with GEANT4 to simulate the detector response to ground particles. The results show that, with the assumption of a bi-modal proton and iron mix, iron is the dominant component of cosmic rays between 5x10^17 and 10^19 eV, assuming that hadronic interactions are well-described by QGSJET at this energy range.

Accurate measurements of the lateral distribution of extensive air showers produced by cosmic rays of energy greater than $10^{17}$ eV were made in the 1970s by Linsley. At the time, the state of knowledge about the best hadronic interaction model to use to describe such data prevented conclusions about the cosmic ray mass from being drawn. We have used a modern model, {\sc qgsjet98}, to infer the primary mass from these data, using the very careful records left by Linsley. We find that at a median energy of $10^{18}$ eV, the data are well described by an iron-dominated composition (88 $\pm$ 6(stat) $\pm$ 18 (syst))%. We discuss the systematic errors in this estimate that arise from model uncertainties and from the range of energies used in the work of Linsley. These data are used with the permission of the late John Linsley to whom this work is dedicated.

In the course of its operation, the Volcano Ranch array collected data on the lateral distribution of showers produced by cosmic rays at energies above $10^{17}$ {\rm eV}. From these data very precise measurements of the steepness of the lateral distribution function, characterized by the $\eta$ parameter, were made. The current availability of sophisticated hadronic interaction models has prompted a reinterpretation of the measurements. We use the interaction models {\sc qgsjet} and {\sc sibyll} in the {\sc aires} Monte Carlo code to generate showers together with {\sc geant4} to simulate the response of the detectors to ground particles. As part of an effort to estimate the primary mass composition of cosmic rays at this energy range, we present the results of our preliminary analysis of the distribution of $\eta$.