2025-04-08 | | Total: 4
Marine biogeochemistry models are critical for forecasting, as well as estimating ecosystem responses to climate change and human activities. Data assimilation (DA) improves these models by aligning them with real-world observations, but marine biogeochemistry DA faces challenges due to model complexity, strong nonlinearity, and sparse, uncertain observations. Existing DA methods applied to marine biogeochemistry struggle to update unobserved variables effectively, while ensemble-based methods are computationally too expensive for high-complexity marine biogeochemistry models. This study demonstrates how machine learning (ML) can improve marine biogeochemistry DA by learning statistical relationships between observed and unobserved variables. We integrate ML-driven balancing schemes into a 1D prototype of a system used to forecast marine biogeochemistry in the North-West European Shelf seas. ML is applied to predict (i) state-dependent correlations from free-run ensembles and (ii), in an ``end-to-end'' fashion, analysis increments from an Ensemble Kalman Filter. Our results show that ML significantly enhances updates for previously not-updated variables when compared to univariate schemes akin to those used operationally. Furthermore, ML models exhibit moderate transferability to new locations, a crucial step toward scaling these methods to 3D operational systems. We conclude that ML offers a clear pathway to overcome current computational bottlenecks in marine biogeochemistry DA and that refining transferability, optimizing training data sampling, and evaluating scalability for large-scale marine forecasting, should be future research priorities.
We present a unified variational treatment of the magnetic dipole matrix elements, Einstein coefficients and line strength for general open-shell diatomic molecules in the general purpose diatomic code Duo. Building on previous work in which similar expressions for the electric quadrupole transitions were developed, we also present a complete ab initio spectroscopic model for the infrared, electric dipole-forbidden, spectrum of the 16O_2 molecule. The model covers seven states, namely the X ^3\Sigma^-_g, a ^1\Delta_g, b ^1\Sigma^+_g, I ^1\Pi_g, II ^1\Pi_g, I ^3\Pi_g and II ^3\Pi_g states, for which 7 potential energy, 6 electronic angular momentum, 7 spin-orbit, and 14 quadrupole moment curves are calculated using ic-MRCI theory and an aug-cc-pV5Z basis set. These curves are diabatised to remove avoided crossings between the excited \Pi states, and the resultant properties are used to produce a line list for higher-order transitions of astrophysical interest.
Decadal-scale interactions between the tropical Atlantic and Pacific Oceans play a crucial role in global climate variability through bidirectional teleconnections. Current climate models show persistent biases in representing these basin interactions, particularly in simulating the Atlantic's influence on Pacific climate. Using historical simulations from 27 CMIP6 models and two configurations of the Energy Exascale Earth System Model (E3SM) during 1950-2015, we systematically evaluate tropical Atlantic-Pacific teleconnections through both Walker circulation and extratropical wave responses. Most models exhibit Pacific-dominated teleconnections, contradicting observational evidence of Atlantic control on Pacific variability during the past 40 years. By developing a performance metric that combines tropical circulation patterns and extratropical wave propagation, we identify two distinct model behaviors: high-skill models capture the bidirectional Atlantic-Pacific teleconnections with a secondary symptom of systematic 20-degree westward shifts in convective centers, while low-skill models display amplified Pacific dominance through reversed Walker circulation responses warming in both tropical basins. Comparative analysis between standard E3SMv2 and its multi-scale modeling framework configuration demonstrates that implementing more sophisticated cloud-scale processes alone, with limited model tuning, cannot resolve these teleconnection biases. Our results identify four CMIP6 models and E3SMv2 that effectively reproduce observed teleconnection pathways, offering a comprehensive diagnostic framework for evaluating decadal Atlantic-Pacific interactions in climate models.
A hothouse climate may develop throughout Earth's history and its warming future and on potentially habitable exoplanets near the inner edge of the habitable zone. Previous studies suggested that near-surface atmospheric inversion (NAIV) with planetary boundary air temperature being higher than the air temperature adjacent to the surface, is a pronounced phenomenon in hothouse climates. However, the underlying mechanisms are unclear. Here we show that lower-tropospheric radiative heating is necessary but not independently sufficient in forming the NAIV. Instead, the dynamic heating induced by large-scale subsidence is essential. With the prescribed reasonable large-scale subsidence, NAIV appears in small-domain cloud-resolving simulations, which was not observed in previous studies. Surface evaporative cooling also contributes to the formation of the NAIV. Besides NAIV, we find that surface inversion (SIV) with the air adjacent to the surface being warmer than the underlying sea surface is also a distinct phenomenon in hothouse climates. SIV is caused by strong surface evaporative cooling and large atmospheric shortwave absorption. These two types of inversion strongly stabilize the atmosphere, weaken atmospheric circulation, dry the free troposphere, and suppress the hydrological cycle.