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...

Full description

Bibliographic Details
Main Authors: J.-S. Gagnon, S. Lovejoy, D. Schertzer
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&times;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 &le;2 holds to within &plusmn;45% down to scales &asymp;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&le;&alpha;&le;2, 0 means monofractal) and <i>C<sub>1</sub></i> is the degree of sparseness of the surface (0&le;<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&times;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 &le;2 holds to within &plusmn;45% down to scales &asymp;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&le;&alpha;&le;2, 0 means monofractal) and <i>C<sub>1</sub></i> is the degree of sparseness of the surface (0&le;<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