One important component of deep crustal reflection seismic data in the absence of drill-hole data and surface-outcrop constraints is classifying and quantifying reflectivity patterns. One approach to this component uses a recently developed data-decomposition technique, seismic skeletonization. Skeletonized coherent events and their attributes are identified and stored in a relational database, allowing easy visualization and parameterization of the reflected wavefield. Because one useful attribute, the instantaneous frequency, is difficult to derive within the current framework of skeletonization, timefrequency analysis and a new method, empirical mode skeletonization, are used to derive it. Other attributes related to timefrequency analysis that can be derived from the methods can be used for shallow and deep reflection seismic interpretation and can supplement the seismic attributes accrued from seismic skeletonization. Bright reflections observed from below the sedimentary basin in the Southern Alberta Lithosphere Transect have recently been interpreted to be caused by highly reflective sills. Timefrequency analysis of one of these reflections shows the lateral variation of energy with instantaneous frequency for any given time and the lateral variation of energy with time for any instantaneous frequency. Results from empirical mode skeletonization for the same segment of data illustrate the differences in the instantaneous frequencies among the intrinsic modes of the data. Thus, timefrequency distribution of amplitude or energy for any signal may be a good indicator of compositional differences that can vary from one location to another.