We apply a combined study in order to investigate the dynamics of cosmological models incorporating nonlinear electrodynamics (NLED). The study is based on the simultaneous investigation of such fundamental aspects as stability and causality, complemented with a dynamical systems investigation of the involved models, as well as Bayesian inference for parameter estimation. We explore two specific NLED models: the power-law and the rational Lagrangian. We present the theoretical framework of NLED coupled with general relativity, followed by an analysis of the stability and causality of the various NLED Lagrangians. We then perform a detailed dynamical analysis to identify the ranges where these models are stable and causal. Our results show that the power-law Lagrangian model transitions through various cosmological phases, evolving from a Maxwell radiation-dominated state to a matter-dominated state. For the rational Lagrangian model, including the Maxwell term, stable and causal behavior is observed within specific parameter ranges, with critical points indicating the evolutionary pathways of the universe. To validate our theoretical findings, we perform Bayesian parameter estimation using a comprehensive set of observational data, including cosmic chronometers, baryon acoustic oscillation (BAO) measurements, and supernovae type Ia (SNeIa). The estimated parameters for both models align with the expected values for the current universe, particularly the matter density ${\Omega }_m$ and the Hubble parameter h. However, the parameters of the models are not tightly constrained within the prior ranges. Our combined studies approach rules out the mentioned models as an appropriate description of the cosmos. Our results highlight the need for further refinement and exploration of NLED-based cosmological models to fully integrate them into the standard cosmological framework.

This study examines the spatio-temporal evolution of precipitation in the State of Guanajuato, Mexico, from 1981 to 2016 by analyzing monthly series from 65 meteorological stations. A rigorous data quality protocol was implemented, selecting stations with more than 30 years of continuous data and less than 10% missing values. Multiple Imputation by Chained Equations (MICE) with Predictive Mean Matching was applied to handle missing data, preserving the statistical properties of the time series as validated by Kolmogorov–Smirnov tests ($p=1.000$ for all stations). Homogeneity was assessed using Pettitt, SNHT, Buishand, and von Neumann tests, classifying 60 stations (93.8%) as useful, 3 (4.7%) as doubtful, and 2 (3.1%) as suspicious for monthly analysis. Breakpoints were predominantly clustered around periods of instrumental changes (2000–2003 and 2011–2014), underscoring the necessity of homogenization prior to trend analysis. The Trend-Free Pre-Whitening Mann–Kendall (TFPW-MK) test was applied to account for significant first-order autocorrelation (${\rho }_1 > 0.3$) present in all series. The analysis revealed no statistically significant monotonic trends in monthly precipitation at any of the 65 stations ($\alpha =0.05$). While 75.4% of the stations showed slight non-significant increasing tendencies (Kendall’s $\tau$ range: 0.0016 to 0.0520) and 24.6% showed non-significant decreasing tendencies ($\tau$ range: −0.0377 to −0.0008), Sen’s slope estimates were negligible (range: −0.0029 to 0.0111 mm/year) and statistically indistinguishable from zero. No discernible spatial patterns or correlation between trend magnitude and altitude ($\rho =-0.022$, $p>0.05$) were found, indicating region-wide precipitation stability during the study period. The integration of advanced imputation, multi-test homogenization, and robust trend detection provides a comprehensive framework for hydroclimatic analysis in semi-arid regions. These findings suggest that Guanajuato’s severe water crisis cannot be attributed to declining precipitation but rather to anthropogenic factors, primarily unsustainable groundwater extraction for agriculture.

In this paper we perform a thorough dynamical systems analysis of the cubic galileon model non-minimally coupled to the dark matter. Three well-known classes of interacting models are considered where the energy exchange between the dark components is a function of the dark matter density and of the dark energy density: $Q_1=3\alpha H\rho_m$, $Q_2=3\beta\rho_m\dot{\phi}$ and $Q_3=3 \epsilon H \rho_\phi$, respectively. We are able to show the global asymptotic dynamics of the model for the exponential potential in a homogeneous and isotropic background. The cosmological implications of the proposed scenarios are explored and it is found that, in addition to the appearance of new equilibrium configurations that do not appear neither in the non-interacting cubic galileon model nor in the interacting quintessence model, there is a significant impact of the non-minimal coupling through modification of the stability properties of the critical points. The resulting cosmological scenario provides a bigbang origin of the cosmic expansion, an early transient stage of inflationary expansion, as well as matter-scaling late time stable state of the universe, among other solutions of lesser cosmological interest. This work extends previous studies of coupled dark energy to a broader class of gravitational theories.