The presence and configuration of subglacial water and sediment maintain the fast flow of arterial ice streams in Antarctica and airborne ice-penetrating radar data represent a potential resource of information about the ice-bed interface. In this thesis an original contribution to the exploration of Antarctic subglacial environments is made through the analysis of airborne surveys from Evans, Institute and Möller Ice Streams, West Antarctica. The primary approach employed is the derivation of bed-returned power (BRP), a proxy for ice-bed reflectivity, which is strongly influenced by the presence of liquid water. Estimating radar englacial attenuation (EA) accurately is a critical part of BRP analysis and a modelled approach is primarily used. BRP is derived across Evans Ice Stream and shows large-scale patterns relating to hypothesised hydrological and geological contrasts at the ice-bed. These results are developed to investigate the influence of: (1) adopted EA correction; (2) the influence of assigned ice dielectric properties in modelled EA; (3) subglacial roughness and (4) the spatial scale over which BRP is derived. Some areas of high basal drag can be detected with BRP analysis, indicating that variations in subglacial hydrology are responsible for their existence. The widely-used empirical method of estimating EA by relating ice thickness to uncorrected BRP is shown to be unreliable where ice properties change along a transect. Monte Carlo error analysis of modelled EA shows that poorly constrained ice dielectric properties also result in significant BRP uncertainty. BRP beneath Institute and Möller Ice Streams is derived on catchment- and local-scales over hypothesised subglacial features. Bungenstock Ice Rise is marked with a clear BRP signal but the locations of "active" lakes, as delineated by satellite altimetry, do not. The sensitivity of idealised flow paths to surface change and grid size are investigated. Potential future research directions regarding BRP analysis are discussed.