---
title: 'Gaia-Enceladus/Sausage: Milky Way Merger Remnant'
url: https://www.emergentmind.com/topics/gaia-enceladus-sausage-ges-4d4a66cf-ad59-4db2-add2-9c3f7ac0d0fa
type: topic
---

# Gaia-Enceladus/Sausage: Milky Way Merger Remnant

Gaia-Enceladus/Sausage, also written Gaia-Sausage-Enceladus and often abbreviated GSE or GS/E, is the accretion remnant associated with the Milky Way’s last major or last significant merger. In current usage, the term denotes both the disrupted progenitor dwarf galaxy and the present-day stellar debris that dominates much of the inner stellar halo. Observationally, the system is identified through strongly radial kinematics, low net Galactic rotation, and distinctive abundance patterns; dynamically, it is widely placed at a lookback time of roughly \(8\)–\(11\) Gyr and linked to a major restructuring of the Galaxy’s halo and disk [2606.04462].

## 1. Definition and historical placement

GSE was first recognized in Gaia-era phase-space analyses as a nearly radially anisotropic component of the inner halo. Its canonical signature is a large radial-to-tangential velocity-dispersion ratio, usually summarized by the anisotropy parameter
\[
\beta \equiv 1 - \frac{\sigma_\phi^2 + \sigma_\theta^2}{2\sigma_r^2},
\]
with values around \(\beta \approx 0.8\) in the Solar neighborhood in early characterizations, and often \(\beta \approx 0.8\)–\(0.9\) for GSE-like debris in later summaries. A second common diagnostic is the orbital eccentricity
\[
e \equiv \frac{r_{\rm apo} - r_{\rm peri}}{r_{\rm apo} + r_{\rm peri}},
\]
with GSE selections frequently requiring \(e \gtrsim 0.7\) [2408.02723].

The system is usually interpreted as the remnant of the last major merger experienced by the Milky Way, and several studies place the event between \(8\) and \(11\) Gyr ago. In this picture, GSE contributed a large fraction of the inner halo, produced the characteristic “sausage” morphology in velocity space, and deposited stars on low-\(L_z\), high-\(J_R\), highly eccentric orbits. Orbit-superposition reconstructions further indicate that the GSE-like component is nearly spherical, with \(q \sim 0.9\), and exhibits an outer ridge at \(r \sim 20\) kpc, interpreted as an apocentric pile-up of the progenitor’s orbit [2108.10525].

## 2. Phase-space selection and present-day debris field

No single observational definition of GSE is universal. Common selections use eccentricity, energy–angular momentum space, radial-action–angular-momentum space, or combined chemical tagging. A widely used action-space box requires
\[
-500 \le L_z \le 500\ {\rm kpc\,km\,s^{-1}}
\]
and
\[
30 \le \sqrt{J_R} \le 55\ (\mathrm{kpc\,km\,s^{-1}})^{1/2},
\]
with a closely related “clean” SkyMapper–Gaia RVS selection adopting \(30 \le \sqrt{J_R} \le 50\) and the same \(L_z\) limits. That latter sample yielded \(N=679\) stars and a narrow MDF centered at \({\rm [Fe/H]}=-1.17\) dex with dispersion \(0.34\) dex, and was argued to minimize contamination relative to velocity-space or \(E\)–\(L_z\) selections [2003.11039].

Selection dependence remains a central methodological issue. A systematic comparison of five observational strategies—eccentricity, \(E\)–\(L_z\), \(J_r\)–\(L_z\), action diamond, and \([{\rm Mg/Mn}]\)–\([{\rm Al/Fe}]\)—found samples ranging from \(144\) to \(1{,}279\) stars, with markedly different contamination and completeness. In the Auriga benchmark, the \(J_r\)–\(L_z\) method was best for sample purity, while the eccentricity method was best for completeness. The same study concluded that GSE selection cannot be treated as interchangeable across science cases [2306.00770].

The debris field is not confined to the local halo. HAC and VOD have orbital eccentricity distributions comparable to GSE, strongly peaked MDFs, and indistinguishable abundance distributions in the analyzed chemical spaces, supporting their interpretation as unmixed GSE debris in two major halo overdensities [2207.11869]. At smaller Galactocentric radii, RR Lyrae mapping of the inner-central halo indicates a much lower GSE fraction than in the Solar neighborhood: only \(6\)–\(9 \pm 2\%\) of the inner-central halo RR Lyrae population is consistent with a GES origin, and very few RR Lyrae stars with bulge-confined orbits are expected to have originated from GES [2507.11741].

## 3. Metallicity, abundance ratios, and age structure

Published GSE metallicity distributions are similar in broad form but not identical in detail. A self-consistent APOGEE DR17 \(+\) Gaia EDR3 characterization reported a median \({\rm [Fe/H]}=-1.22\) dex with \(0.23\) dex dispersion for stars selected by novel chemodynamical criteria [2206.10505]. Other APOGEE-based selections found a mean \({\rm [Fe/H]}\sim -1.15\) with an approximately Gaussian width of \(\sim 0.3\) dex, while high-eccentricity SEGUE-based fits for \(e>0.7\) yielded \(\mu_{\rm [Fe/H]}=-1.39\pm0.01\) and \(\sigma_{\rm [Fe/H]}=0.20\) dex in one comparison sample [2105.12141].

In abundance space, GSE is typically described as a single \([{\rm Mg/Fe}]\) sequence with the onset of SN Ia enrichment marked by an \(\alpha\)-knee. Different selections place that knee at \({\rm [Fe/H]}\sim -1.4\) or \({\rm [Fe/H]}\simeq -1.2\). APOGEE analyses also emphasize uniformly low \([{\rm Al/Fe}]\sim -0.25\) dex relative to the Milky Way disc, and a tight accreted locus at high \([{\rm Mg/Mn}]\) and low \([{\rm Al/Fe}]\); one empirical definition of the chemically unevolved accreted locus is \([{\rm Al/Fe}]<-0.1\) and \([{\rm Mg/Mn}]>+0.35\) [2301.01302]. At the same time, GSE stars show an excess of \([{\rm Al/Fe}]\), \([{\rm Mg/Mn}]\), and \([{\rm Mg/Fe}]\) in comparison to surviving Milky Way dwarf satellites, a contrast interpreted as a difference in star-formation efficiencies and timescales [2206.10505].

Neutron-capture abundances add a further temporal constraint. A SAGA-based compilation found that the joint behavior of \([{\rm Eu/Mg}]\), \([{\rm Ba/Mg}]\), and \([{\rm Eu/Ba}]\) indicates a prolonged period of slow star formation lasting over \(2\) Gyr, followed by abrupt quenching, with no super-solar \([{\rm Fe/Mg}]\) and a truncation near \({\rm [Fe/H]}\approx -0.5\). In that interpretation, GSE lies between Sculptor-like short-duration and Fornax-like very extended star-formation histories [2505.06606].

Direct age inference from Gaia CMD fitting yields a more structured chronology. ChronoGal recovered three main populations and a fourth smaller one arranged along an almost linear age–\([{\rm M/H}]\) relation. The three oldest populations correspond to the bulk of star formation lasting for at least \(\sim 3\)–\(4\) Gyr and ending about \(10\) Gyr ago, with metallicities from \(-1.7\) to \(-0.8\). These were grouped into an “isolation epoch” from \(\sim 13.5\) to \(\sim 12.0\) Gyr and a “merger-event epoch” from \(\sim 12.0\) to \(\sim 10.5\) Gyr. The fourth population, at \(\sim 8.5\) Gyr and \({\rm [M/H]}\sim -0.4\), has an unclear link to GSE [2502.20439].

## 4. Progenitor mass, geometry, globular clusters, and internal gradients

Mass estimates for the progenitor vary substantially with methodology. Redshift-dependent mass–metallicity inversions using \({\rm [Fe/H]}\approx -1.17\) and a mean stellar age of \(10\) Gyr gave \(M_\star \simeq 10^{8.85}-10^{9.85}\,M_\odot\) in one study [2003.11039]. Related estimates gave \(M_{\star,\rm GSE}\sim10^{8.5-9.5}\,M_\odot\), and a benchmark across several observational selections yielded an average \(M_\star\) of \(\sim 4\times10^8\,M_\odot\) from the \(z=0\) mass–metallicity relation [2306.00770].

Direct density-modeling of high-purity samples produces lower stellar masses. A distribution-function-corrected fit to APOGEE DR16 \(+\) Gaia found that GS/E has a shallow density profile in the inner Galaxy, a break between \(15\) and \(25\) kpc, and a triaxial shape with axis ratios \(1:0.55:0.45\). The major axis is oriented about \(80^\circ\) from the Sun–Galactic centre line and \(16^\circ\) above the plane. The inferred stellar mass was
\[
1.45^{+0.92}_{-0.51}\,{\rm(stat.)}^{+0.13}_{-0.37}\,{\rm(sys.)}\times10^8\,M_\odot,
\]
implying that GS/E could make up as little as \(15\)–\(25\%\) of the Milky Way stellar halo [2306.03084]. A later re-analysis using a radially varying anisotropy profile increased that mass estimate to
\[
2.29^{+0.95}_{-0.63}\times10^8\,M_\odot,
\]
while retaining GS/E on the low end of the literature range [2509.04557].

The globular-cluster connection is especially important. In a joint analysis of stars and probable GCs, APOGEE observations showed no evidence for atypical \({\rm [Fe/H]}\) spreads among probable GSE globular clusters, with the exception of \(\omega\) Centauri. Under the assumption that \(\omega\) Cen is a stripped nuclear star cluster, the stellar mass of its progenitor was estimated as \(M_\star \approx 1.3\times10^9\,M_\odot\), described as well within literature expectations for GSE, leading to the proposal that GSE is the best available candidate for the original host galaxy of \(\omega\) Cen [2206.10505].

Internal chemical structure is another active line of reconstruction. Auriga analogues of GES-like progenitors show negative stellar metallicity gradients at infall, ranging from \(-0.09\) to \(-0.03\) dex/kpc against radius, which become much shallower by \(z=0\), with percentage changes of \(78\)–\(98\%\) in radial metallicity gradient. This suggests that present-day debris may retain only a blurred version of the progenitor’s original chemical cartography [2509.24705].

## 5. Dynamical and chemo-dynamical impact on the Milky Way

Simulations consistently assign a major dynamical role to the GSE merger, but the detailed mechanism differs by observable. In Auriga cosmological MHD simulations, the GES merger is likely to have been gas-rich and to contribute \(10\)–\(50\%\) of the gas to a merger-induced centrally concentrated starburst. That event rapidly formed a compact, rotationally supported thick disc and heated the proto-disc, with local \(\sigma_z\) increases of \(\Delta \sigma_z \simeq 20\)–\(60\ {\rm km\,s^{-1}}\). The same simulations linked the merger mass and timing to the fraction and age distribution of local counter-rotating stars [2001.06009].

A separate Auriga analysis examined bar formation. In the Au-18 realization, a GES-like merger with peak stellar mass \(M_{\star,\rm sat}\approx5\times10^8\,M_\odot\) reached first pericentre at \(t_{\rm lookback}\approx8.9\) Gyr and final coalescence at \(t_{\rm lookback}\approx8.5\) Gyr. Bar seeds appeared immediately after first pericentre, and a permanent bar was in place by \(t_{\rm lookback}\approx7.7\) Gyr, whereas removing the GES-18 progenitor delayed stable bar formation by \(\gtrsim 2\) Gyr [2312.02318].

The Galactic warp has likewise been modeled as a GSE consequence. A GIZMO gas-rich merger with host virial mass \(5\times10^{11}\,M_\odot\), satellite virial mass \(2\times10^{11}\,M_\odot\), \(e=0.75\), and retrograde bias produced an S-shaped warp in both the stellar and H I disks after the merger. The fitted onset radius was \(R_w\approx11\) kpc, the best-fit slope was \(a\approx0.17\,{\rm kpc\,kpc^{-1}}\), and the young-star warp amplitude at \(R\approx16\) kpc reached \(2.0\)–\(2.5\) kpc. In that simulation the warp persisted from \(\sim3.25\) Gyr post-merger to \(9\) Gyr, i.e. for \(>5\) Gyr [2409.03264].

Chemo-dynamical models further propose that the retrograde or radial nature of the merger trajectory matters for disk abundance bimodality. In a numerical GCE framework with radial gas flows, a retrograde GSE merger near \(z\sim2\) was found to drive a major sinking event, lower gas surface density across much of the disk, and accelerate the decline in \([\alpha/{\rm Fe}]\) through continuing SN Ia enrichment. In that scenario, retrograde trajectories increase the number of low-\([\alpha/{\rm Fe}]\) stellar populations relative to prograde trajectories and can help generate the observed bimodality when a preexisting high-\(\alpha\) population is already substantial [2510.08688].

## 6. Debates, composite scenarios, and current synthesis

The classical picture treats GSE as the debris of a single ancient merger, but recent work has raised two distinct complications: multiple stripping episodes within one progenitor, and genuinely composite origins involving more than one progenitor.

Evidence for a multi-passage history comes from APOGEE DR17 analyses that separate GES into a low-energy population of \(325\) field stars and a high-energy population of \(445\) field stars. The low-energy component is chemically more enriched, the high-energy component more metal-poor, and the difference is reproduced by inside-out chemical-evolution models with a former metallicity gradient of \(d{\rm [Fe/H]}/dr \simeq -0.08\ {\rm dex\,kpc^{-1}}\). This work explicitly interprets GES as a multi-passage event through the Milky Way disc, with a first passage \(11\)–\(12\) Gyr ago stripping outer stars and a second around \(10\) Gyr ago stripping inner stars [2602.00228].

A stronger revision appears in a DESI-based unsupervised clustering study using GS\(^3\) Hunter. That analysis identified \(17\) halo structures and argued that the GSE region itself contains four distinct substructures: GSE-GSH1 (\(12\) Gyr), GSE-GSH2 (\(10\) Gyr), GSE-GSH3 (\(8\) Gyr), and GSE-GSH4 (\(7\) Gyr). Although all four broadly match the overall GSE phase-space distribution and abundance patterns, they display distinct orbital actions and chemical abundances. The authors therefore argue that GSE is not the remnant of a single accretion event, but rather a composite structure assembled through multiple, sequential merger episodes during the early Milky Way [2606.04462].

Cosmological simulations also weaken any presumption that high radial anisotropy implies a single progenitor. In IllustrisTNG50, GSE-like debris was found in \(32\) of \(98\) Milky Way analogues, and two-merger GSEs accounted for a third of those cases. Kinematics alone could not robustly distinguish single-merger from two-merger GSEs; chemistry and star-formation histories were more informative, because single-merger GSEs tended to be more metal-rich and younger, with median infall times \(5.9^{+3.3}_{-2.0}\) Gyr ago, while the galaxies in two-merger GSEs had median infall times around \(10.7^{+1.2}_{-3.7}\) Gyr ago [2408.02723].

A further complication is overlap with other accreted structures. In the APOGEE DR17 \(+\) Gaia EDR3 reconstruction, stars assigned to Sequoia essentially overlap the GSE footprint in all analyzed chemical-abundance spaces, but present lower metallicities [2206.10505]. This makes chemical separation non-trivial and reinforces the broader conclusion that neither kinematics alone nor chemistry alone is sufficient in every regime.

Taken together, the current literature supports a stable core picture—an old, radially biased, low-\(L_z\), inner-halo-dominant accretion remnant with a major role in Milky Way evolution—while leaving open whether the observed GSE locus is the debris of one progenitor, one progenitor stripped in multiple passages, or a chemically and chronologically composite structure. The newer results suggest that “Gaia-Enceladus/Sausage” may function less as a single-object label than as the name for the dominant region of phase space in which several early accretion histories overlap.

Source: https://www.emergentmind.com/topics/gaia-enceladus-sausage-ges-4d4a66cf-ad59-4db2-add2-9c3f7ac0d0fa