Multifractal earth topography
This paper shows how modern ideas of scaling can be used to model topography with various morphologies and also to accurately characterize topography over wide ranges of scales. Our argument is divided in two parts. We first survey the main topographic models and show that they are based on convolut...
Main Authors: | , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Copernicus Publications
2006-01-01
|
Series: | Nonlinear Processes in Geophysics |
Online Access: | http://www.nonlin-processes-geophys.net/13/541/2006/npg-13-541-2006.pdf |
_version_ | 1818928458809475072 |
---|---|
author | J.-S. Gagnon S. Lovejoy S. Lovejoy D. Schertzer |
author_facet | J.-S. Gagnon S. Lovejoy S. Lovejoy D. Schertzer |
author_sort | J.-S. Gagnon |
collection | DOAJ |
description | This paper shows how modern ideas of scaling can be used to model topography with various morphologies and also to accurately characterize topography over wide ranges of scales. Our argument is divided in two parts. We first survey the main topographic models and show that they are based on convolutions of basic structures (singularities) with noises. Focusing on models with large numbers of degrees of freedom (fractional Brownian motion (fBm), fractional Levy motion (fLm), multifractal fractionally integrated flux (FIF) model), we show that they are distinguished by the type of underlying noise. In addition, realistic models require anisotropic singularities; we show how to generalize the basic isotropic (self-similar) models to anisotropic ones. Using numerical simulations, we display the subtle interplay between statistics, singularity structure and resulting topographic morphology. We show how the existence of anisotropic singularities with highly variable statistics can lead to unwarranted conclusions about scale breaking. <P> We then analyze topographic transects from four Digital Elevation Models (DEMs) which collectively span scales from planetary down to 50 cm (4 orders of magnitude larger than in previous studies) and contain more than 2×10<sup>8</sup> pixels (a hundred times more data than in previous studies). We use power spectra and multiscaling analysis tools to study the global properties of topography. We show that the isotropic scaling for moments of order ≤2 holds to within ±45% down to scales ≈40 m. We also show that the multifractal FIF is easily compatible with the data, while the monofractal fBm and fLm are not. We estimate the universal parameters (α, <i>C<sub>1</sub></i>) characterizing the underlying FIF noise to be (1.79, 0.12), where α is the degree of multifractality (0≤α≤2, 0 means monofractal) and <i>C<sub>1</sub></i> is the degree of sparseness of the surface (0≤<i>C<sub>1</sub></i>, 0 means space filling). In the same way, we investigate the variation of multifractal parameters between continents, oceans and continental margins. Our analyses show that no significant variation is found for (α, <i>C<sub>1</sub></i>) and that the third parameter <i>H</i>, which is a degree of smoothing (higher <i>H</i> means smoother), is variable: our estimates are <i>H</i>=0.46, 0.66, 0.77 for bathymetry, continents and continental margins. An application we developped here is to use (α, <i>C<sub>1</sub></i>) values to correct standard spectra of DEMs for multifractal resolution effects. |
first_indexed | 2024-12-20T03:29:14Z |
format | Article |
id | doaj.art-8d65d47cc88041f59a27d76b99f675d1 |
institution | Directory Open Access Journal |
issn | 1023-5809 1607-7946 |
language | English |
last_indexed | 2024-12-20T03:29:14Z |
publishDate | 2006-01-01 |
publisher | Copernicus Publications |
record_format | Article |
series | Nonlinear Processes in Geophysics |
spelling | doaj.art-8d65d47cc88041f59a27d76b99f675d12022-12-21T19:55:01ZengCopernicus PublicationsNonlinear Processes in Geophysics1023-58091607-79462006-01-01135541570Multifractal earth topographyJ.-S. GagnonS. LovejoyS. LovejoyD. SchertzerThis paper shows how modern ideas of scaling can be used to model topography with various morphologies and also to accurately characterize topography over wide ranges of scales. Our argument is divided in two parts. We first survey the main topographic models and show that they are based on convolutions of basic structures (singularities) with noises. Focusing on models with large numbers of degrees of freedom (fractional Brownian motion (fBm), fractional Levy motion (fLm), multifractal fractionally integrated flux (FIF) model), we show that they are distinguished by the type of underlying noise. In addition, realistic models require anisotropic singularities; we show how to generalize the basic isotropic (self-similar) models to anisotropic ones. Using numerical simulations, we display the subtle interplay between statistics, singularity structure and resulting topographic morphology. We show how the existence of anisotropic singularities with highly variable statistics can lead to unwarranted conclusions about scale breaking. <P> We then analyze topographic transects from four Digital Elevation Models (DEMs) which collectively span scales from planetary down to 50 cm (4 orders of magnitude larger than in previous studies) and contain more than 2×10<sup>8</sup> pixels (a hundred times more data than in previous studies). We use power spectra and multiscaling analysis tools to study the global properties of topography. We show that the isotropic scaling for moments of order ≤2 holds to within ±45% down to scales ≈40 m. We also show that the multifractal FIF is easily compatible with the data, while the monofractal fBm and fLm are not. We estimate the universal parameters (α, <i>C<sub>1</sub></i>) characterizing the underlying FIF noise to be (1.79, 0.12), where α is the degree of multifractality (0≤α≤2, 0 means monofractal) and <i>C<sub>1</sub></i> is the degree of sparseness of the surface (0≤<i>C<sub>1</sub></i>, 0 means space filling). In the same way, we investigate the variation of multifractal parameters between continents, oceans and continental margins. Our analyses show that no significant variation is found for (α, <i>C<sub>1</sub></i>) and that the third parameter <i>H</i>, which is a degree of smoothing (higher <i>H</i> means smoother), is variable: our estimates are <i>H</i>=0.46, 0.66, 0.77 for bathymetry, continents and continental margins. An application we developped here is to use (α, <i>C<sub>1</sub></i>) values to correct standard spectra of DEMs for multifractal resolution effects.http://www.nonlin-processes-geophys.net/13/541/2006/npg-13-541-2006.pdf |
spellingShingle | J.-S. Gagnon S. Lovejoy S. Lovejoy D. Schertzer Multifractal earth topography Nonlinear Processes in Geophysics |
title | Multifractal earth topography |
title_full | Multifractal earth topography |
title_fullStr | Multifractal earth topography |
title_full_unstemmed | Multifractal earth topography |
title_short | Multifractal earth topography |
title_sort | multifractal earth topography |
url | http://www.nonlin-processes-geophys.net/13/541/2006/npg-13-541-2006.pdf |
work_keys_str_mv | AT jsgagnon multifractalearthtopography AT slovejoy multifractalearthtopography AT slovejoy multifractalearthtopography AT dschertzer multifractalearthtopography |