**Scientific analysis based on the primary source:** Mildner, S. (2025/2026). *A new interpretation of Ptolemy's Germania Magna: Employing computer-assisted image distortion of a medieval map by Donnus Nicolaus Germanus to examine post-glacial geodynamics in Europe*. EarthArXiv (Preprint). https://doi.org/10.31223/X5313T
---
## 1. Introduction and Research Interest
The historical geography of *Germania Magna* – the territory east of the Rhine and north of the Danube described by Claudius Ptolemy in the *Geographike Hyphegesis* (ca. 150 AD) – constitutes one of the methodologically most demanding fields of classical studies and geodetic research. The currently paradigmatically influential reference model, the statistical-geodetic rectification of the TU Berlin group (Karlsen et al., 2011), explains deviations between Ptolemaic coordinates and modern topography primarily as measurement errors of ancient instruments or as transmission artefacts.
Sven Mildner of Dresden opposes this concept with a fundamentally different approach. The **primary explanatory principle** of his model is not the correction of errors in the ancient map itself, but the recognition that the **northern reference line** – the coastline of the *Oceanus Germanicus* – lay approximately **120 km further south** during antiquity than today. Since medieval cartographers such as Donnus Nicolaus Germanus were unaware of this shift, they projected Ptolemaic coordinates onto an already geographically transformed landscape. The result was a systematic northward stretching of the map image, which inevitably produced a proportional eastward displacement of all eastern coordinates – thereby shifting the Ptolemaic *Vistula Fluvius* from its original Lusatian context all the way to the Polish Vistula. Geodynamic processes (reactivation of the Caledonian Deformation Front [CDF], lateral extrusion, block rotations) represent a **secondary**, here quantitatively investigated component.
The present paper formulates Mildner's rectification model in explicit mathematical terms, derives an affine coordinate transformation from the three invariant river-mouth anchor points, and conducts a systematic **residual analysis** for all gazetteer points. The resulting residual patterns are statistically tested and geodynamically interpreted.
---
## 2. The Cartometric Transformation Model
### 2.1 Scaling of the Ptolemaic Degree of Longitude
The core element of Mildner's rectification is an empirically determined, spatially fixed scaling factor $k$ for the Ptolemaic degree of longitude. This is derived from the physical distance between the mouths of two invariant reference rivers – the Rhenus Fluvius (central mouth, $\lambda_P = 27{.}00°$) and the Albis Fluvius ($\lambda_P = 31{.}00°$), as well as the Vistula Fluvius ($\lambda_P = 45{.}00°$; Mildner identification: Oderberg, at the mouth of the reconstructed "United Vistula" into the *Oceanus Germanicus*).
The two independent baseline estimates are (Mildner, 2025/2026, section on degree-length determination):
$$k_1 = \frac{d_{\text{Rh–El}}}{\Delta\lambda_{P,\,\text{Rh–El}}} = \frac{\approx 115\,\text{km}}{4°} \approx 28{.}75\,\frac{\text{km}}{°}$$
$$k_2 = \frac{d_{\text{Rh–Vi}}}{\Delta\lambda_{P,\,\text{Rh–Vi}}} = \frac{\approx 490\,\text{km}}{18°} \approx 27{.}22\,\frac{\text{km}}{°}$$
The weighted mean (weighted by baseline length in Ptolemaic degrees: $4°$ and $18°$ respectively) yields:
$$k = \frac{4 \cdot k_1 + 18 \cdot k_2}{22} = \frac{115 + 490}{22} \approx 27{.}5\,\frac{\text{km}}{°}$$
For the back-transformation into geographic degrees of longitude at mean latitude $\bar{\phi} \approx 52{.}5°\,\text{N}$:
$$\frac{\Delta\lambda_{\text{mod}}}{{\Delta\lambda_P}} = \frac{k}{111{.}3\,\text{km/°} \cdot \cos\bar{\phi}} = \frac{27{.}5}{111{.}3 \times 0{.}609} = 0{.}406\,\frac{°_{\text{geogr.}}}{°_{\text{Ptol.}}}$$
### 2.2 Affine Coordinate Transformation Model
The complete coordinate transformation from Ptolemaic to modern geographic coordinates is modelled as an **affine mapping**:
$$\lambda_{\text{mod}} = a_1 + a_2 \cdot \lambda_P + a_3 \cdot \phi_P$$
$$\phi_{\text{mod}} = b_1 + b_2 \cdot \lambda_P + b_3 \cdot \phi_P$$
with six coefficients to be determined. The minimisation functional for the least-squares adjustment over all $n$ gazetteer points is:
$$S = \sum_{i=1}^{n} w_i \left[\left(\lambda_{\text{mod},i} - \hat{\lambda}_{\text{mod},i}\right)^2 + \left(\phi_{\text{mod},i} - \hat{\phi}_{\text{mod},i}\right)^2\right] \rightarrow \min$$
Since exactly three invariant anchor points (Rhine, Elbe, and Vistula mouths) are available for calibration, the system of equations for both the longitude and the latitude transformation is exactly determined. The three anchor points define the transformation parameters completely.
### 2.3 Solution of the System of Equations
**Calibration points** (coordinate values in decimal degrees):
| Point | $\lambda_P$ | $\phi_P$ | $\lambda_{\text{mod}}$ | $\phi_{\text{mod}}$ |
|---|---|---|---|---|
| Rhenus Fl. (central mouth) | 27.00 | 53.167 | 6.750 | 52.250 |
| Albis Fl. (mouth) | 31.00 | 56.250 | 8.583 | 53.183 |
| Vistula Fl. (mouth/Oderberg) | 45.00 | 56.000 | 14.150 | 52.867 |
The linear system of equations for the longitude transformation reads in matrix form:
$$\begin{pmatrix} 1 & 27{.}00 & 53{.}17 \\ 1 & 31{.}00 & 56{.}25 \\ 1 & 45{.}00 & 56{.}00 \end{pmatrix} \begin{pmatrix} a_1 \\ a_2 \\ a_3 \end{pmatrix} = \begin{pmatrix} 6{.}750 \\ 8{.}583 \\ 14{.}150 \end{pmatrix}$$
The solution yields the following **transformation parameters**:
$$\boxed{\lambda_{\text{mod}} = -8{.}114 + 0{.}3989\,\lambda_P + 0{.}0770\,\phi_P}$$
$$\boxed{\phi_{\text{mod}} = +35{.}458 - 0{.}0167\,\lambda_P + 0{.}3244\,\phi_P}$$
The **longitude scaling parameter** $a_2 = 0{.}3989$ corresponds in ground kilometres at $\bar{\phi} = 52{.}5°\,\text{N}$:
$$k_{\lambda} = a_2 \cdot 111{.}3 \cdot \cos(52{.}5°) = 0{.}3989 \times 67{.}7 = 27{.}0\,\frac{\text{km}}{°_P}$$
This confirms Mildner's stated value of $\approx 28$ km per Ptolemaic degree of longitude with a deviation of only 3.5 %. The **latitude scaling parameter** $b_3 = 0{.}3244$ demonstrates the characteristic strong latitude compression of the Ptolemaic system at higher latitudes:
$$k_{\phi} = b_3 \cdot 111{.}3 = 0{.}3244 \times 111{.}3 = 36{.}1\,\frac{\text{km}}{°_P}$$
The pronounced asymmetry between longitude ($\approx 27$ km/°) and latitude scaling ($\approx 36$ km/°) reflects the systematic latitude distortion in the Ptolemaic coordinate system for the northern European region – the degrees of latitude are more strongly compressed at $\phi > 52°$ in the Ptolemaic system than the degrees of longitude.
---
## 3. Residual Analysis of the Gazetteer
### 3.1 Methodology
For all 22 non-calibration points in the gazetteer (river sources, additional river mouths, settlements, mountains), the **prediction residuals** are calculated as the difference between the affinely transformed prediction $(\hat{\lambda}_{\text{mod}},\, \hat{\phi}_{\text{mod}})$ and Mildner's identification $(\lambda_{\text{Mild}},\, \phi_{\text{Mild}})$:
$$\Delta\lambda = \hat{\lambda}_{\text{mod}} - \lambda_{\text{Mild}}, \qquad \Delta\phi = \hat{\phi}_{\text{mod}} - \phi_{\text{Mild}}$$
The conversion to ground kilometres is performed at the mean latitude of each point $\bar{\phi}$:
$$\Delta\lambda_{\text{km}} = \Delta\lambda \cdot 111{.}3 \cdot \cos\bar{\phi}, \qquad \Delta\phi_{\text{km}} = \Delta\phi \cdot 111{.}3$$
The scalar total residual vector is calculated as the Euclidean norm:
$$r_i = \sqrt{\Delta\lambda_{\text{km},i}^2 + \Delta\phi_{\text{km},i}^2}$$
### 3.2 Results of the Residual Analysis
**Table 1:** Residual analysis of all gazetteer points. A negative $\Delta\lambda_{\text{km}}$ indicates that the transformed prediction lies west of the Mildner identification (i.e., the identification site lies further east than the linear model predicts).
| No. | Ptolemaic Name | Identification (Mildner) | $\lambda_P$ | $\phi_P$ | $\hat\lambda$ | $\hat\phi$ | $\lambda_M$ | $\phi_M$ | $\Delta\lambda_{\text{km}}$ | $\Delta\phi_{\text{km}}$ | $r$ [km] | Group |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| **K1** | Rhenus Fl. (mouth) | Hengelo/Enschede | 27.00 | 53.17 | 6.750 | 52.250 | 6.750 | 52.250 | 0.0 | 0.0 | **0.0** | Cal. |
| **K2** | Albis Fl. (mouth) | NW of Bremen | 31.00 | 56.25 | 8.583 | 53.183 | 8.583 | 53.183 | 0.0 | 0.0 | **0.0** | Cal. |
| **K3** | Vistula Fl. (mouth) | Oderberg | 45.00 | 56.00 | 14.150 | 52.867 | 14.150 | 52.867 | 0.0 | 0.0 | **0.0** | Cal. |
| **F1** | Vistula main source | Königswartha | 44.00 | 52.50 | 13.481 | 51.751 | 14.367 | 51.283 | −61.5 | +52.1 | 80.5 | Lusatia |
| **F2** | Vistula western source | Königsbrück/Pulsnitz | 40.17 | 52.67 | 11.964 | 51.872 | 13.883 | 51.267 | −133.0 | +67.4 | 149.2 | Lusatia |
| **F3** | Chalusus Fl. (mouth) | Havelberg | 37.00 | 56.00 | 10.957 | 53.006 | 12.067 | 52.817 | −74.5 | +21.0 | 77.4 | Coast |
| **F4** | Suebus Fl. (mouth) | Neuruppin/Fehrbellin | 39.50 | 56.00 | 11.955 | 52.964 | 12.900 | 52.817 | −63.4 | +16.4 | 65.5 | Coast |
| **F5** | Viadua Fl. (mouth) | Finowfurt/Marienwerder | 42.50 | 56.00 | 13.151 | 52.914 | 13.633 | 52.850 | −32.3 | +7.1 | 33.1 | Coast |
| **S1** | Agrippinensis* | Cologne (Altstadt) | 27.67 | 51.17 | 6.860 | 51.597 | 6.958 | 50.941 | −6.9 | +73.0 | 73.3 | Gallia Belg. |
| **S2** | Aliso* | Haltern am See | 28.00 | 51.50 | 7.021 | 51.696 | 7.324 | 51.712 | −20.9 | −1.8 | 21.0 | Gallia Belg. |
| **S3** | Budorigum | Doberlug-Kirchhain | 41.00 | 52.67 | 12.296 | 51.858 | 13.554 | 51.616 | −87.1 | +26.9 | **91.2** | Elster-Cl. |
| **S4** | Calisia | Calau | 43.75 | 52.83 | 13.406 | 51.866 | 13.960 | 51.743 | −38.2 | +13.7 | 40.6 | Lusatia-E |
| **S5** | Limis Lucus | Baruth/Mark | 41.00 | 53.50 | 12.361 | 52.128 | 13.503 | 52.045 | −78.2 | +9.2 | **78.7** | Elster-Cl. |
| **S6** | Lugidunum | Falkenberg/Elster | 39.50 | 52.50 | 11.686 | 51.826 | 13.269 | 51.606 | −109.5 | +24.5 | **112.2** | Elster-Cl. |
| **S7** | Stragona | Herzberg/Elster | 39.67 | 52.33 | 11.740 | 51.773 | 13.200 | 51.667 | −101.0 | +11.8 | **101.7** | Elster-Cl. |
| **S8** | Treva | Bremen | 33.00 | 55.67 | 9.336 | 52.965 | 8.939 | 53.036 | +26.6 | −7.9 | 27.8 | Coast-W |
| **S9** | Lirimiris | Bispingen/Soltau | 34.50 | 55.50 | 9.922 | 52.886 | 10.010 | 53.107 | −5.9 | −24.6 | 25.3 | Coast-W |
| **G1** | Asciburgius Mons NW | Fläming/E. of Magdeburg | 39.00 | 54.00 | 11.601 | 52.325 | 12.283 | 52.167 | −46.5 | +17.6 | 49.7 | Fläming |
| **G2** | Asciburgius Mons SE | Calauer Schweiz/Senftenberg | 44.00 | 52.50 | 13.481 | 51.751 | 13.933 | 51.583 | −31.3 | +18.7 | 36.3 | Fläming |
| **G3** | Melibocus Mons W | Harz/Weser-Leine Highlands | 33.00 | 52.50 | 9.093 | 51.935 | 9.433 | 52.000 | −23.3 | −7.2 | 24.4 | Harz |
| **G4** | Melibocus Mons E | Harz/Eisleben | 37.00 | 52.50 | 10.688 | 51.868 | 11.283 | 51.517 | −41.2 | +39.1 | 55.2 | Harz |
| **G5** | Sudete Mons W | Thuringian Forest/Kassel | 34.00 | 50.00 | 9.299 | 51.110 | 9.583 | 51.233 | −19.8 | −13.7 | 23.9 | Thuringia |
| **G6** | Sudete Mons E | Thuringian Slate Mts. | 40.00 | 50.00 | 11.692 | 51.010 | 11.533 | 50.367 | +11.2 | +71.6 | 72.5 | Thuringia |
| **G7** | Sarmate Mons N | Lusatian Highlands | 43.50 | 50.50 | 13.127 | 51.113 | 14.467 | 51.133 | −93.9 | −2.2 | 94.0 | Lusatia |
| **G8** | Abnobae Mons W | Taunus (Hoher Taunus) | 31.00 | 49.00 | 8.025 | 50.836 | 7.867 | 50.017 | +11.2 | +91.2 | 91.9 | South |
*Agrippinensis and Aliso derive from Chapter 8 of the *Geographike Hyphegesis* (Gallia Belgica) and were recorded under a different measurement system; their latitude residuals are therefore to be evaluated separately.*
---
## 4. Statistical Evaluation and Group Analysis
### 4.1 Regional RMSE Analysis
From the residual magnitudes, the root mean square error (RMSE) is calculated for the homogeneous subgroups:
$$\text{RMSE}_G = \sqrt{\frac{1}{n_G}\sum_{i \in G} r_i^2}$$
| Group | Points | $n$ | RMSE [km] | Mean $\Delta\lambda_{\text{km}}$ |
|---|---|---|---|---|
| Calibration (river mouths) | K1–K3 | 3 | **0.0** | 0.0 |
| **Elster Cluster** (Budorigum, Limis Lucus, Lugidunum, Stragona) | S3, S5, S6, S7 | 4 | **96.8** | −94.0 |
| Coastal settlements W (Treva, Lirimiris) | S8, S9 | 2 | **26.6** | +10.4 |
| Coastal rivers (Chalusus, Suebus, Viadua) | F3–F5 | 3 | **59.5** | −56.7 |
| Fläming (Asciburgius Mons NW/SE) | G1, G2 | 2 | **43.0** | −38.9 |
| Harz (Melibocus Mons) | G3, G4 | 2 | **42.7** | −32.3 |
| Thuringian Forest (Sudete Mons) | G5, G6 | 2 | **54.0** | −4.3 |
| Lusatia/Sarmate | G7, F1 | 2 | **87.3** | −77.7 |
| Gallia Belgica | S1, S2 | 2 | **47.2** | −13.9 |
The most striking result is the factor of **3.6** between the RMSE of the Elster Cluster (96.8 km) and the RMSE of the coastal settlements (26.6 km). This discrepancy is incompatible with spatially uniform measurement errors and requires a geodynamic explanation.
### 4.2 t-Test for Systematic Eastward Offset of the Elster Cluster
For the four Elster Cluster points, the $\Delta\lambda$ residuals in degrees are analysed:
$$\Delta\lambda \in \{-1{.}258°;\; -1{.}142°;\; -1{.}583°;\; -1{.}460°\}$$
$$\overline{\Delta\lambda} = \frac{-1{.}258 - 1{.}142 - 1{.}583 - 1{.}460}{4} = -1{.}361°$$
$$s_{\Delta\lambda} = \sqrt{\frac{\sum_{i=1}^{4}(\Delta\lambda_i - \overline{\Delta\lambda})^2}{n-1}} = \sqrt{\frac{0{.}0106 + 0{.}0480 + 0{.}0493 + 0{.}0098}{3}} = 0{.}198°$$
The one-tailed t-test for $H_0: \mu_{\Delta\lambda} = 0$ (no systematic offset) against $H_1: \mu_{\Delta\lambda} < 0$ (systematic western position of the prediction, i.e. eastward offset of the identification) yields:
$$t = \frac{\overline{\Delta\lambda}}{s_{\Delta\lambda}/\sqrt{n}} = \frac{-1{.}361}{0{.}198/\sqrt{4}} = \frac{-1{.}361}{0{.}099} = -13{.}7$$
For $df = n - 1 = 3$, the critical t-value at $\alpha = 0{.}001$ (one-tailed) is $t_{\text{crit}} = -7{.}45$. Since $|t| = 13{.}7 > 7{.}45$, $H_0$ is rejected at the **0.1% significance level** ($p < 0{.}001$). The mean offset of:
$$\overline{\Delta\lambda}_{\text{km}} = -1{.}361° \times 111{.}3 \times \cos(52°) = -1{.}361 \times 68{.}4 = -93{.}1\,\text{km}$$
is therefore **statistically highly significant**. The negative direction of the offset (prediction lies west of the identification) means: the entire Elster/Fläming/Lusatia block is systematically approximately **93 km further east** in Mildner's model than the linear coastal transformation predicts.
### 4.3 Moran's I – Qualitative Spatial Autocorrelation
With only $n = 22$ non-calibration points, a formal Moran's I test is statistically not very informative. Qualitatively, however, the residual structure shows clear positive spatial autocorrelation: the four Elster Cluster points (S3, S5, S6, S7), which are spatially adjacent in the Elster/Lusatia corridor, all display strongly negative $\Delta\lambda_{\text{km}}$ values (−78 to −110 km), while the geographically more distant Treva (Bremen, +27 km) and Lirimiris (Bispingen, −6 km) show markedly smaller residuals and, in the case of Treva, even oppositely directed ones. This pattern – spatially coherent signal at geographically proximate points, smaller residuals at more distant regions – is the classic signature of a **geodynamically localised block offset**, not of uniform measurement imprecision.
---
## 5. Geodynamic Interpretation of the Residual Patterns
### 5.1 The Elster Cluster: Signature of the Tectonic Eastward Offset
The statistically highly significant eastward offset of the Elster/Fläming/Lusatia Cluster ($\overline{\Delta} \approx 93$ km eastward, $p < 0{.}001$) cannot be explained by ancient measurement errors or random identification uncertainty. Four geographically proximate points, independently identified through different methodologies, display the same sign and a similar magnitude of the offset. This is consistent with the hypothesis of a **regional, coherent crustal block offset**.
Within the framework of Mildner's model (Mildner, 2025/2026), this offset is explained by transpressively reactivated CDF tectonics: the reactivation of the Caledonian Deformation Front generated a NNW-SSE directed compressional regime. As a consequence of the **lateral extrusion** of sedimentary masses from the northwest (towards the Cimbrian Peninsula), accretionary wedges advanced against the western flank of the Fläming. Since the Lusatian crustal block in the southeast (Senftenberg area) was anchored as a **fixed rotation point** (pivot), no linear evasion occurred; instead a **dextral rotation** of the Fläming massif took place. The western limb was displaced eastward, resulting in the ~90 km eastward offset of today's Elster-Lusatia complex as postulated in the model.
Geophysically, Deutschmann et al. (2018) support this mechanism: the polyphase fault history west of Rügen within the TESZ documents six tectonic reactivation phases from the Caledonian collision to the Late Cretaceous–Palaeogene inversion tectonics. Lyngsie & Thybo (2007) document a 150 km wide overthrust zone of Avalonian crust over the Baltica lower crustal shield – the deep-crustal geometry that mechanically enables such a crustal block offset in geologically younger time.
### 5.2 The Geochemical Verification (Anthracite near Budorigum/Doberlug-Kirchhain)
The residual of Budorigum (S3) with $r = 91{.}2$ km is of particular significance. This point lies precisely at the tectonic hinge zone postulated in the model. The existence of **near-surface anthracite deposits** at Doberlug-Kirchhain (standard deep-burial anthracite forms only above $T > 150\,°C$ at depths of several kilometres) is interpreted in the Mildner model as a direct proxy for **stress metamorphism** at the upper boundary of the shear zone. The convergence of two completely independent lines of evidence – the cartometric residual vector (from the coordinate transformation) and the petrographic anomaly (from field geology) – at the identical geographic location is methodologically extraordinarily strong and cannot be explained by coincidental agreement.
### 5.3 Coast-Proximate Points and Southern Outliers
The small residuals of the coastal settlements Treva (Bremen, $r = 27{.}8$ km) and Lirimiris (Bispingen, $r = 25{.}3$ km) are **methodologically expected**: since the calibration is based on the three river-mouth points (likewise near the coast), the transformation fit is optimal in the coastal region. The small residuals validate the internal consistency of the calibration framework.
The southern Abnobae Mons (Taunus, G8) displays a large $\Delta\phi_{\text{km}} = +91$ km, reflecting the nonlinear latitude behaviour of the Ptolemaic coordinate system in very southern regions of *Germania Magna*. The affine calibration on the northern coastal line points is considerably less precise for southern projections – which Mildner employs as an argument for the three-dimensional character of Ptolemaic latitude distortion.
---
## 6. Methodological Defence against Criticism
### 6.1 The Rubber-Sheeting Argument and its Refutation
Critics argue that through arbitrary digital map distortion, infinitely many alternative fits could be generated (the so-called rubber-sheeting accusation). This argument fundamentally misses the difference between uncontrolled topological morphing and the **strictly regulated morphometric model** presented here. The transformation is constrained by an **inviolable matrix of simultaneous constraints**:
1. **Geometric scaling rigidity:** $k = 28$ km/° derived from the empirically measurable Rhine-Elbe baseline – not locally optimisable.
2. **Hydrographic topological constraint:** An alternative model must find a river system in which two separate major source streams travel more than 50% of their northward course *south* of a specific mountain range and converge *east* of it – a topological requirement exactly fulfilled in the Lusatian Schwarze Elster/Spree system, but not in the Polish Vistula system.
3. **Cartographic curvature constraint:** The purely graphic bend of the Asciburgius Mons on the Germanus map must correspond to a geologically verified tectonic hinge zone – a constraint equally non-freely selectable.
4. **Geochemical anchor point:** The coordinate-geometrically positioned identification Budorigum = Doberlug-Kirchhain falls precisely on a stress metamorphism anomaly independently measurable from cartometry.
No AI-based or purely statistically-geodetic rectification algorithm can satisfy these four simultaneous constraints whilst preserving topological reality, because such algorithms minimise the Euclidean error across all data points – which inevitably pulls the Vistula back onto the Vistula and leaves the hydrographic paradox of Central Bohemia unresolved (Mildner, 2025/2026, section 9.0.1).
### 6.2 Criticism: Archaeological Finds Refute the Model
Mildner's hypothesis does not generally deny the existence of pre-catastrophic settlement traces. Rather, the model argues with nuance:
- In the marginal zones of the impact area, where geodynamic processes operated primarily in deeper crustal layers, surface structures could partially survive.
- The dating of stratigraphic layers is methodologically limited by zircon age inheritance, ¹⁴C resetting through CO₂ input from secondary volcanism, and OSL inaccuracies under turbulent sedimentation.
- Volkmann's (2014) archaeological findings document non-linear settlement breaks within a few decades – precisely *inconsistent* with gradual transformation models, but consistent with a geodynamically catastrophic explanation.
### 6.3 Falsifiability and Scientific Status
The model is explicitly falsifiable. Concrete test points include:
- Targeted archaeological deep drilling at the newly calculated coordinates (e.g., Budorigum/Doberlug-Kirchhain, Lugidunum/Falkenberg/Elster): if no traces of a significant settlement from the 1st–5th century AD are found, the model is refutable.
- Micromorphological analysis of Dark Earth horizons against the criteria P1–P6 of the Event-Dark-Earth test concept (Mildner, 2026).
- Independent dating of the matrix phase of the Český Kráter (Rajlich, 1992) through high-precision OSL measurements of the amorphous binder phase – which would either support or refute the age inheritance hypothesis.
---
## 7. Conclusions
The present analysis has subjected Mildner's geodynamic rectification model of *Germania Magna* to a complete, coordinate-based **residual analysis of the gazetteer** for the first time. The key results are:
1. **Scaling consistency:** The longitude scaling parameter $k = 27{.}0$ km/° derived from the affine calibration agrees with Mildner's postulated $\approx 28$ km/° within measurement precision. The internal consistency of the model is thus cartometrically verified.
2. **Statistically significant eastward offset structure:** The Elster/Fläming/Lusatia Cluster displays a highly significant eastward offset of $\overline{\Delta} = -93{.}1$ km ($t = -13{.}7$, $p < 0{.}001$, $df = 3$). This is incompatible with uniform measurement errors and requires a geodynamic explanation: the lateral, CDF-reactivation-induced eastward offset of the Fläming-Lusatia block.
3. **Spatially autocorrelated residual structure:** Coastal points (RMSE $\approx 27$ km) show qualitatively good model fit; the Elster inland cluster shows three times higher residual magnitudes (RMSE $\approx 97$ km). This spatial structure is consistent with a **geodynamically controlled, non-uniform deformation** of the Ptolemaic landscape since antiquity.
4. **Geochemical convergence:** The convergence of the cartometric identification (Budorigum = Doberlug-Kirchhain) with the measured anthracite stress metamorphism anomaly at precisely this location provides a methodologically independent verification of the model.
5. **Methodological superiority:** The model is not a rubber-sheeting approach. It is rigidly constrained by four simultaneous, independent constraints (scaling, hydrography, curvature morphology, geochemistry) and is thereby fundamentally distinguished from arbitrary map distortion.
The analysis demonstrates that Mildner's rectification approach is not only cartographically coherent, but statistically significant and geodynamically founded. It merits a systematic empirical examination through targeted archaeological deep prospection and high-resolution geophysical mapping of the identified regions.
---
## References (APA)
Abbott, D. H., Breger, D., Biscaye, P. E., Barron, J. A., Juhl, R. A., & McCafferty, P. (2014). What caused terrestrial dust loading and climate downturns between A.D. 533 and 540? In G. Keller & A. C. Kerr (Eds.), *Volcanism, impacts, and mass extinctions: Causes and effects* (GSA Special Papers, Vol. 505, pp. 421–438). Geological Society of America. https://doi.org/10.1130/2014.2505(23)
Deutschmann, A., Meschede, M., & Obst, K. (2018). Fault system evolution in the Baltic Sea area west of Rügen, NE Germany. *Geological Society, London, Special Publications, 469*, 83–98. https://doi.org/10.1144/SP469.24
Gaberz, S. (2014). *Dark Earth – die schwarze Schicht* [Diploma thesis, Karl-Franzens-Universität Graz]. UNIGRAZonline. https://unipub.uni-graz.at/obvugrhs/content/titleinfo/239930
Geersen, J., Bradtmöller, M., Schneider von Deimling, J., Feldens, P., Auer, J., Held, P., Lohrberg, A., Supka, R., Hoffmann, J. J. L., Eriksen, B. V., Rabbel, W., Karlsen, H.-J., Krastel, S., Brandt, D., Heuskin, D., & Lübke, H. (2024). A submerged Stone Age hunting architecture from the Western Baltic Sea. *Proceedings of the National Academy of Sciences of the United States of America, 121*, e2312008121. https://doi.org/10.1073/pnas.2312008121
Götze, J., Möckel, R., Pan, Y., & Müller, A. (2024). Geochemistry and formation of agate-bearing lithophysae in Lower Permian volcanics of the NW-Saxonian Basin (Germany). *Mineralogy and Petrology, 118*, 23–40. https://doi.org/10.1007/s00710-023-00841-2
Karlsen, H.-J., Marx, C., & Lelgemann, D. (2011). Germania magna – ein neuer Blick auf eine alte Karte: entzerrte geographische Daten des Ptolemaios für die antiken Orte zwischen Rhein und Weichsel. *Germania, 89*, 115–155. https://doi.org/10.11588/ger.2011.96480
Lambeck, K., Smither, C., & Johnston, P. (1998). Sea-level change, glacial rebound and mantle viscosity for northern Europe. *Geophysical Journal International, 134*(1), 102–144. https://doi.org/10.1046/j.1365-246x.1998.00541.x
Lyngsie, S. B., & Thybo, H. (2007). A new tectonic model for the Laurentia–Avalonia–Baltica sutures in the North Sea: A case study along MONA LISA profile 3. *Tectonophysics, 429*, 201–227. https://doi.org/10.1016/j.tecto.2006.09.017
Mildner, S. (2025/2026). *A new interpretation of Ptolemy's Germania Magna: Employing computer-assisted image distortion of a medieval map by Donnus Nicolaus Germanus to examine post-glacial geodynamics in Europe* [EarthArXiv preprint]. https://doi.org/10.31223/X5313T
Nielsen, S., Stephenson, R., & Thomsen, E. (2007). Dynamics of Mid-Palaeocene North Atlantic rifting linked with European intra-plate deformations. *Nature, 450*, 1071–1074. https://doi.org/10.1038/nature06379
PAGES 2k Consortium. (2013). Continental-scale temperature variability during the past two millennia. *Nature Geoscience, 6*, 339–346. https://doi.org/10.1038/ngeo1797
Rajlich, P. (1992). *Bohemian circular structure, Czechoslovakia: Search for the impact evidence* [Conference paper]. Geoterra.
Seyhan, İ. (1971/2023). Die Entstehung des vulkanischen Kaolins und das Andesit-Problem. *Bulletin of the Mineral Research and Exploration, 76*.
Thybo, H., & Nielsen, C. A. (2012). Seismic velocity structure of crustal intrusions in the Danish Basin. *Tectonophysics, 572*, 64–75. https://doi.org/10.1016/j.tecto.2011.11.019
Volkmann, A. (2014). Neues zur „Odergermanischen Gruppe": Das innere Barbaricum an der unteren Oder im 5.–6. Jh. AD. Universitätsbibliothek Heidelberg. http://nbn-resolving.de/urn:nbn:de:bsz:16-heidok-159188
Wiktionary. (2023). *ustulate*. https://en.wiktionary.org/w/index.php?title=ustulate&oldid=76097086
