The distribution of data records for the maximum horizontal stress orientation SHmax in the Earth’s crust is sparse and very unequally. In order to analyse the stress pattern and its wavelength or to predict the mean SHmax orientation on a regular grid, statistical interpolation as conducted e.g. by Coblentz and Richardson (1995), Müller et al. (2003), Heidbach and Höhne (2008), Heidbach et al. (2010) or Reiter et al. (2014) is necessary. Based on their work we wrote the Matlab® script Stress2Grid that provides several features to analyse the mean SHmax pattern. The script facilitates and speeds up this analysis and extends the functionality compared to aforementioned publications. The script is complemented by a number of example and input files as described in the WSM Technical Report (Ziegler and Heidbach, 2017, http://doi.org/10.2312/wsm.2017.002).
The script provides two different concepts to calculate the mean SHmax orientation on a regular grid. The first is using a fixed search radius around the grid point and computes the mean SHmax orientation if sufficient data records are within the search radius. The larger the search radius the larger is the filtered wavelength of the stress pattern. The second approach is using variable search radii and determines the search radius for which the variance of the mean SHmax orientation is below a given threshold. This approach delivers mean SHmax orientations with a user-defined degree of reliability. It resolves local stress perturbations and is not available in areas with conflicting information that result in a large variance. Furthermore, the script can also estimate the deviation between plate motion direction and the mean SHmax orientation.
Authors
Ziegler, Moritz;GFZ German Research Centre for Geosciences, Potsdam, Germany
Heidbach, Oliver;GFZ German Research Centre for Geosciences, Potsdam, Germany
Contact
Heidbach, Oliver; GFZ German Research Centre for Geosciences, Potsdam, Germany;
Contributors
GFZ German Rersearch Centre for Geosciences
Keywords
World Stress Map, Stress orientation data analysis tool for Matlab
affiliation (affiliationIdentifier=0000-0002-2288-2820 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
affiliation (affiliationIdentifier=0000-0001-8009-5422 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
titles
title: Matlab script Stress2Grid
publisher: GFZ Data Services
publicationYear: 2017
subjects
subject: World Stress Map
subject: Stress orientation data analysis tool for Matlab
affiliation (affiliationIdentifier=0000-0002-2288-2820 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
affiliation (affiliationIdentifier=0000-0001-8009-5422 affiliationIdentifierScheme=ORCID): GFZ German Research Centre for Geosciences, Potsdam, Germany
contributor (contributorType=HostingInstitution)
contributorName: GFZ German Rersearch Centre for Geosciences
affiliation (affiliationIdentifier= affiliationIdentifierScheme=): GFZ German Research Centre for Geosciences, Potsdam, Germany
contributor (contributorType=ContactPerson)
contributorName: Heidbach, Oliver
affiliation (affiliationIdentifier= affiliationIdentifierScheme=): GFZ German Research Centre for Geosciences, Potsdam, Germany
CharacterString: The distribution of data records for the maximum horizontal stress orientation SHmax in the Earth’s crust is sparse and very unequally. In order to analyse the stress pattern and its wavelength or to predict the mean SHmax orientation on a regular grid, statistical interpolation as conducted e.g. by Coblentz and Richardson (1995), Müller et al. (2003), Heidbach and Höhne (2008), Heidbach et al. (2010) or Reiter et al. (2014) is necessary. Based on their work we wrote the Matlab® script Stress2Grid that provides several features to analyse the mean SHmax pattern. The script facilitates and speeds up this analysis and extends the functionality compared to aforementioned publications. The script is complemented by a number of example and input files as described in the WSM Technical Report (Ziegler and Heidbach, 2017, http://doi.org/10.2312/wsm.2017.002).
The script provides two different concepts to calculate the mean SHmax orientation on a regular grid. The first is using a fixed search radius around the grid point and computes the mean SHmax orientation if sufficient data records are within the search radius. The larger the search radius the larger is the filtered wavelength of the stress pattern. The second approach is using variable search radii and determines the search radius for which the variance of the mean SHmax orientation is below a given threshold. This approach delivers mean SHmax orientations with a user-defined degree of reliability. It resolves local stress perturbations and is not available in areas with conflicting information that result in a large variance. Furthermore, the script can also estimate the deviation between plate motion direction and the mean SHmax orientation.
CI_OnLineFunctionCode (codeList=http://www.isotc211.org/2005/resources/Codelist/gmxCodelists.xml#CI_OnLineFunctionCode codeListValue=http://www.isotc211.org/2005/resources/Codelist/gmxCodelists.xml#CI_OnLineFunctionCode_information): information
Released
Software
Abstract
;GFZ German Research Centre for Geosciences, Potsdam, Germany
Files