% ============================
% figures.tex
% ============================

\section*{Figures}

\begin{figure}[htbp]
\centering
\includegraphics[width=\textwidth]{figures/fig_composite_return_diagnostics.png}
\caption{
Composite return diagnostics comparing the continuous and event-driven TPW return models.
Top panel shows the model-derived global reorganization envelope (event-rate proxy).
Subsequent panels show mean absolute effective elevation adjustment rates
($|\Delta Z_{\mathrm{eff}} / \Delta t|$) at Early Homo and Early Civilization sites,
followed by the variance of site-level stability metrics.
All quantities are evaluated on a common temporal axis (ka BP).
}
\label{fig:composite_return_diagnostics}
\end{figure}

\begin{figure}[htbp]
\centering
\includegraphics[width=0.9\textwidth]{figures/fig_sites_dzdt_timeseries.png}
\caption{
Mean absolute effective elevation adjustment rate
($|\Delta Z_{\mathrm{eff}} / \Delta t|$) through time at occupied sites.
Early Homo and Early Civilization site ensembles are shown separately.
Sharp peaks correspond to discrete return steps in the event-driven model,
while the continuous model produces a smoother envelope.
}
\label{fig:sites_dzdt_timeseries}
\end{figure}

\begin{figure}[htbp]
\centering
\includegraphics[width=0.9\textwidth]{figures/fig_model_event_rate_envelope.png}
\caption{
Comparison of the continuous return model and the event-driven return model
expressed as normalized global reorganization envelopes.
The event-driven model concentrates adjustment into discrete episodes,
while the continuous model distributes adjustment over time via viscoelastic relaxation.
}
\label{fig:model_event_rate_envelope}
\end{figure}

% ============================
% Temporal null tests
% ============================

\begin{figure}[htbp]
\centering
\includegraphics[width=0.48\textwidth]{figures/fig_temporal_null_civ_variance.png}
\includegraphics[width=0.48\textwidth]{figures/fig_temporal_null_civ_dzdt.png}
\caption{
Temporal null distributions for Early Civilization sites.
Left: variance of site-level stability metric under temporal phase randomization.
Right: mean absolute effective elevation adjustment rate
($|\Delta Z_{\mathrm{eff}} / \Delta t|$).
Dashed vertical lines indicate observed values from the return model.
}
\label{fig:temporal_null_civ}
\end{figure}

\begin{figure}[htbp]
\centering
\includegraphics[width=0.48\textwidth]{figures/fig_temporal_null_homo_variance.png}
\includegraphics[width=0.48\textwidth]{figures/fig_temporal_null_homo_dzdt.png}
\caption{
Temporal null distributions for Early Homo sites.
Left: variance of site-level stability metric under temporal phase randomization.
Right: mean absolute effective elevation adjustment rate
($|\Delta Z_{\mathrm{eff}} / \Delta t|$).
Dashed vertical lines indicate observed values from the return model.
}
\label{fig:temporal_null_homo}
\end{figure}

% ============================
% Spatial null tests
% ============================

\begin{figure}[htbp]
\centering
\includegraphics[width=0.48\textwidth]{figures/fig_null_civ_variance.png}
\includegraphics[width=0.48\textwidth]{figures/fig_null_civ_dzdt.png}
\caption{
Spatial null distributions for Early Civilization sites.
Left: variance of site-level stability metric under spatial randomization.
Right: mean absolute effective elevation adjustment rate
($|\Delta Z_{\mathrm{eff}} / \Delta t|$).
Observed values are shown as dashed vertical lines.
}
\label{fig:spatial_null_civ}
\end{figure}

\begin{figure}[htbp]
\centering
\includegraphics[width=0.48\textwidth]{figures/fig_null_homo_variance.png}
\includegraphics[width=0.48\textwidth]{figures/fig_null_homo_dzdt.png}
\caption{
Spatial null distributions for Early Homo sites.
Left: variance of site-level stability metric under spatial randomization.
Right: mean absolute effective elevation adjustment rate
($|\Delta Z_{\mathrm{eff}} / \Delta t|$).
Observed values are shown as dashed vertical lines.
}
\label{fig:spatial_null_homo}
\end{figure}

% ============================
% End figures
% ============================