Unit-level and area-level small area estimation under heteroscedasticity using digital aerial photogrammetry data

In many applications, estimates are required for small sub-populations with so few (or no) sample plots that direct estimators that do not utilize auxiliary variables (e.g. remotely sensed data) are not applicable or result in low precision. This problem is overcome in small area estimation (SAE) by linking the variable of interest to auxiliary variables using a model. Two types of models can be distinguished based on the scale on which they operate: i) Unit-level models are applied in the well-known area-based approach (ABA) and are commonly used in forest inventories supported by fine-resolution 3D remote sensing data such as airborne laser scanning (ALS) or digital aerial photogrammetry (AP); ii) Area-level models, where the response is a direct estimate based on a sample within the domain and the explanatory variables are aggregated auxiliary variables, are less frequently applied. Estimators associated with these two model types can make use of sample plots within domains if available and reduce to so-called synthetic estimators in domains where no sample plots are available. We used both model types and their associated model-based estimators in the same study area with AP data as auxiliary variables. Heteroscedasticity, i.e. for continuous dependent variables typically an increasing dispersion of residuals with increasing predictions, is often observed in models linking field- and remotely sensed data. This violates the model assumption that the distribution of the residual errors is constant. Complying with model assumptions is required for model-based methods to result in reliable estimates. Addressing heteroscedasticity in models had considerable impacts on standard errors. When complying with model assumptions, the precision of estimates based on unit-level models was, on average, considerably greater (29%-31% smaller standard errors) than those based on area-level models. Area-level models may nonetheless be attractive because they allow the use of sampling designs that do not easily link to remotely sensed data, such as variable radius plots.