Dissolved Oxygen (DO) plays a key role in reactive processes and microbial dynamics in the critical zone. While the general view is that oxygen is rapidly depleted in soils and that deeper compartments are anoxic, recent observations showed that fractures can provide rapid pathways for deep oxygen penetration, triggering unexpected biogeochemical processes. As it is transported in the subsurface, DO reacts with electron donors, such as $Feˆ{2+}$ coming from mineral dissolution, hence influencing rockweathering. Yet, little is known about the factors controlling the spatial heterogeneity and distribution of oxygen with depth. Here we present analytical expressions describing the coupled evolution of DO and $Feˆ{2+}$ as a function of fluid travel time in crystalline rocks. Our model, validated with reactive transport simulations, predicts a linear decay of DO with time, followed by a rapid non-linear increase of $Feˆ{2+}$ concentrations up to an equilibrium state. Relative effects of the reducing capacity of the bedrock and of transport velocity are quantified through a Damkohler number, capturing key hydrological and geological