Metallicity Gradient–Age Relation (MGAR) in Galaxies

• MGAR is the relation linking metallicity gradients to stellar ages across different galaxy types, serving as a probe of formation and evolution.
• Studies show that in dwarfs and the Milky Way, older stellar populations display steeper negative metallicity gradients, with specific correlations quantified in simulations and surveys.
• Research employs IFU surveys, cosmological simulations, and spectral analyses to highlight the roles of feedback-driven migration, disk settling, and environmental quenching.

Searching arXiv for recent and foundational papers on metallicity gradient–age relations across dwarf galaxies, the Milky Way, and IFU surveys. The Metallicity Gradient–Age Relation (MGAR) is the family of relations that connect a metallicity gradient to stellar age, an age distribution, or the age structure of a population. In the current literature, the term is used in several closely related senses: in dwarf-galaxy work it usually denotes a correlation between the strength of the stellar metallicity gradient and a galaxy-wide median or light-weighted stellar age; in Milky Way work it usually denotes the age dependence of radial or vertical metallicity gradients measured for mono-age populations; and in spatially resolved IFU studies of external galaxies it often refers to tests for a direct correlation between metallicity gradients and age gradients within galaxies. Across these usages, MGAR functions as a probe of feedback-driven migration, disk flaring, upside-down disk settling, radial migration, baryon cycling, and environmental quenching (Mercado et al., 2020, Ciucă et al., 2017, Zheng et al., 2016).

1. Definitions, observables, and measurement conventions

The quantity called a “metallicity gradient” is not uniform across the literature. In dwarf-galaxy studies it is commonly the slope of a linear fit to median stellar metallicity versus projected radius, often within $2R_{1/2}$, with $\gamma_z$ reported in dex per $R_{1/2}$; in nearby-dwarf IFU work it is written as $\nabla[Z/H]=d[Z/H]/d(R/R_e)$ in dex per effective radius; in Milky Way vertical studies it is $d[\mathrm{M}/\mathrm{H}]/d|Z|$ or $d[\mathrm{Fe/H}]/d|Z|$ in dex kpc$^{-1}$; and in Milky Way radial studies it is usually $d[\mathrm{Fe/H}]/dR_g$, with $R_g$ the guiding-center radius, again in dex kpc$^{-1}$. IFU studies of external galaxies also report $\gamma_z$0 or $\gamma_z$1 per $\gamma_z$2, while age is represented either by global median stellar age, light-weighted age, mass-weighted age, or discrete mono-age bins (Mercado et al., 2020, Li et al., 15 Jun 2026, Ciucă et al., 2017, Vickers et al., 2021, Zheng et al., 2016).

Context	Gradient definition	Age variable
Dwarf galaxies	$\gamma_z$3, $\gamma_z$4	$\gamma_z$5, light-weighted age
MW vertical structure	$\gamma_z$6, $\gamma_z$7	mono-age bins, $\gamma_z$8 proxy
MW radial structure	$\gamma_z$9, $R_{1/2}$0	age bins
IFU stellar populations	$R_{1/2}$1, $R_{1/2}$2 per $R_{1/2}$3	age gradients, light/mass-weighted age

The age variable is equally heterogeneous. Some studies use galaxy-wide median stellar age $R_{1/2}$4, others use light-weighted stellar age, and Galactic archaeology studies isolate narrow age bins or use mono-temperature populations as an age proxy for young stars. This suggests that direct numerical comparison across papers requires care, because different MGARs need not trace the same physical quantity even when the phrase “metallicity gradient–age relation” is used.

2. Dwarf galaxies: the clearest direct MGAR

The strongest direct evidence for an MGAR comes from dwarf galaxies. In FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies with stellar masses between $R_{1/2}$5 and $R_{1/2}$6, combined with observational data for 10 Local Group dwarfs, stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is tightly correlated with galaxy-wide median stellar age: older galaxies have stronger negative gradients and younger galaxies have flatter gradients. The best-fit relations are

$R_{1/2}$7

for simulated galaxies and

$R_{1/2}$8

for observed Local Group dwarfs, with a Spearman correlation coefficient $R_{1/2}$9 in the simulation sample. The same work reports that gas-phase metallicity gradients are uniformly weak and show no correlation with stellar metallicity gradients, and it identifies no clear dependence of gradient strength on galaxy mass, size, rotation $\nabla[Z/H]=d[Z/H]/d(R/R_e)$0, or morphology, with age as the primary correlator (Mercado et al., 2020).

That framework assigns the dwarf-galaxy MGAR to two competing internal mechanisms. The first is the steady “puffing” of old, metal-poor stars by feedback-driven potential fluctuations, which moves older populations outward and establishes a negative stellar metallicity gradient. The second is the accretion of extended, metal-rich gas at late times, which fuels extended metal-rich star formation and washes out pre-existing gradients. In that picture, older galaxies are those in which puffing dominates, whereas younger galaxies are those in which late-time, spatially extended star formation has flattened the gradient (Mercado et al., 2020).

A larger observational IFU study of 90 nearby low-mass galaxies with VLT/MUSE, spanning $\nabla[Z/H]=d[Z/H]/d(R/R_e)$1 to $\nabla[Z/H]=d[Z/H]/d(R/R_e)$2, recovered a robust negative correlation between $\nabla[Z/H]=d[Z/H]/d(R/R_e)$3 and light-weighted stellar age measured out to $\nabla[Z/H]=d[Z/H]/d(R/R_e)$4: older dwarf galaxies have steeper, more negative gradients. For light-weighted age the Pearson coefficient is $\nabla[Z/H]=d[Z/H]/d(R/R_e)$5, with best-fit slope $\nabla[Z/H]=d[Z/H]/d(R/R_e)$6 at $\nabla[Z/H]=d[Z/H]/d(R/R_e)$7; the correlation persists after controlling for stellar mass, effective radius, surface brightness, and large-scale environment, and is strongest in the intermediate-mass regime $\nabla[Z/H]=d[Z/H]/d(R/R_e)$8. The same study, however, finds that the H I deficiency parameter is the second-strongest correlate after stellar age, with $\nabla[Z/H]=d[Z/H]/d(R/R_e)$9 and $d[\mathrm{M}/\mathrm{H}]/d|Z|$0, such that galaxies with higher H I deficiency tend to have more negative gradients. It therefore argues for a synergy of feedback-driven dynamical heating and environment-driven outside-in quenching, while noting that only the latter has robust observational support (Li et al., 15 Jun 2026).

The Small Magellanic Cloud supplies a related but conceptually distinct case. A Ca II triplet metallicity study of 3037 field red giant stars finds a median metallicity of $d[\mathrm{M}/\mathrm{H}]/d|Z|$1 and a radial gradient of $d[\mathrm{M}/\mathrm{H}]/d|Z|$2 dex deg$d[\mathrm{M}/\mathrm{H}]/d|Z|$3 over the inner $d[\mathrm{M}/\mathrm{H}]/d|Z|$4. That gradient is interpreted as the result of an increasing fraction of young stars with decreasing galactocentric radius, coupled with a uniform global age–metallicity relation. In other words, the abundance gradient is attributed to the spatial mixing of populations with different ages, not to a spatial change in the age–metallicity relation itself (Dobbie et al., 2014).

3. Vertical MGAR in the Milky Way

In the Milky Way, one major form of MGAR is the dependence of the vertical metallicity gradient on stellar age. Using 18,435 dwarf stars from TGAS and RAVE DR5, divided into five mono-age populations between 0 and 11 Gyr, a hierarchical regression analysis finds no vertical metallicity gradient in the youngest population and increasingly steeper negative vertical metallicity gradients in older populations. The measured $d[\mathrm{M}/\mathrm{H}]/d|Z|$5 values are as follows (Ciucă et al., 2017):

Age bin (Gyr)	Unweighted gradient	Weighted gradient
0–2	$d[\mathrm{M}/\mathrm{H}]/d|Z|$6	$d[\mathrm{M}/\mathrm{H}]/d|Z|$7
2–4	$d[\mathrm{M}/\mathrm{H}]/d|Z|$8	$d[\mathrm{M}/\mathrm{H}]/d|Z|$9
4–6	$d[\mathrm{Fe/H}]/d|Z|$0	$d[\mathrm{Fe/H}]/d|Z|$1
6–8	$d[\mathrm{Fe/H}]/d|Z|$2	$d[\mathrm{Fe/H}]/d|Z|$3
8–11	$d[\mathrm{Fe/H}]/d|Z|$4	$d[\mathrm{Fe/H}]/d|Z|$5

The same analysis finds that the metallicity at the disc plane remains almost constant between 2 and 8 Gyr and becomes significantly lower for the $d[\mathrm{Fe/H}]/d|Z|$6 Gyr population, while the intrinsic metallicity dispersion increases steadily with age. The interpretation advanced there is a flaring thin star-forming disc, in which mono-age populations form without vertical gradients and later acquire a negative vertical MGAR through radial mixing in a flaring geometry (Ciucă et al., 2017).

A 2025 LAMOST analysis of mono-temperature populations reaches a complementary conclusion for the youngest Galactic populations. After correcting $d[\mathrm{Fe/H}]/d|Z|$7 for $d[\mathrm{Fe/H}]/d|Z|$8-dependent systematics and subtracting the radial metallicity gradient, it finds that the vertical metallicity gradient approaches zero as stellar effective temperature increases, or equivalently as age decreases, across various Galactocentric distances. In that work, broad-age G dwarfs show steeply negative gradients from $d[\mathrm{Fe/H}]/d|Z|$9 to $^{-1}$0 dex kpc$^{-1}$1, whereas the hottest OBA stars have slopes of approximately $^{-1}$2 dex kpc$^{-1}$3. Independent validation with 295 open clusters younger than 3 Gyr and 976 classical Cepheids shows that, within a given narrow age range, the vertical metallicity gradient is effectively zero. This is taken as support for upside-down disk formation, in which the youngest and most metal-rich stars dominate the midplane while older and more metal-poor stars formed at larger vertical heights and currently tend to be found there (Long et al., 2 Apr 2025).

Taken together, these studies imply that the vertical MGAR depends strongly on whether one measures genuinely mono-age populations or mixed-age samples. The steep negative vertical gradient of a broad stellar sample need not represent the birth structure of any single cohort; it can emerge from superposing cohorts with different ages, scale heights, and metallicities.

4. Radial MGAR in the Milky Way

A second Galactic usage of MGAR concerns the age dependence of the radial metallicity gradient. Using $^{-1}$4 million LAMOST and Gaia stars, with a working sample of 1.3 million, one study measures $^{-1}$5 in age bins and finds that thin-disk selections yield radial metallicity gradients that grow shallower for the oldest stars. For the youngest thin-disk stars the gradient is approximately $^{-1}$6 dex kpc$^{-1}$7, and the temporal flattening is approximately $^{-1}$8 dex kpc$^{-1}$9 Gyr$d[\mathrm{Fe/H}]/dR_g$0, so that gradients approach zero at ages $d[\mathrm{Fe/H}]/dR_g$1–10 Gyr. In contrast, thick-disk selections reveal a slightly positive radial metallicity gradient of $d[\mathrm{Fe/H}]/dR_g$2 dex kpc$d[\mathrm{Fe/H}]/dR_g$3, similar in magnitude at all ages. The same study argues that the trend is strongest for very small orbital $d[\mathrm{Fe/H}]/dR_g$4, consistent with a provenance bias in which churning is most efficient for vertically cold orbits (Vickers et al., 2021).

A later re-analysis with LAMOST DR8 and asteroseismologically calibrated ages reports the same basic pattern while modifying the quantitative rate. In that study the thin disk again shows a steadily flattening MGAR, with the youngest thin-disk stars at $d[\mathrm{Fe/H}]/dR_g$5 dex kpc$d[\mathrm{Fe/H}]/dR_g$6, the oldest thin-disk stars nearly flat, and a flattening rate of $d[\mathrm{Fe/H}]/dR_g$7 dex kpc$d[\mathrm{Fe/H}]/dR_g$8 Gyr$d[\mathrm{Fe/H}]/dR_g$9. The thick disk is found to have a mildly positive gradient of $R_g$0 dex kpc$R_g$1, roughly constant with age. Two explanations are left open: large-scale radial migration induced by transient spiral arms, or a time-dependent steepening of the interstellar-medium metallicity gradient (Chen et al., 28 Aug 2025).

Older Gaia-ESO results already indicated that radial metallicity gradients vary with age and height above the plane. For young stars with ages $R_g$2 Gyr and $R_g$3 pc, $R_g$4 dex kpc$R_g$5; for the vertical bin $R_g$6 pc, the gradient is $R_g$7 dex kpc$R_g$8; and for stars older than 12 Gyr the measured gradient is $R_g$9, formally positive but with large uncertainty. In the same sample, $^{-1}$0 gradients become more negative with increasing age and height (Bergemann et al., 2014).

Other Milky Way work emphasizes the youngest populations as tracers of the present-day gradient. Using APOGEE-Gaia DR2 and the youngest mean-age bin at each radial zone in the midplane, one study measures a present-day metallicity gradient of $^{-1}$1 dex kpc$^{-1}$2, in agreement with Cepheids and young field stars (Feuillet et al., 2019).

Not all Galactic analyses recover a strongly age-dependent radial MGAR. An APOGEE red-giant study with machine-learning ages finds that stars up to $^{-1}$3 Gyr follow a similar relation between metallicity and Galactocentric radius, with a slope in $^{-1}$4 generally near $^{-1}$5 dex kpc$^{-1}$6 and only minor fluctuations among mono-age bins. That result motivates an equilibrium interpretation in which the gas-phase gradient reaches a nearly constant normalization early in the disk lifetime, rather than a picture in which old stars at fixed radius are systematically more metal-poor (Johnson et al., 2024).

5. External-galaxy surveys and the non-universality of direct gradient–gradient coupling

Outside the Local Group dwarfs and the Milky Way, the direct coupling between age and metallicity gradients is not universal. A long-slit spectroscopic study of 40 early-type galaxies spanning $^{-1}$7 to $^{-1}$8 finds that the mean metallicity gradients are approximately $^{-1}$9 and the mean age gradients $\gamma_z$00 dex per decade of radius, with wide object-to-object spread. However, it finds no evidence for a correlation between metallicity gradient and luminosity, velocity dispersion, central age, or age gradient, and likewise no correlation between age gradient and other parameters in bright early-type galaxies. Only faint early-types with $\gamma_z$01 mag show a correlation between age gradient and luminosity, in the sense that the age gradient becomes more positive for fainter galaxies (Koleva et al., 2011).

A MaNGA analysis of 1105 galaxies reaches a similar conclusion at larger sample size. Using STARLIGHT on IFU spectroscopy, it finds mean age and metallicity gradients close to zero but slightly negative, with mean mass-weighted metallicity gradients of $\gamma_z$02 dex $\gamma_z$03 for disks and $\gamma_z$04 dex $\gamma_z$05 for ellipticals, and mean mass-weighted age gradients of $\gamma_z$06 dex $\gamma_z$07 for disks and $\gamma_z$08 dex $\gamma_z$09 for ellipticals. Yet it does not identify a significant statistical MGAR and finds weak or no statistically significant dependence of these gradients on large-scale structure type or local density (Zheng et al., 2016).

A second MaNGA study of $\gamma_z$10 galaxies emphasizes morphology. Early-type galaxies have flat age gradients and negative metallicity gradients that steepen modestly with mass, whereas late-type galaxies have negative age and negative metallicity gradients, with age gradients steepening strongly with increasing galaxy mass and metallicity gradients remaining nearly constant with mass. That work explicitly reports no strong direct correlation within a given galaxy between age and metallicity gradients, implying that the two gradients respond to different physical drivers in early- and late-type systems (Parikh et al., 2021).

These surveys therefore suggest that a direct, one-parameter MGAR is not generic across the broader galaxy population. In many massive-galaxy samples, age and metallicity gradients are only weakly coupled or effectively decoupled.

6. Physical interpretations, competing frameworks, and open issues

Several mechanisms recur in MGAR studies. In dwarfs, feedback-driven potential fluctuations can “puff” old, metal-poor stars to large radii, while late-time accretion of extended, metal-rich gas can flatten or erase older gradients; this produces an age–gradient sequence without requiring environmental processes (Mercado et al., 2020). More recent nearby-dwarf IFU work accepts a role for feedback-driven migration but places stronger observational weight on environment-driven outside-in quenching, because H I deficiency correlates significantly with $\gamma_z$11 and many dwarfs show positive age gradients, with older outer regions (Li et al., 15 Jun 2026).

In the Milky Way thin disk, one explanatory line centers on churning. The steady flattening of the radial metallicity gradient with age, especially for vertically cold orbits, is interpreted as a hallmark feature of radial migration (Vickers et al., 2021). A later study preserves that possibility but also admits an alternative in which the observed MGAR mainly reflects the time dependence of the ISM metallicity gradient itself rather than post-birth migration (Chen et al., 28 Aug 2025).

Vertical MGAR studies emphasize disk geometry and settling. The 2017 mono-age analysis argues for a flaring thin star-forming disk, whereas the 2025 LAMOST work argues that young mono-age populations have essentially zero vertical gradient and that steep gradients arise from age mixing in an upside-down disk-formation sequence (Ciucă et al., 2017, Long et al., 2 Apr 2025). A related chemodynamical simulation of disk stars shows how an age–metallicity relation, an age–velocity-dispersion relation, and disk flaring can reverse the sign of the radial metallicity gradient at large heights above the plane: the simulated gradient changes from $\gamma_z$12 dex kpc$\gamma_z$13 at $\gamma_z$14–1 kpc to $\gamma_z$15 dex kpc$\gamma_z$16 at $\gamma_z$17–3 kpc because younger, more metal-rich, lower-dispersion stars reach large $\gamma_z$18 more readily in the outer disk (Rahimi et al., 2013).

A different framework challenges strong age evolution of the Milky Way radial MGAR altogether. The equilibrium model proposed for APOGEE red giants treats the radial metallicity gradient as the outcome of a local balance among metal production, accretion, and outflow, with a nearly age-invariant normalization out to $\gamma_z$19 Gyr. In that model,

$\gamma_z$20

and the radial slope is

$\gamma_z$21

Major perturbations such as mergers can distort the gradient temporarily, but the system relaxes back toward equilibrium on $\gamma_z$22Gyr timescales (Johnson et al., 2024).

The literature on star-forming regions extends the same logic to assembly time. In MaNGA and EAGLE, strong gas-phase metallicity gradients are associated with recent accretion or merger activity, and galaxies gradually evolve back toward weak gradients and the median mass–metallicity relation over $\gamma_z$23 Gyr. This suggests that, in gas as well as stars, metallicity gradients can act as a fossil record of recent assembly events (Jara-Ferreira et al., 2024).

Taken together, these results suggest that MGAR is not a single universal law. Its sign, slope, and physical interpretation depend on whether the tracer is stellar or gas-phase, whether the geometry is radial or vertical, whether age is defined for a galaxy as a whole or for mono-age cohorts, and whether the dominant mechanism is feedback-driven migration, disk settling, equilibrium enrichment, or environmental stripping. The term is therefore best understood as a framework for relating chemical structure to temporal structure, rather than as one invariant correlation across all galaxy classes.