The impact of COVID-19 on relative changes in aggregated mobility using mobile-phone data

Abstract: Evaluating relative changes leads to additional insights which would remain hidden when only evaluating absolute changes. We analyze a dataset describing mobility of mobile phones in Austria before, during COVID-19 lock-down measures until recent. By applying compositional data analysis we show that formerly hidden information becomes available: we see that the elderly population groups increase relative mobility and that the younger groups especially on weekends also do not decrease their mobility as much as the others.

\section{Introduction}

In order to assist authors in the process of preparing a manuscript for a journal, the Taylor \& Francis `\textsf{Interact}' layout style has been implemented as a \LaTeXe\ class file based on the \texttt{article} document class. A \textsc{Bib}\TeX\ bibliography style file and a sample bibliography are also provided in order to assist with the formatting of your references.

Commands that differ from or are provided in addition to standard \LaTeXe\ are described in this document, which is \emph{not} a substitute for a \LaTeXe\ tutorial.

The \texttt{interacttfssample.tex} file can be used as a template for a manuscript by cutting, pasting, inserting and deleting text as appropriate, using the preamble and the \LaTeX\ environments provided (e.g.\ \verb"\begin{abstract}", \verb"\begin{keywords}").

\subsection{The \textsf{Interact} class file}\label{class}

The \texttt{interact} class file preserves the standard \LaTeXe\ interface such that any document that can be produced using \texttt{article.cls} can also be produced with minimal alteration using the \texttt{interact} class file as described in this document.

If your article is accepted for publication it will be typeset as the journal requires in Minion Pro and/or Myriad Pro. Since most authors will not have these fonts installed, the page make-up is liable to alter slightly with the change of font. Also, the \texttt{interact} class file produces only single-column format, which is preferred for peer review and will be converted to two-column format by the typesetter if necessary during preparation of the proofs. Please therefore do not try to match the typeset format exactly, but use the standard \LaTeX\ fonts instead and ignore details such as slightly long lines of text or figures/tables not appearing in exact synchronization with their citations in the text: these details will be dealt with by the typesetter. Similarly, it is unnecessary to spend time addressing warnings in the log file -- if your .tex file compiles to produce a PDF document that correctly shows how you wish your paper to appear, such warnings will not prevent your source files being imported into the typesetter's program.

\subsection{Submission of manuscripts prepared using \emph{\LaTeX}}

Manuscripts for possible publication should be submitted to the Editors for review as directed in the journal's Instructions for Authors, and in accordance with any technical instructions provided in the journal's ScholarOne Manuscripts or Editorial Manager site. Your \LaTeX\ source file(s), the class file and any graphics files will be required in addition to the final PDF version when final, revised versions of accepted manuscripts are submitted.

Please ensure that any author-defined macros used in your article are gathered together in the preamble of your .tex file, i.e.\ before the \verb"\begin{document}" command. Note that if serious problems are encountered in the coding of a document (missing author-defined macros, for example), the typesetter may resort to rekeying it.

\section{Using the \texttt{interact} class file}

For convenience, simply copy the \texttt{interact.cls} file into the same directory as your manuscript files (you do not need to install it in your \TeX\ distribution). In order to use the \texttt{interact} document class, replace the command \verb"\documentclass{article}" at the beginning of your document with the command \verb"\documentclass{interact}".

The following document-class options should \emph{not} be used with the \texttt{interact} class file:
\begin{itemize}
  \item \texttt{10pt}, \texttt{11pt}, \texttt{12pt} -- unavailable;
  \item \texttt{oneside}, \texttt{twoside} -- not necessary, \texttt{oneside} is the default;
  \item \texttt{leqno}, \texttt{titlepage} -- should not be used;
  \item \texttt{twocolumn} -- should not be used (see Subsection~\ref{class});
  \item \texttt{onecolumn} -- not necessary as it is the default style.
\end{itemize}
To prepare a manuscript for a journal that is printed in A4 (two column) format, use the \verb"largeformat" document-class option provided by \texttt{interact.cls}; otherwise the class file produces pages sized for B5 (single column) format by default. The \texttt{geometry} package should not be used to make any further adjustments to the page dimensions.

\section{Additional features of the \texttt{interact} class file}

\subsection{Title, authors' names and affiliations, abstracts and article types}

The title should be generated at the beginning of your article using the \verb"\maketitle" command.
In the final version the author name(s) and affiliation(s) must be followed immediately by \verb"\maketitle" as shown below in order for them to be displayed in your PDF document.
To prepare an anonymous version for double-blind peer review, you can put the \verb"\maketitle" between the \verb"\title" and the \verb"\author" in order to hide the author name(s) and affiliation(s) temporarily.
Next you should include the abstract if your article has one, enclosed within an \texttt{abstract} environment.
The \verb"\articletype" command is also provided as an \emph{optional} element which should \emph{only} be included if your article actually needs it.
For example, the titles for this document begin as follows:
\begin{verbatim}
\articletype{ARTICLE TEMPLATE}

\title{Taylor \& Francis \LaTeX\ template for authors (\textsf{Interact}
layout + reference style S)}

\author{
\name{A.~N. Author\textsuperscript{a}\thanks{CONTACT A.~N. Author.
Email: latex.helpdesk@tandf.co.uk} and John Smith\textsuperscript{b}}
\affil{\textsuperscript{a}Taylor \& Francis, 4 Park Square, Milton
Park, Abingdon, UK; \textsuperscript{b}Institut f\"{u}r Informatik,
Albert-Ludwigs-Universit\"{a}t, Freiburg, Germany} }

\maketitle

\begin{abstract}
This template is for authors who are preparing a manuscript for a
Taylor \& Francis journal using the \LaTeX\ document preparation system
and the \texttt{interact} class file, which is available via selected
journals' home pages on the Taylor \& Francis website.
\end{abstract}
\end{verbatim}

An additional abstract in another language (preceded by a translation of the article title) may be included within the \verb"abstract" environment if required.

A graphical abstract may also be included if required. Within the \verb"abstract" environment you can include the code
\begin{verbatim}
\\\resizebox{25pc}{!}{\includegraphics{abstract.eps}}
\end{verbatim}
where the graphical abstract is to appear, where \verb"abstract.eps" is the name of the file containing the graphic (note that \verb"25pc" is the recommended maximum width, expressed in pica, for the graphical abstract in your manuscript).

\subsection{Abbreviations}

A list of abbreviations may be included if required, enclosed within an \texttt{abbreviations} environment, i.e.\ \verb"\begin{abbreviations}"\ldots\verb"\end{abbreviations}", immediately following the \verb"abstract" environment.

\subsection{Keywords}

A list of keywords may be included if required, enclosed within a \texttt{keywords} environment, i.e.\ \verb"\begin{keywords}"\ldots\verb"\end{keywords}". Additional keywords in other languages (preceded by a translation of the word `keywords') may also be included within the \verb"keywords" environment if required.

\subsection{Subject classification codes}

AMS, JEL or PACS classification codes may be included if required. The \texttt{interact} class file provides an \texttt{amscode} environment, i.e.\ \verb"\begin{amscode}"\ldots\verb"\end{amscode}", a \texttt{jelcode} environment, i.e.\ \verb"\begin{jelcode}"\ldots\verb"\end{jelcode}", and a \texttt{pacscode} environment, i.e.\ \verb"\begin{pacscode}"\ldots\verb"\end{pacscode}" to assist with this.

\subsection{Additional footnotes to the title or authors' names}

The \verb"\thanks" command may be used to create additional footnotes to the title or authors' names if required. Footnote symbols for this purpose should be used in the order
$^\ast$~(coded as \verb"$^\ast$"), $\dagger$~(\verb"$\dagger$"), $\ddagger$~(\verb"$\ddagger$"), $\S$~(\verb"$\S$"), $\P$~(\verb"$\P$"), $\|$~(\verb"$\|$"),
$\dagger\dagger$~(\verb"$\dagger\dagger$"), $\ddagger\ddagger$~(\verb"$\ddagger\ddagger$"), $\S\S$~(\verb"$\S\S$"), $\P\P$~(\verb"$\P\P$").

Note that any \verb"footnote"s to the main text will automatically be assigned the superscript symbols 1, 2, 3, etc. by the class file.\footnote{If preferred, the \texttt{endnotes} package may be used to set the notes at the end of your text, before the bibliography. The symbols will be changed to match the style of the journal if necessary by the typesetter.}

\section{Some guidelines for using the standard features of \LaTeX}

\subsection{Sections}

The \textsf{Interact} layout style allows for five levels of section heading, all of which are provided in the \texttt{interact} class file using the standard \LaTeX\ commands \verb"\section", \verb"\subsection", \verb"\subsubsection", \verb"\paragraph" and \verb"\subparagraph". Numbering will be automatically generated for all these headings by default.

\subsection{Lists}

Numbered lists are produced using the \texttt{enumerate} environment, which will number each list item with arabic numerals by default. For example,
\begin{enumerate}
  \item first item
  \item second item
  \item third item
\end{enumerate}
was produced by
\begin{verbatim}
\begin{enumerate}
  \item first item
  \item second item
  \item third item
\end{enumerate}
\end{verbatim}
Alternative numbering styles can be achieved by inserting an optional argument in square brackets to each \verb"item", e.g.\ \verb"\item[(i)] first item"\, to create a list numbered with roman numerals at level one.

Bulleted lists are produced using the \texttt{itemize} environment. For example,
\begin{itemize}
  \item First bulleted item
  \item Second bulleted item
  \item Third bulleted item
\end{itemize}
was produced by
\begin{verbatim}
\begin{itemize}
  \item First bulleted item
  \item Second bulleted item
  \item Third bulleted item
\end{itemize}
\end{verbatim}

\subsection{Figures}

The \texttt{interact} class file will deal with positioning your figures in the same way as standard \LaTeX. It should not normally be necessary to use the optional \texttt{[htb]} location specifiers of the \texttt{figure} environment in your manuscript; you may, however, find the \verb"[p]" placement option or the \verb"endfloat" package useful if a journal insists on the need to separate figures from the text.

Figure captions appear below the figures themselves, therefore the \verb"\caption" command should appear after the body of the figure. For example, Figure~\ref{sample-figure} with caption and sub-captions is produced using the following commands:
\begin{verbatim}
\begin{figure}
\centering
\subfloat[An example of an individual figure sub-caption.]{%
\resizebox*{5cm}{!}{\includegraphics{graph1.eps}}}\hspace{5pt}
\subfloat[A slightly shorter sub-caption.]{%
\resizebox*{5cm}{!}{\includegraphics{graph2.eps}}}
\caption{Example of a two-part figure with individual sub-captions
 showing that captions are flush left and justified if greater
 than one line of text.} \label{sample-figure}
\end{figure}
\end{verbatim}
\begin{figure}
\centering
\subfloat[An example of an individual figure sub-caption.]{%
\resizebox*{5cm}{!}{\includegraphics{graph1.eps}}}\hspace{5pt}
\subfloat[A slightly shorter sub-caption.]{%
\resizebox*{5cm}{!}{\includegraphics{graph2.eps}}}
\caption{Example of a two-part figure with individual sub-captions
 showing that captions are flush left and justified if greater
 than one line of text.} \label{sample-figure}
\end{figure}

To ensure that figures are correctly numbered automatically, the \verb"\label" command should be included just after the \verb"\caption" command, or in its argument.

The \verb"\subfloat" command requires \verb"subfig.sty", which is called in the preamble of the \texttt{interacttfssample.tex} file (to allow your choice of an alternative package if preferred) and included in the \textsf{Interact} \LaTeX\ bundle for convenience. Please supply any additional figure macros used with your article in the preamble of your .tex file.

The source files of any figures will be required when the final, revised version of a manuscript is submitted. Authors should ensure that these are suitable (in terms of lettering size, etc.) for the reductions they envisage.

The \texttt{epstopdf} package can be used to incorporate encapsulated PostScript (.eps) illustrations when using PDF\LaTeX, etc. Please provide the original .eps source files rather than the generated PDF images of those illustrations for production purposes.

\subsection{Tables}

The \texttt{interact} class file will deal with positioning your tables in the same way as standard \LaTeX. It should not normally be necessary to use the optional \texttt{[htb]} location specifiers of the \texttt{table} environment in your manuscript; you may, however, find the \verb"[p]" placement option or the \verb"endfloat" package useful if a journal insists on the need to separate tables from the text.

The \texttt{tabular} environment can be used as shown to create tables with single horizontal rules at the head, foot and elsewhere as appropriate. The captions appear above the tables in the \textsf{Interact} style, therefore the \verb"\tbl" command should be used before the body of the table. For example, Table~\ref{sample-table} is produced using the following commands:
\begin{table}
\tbl{Example of a table showing that its caption is as wide as
 the table itself and justified.}
{\begin{tabular}{lcccccc} \toprule
 & \multicolumn{2}{l}{Type} \\ \cmidrule{2-7}
 Class & One & Two & Three & Four & Five & Six \\ \midrule
 Alpha\textsuperscript{a} & A1 & A2 & A3 & A4 & A5 & A6 \\
 Beta & B2 & B2 & B3 & B4 & B5 & B6 \\
 Gamma & C2 & C2 & C3 & C4 & C5 & C6 \\ \bottomrule
\end{tabular}}
\tabnote{\textsuperscript{a}This footnote shows how to include
 footnotes to a table if required.}
\label{sample-table}
\end{table}
\begin{verbatim}
\begin{table}
\tbl{Example of a table showing that its caption is as wide as
 the table itself and justified.}
{\begin{tabular}{lcccccc} \toprule
 & \multicolumn{2}{l}{Type} \\ \cmidrule{2-7}
 Class & One & Two & Three & Four & Five & Six \\ \midrule
 Alpha\textsuperscript{a} & A1 & A2 & A3 & A4 & A5 & A6 \\
 Beta & B2 & B2 & B3 & B4 & B5 & B6 \\
 Gamma & C2 & C2 & C3 & C4 & C5 & C6 \\ \bottomrule
\end{tabular}}
\tabnote{\textsuperscript{a}This footnote shows how to include
 footnotes to a table if required.}
\label{sample-table}
\end{table}
\end{verbatim}

To ensure that tables are correctly numbered automatically, the \verb"\label" command should be included just before \verb"\end{table}".

The \verb"\toprule", \verb"\midrule", \verb"\bottomrule" and \verb"\cmidrule" commands are those used by \verb"booktabs.sty", which is called by the \texttt{interact} class file and included in the \textsf{Interact} \LaTeX\ bundle for convenience. Tables produced using the standard commands of the \texttt{tabular} environment are also compatible with the \texttt{interact} class file.

\subsection{Landscape pages}

If a figure or table is too wide to fit the page it will need to be rotated, along with its caption, through 90$^{\circ}$ anticlockwise. Landscape figures and tables can be produced using the \verb"rotating" package, which is called by the \texttt{interact} class file. The following commands (for example) can be used to produce such pages.
\begin{verbatim}
\setcounter{figure}{1}
\begin{sidewaysfigure}
\centerline{\epsfbox{figname.eps}}
\caption{Example landscape figure caption.}
\label{landfig}
\end{sidewaysfigure}
\end{verbatim}
\begin{verbatim}
\setcounter{table}{1}
\begin{sidewaystable}
 \tbl{Example landscape table caption.}
  {\begin{tabular}{@{}llllcll}
    .
    .
    .
  \end{tabular}}\label{landtab}
\end{sidewaystable}
\end{verbatim}
Before any such float environment, use the \verb"\setcounter" command as above to fix the numbering of the caption (the value of the counter being the number given to the preceding figure or table). Subsequent captions will then be automatically renumbered accordingly. The \verb"\epsfbox" command requires \verb"epsfig.sty", which is called by the \texttt{interact} class file and is also included in the \textsf{Interact} \LaTeX\ bundle for convenience.

Please note that if the \verb"endfloat" package is used, one or both of the commands
\begin{verbatim}
\DeclareDelayedFloatFlavor{sidewaysfigure}{figure}
\DeclareDelayedFloatFlavor{sidewaystable}{table}
\end{verbatim}
will need to be included in the preamble of your .tex file, after the \verb"endfloat" package is loaded, in order to process any landscape figures and/or tables correctly.

\subsection{Theorem-like structures}

A predefined \verb"proof" environment is provided by the \texttt{amsthm} package (which is called by the \texttt{interact} class file), as follows:
\begin{proof}
More recent algorithms for solving the semidefinite programming relaxation are particularly efficient, because they explore the structure of the MAX-CUT problem.
\end{proof}
\noindent This was produced by simply typing:
\begin{verbatim}
\begin{proof}
More recent algorithms for solving the semidefinite programming
relaxation are particularly efficient, because they explore the
structure of the MAX-CUT problem.
\end{proof}
\end{verbatim}
Other theorem-like environments (theorem, definition, remark, etc.) need to be defined as required, e.g.\ using \verb"\newtheorem{theorem}{Theorem}" in the preamble of your .tex file (see the preamble of \verb"interacttfssample.tex" for more examples). You can define the numbering scheme for these structures however suits your article best. Please note that the format of the text in these environments may be changed if necessary to match the style of individual journals by the typesetter during preparation of the proofs.

\subsection{Mathematics}

\subsubsection{Displayed mathematics}

The \texttt{interact} class file will set displayed mathematical formulas centred on the page without equation numbers if you use the \texttt{displaymath} environment or the equivalent \verb"\[...\]" construction. For example, the equation
\[
 \hat{\theta}_{w_i} = \hat{\theta}(s(t,\mathcal{U}_{w_i}))
\]
was typeset using the commands
\begin{verbatim}
\[
 \hat{\theta}_{w_i} = \hat{\theta}(s(t,\mathcal{U}_{w_i}))
\]
\end{verbatim}

For those of your equations that you wish to be automatically numbered sequentially throughout the text for future reference, use the \texttt{equation} environment, e.g.
\begin{equation}
 \hat{\theta}_{w_i} = \hat{\theta}(s(t,\mathcal{U}_{w_i}))
\end{equation}
was typeset using the commands
\begin{verbatim}
\begin{equation}
 \hat{\theta}_{w_i} = \hat{\theta}(s(t,\mathcal{U}_{w_i}))
\end{equation}
\end{verbatim}

Part numbers for sets of equations may be generated using the \texttt{subequations} environment, e.g.
\begin{subequations} \label{subeqnexample}
\begin{equation}
     \varepsilon \rho w_{tt}(s,t) = N[w_{s}(s,t),w_{st}(s,t)]_{s},
     \label{subeqnparta}
\end{equation}
\begin{equation}
     w_{tt}(1,t)+N[w_{s}(1,t),w_{st}(1,t)] = 0,   \label{subeqnpartb}
\end{equation}
\end{subequations}
which was typeset using the commands
\begin{verbatim}
\begin{subequations} \label{subeqnexample}
\begin{equation}
     \varepsilon \rho w_{tt}(s,t) = N[w_{s}(s,t),w_{st}(s,t)]_{s},
     \label{subeqnparta}
\end{equation}
\begin{equation}
     w_{tt}(1,t)+N[w_{s}(1,t),w_{st}(1,t)] = 0,   \label{subeqnpartb}
\end{equation}
\end{subequations}
\end{verbatim}
This is made possible by the \texttt{amsmath} package, which is called by the class file. If you put a \verb"\label" just after the \verb"\begin{subequations}" command, references can be made to the collection of equations, i.e.\ `(\ref{subeqnexample})' in the example above. Or, as the example also shows, you can label and refer to each equation individually -- i.e.\ `(\ref{subeqnparta})' and `(\ref{subeqnpartb})'.

Displayed mathematics should be given end-of-line punctuation appropriate to the running text sentence of which it forms a part, if required.

\subsubsection{Math fonts}

\paragraph{Superscripts and subscripts}
Superscripts and subscripts will automatically come out in the correct size in a math environment (i.e.\ enclosed within \verb"\(...\)" or \verb"$...$" commands in running text, or within \verb"\[...\]" or the \texttt{equation} environment for displayed equations). Sub/superscripts that are physical variables should be italic, whereas those that are labels should be roman (e.g.\ $C_p$, $T_\mathrm{eff}$). If the subscripts or superscripts need to be other than italic, they must be coded individually.

\paragraph{Upright Greek characters and the upright partial derivative sign}
Upright lowercase Greek characters can be obtained by inserting the letter `u' in the control code for the character, e.g.\ \verb"\umu" and \verb"\upi" produce $\umu$ (used, for example, in the symbol for the unit microns -- $\umu\mathrm{m}$) and $\upi$ (the ratio of the circumference of a circle to its diameter). Similarly, the control code for the upright partial derivative $\upartial$ is \verb"\upartial". Bold lowercase as well as uppercase Greek characters can be obtained by \verb"{\bm \gamma}", for example, which gives ${\bm \gamma}$, and \verb"{\bm \Gamma}", which gives ${\bm \Gamma}$.

\section*{Acknowledgement(s)}

An unnumbered section, e.g.\ \verb"\section*{Acknowledgements}", may be used for thanks, etc.\ if required and included \emph{in the non-anonymous version} before any Notes or References.

\section*{Disclosure statement}

An unnumbered section, e.g.\ \verb"\section*{Disclosure statement}", may be used to declare any potential conflict of interest and included \emph{in the non-anonymous version} before any Notes or References, after any Acknowledgements and before any Funding information.

\section*{Funding}

An unnumbered section, e.g.\ \verb"\section*{Funding}", may be used for grant details, etc.\ if required and included \emph{in the non-anonymous version} before any Notes or References.

\section*{Notes on contributor(s)}

An unnumbered section, e.g.\ \verb"\section*{Notes on contributors}", may be included \emph{in the non-anonymous version} if required. A photograph may be added if requested.

\section*{Nomenclature/Notation}

An unnumbered section, e.g.\ \verb"\section*{Nomenclature}" (or \verb"\section*{Notation}"), may be included if required, before any Notes or References.

\section*{Notes}

An unnumbered `Notes' section may be included before the References (if using the \verb"endnotes" package, use the command \verb"\theendnotes" where the notes are to appear, instead of creating a \verb"\section*").

\section{References}

\subsection{References cited in the text}

References should be cited in the text by numbers in square brackets based on the order in which they appear in an alphabetical list of references at the end of the document (not the order of citation), so the first reference cited in the text might be [23]. For example, these may take the forms [32], [5,\,6,\,14], [21--55] (\emph{not} [21]--[55]). For further details on this reference style, see the Instructions for Authors on the Taylor \& Francis website.

Each bibliographic entry has a key, which is assigned by the author and is used to refer to that entry in the text. In this document, the key \verb"Ali95" in the citation form \verb"\cite{Ali95}" produces `\cite{Ali95}', and the keys \verb"Bow76" and \verb"BGP02" in the citation form \verb"\cite{Bow76,BGP02}" produce `\cite{Bow76,BGP02}'. The citation for a range of bibliographic entries (e.g. `\cite{GMW81,Hor96,Con96,Har97,FGK03,Fle80,Ste98,Ell98,Str97,Coo03,Kih01,Hol03,Hai01,Hag03,Hov03,Mul00,GHGsoft,Pow00}') will automatically be produced by \verb"\cite{GMW81,Hor96,Con96,Har97,FGK03,Fle80,Ste98,Ell98,Str97,Coo03,Kih01," \verb"Hol03,Hai01,Hag03,Hov03,Mul00,GHGsoft,Pow00}".
Optional notes may be included at the beginning and/or end of a citation by the use of square brackets, e.g. \verb"\cite[cf.][]{Tay03}" produces `\cite[cf.][]{Tay03}', and \verb"\cite[see][and references therein]{Wei95}" produces `\cite[see][and references therein]{Wei95}'.

\subsection{The list of references}

To produce the list of references, the bibliographic data about each reference item should be listed in \texttt{thebibliography} environment in alphabetical order. References with the same author or group of authors are further sorted chronologically, beginning with the earliest. The following list shows some sample references prepared in Taylor \& Francis' Reference Style S.

\section{Introduction}
\label{sec:intro}
The usage data of mobile-phones is used in a variety of different areas,
such as during the COVID-19 pandemic \cite{Pepe2020, JaysonS2020, Gao2020, Jeffrey2020, Yabe2020,Vollmer2020,Xu2020, JRC2020b, JRC2020, JRC2020a, Heuzroth}, customer segmentation \cite{Aheleroff2011}, identification of personality traits and lifestyle \cite{Chittaranjan2011, Hillebrand2020}, the analysis of large social networks \cite{Aksu2019, Al-Molhem2019, Aledavood2018}, hotspot detection \cite{Nika2016}, prediction of movement \cite{Dao2019}, mode of transport identification \cite{Zhao2020}, credit scoring \cite{Liu2018}, disaster recovery \cite{Andrade, Marzuoli}, analysis of sleeping behavior of the population \cite{Monsivais2017a}, migration \cite{Isaacman2018} and land usage classification \cite{Chen2019a, Lenormand2015}.

The location of a mobile-phone is known for the \gls{mno} of the \gls{gsm} network.
We have partnered with an \gls{mno} in Austria to get access to such anonymized data.
We have defined an aggregation method to understand the overall aggregated mobility of the whole population.
Our data set as well as the aggregation, anonymization approach and the various phases of the lock-down are outlined in detail in \cite{heiler2020covidbasic}.

With the outbreak of COVID-19 and the subsequent lock-down in Austria, the mobility
behavior of the population has changed significantly.
This is reflected in the mobility data derived from mobile-phone information \cite{heiler2020covidbasic}.
To measure mobility, an appropriate measure needs to be established which reflects the mobility.
One possibility is the \gls{rog} \cite{Gooch2011}.
It is formally defined below, and refers to the time-weighted distance of the movement locations to the main location.
We compute it on a daily level.
Its unit is meters, and the values
are strictly positive.
In this work the analyze the aggregated (median) \gls{rog} of the whole population of Austria for various groups as a time series.
The groups are defined by gender- or age groups.

Traditionally, a comparison is made in terms
of \textit{absolute information}, i.e., the \gls{rog} time series
values of the different groups are analyzed in their unit of meters.
We have conducted such an analysis \cite{heilerAndReish2020CovidGender} which focuses on gender differences.

An alternative is to compare \textit{relative information}, for example the
\gls{rog} of the males with respect to females, or in terms
of the ratio males to females.
This leads to a dimensionless time series, and to a different aspect of data analysis which emphasizes the differences between the individual groups.
A joint increase or decrease in both groups may not lead to a big change of the ratio.
On the other hand, the ratio will change if the values of one group increase, and at the same time they decrease in the other group, or vice versa.
Here again, the relative change rather than the absolute change is important.
For example, if the \gls{rog} changes from 1000m to 2000m in one group, and
from 2000m to 1000m in the other group, the ratio would change from 1/2 to 2.
The same change could be observed if the absolute values in both groups would be bigger by a factor of 10.
Thus, absolute values are no longer
relevant in this consideration, because a multiplication by any positive
constant leads to the same ratio.
This is still trivial in case of comparing two groups, but it is no longer
straightforward when relative information of several groups, such as age classes,
should be compared.
\emph{Compositional data analysis} is devoted to this problem
of analyzing relative information \cite{Aitchison86,pawlowsky15,filzmoser2018applied}.
In fact, compositional data analysis is frequently used in geosciences, but also more and more in
other fields such as
biology \cite{espinoza2020applications},
bioinformatics \cite{quinn2018understanding},
economics \cite{trinh2019relations},
marketing \cite{joueid2018marketing},
medicine \cite{dumuid2020compositional}, etc.

We contribute compositional analysis of the movement data during the COVID-19 analysis which comes to the conclusion that special groups (elderly and young cohorts during weekends) need an additional caring treatment to improve their security.

This work is structured as follows.
In Section~\ref{sec:coda} we give a
brief mathematical introduction to compositional data
analysis.
Section~\ref{sec:data} provides more details about the mobile-phone data used and about the quantities derived.
In addition to the mobility measured as the \gls{rog}
we will furthermore investigate the call duration per day, again aggregated by the median.
Section~\ref{sec:results} presents comparisons of the analysis based on absolute and on relative information, and the final Section~\ref{sec:conclusions} summarizes the findings

\section{Materials and methods}
\subsection{Compositional data analysis}
\label{sec:coda}
From the point of view of compositional data analysis, a composition is defined as multivariate information, consisting of strictly positive values, where the absolute numbers as such are not of interest, and only relative
information is relevant for the analysis
\cite{filzmoser2018applied}.
A composition can be given for example by the median \gls{rog} values of different age categories for a certain day, and every age category is denoted as a compositional part.
We use the notation $x_1, \ldots ,x_D$ for the compositional parts of $D$ categories,
and the composition is written as the (column) vector $\mathbf{x}=(x_1,\ldots ,x_D)'$.
For every day recorded in our data base we will observe such a composition, which
in fact leads to a multivariate compositional time series.
The interest is in relative information in terms of the ratios, and thus all pairs $x_j/x_k$, for $j,k= 1,\ldots ,D$, should be considered in the analysis.
Obviously, the pairs for $j=k$ are not relevant, and pairs of the reverse ratio $x_k/x_j$ do not contain potentially new information.
This motivates to consider the logarithm of the ratios, $\ln(x_j/x_k)$, so-called log-ratios.
The reverse ratios have a different sign, and thus do not need to be considered, and their variance is the same as for the original ratio.
Moreover, the distributions of log-ratios tend to be more symmetric than without a logarithm \cite{pawlowsky15}.

Still, the resulting $D(D-1)$ pairs $\ln(x_j/x_k)$, for $k>j$, only live in a subspace of dimension $\leq D-1$ \cite{filzmoser2018applied}, and thus it is natural to aggregate this information.
Consider an aggregation
\begin{equation}
\label{eq:clr1}
y_1=\frac{1}{D}\left(\ln\frac{x_1}{x_2}+ \ldots + \ln \frac{x_1}{x_D}\right)=
\ln\frac{x_1}{g(\mathbf{x})} \ ,
\end{equation}
where
$$
g(\mathbf{x})=\sqrt[D]{\prod_{j=1}^D x_j}
$$
is the geometric mean of the composition $\mathbf{x}$.
Then, $y_1$ represents all relative information about the part $x_1$ to the other
parts in the composition in a form of an average of the log-ratios.
This leads to the definition of so-called \gls{clr} coefficients \cite{Aitchison86}
\begin{equation}
\label{eq:clr}
\mathbf{y}=(y_1,\ldots ,y_D)'
\quad \mbox{ with } \quad
y_j=\ln \frac{x_j}{g(\mathbf{x})} \ .
\end{equation}
The vector $\mathbf{y}$ contains all relative information about $\mathbf{x}$ in the above sense.
It consists of $D$ components $y_j$ which are associated with the relative information about the corresponding part $x_j$.
However, it turns out that $y_1+\ldots +y_D=0$, and thus a representation of data in terms of \gls{clr} coefficients leads to singularity \cite{filzmoser2018applied}.
Although there are ways to circumvent this issue \cite{filzmoser2018applied}, we will proceed with \gls{clr} coefficients for the following analysis for simplicity.

Consider now a multivariate compositional time series
$\mathbf{x}_{t}=(x_{t1},\ldots ,x_{tD})'$, for the time points $t=1,\ldots ,T$,
and the observations $x_{tj}$ for each part $j\in \{1,\ldots ,D\}$.
The time series expressed in \gls{clr} coefficients is
$\mathbf{y}_t=(y_{t1},\ldots ,y_{tD})'$, with $y_{tj}=\ln(x_{tj}/g(\mathbf{x}_t))$,
with the geometric mean $g(\mathbf{x}_t)=(\prod_{j=1}^D x_{tj})^{1/D}$ per time point.
Since this data representation only reflects relative information of the time series, an additional visualization of the absolute time series values can be interesting
to get a more complete picture.

The \gls{clr} coefficients result in multivariate data that can be analyzed with the traditional multivariate statistical methods \cite{filzmoser2018applied}.
A prominent way to represent the information in a lower-dimensional space is
to use \gls{pca}.
Since PCA is sensitive to data outliers or inhomogeneous data, robust versions have been proposed, also in the compositional data analysis framework \cite{filzmoser2018applied}.
The resulting loadings and scores are commonly represented in a biplot
    to get an overview of the multivariate data \cite{AG02}.

\subsection{Mobile-phone data}
\label{sec:data}

In this work we analyze two measures obtained from the mobile phone data, the call duration and the radius of gyration \gls{rog}.
While the meaning of the former is straightforward, the latter needs
to be defined.

Consider an individual $i:=i(t)$ at a certain day
$t\in\{1,\ldots ,T\}$.
For reasons of data privacy, the index of the individual will change every day.
The  data is made available to the researchers already anonymized with a daily changing key.

Furthermore, the current location of the individual's mobile phone is available at the time points $t_{\tau}=t+\tau_t$, for a number of  time points $T_t$ per day, where $\tau_t \in[0,1)$.
The corresponding $x$- and $y$-coordinates are denoted by $(\xi_{it_{\tau}}, \eta_{it_{\tau}})$.
With this information, the stay duration $l_{it_{\tau}}$ for individual $i$ at time point $t_{\tau}$ can be computed, which is used to calculate a weighted average
$(\bar{\xi}_{it},\bar{\eta}_{it})=
\left(
\frac{\sum_{t_{\tau}}l_{it_{\tau}}\xi_{it_{\tau}}}
{\sum_{t_{\tau}}l_{it_{\tau}}},
\frac{\sum_{t_{\tau}}l_{it_{\tau}}\eta_{it_{\tau}}}
{\sum_{t_{\tau}}l_{it_{\tau}}}
\right)$
for individual $i$ for day $t$.
These coordinates are in the middle of the area covered by  all the locations which were visited during the day and are dominated by the two most prominently used (longest used) locations: home and work location.

Denote $d^2_{it_{\tau}}=(\bar{\xi}_{it}-\xi_{it_{\tau}})^2
+(\bar{\eta}_{it}-\eta_{it_{\tau}})^2$ as the squared Euclidean distance
between the coordinates $(\bar{\xi}_{it},\bar{\eta}_{it})$ and
$(\xi_{it_{\tau}}, \eta_{it_{\tau}})$.
This requires the coordinate system to be local in order to obtain valid, i.e. less distorted results.
Otherwise, a Haversine distance could be used instead in case of epsg:4326 WGS-84 projection of the coordinates.

The \gls{rog} for individual $i$
and day $t$ is then defined as
\begin{equation}
R_{it} = \sqrt{\frac{\sum_{t_{\tau}} d^2_{it_{\tau}}}
{\sum_{t_{\tau}} l_{it_{\tau}}}} ,
\label{eq:rog}
\end{equation}
and it thus represents a distance to the center of all the places of stay during that day $t$ weighted by the lengths of the stay duration at the different places.

Details are available and especially a description of how the large quantity of data was handled is available in \cite{heiler2020covidbasic}.

Using additional metadata, an individual $i$ can be assigned to a gender
group (female, male), to an age group (here we consider the age groups in 15 year intervals:
15-29, 30-44, 45-59, 60-74, and 75$+$), and to an Austrian district of the daily night location to derive the groups.
Since the distributions are generally very right-skewed, we work with the median per group and day in the following and also ensure k-anonymity for each one.

The resulting time series can be directly investigated in terms of  their absolute information, and they can be compared to an analysis based on relative information.

\section{Results}
\label{sec:results}

\subsection{Mobility measured by \gls{rog}}
The results reported in this section refer to the median values of
the \gls{rog} per group.
To begin with, Figure~\ref{fig:plot_r50abs} shows the absolute values for the females (top) and males (bottom) for different age groups.
The legend indicates the considered age groups: 15 for age 15-29, 30 for age 30-44, 45 for age 45-59, 60 for age 60-74, and 75 for age elder than 75. For all of the following time series  plots, the vertical dashed lines indicate the date March 16\textsuperscript{th}, 2020, when the restrictions came into action, and the date April 6\textsuperscript{th}, 2020, when they were relaxed.
The data considered here are from the period February 1\textsuperscript{st} until August 9\textsuperscript{th}, 2020.
The lock-down is clearly visible in the plots by an abrupt decay of the median \gls{rog} values in all age classes for both genders. After the lock-down, the order of the values still remains the same, from the eldest group with the smallest values, and the youngest group with the highest values, but it is on a much smaller level. The level then increased more or less systematically until the middle of June. Afterwards, the level
is not changing a lot, it is lower than at the beginning, and weekly patterns are clearly visible.
Note that these weekly time series patterns that are very regular at the beginning are getting somehow distorted, partially also due to holidays (April 13\textsuperscript{th}, May 1\textsuperscript{st}, May 21\textsuperscript{th}, June 1\textsuperscript{st}, June 11\textsuperscript{th}), and they never get back to this regularity.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/plot_r50absaug09.pdf}}
\caption{Median \gls{rog} values for different age groups over time
for females (top) and males (bottom) in different age groups.}
\label{fig:plot_r50abs}
\end{figure}

Figure~\ref{fig:plot_r50rel} focuses on the relative information
contained in the median \gls{rog} values.
We consider the female age groups and the male age groups separately as two compositions.
The plots show the corresponding \gls{clr} coefficients for females (top) and males (bottom). While in Figure~\ref{fig:plot_r50abs} we have essentially
seen a decline of all values at the beginning of the lock-down phase,
followed by an increase, we did not pay attention how differently
the age groups declined and increased. This is the purpose of
the relative view in Figure~~\ref{fig:plot_r50rel}, where we mainly
investigate the developments of the age groups to each other.

In both plots of Figure~~\ref{fig:plot_r50rel} we can see roughly the same pattern after the lock-down: the biggest relative changes are visible for the youngest and the oldest age group, but they go into different directions.
While group 15 had the biggest decline, group 75$+$ increased the values relative to the other age groups.
This seems to be counter-intuitive, but it can be explained by the fact that the geometric mean also went down significantly, and the ratio of the values of group 75$+$ to the geometric mean then even increased after the lock-down.
Another interesting phenomenon is that the groups 60 and 75$+$ show the biggest increase in mobility (in a relative sense) during the weekends in this lock-down period. Although on a different level, the values from July
show a similar structure to those from February. It is interesting to note
that the youngest age group 15 shows a somehow mirrored weekly pattern
compared to the elder age groups. This is not visible when looking at the
absolute values in Figure~\ref{fig:plot_r50abs}.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/plot_r50relaug09.pdf}}
\caption{\gls{clr} coefficients of median \gls{rog} values for the female (top) and the male (bottom) composition.
For the legend see Figure~\ref{fig:plot_r50abs}.}
\label{fig:plot_r50rel}
\end{figure}

Relative information could also be understood in terms of data proportions.
In particular, one could compute the proportion of a group on the
total per time point, which in fact corresponds to normalizing the
data per time point to a value of 1.
Such a proportional presentation
is shown in Figure~\ref{fig:plot_r50prop} for the \gls{rog}
values of the female age groups.
Obviously, the information contained in
this representation is different from \gls{clr} coefficients which focus on
log-ratio information.
One can hardly see any differences between the
lock-down period and the remaining period, and thus this kind of
``relative view'' is not valuable for the analysis.
\begin{figure}[htbp]
\centerline{\includegraphics[width=\linewidth]{img/plot_r50propaug09.pdf}}
\caption{Proportional presentation of the median \gls{rog} values for the
female age groups. For each time point, the data are normalized to a value
of 1.}
\label{fig:plot_r50prop}
\end{figure}

The median \gls{rog} values for the female and male age groups are analyzed
in the following with PCA. Here, the method ROBPCA \cite{hubert-ROBPCA:2005}
is taken, a robust version of PCA which downweights outlying observations.
Figure~\ref{fig:biplot_r50} shows the biplot of the first two principal components (PCs)
for the clr coefficients. The coloring is according to the time phases:
green before the lock-down, pink during the lock-down period, purple
after lock-down until mid of June, and
light-blue after this period. The left biplot for the females
identifies these four periods as clear clusters, while there is more
overlap visible in the right biplot for the males.
For the females, the direction of the first PC (71\% explained
variance) shows a transition
of the relative \gls{rog} values from the young generation (f15, f30)
before lock-down to the old (f75) one during lock-down, and then back to the
center. Thus, younger and elder females show a contrasting behavior in this
time period, which was already observed in Figure~~\ref{fig:plot_r50rel} (top panel).
The second PC (21\% explained variance) shows also differences between
the time periods, but it also reveals weekend effects. Especially
on Sundays, the mobility for group f15 was
bigger before and after lock-down, but it moved to group f75 during the
lock-down phase.

The data structure in the biplot for the males (right plot) looks a bit
different, but leads to similar conclusions. PC1
explains 69\% and PC2
25\% of the variance. Groups m75 and m15 have a similarly diverging
behavior of Sunday mobility as observed for the females.
The weekdays of the lock-down phase are in the center of the distribution,
while for the females they were clearly moved towards group f75.
On the other hand, the weekdays in the first time period (February 1\textsuperscript{st} - March 15\textsuperscript{th})
are better distinguishable from the weekdays of the last period (June 15 -
August 9); a possible explanation is the fact that the working male
population changed the mobility behavior more significantly than
that of females due to home office.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/biplot_r50aug09.pdf}}
\caption{Biplots of the \gls{clr} coefficients of the median \gls{rog} values
for females (left) and male (right) age groups. Green color for period before the lock-down, pink for lock-down period, purple after lock-down until mid of June, and light-blue after this period.}
    \label{fig:biplot_r50}
\end{figure}

A quite contrasting view is revealed in Figure~\ref{fig:biplot_r50abs}, which
shows the robust PCA results for the absolute
values of \gls{rog}, for females (left) and males (right). In both
analyses, PC1 explains 98\% of the variability, and this direction
essentially reflects the big change of the \gls{rog} over this time
period. Otherwise, there is not much information left in these analyses,
which reflects the limited usefulness of absolute information if the task
is to compare age groups.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/biplot_r50absaug09.pdf}}
\caption{Biplots of the (absolute) median \gls{rog} values
for females (left) and male (right) age groups. Green color for period before the lock-down, pink for lock-down period, purple after lock-down until mid of June, and light-blue after this period.
}
    \label{fig:biplot_r50abs}
\end{figure}

\subsection{Interaction measured by call duration}

Figure~\ref{fig:plot_mo50} investigates the median call duration, reported
in seconds, again for the two genders and the age groups. The absolute values
are shown in the upper plot jointly for males and females. Here we observe the
reverse ordering of the age groups compared to the plots for the \gls{rog}
values: the lowest values are for the youngest group, and the biggest for the
oldest group. The values of the females are systematically higher than those
of the males. It is interesting to see that the call durations already started
to increase one week before the lock-down. While the \gls{rog} time series had
their peaks during the weekend, we have the opposite here. This pattern, however,
seems to change after the lock-down for group f75 (uppermost line),
and it went back to \emph{normality} only later on.

The bottom plot of Figure~\ref{fig:plot_mo50} presents the \gls{clr} coefficients, which
are separately calculated for females and males, but presented here jointly for
easy comparison. Although the absolute values of the youngest age group
also increased with the lock-down, the increase was smaller compared to the
other groups, which is reflected by decreasing \gls{clr} coefficients. The pattern
of f15 and m15 has also an interesting structure: Before the lock-down,
the groups had quite different behavior within their gender-group, but during
the lock-down phase they became quite similar. From June on, they show again a
similar behavior as at the beginning. Another interesting phenomenon can be
seen after the lock-down: the two oldest groups show a contrary behavior to
the other groups during the weekends. Their decline in call duration during
the weekends was much smaller than that of the other age groups.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/plot_mo50aug09.pdf}}
\caption{Median values of call durations per gender and age group over time (top),
and clr representations separately for female and male age groups (bottom).}
\label{fig:plot_mo50}
\end{figure}

Figure~\ref{fig:biplot_mo50} presents biplots of a robust PCA for the \gls{clr}
coefficients for the female (left) and male (right) age groups. The coloring
is taken as in the previous biplots, green before lock-down, pink during,
purple after lock-down, and light-blue from June 15\textsuperscript{th} onwards.
PC1 explains 72\% of the variability for the
females, 54\% for males, and PC1 and PC2 together explain about 98\% variance
in both cases. The different groups which are visible in the biplots are
essentially weekend-effects or affects due to the lock-down.
These grouping effects are essentially caused by the youngest and
oldest age groups.
When comparing the first observed time period with the last one,
we can find quite clear differences in the corresponding PCA scores.
These differences are essentially caused by the changing contrasting
behavior between the youngest group and the elder groups; groups
f75 and m75 (and also m30) do not seem to contribute to this difference.
A possible explanation is the exploration of alternative methods for
communication, especially for the elder groups.
\begin{figure
\centerline{\includegraphics[width=\linewidth]{img/biplot_mo50aug09.pdf}}
\caption{Biplots of the \gls{clr} coefficients of the median call duration values
for females (left) and male (right) age groups. Green color for the period before the lock-down, pink for lock-down period, light-blue after this period.}
    \label{fig:biplot_mo50}
\end{figure}

\subsection{Interactions between source and destination}

It can be recorded who is actively calling a person, and who is receiving a call.
The former person is called \emph{source}, and the latter \emph{destination}.
Here we investigate the median \gls{rog} values for the different age groups of the females and males.
However, the data set is more complex now, because a person from a certain age group can be source, while the destination can originate from a different age group.
Moreover, both source and destination will have
specific median \gls{rog} values.

Figure~\ref{fig:plot_45_75} illustrates these data for four specific cases:
source f45 (f45$\_$src) with destination f75 (f75$\_$dst), and source
f75 (f75$\_$src) with destination f45 (f45$\_$dst). In both cases, the
median \gls{rog} values can be taken from the source group or from the
destination group, see also
figure legend. Throughout the whole period (here from February 1\textsuperscript{st} - July 26\textsuperscript{th}), the median \gls{rog} values from
the source groups (solid lines) have slightly higher
values than those of the destination groups (dashed lines) for the same
age classes, which can be expected because people from the source groups might
call from a place outside their usual environment. While the lines
are on a similar level at the beginning and at the end of the considered
period, the weekly periodicity changes, probably caused by the
summer holidays.
\begin{figure}[htbp]
\centerline{\includegraphics[width=\linewidth]{img/plot_45_75aug09.pdf}}
\caption{Median \gls{rog} values for the female age groups $f45$ and $f75$,
depending whether they actively call (src) or they passively receive
the call (dst). For example, line f75$\_$src -- f45$\_$dst refers to the
median \gls{rog} values for females in age group f75, actively calling
females in age group f45.}
    \label{fig:plot_45_75}
\end{figure}

In the following analyses we are interested in the similarity of the relative
\gls{rog} values in terms of correlations, before
lock-down (February 1\textsuperscript{st} -- March 15\textsuperscript{th}) and after (March 16\textsuperscript{th} -- May 31\textsuperscript{th}). In order to
investigate relative information, the \gls{clr} coefficients are computed for
a composition with all 25 age combinations of the source-destination groups
and all 25 age combinations of the destination-source groups, separately
for females and males.
Figure~\ref{fig:plotcorf} shows the resulting
correlation matrix for the females as a heat map, left for time points before the
lock-down, and right after lock-down.
The row and column labels are referring
to the group numbers.
For example, src1-3 refers to the time series
f15$\_$src -- f30$\_$dst, or dst5-1 is the series f75$\_$dst -- f15$\_$src.
The heatmaps show that the correlation structure before and after lock-down
has clearly changed.
Afterwards, there are more blocks with higher (absolute)
correlations, and thus more similarity or dissimilarity between certain
age groups.
In general, there is a more pronounced difference after lock-down
in the mobility behavior between the younger and the elder age groups.
\begin{figure
\centerline{\includegraphics[width=0.48\linewidth]{img/plotcorf1.pdf}
\includegraphics[width=0.48\linewidth]{img/plotcorf2.pdf}}
\caption{Correlations of the \gls{clr} coefficients for median \gls{rog} values
for the female age groups (1 to 5, referring to 15 to 75$+$),
when they are actively calling (src) or passively receiving a call (dst),
recorded before March 16\textsuperscript{th}, 2020 (left), and afterwards (right).}
    \label{fig:plotcorf}
\end{figure}

\subsection{Incorporating spatial location}
The mobile phone data also contain information about the location, in
our case about the Austrian political district in which the phone has been used.
The Austrian regions had different restrictions during the
lock-down phase, and in particular people from all districts in Tirol
had the strongest movement restrictions.
Thus, in Figure~\ref{fig:plot_f50AB} we compare the median
\gls{rog} values for Kitzb\"uhel, a district in Tirol, and Zell am See,
which is also a rural district but located in Salzburg. The absolute values
of the female age groups are shown in the upper plots, while the \gls{clr} coefficients are presented in the
lower plots.
Since the same scale is used along the vertical axes, one can clearly see the
difference in mobility during the lock-down period in Kitzb\"uhel and
Zell am See, and this is also visible in the \gls{clr} coefficients.
For Kitzb\"uhel, there is much smaller variability of the values
during lock-down, and also the relative differences between the age groups
become much smaller.
The change in the relative differences is not so
pronounced for Zell am See.
This means that also from a relative point of view,
the data structure changes completely in Kitzb\"uhel due to the restrictions.
\begin{figure
\includegraphics[width=\linewidth]{img/plot_f50ABaug09.pdf}
\caption{Median \gls{rog} values for Kitzb\"uhel (left) and Zell am See (right)
as absolute (top) and relative (bottom) information.}
    \label{fig:plot_f50AB}
\end{figure}

Figure~\ref{fig:plotm30AB} focuses on the male age group m30, and compares the
composition of all districts in Tirol with that of all districts in Salzburg.
The dashed lines refer to the district capitals (Innsbruck and Salzburg,
respectively).
These districts behave differently compared to the other districts
which are rural with many people commuting to their work place.
The values of the districts in Tirol (except Innsbruck) get closer to each other
after lock-down, and they start to diverge only in the middle of April.
This may be explained by a similar mobility behavior of the m30 group within
this period, probably caused by home-office or reduced working time. This seems
different in districts of Salzburg, where the \gls{clr} coefficients show more
variability after lock-down.
\begin{figure
\includegraphics[width=\linewidth]{img/plot_m30ABaug09.pdf}
\caption{\gls{clr} coefficients of median \gls{rog} values for the male group m30 in
all districts of Tirol
(top) and Salzburg (bottom).
The capitals are shown as dashed lines.}
    \label{fig:plotm30AB}
\end{figure}

\section{Discussion}
\label{sec:conclusions}

By analyzing relative changes using compositional data analysis methodologies formerly hidden insights can be identified.
In this work of analyzing mobility during the COVID-19 lock-down measures we see that certain age-groups of the population (elderly, young during weekends) do restrict mobility less than other members of the population.
Especially for the elderly which are at high risk of infections potentially additional reminders should be sent to adhere to the interventions.
Similarly, for the young groups on weekends additional reminders to use mouth nose protection could be useful.

\section*{Acknowledgements}
This work was funded by the WWTF under project COV COV20-035.

\section*{Declaration of interest}
We have no competing interests to declare.
\printbibliography

\end{document}

Physics and Society (physics.soc-ph); Applications (stat.AP); Computation (stat.CO)


Magnetic anisotropy switching in (Ga,Mn)As with increasing hole concentration

Abstract: We study a possible mechanism of the switching of the magnetic easy axis as a function of hole concentration in (Ga,Mn)As epilayers. In-plane uniaxial magnetic anisotropy along [110] is found to exceed intrinsic cubic magnetocrystalline anisotropy above a hole concentration of p = 1.5 * 10^21 cm^-3 at 4 K. This anisotropy switching can also be realized by post-growth annealing, and the temperature-dependent ac susceptibility is significantly changed with increasing annealing time. On the basis of our recent scenario [Phys. Rev. Lett. 94, 147203 (2005); Phys. Rev. B 73, 155204 (2006).], we deduce that the growth of highly hole-concentrated cluster regions with [110] uniaxial anisotropy is likely the predominant cause of the enhancement in [110] uniaxial anisotropy at the high hole concentration regime. We can clearly rule out anisotropic lattice strain as a possible origin of the switching of the magnetic anisotropy.

\section{INTRODUCTION}

Ferromagnetism in III-V magnetic alloy semiconductors has been studied in many theoretical and experimental approaches,\cite{Ohno2,Matsukura,Dietl2,Abolfath,Keavney} invoked to $p-d$ exchange interaction between hole carriers and magnetic impurities. One of the most fascinating magnetic features in (Ga,Mn)As epilayers is the temperature variation of the magnetic anisotropy from $\left\langle 100 \right\rangle$ to [110] with increasing temperature,\cite{Welp,Liu,Hamaya2,Sawicki} where the origin of [110] uniaxial anisotropy remains still unclear. Also, the mechanism of the anisotropy switching from $\left\langle 100 \right\rangle$ to [110] with increasing temperature had not been elucidated, although theoretical works\cite{Dietl2,Abolfath} have made an attempt to correlate the magnetic anisotropy with the shape of the Fermi surface of valence hole subbands.

In regard to this point, we have recently proposed another two-phase scenario to explain the magnetic properties, where two ferromagnetic phases of [110] uniaxial anisotropic clusters and $\left\langle 100 \right\rangle$ matrix coexist in a (Ga,Mn)As epilayer:\cite{HamayaPRL,HamayaPRB} [110] uniaxial anisotropic clusters have relatively high hole concentration in the epilayer. Using this model, we can clearly explain the origin of the temperature-dependent anisotropy switching: with increasing temperature, the ferromagnetic $\left\langle 100 \right\rangle$ matrix phase becomes paramagnetic, and, the [110] uniaxial anisotropic clusters dominate the magnetism of (Ga,Mn)As epilayers near the Curie temperature.\cite{HamayaPRL,HamayaPRB} This cluster/matrix scenario is based on the fact that the low-temperature grown (Ga,Mn)As epilayers have point defects such as Mn interstitials (Mn$_\mathrm{I}$), \cite{Yu,Blinowski} and the Mn$_\mathrm{I}$ donors are diffused by post-growth low-temperature annealing.\cite{Ku,Edmonds} Also, the presence of energetically favorable Mn clusters, e.g., Mn$_\mathrm{Ga}$-Mn$_\mathrm{Ga}$ and Mn$_\mathrm{Ga}$-Mn$_\mathrm{I}$-Mn$_\mathrm{Ga}$, is proposed by theoretical studies \cite{Mahadevan,Raebiger} and direct observation using scanning tunneling microscopy.\cite{Sullivan} Recently, Sato {\it et al.} theoretically suggested the presence of spinodal decomposition phases in magnetic semiconductors such as (Ga,Mn)As and (Ga,Mn)N, where the phases bring about a high Curie temperature.\cite{Sato}

On the other hand, the effect of hole concentration ($p$) on the magnetic anisotropy has also been examined experimentally. For example, it has been reported that $\left\langle 100 \right\rangle$ magnetocrystalline anisotropy can be controlled by changing $p$ in micro-structured (Ga,Mn)As wires\cite{HamayaJJAP} and modulation-doped (Ga,Mn)As/(Ga,Al)As heterostructures,\cite{Liu2} indicating that $p$ is a predominant parameter in determining the cubic magnetocrystalline anisotropy. In terms of in-plane uniaxial anisotropy along [110], we have demonstrated that the [110] uniaxial anisotropy is enhanced in an annealed (Ga,Mn)As epilayer compared with the corresponding as-grown epilayer.\cite{Kato} Also, Sawicki {\it et al.}\cite{Sawicki2} suggested that the in-plane uniaxial anisotropy depends on $p$ using the $p-d$ Zener model on the assumption of a small trigonal distortion. Since such trigonal distortion is likely to be attributed to the distribution of Mn ions (Mn clusters), this assumption seems to be compatible with our scenario: the Mn clusters cause inhomogeneity of $p$ in an epilayer and the magnetic anisotropy. To obtain a common view about the $p$ dependence of the magnetic anisotropy, we hence need to reconsider the magnetic anisotropy problem as a function of $p$ on the basis of the cluster scenario.\cite{HamayaPRL,HamayaPRB}
\begin{figure}[t]
\includegraphics[width=8cm]{Fig1.eps}
\caption{$M-H$ curves of (Ga,Mn)As epilayers at 4 K. Magnetic field is applied parallel to either
$\left\langle 100 \right\rangle$, [110], or [1$\overline{1}$0]. Hole concentrations ($p$) of three samples are
(a) 3$\times$$10^{20}$ cm$^{-3}$, (b) 7$\times$$10^{20}$ cm$^{-3}$, and (c) 1.5$\times$$10^{21}$ cm$^{-3}$.}
\end{figure}

In this work, we show a clear manifestation of hole-concentration-dependent cubic and uniaxial anisotropy in (Ga,Mn)As epilayers, and find switching of the magnetic easy axis from $\left\langle 100 \right\rangle$ to [110] in the high hole concentration regime ($p\gtrsim$10$^{21}$ cm$^{-3}$) at 4 K. Comparison between ac susceptibility data of two separate epilayers with $p\sim$10$^{20}$ cm$^{-3}$ and $p\sim$10$^{21}$ cm$^{-3}$ reveals that [110] magnetic anisotropy arises from highly hole-concentrated (Ga,Mn)As clusters in the (Ga,Mn)As epilayer and such clusters are becoming predominant with increasing hole concentration.

\section{SAMPLES AND EXPERIMENT}
(Ga,Mn)As epilayers with a thickness of 100 or 200 nm were grown on a semi-insulating GaAs (001) substrate
using low-temperature molecular beam epitaxy at 190 $-$ 240$^{\circ}$C. Prior to growing the (Ga,Mn)As, a 400-nm-thick GaAs buffer layer was grown in an As-stabilized condition at 590$^{\circ}$C. The Mn contents $x$ in the Ga$_{1-x}$Mn$_{x}$As epilayers were estimated to range from 0.024 to 0.088. After growth, some of the epilayers were annealed at 250$^{\circ}$C.\cite{Hayashi,Ku,Edmonds} Structural characterizations were performed by high resolution x-ray diffraction (HRXRD) using a Philips MRD. Magnetization was measured with a superconducting quantum interference device (SQUID)
magnetometer. Ac susceptibility measurements were carried out using a physical property measurement system (PPMS).
 The hole carrier concentrations ($p$) which correspond to the ionized Mn acceptor concentration were measured with an electrochemical capacitance-voltage method at room temperature.\cite{Moriya} The samples we use show typical hole-mediated ferromagnetic properties; Curie temperature ($T_{c}$) and hole concentration ($p$) range from 55 K to 135 K and from 3$\times$$10^{20}$ cm$^{-3}$ to 1.5$\times$$10^{21}$ cm$^{-3}$, respectively, being  similar to those used in previous work reported by Ku {\it et al}.\cite{Ku}

\section{RESULTS AND DISCUSSION}
Figure 1 shows the field-dependent magnetization ($M-H$ curve) of (Ga,Mn)As epilayers
with different hole concentrations of (a) $p=$ 3$\times$$10^{20}$ cm$^{-3}$, (b) $p=$ 7$\times$$10^{20}$ cm$^{-3}$, and (c)  $p=$ 1.5$\times$$10^{21}$ cm$^{-3}$ at various field orientations at 4 K.
For sample (a), no differences in the $M-H$ curves are seen for field directions along [110] and [1$\overline{1}$0] (Fig. 1(a)). In contrast, for samples (b) and (c), the differences between [110] and [1$\overline{1}$0] become clearer with increasing hole concentration (Figs. 1(b) and (c)), which is evidence of [110] uniaxial anisotropy. Also note that the remanent magnetization for [110] is larger than those for $\left\langle 100 \right\rangle$ in Fig. 1(c), indicating that the magnetic easy axis alters from [100] to a direction near [110] in the high hole concentration regime $p\gtrsim$$10^{21}$ cm$^{-3}$.

To deduce the magnetic anisotropy constants, we assume that the $M-H$ curves along the hard axis can be modeled on coherent rotation.\cite{Welp} The magnetostatic energy of in-plane magnetized (Ga,Mn)As is generally given by $E =$ $K_{u}$sin$^{2}$($\varphi - $45$^\circ$) $+$ ($K_{1}$/4)sin$^{2}2\varphi$ $-$ $MH$cos($\varphi - \theta$), where {\it K$_{u}$} and {\it K$_{1}$} are the in-plane uniaxial and cubic anisotropy constants, respectively, {\it M} is the magnetization, {\it H} is the applied field strength, $\varphi$ is the magnetization direction, and $\theta$ is the applied field direction with respect to GaAs[100].\cite{Welp,Hrabovsky,Hamaya2}
We fit a calculated $M-H$ curve in the field region from 0 Oe to 2 kOe using the conditions $\partial${\it E}/$\partial$$\varphi$ $=$ 0, $\partial$$^2${\it E}/$\partial$$\varphi$$^2$ $>$ 0,
and $\theta =$ 135$^{\circ}$ ([1$\overline{1}$0]) in this formula.\cite{note}
Figure 2 (a) plots {\it K$_{u}$} and {\it K$_{1}$} at 4 K as a function of $p$.
{\it K$_{1}$} decreases from 3$\times$10$^{4}$ to 6$\times$10$^{3}$ (erg/cm$^{3}$) with increasing $p$. In contrast, {\it K$_{u}$} gradually increases from 1$\times$10$^{3}$ to 8$\times$10$^{3}$ (erg/cm$^{3}$) and becomes larger than {\it K$_{1}$} for a sample with $p =$ 1.5$\times$$10^{21}$ cm$^{-3}$, showing significant [110] uniaxial anisotropy. Our experimental data in Fig. 2 (a) clearly show a crossover of the values of {\it K$_{u}$} and {\it K$_{1}$} in the high hole concentration regime $p\gtrsim$$10^{21}$ cm$^{-3}$.
\begin{figure}[t]
\includegraphics[width=8cm]{Fig2.eps}
\caption{(a) {\it K$_{u}$} and {\it K$_{1}$} estimated from $M-H$ curves at 4 K as a function of room-temperature hole concentration. (b) {\it K$_{1}$ } $vs$ Mn content at 4 K. Different symbols are used for as-grown and annealed samples. }
\end{figure}

Although $p$ seems to be the most appropriate parameter determining the magnetic anisotropy, another cause such as Mn content and lattice strain should be checked whether they could represent the experimental behavior. Here, Mn content is estimated from HRXRD measurements.\cite{Shen} As shown in Fig. 2 (b), it is clearly seen that the plots of {\it K$_{1}$} vs Mn content are rather scattered and show no systematic changes as a function of Mn content: even samples with the same Mn content show different {\it K$_{1}$} values depending on whether the sample is as-grown or annealed. Also, the magnetic anisotropy of (Ga,Mn)As was reported to be affected by the shape of valence hole subbands and lattice strain could be a possible cause of the change in the shape of the valence hole subbands.\cite{Dietl2,Abolfath} We have confirmed that all the epilayers used were grown coherently with tetragonal compressive strain and no misfit dislocation, and that no in-plane crystal asymmetry between [110] and [1$\overline{1}$0] is induced by increasing Mn contents or low temperature annealing, ensuring that the influence of the asymmetric lattice strain on the change in the values of {\it K$_{u}$} can be ignored. Welp {\it et al}. also showed that [110] uniaxial anisotropy field did not depend on epilayer thickness and the effects of (Ga,Mn)As surface and (Ga,Mn)As/GaAs(001) interface on [110] uniaxial anisotropy are not significant.\cite{Welp2} These experimental facts support our description that $p$ is most likely the predominant parameter in determining the magnetic anisotropy of (Ga,Mn)As rather than other parameters.
\begin{figure}
\includegraphics[width=8.5cm]{Fig3.eps}
\caption{Zero-field cooled and field cooled magnetizations for samples with various annealing time ((a) as-grown, (b) 1-hour annealed, (c) 16-hour annealed).}
\end{figure}

To gain additional insight into the origin of the change in the magnetic anisotropy with $p$, we examine the magnetic features of as-grown and low-temperature annealed samples (1-hour annealed and 16-hour annealed samples) with a constant Mn content. For the sample with a Mn content of 4.4\%, the [110] uniaxial anisotropy at 4 K and the $p$ were enhanced with increasing annealing time, being consistent with previous studies.\cite{Kato,Sawicki2} Also, no anisotropy switching with increasing temperature was observed for 16-hour annealed sample because of the marked [110] uniaxial anisotropy.\cite{HamayaPRB} For these samples, we measure the zero-field cooled (ZFC) and field cooled (FC) magnetizations. Figure 3 presents the ZFC and FC magnetizations of (a) as-grown sample, (b) 1-hour annealed sample, and (c) 16-hour annealed sample at $h_{dc}$ $=$ 20 Oe along [100]. Clear bifurcations of the ZFC and FC magnetizations are seen for all the samples. If there are no inhomogeneous magnetic phases and the magnetic structure consists of a single-domain as pointed out by Wang {\it et al.},\cite{Wang} we cannot expect the bifurcation of the ZFC and FC magnetizations for the samples.\cite{HamayaPRB} Thus, these bifurcations strongly support the presence of inhomogeneous magnetic phases in (Ga,Mn)As epilayers.\cite{Mahadevan,Raebiger,Sullivan,HamayaPRL,HamayaPRB,Sato}

For our samples, since the single-domain model presented by Wang {\it et al.}\cite{Wang} is not realistic, we hereafter discuss the cluster/matrix model in our recent reports:\cite{HamayaPRL,HamayaPRB} (Ga,Mn)As epilayers that show two different $\left\langle 100 \right\rangle$ cubic anisotropy and [110] uniaxial anisotropy are composed of a matrix phase with relatively low hole concentration and highly-hole concentrated (Ga,Mn)As clusters embedded in the matrix. In this model, the distribution in $p$ is likely due to the inhomogeneous distribution of Mn ions.\cite{Mahadevan,Raebiger,Sullivan,HamayaPRL,HamayaPRB,Sato} In Fig. 4, we also show the temperature-dependent ac susceptibility ($\chi_{2}$"$- T$) curve for the same samples as in Fig. 3. An ac magnetic field of 6 Oe is small enough because features seen in the higher temperature regime are compatible with the magnetization data measured in Fig. 3: the high-temperature features in Fig. 4 are seen at around the temperature where the magnetization disappears in Fig. 3. As shown in our previous study,\cite{HamayaPRL} two peak-to-dip features at $\sim$38 K and $\sim$57 K in the $\chi_{2}$"$- T$ curve are seen for the as-grown sample, arising from the matrix and the clusters with different Curie temperatures.\cite{HamayaPRL,HamayaPRB} Hence, the feature seen in the lower temperature regime ($\sim$38 K) is associated with the Curie temperature of the matrix.\cite{HamayaPRL} For annealed samples, significant changes in the $\chi_{2}$" are found: the sample subjected to 1-hour annealing shows a shift of the low-temperature feature at 38 K towards further lower temperatures, while the high-temperature feature at 57 K shifts towards higher temperatures, causing multiple-dip feature in the $\chi_{2}$"$- T$ curve. After 16-hours annealing, the high-temperature feature in the $\chi_{2}$"$-T$ further changes and shifts towards higher temperature, and the low-temperature feature almost disappears.
\begin{figure}
\includegraphics[width=8.5cm]{Fig4.eps}
\caption{Temperature dependent nonlinear ac susceptibility of the same samples used in dc magnetization measurements (Fig.3) at various ac frequencies.}
\end{figure}

We discuss the above systematic changes in the $\chi_{2}$" $- T$ curves, together with recent studies based on the presence of Mn clusters. Mahadevan {\it et al.}\cite{Mahadevan} theoretically suggested that Mn clusters, i.e., Mn$_\mathrm{Ga}$-Mn$_\mathrm{I}$-Mn$_\mathrm{Ga}$ complex, are formed favorably under As-rich low-temperature growth conditions. Edmonds {\it et al.}\cite{Edmonds} also proposed that the dissociation process of Mn$_\mathrm{I}$-Mn$_\mathrm{Ga}$ pairs is likely promoted in the region where Mn$_\mathrm{Ga}$ concentration is rather high and further low-temperature annealing diffuses the dissociated Mn$_\mathrm{I}$ ions. Since the number of Mn$_\mathrm{I}$ ions compensating holes decreases with annealing in such a region (highly hole-concentrated cluster region), the $p$ should increase with annealing, enhancing the Curie temperature as seen in the shift of the high-temperature feature in the $\chi_{2}$" $- T$ curve towards higher temperature. It is also likely that the volume of the highly hole-concentrated cluster region increases with annealing at the same time. On the other hand, some of the dissociated Mn ions which diffuse in the annealing process may be trapped again in the matrix phase and the dissociated Mn ions compensate holes, resulting in a decrease in the $p$ in the matrix. Therefore, the low-temperature feature in the $\chi_{2}$" $- T$ curve is suppressed with annealing and simultaneously the contribution of the cubic anisotropy which we attribute to the matrix phase is also reduced. The other dissociated Mn ions may segregate at the surface region\cite{Edmonds} and have no contribution to the change in the magnetic properties of the matrix phase. As a result, the increase in the $p$ in the cluster regions is more significant than the decrease in the $p$ in the matrix. These descriptions seem to be consistent with theoretical and experimental results on the Mn distribution in low-temperature grown (Ga,Mn)As epilayers.\cite{Yu,Blinowski,Mahadevan,Raebiger,Sullivan,Ku,Edmonds,HamayaPRL,HamayaPRB}
In this context, the multiple-dip feature seen for 1-hour annealed sample is caused by an increase in the $p$ distribution, giving rise to the inhomogeneous magnetic properties.

From these consideration, we conclude that highly hole-concentrated (Ga,Mn)As clusters formed during the growth are the origin of [110] uniaxial magnetic anisotropy on the basis of our recent picture of the cluster/matrix magnetic structure in (Ga,Mn)As. That is, the switching of the magnetic anisotropy as a function of $p$ originates from the change in the relative volume of highly hole-concentrated (Ga,Mn)As clusters.

\section{CONCLUSION}
In conclusion, we have observed a systematic change in {\it K$_{u}$} and {\it K$_{1}$} as a function of hole concentration ($p$), where the contribution of the [110] uniaxial anisotropy becomes larger than the cubic magnetocrystalline anisotropy at the high hole concentration regime ($p \sim$ 10$^{21}$ cm$^{-3}$) at 4 K. We have revealed that the contribution of [110] uniaxial anisotropy is becoming significant with increasing $p$ in the highly hole-concentrated (Ga,Mn)As clusters and the growth of the clusters. This mechanism is the most probable origin of the switching of magnetic easy axis with increasing $p$.

\vspace{2mm}
Authors thank Prof. H. Munekata of Tokyo Institute of Technology for the collaboration in the study of ferromagnetic semiconductors.
The authors also thank Prof. H. Funakubo of Tokyo Institute of Technology for providing the opportunity to use their facility.

Materials Science (cond-mat.mtrl-sci); Mesoscale and Nanoscale Physics (cond-mat.mes-hall)