The temporal and spatial data analysed in, for example, ecology or climatology, are often hierarchically structured, carrying information in different scales. An important goal of data analysis is then to decompose the observed signal into distinctive hierarchical levels and to determine the size of the features that each level represents. Using differences of smooths, scale space multiresolution analysis decomposes a signal into additive components associated with different levels of scales present in the data. The smoothing levels used to compute the differences are determined by the local minima of the norm of the so-called scale-derivative of the signal. While this procedure accomplishes the first goal, the hierarchical decomposition of the signal, it does not achieve the second goal, the determination of the actual size of the features corresponding to each hierarchical level. Here, we show that the maximum of the scale-derivative norm of an extracted hierarchical component can be used to estimate its characteristic feature size. The feasibility of the method is demonstrated using an artificial image and a time series of a drought index, based on climate reconstructions from long tree ring chronologies.We propose a two-step procedure for the analysis of such hierarchical data. First, the signal is decomposed into components that represent its variation in different levels of hierarchy, and second, the characteristic feature size of each hierarchical level is determined.In the first step, scale space multiresolution analysis (Holmström et al., 2011), an instance of statistical scale space methodology, is applied. Conventional scale space analysis considers smooths of a signal, interpreting each smooth