Approximate multifractal correlation and products of universal multifractal fields, with application to rainfall data

Universal multifractals (UMs) have been widely used to simulate and characterize, with the help of only two physically meaningful parameters, geophysical fields that are extremely variable across a wide range of scales. Such a framework relies on the assumption that the underlying field is generated through a multiplicative cascade process. Derived analysis techniques have been extended to study correlations between two fields not only at a single scale and for a single statistical moment as with the covariance, but across scales and for all moments. Such a framework of joint multifractal analysis is used here as a starting point to develop and test an approach enabling correlations between UM fields to be analysed and approximately simulated. First, the behaviour of two fields consisting of renormalized multiplicative power law combinations of two UMfields is studied. It appears that in the general case the resulting fields can be well approximated byUMfields with known parameters. Limits of this approximation will be quantified and discussed. Techniques to retrieve the UM parameters of the underlying fields as well as the exponents of the combination have been developed and successfully tested on numerical simulations. In a second step tentative correlation indicators are suggested. Finally the suggested approach is implemented to study correlation across scales of detailed rainfall data collected with the help of disdrometers of the Fresnel platform of Ecole des Ponts ParisTech (see available data at https://hmco.enpc.fr/portfolio-archive/taranis-observatory/, last access: 12 March 2020). More precisely, four quantities are used: the rain rate (R), the liquid water content (LWC) and the total drop concentration (N_t) along with the mass weighed diameter (D_m), which are commonly used to characterize the drop size distribution. Correlations across scales are quantified. Their relative strength (very strong between R and LWC, strong between DSD features and R or LWC, almost null between N_t and D_m) is discussed.