More Like this

1. Constraining Earth’s nonlinear mantle viscosity using plate-boundary resolving global inversions

   https://doi.org/10.1073/pnas.2318706121  

   Hu, Jiashun; Rudi, Johann; Gurnis, Michael; Stadler, Georg (July 2024, Proceedings of the National Academy of Sciences)

   Variable viscosity in Earth’s mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 $\pm$0.25 (mean $\pm$SD). A relatively low yield stress of 151 $\pm$19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes. 

   more » « less
2. Latitudinal gradient in the respiration quotient and the implications for ocean oxygen availability

   https://doi.org/10.1073/pnas.2004986117  

   Moreno, Allison R.; Garcia, Catherine A.; Larkin, Alyse A.; Lee, Jenna A.; Wang, Wei-Lei; Moore, J. Keith; Primeau, Francois W.; Martiny, Adam C. (August 2020, Proceedings of the National Academy of Sciences)

   Climate-driven depletion of ocean oxygen strongly impacts the global cycles of carbon and nutrients as well as the survival of many animal species. One of the main uncertainties in predicting changes to marine oxygen levels is the regulation of the biological respiration demand associated with the biological pump. Derived from the Redfield ratio, the molar ratio of oxygen to organic carbon consumed during respiration (i.e., the respiration quotient, $r_{-O2:C}$) is consistently assumed constant but rarely, if ever, measured. Using a prognostic Earth system model, we show that a 0.1 increase in the respiration quotient from 1.0 leads to a 2.3% decline in global oxygen, a large expansion of low-oxygen zones, additional water column denitrification of 38 Tg N/y, and the loss of fixed nitrogen and carbon production in the ocean. We then present direct chemical measurements of $r_{-O2:C}$using a Pacific Ocean meridional transect crossing all major surface biome types. The observed $r_{-O2:C}$has a positive correlation with temperature, and regional mean values differ significantly from Redfield proportions. Finally, an independent global inverse model analysis constrained with nutrients, oxygen, and carbon concentrations supports a positive temperature dependence of $r_{-O2:C}$in exported organic matter. We provide evidence against the common assumption of a static biological link between the respiration of organic carbon and the consumption of oxygen. Furthermore, the model simulations suggest that a changing respiration quotient will impact multiple biogeochemical cycles and that future warming can lead to more intense deoxygenation than previously anticipated. 

   more » « less
3. Multipliers on bi-parameter Haar system Hardy spaces

   https://doi.org/10.1007/s00208-024-02887-9  

   Lechner, R.; Motakis, P.; Müller, P_F_X; Schlumprecht, Th (May 2024, Mathematische Annalen)

   Abstract Let$$(h_I)$$ $\left(h_I\right)$denote the standard Haar system on [0, 1], indexed by$$I\in \mathcal {D}$$ $I\in D$, the set of dyadic intervals and$$h_I\otimes h_J$$ $h_I\otimes h_J$denote the tensor product$$(s,t)\mapsto h_I(s) h_J(t)$$ $(s,t)\mapsto h_I\left(s\right)h_J\left(t\right)$,$$I,J\in \mathcal {D}$$ $I,J\in D$. We consider a class of two-parameter function spaces which are completions of the linear span$$\mathcal {V}(\delta ^2)$$ $V\left({\delta }^2\right)$of$$h_I\otimes h_J$$ $h_I\otimes h_J$,$$I,J\in \mathcal {D}$$ $I,J\in D$. This class contains all the spaces of the formX(Y), whereXandYare either the Lebesgue spaces$$L^p[0,1]$$ $L^p[0,1]$or the Hardy spaces$$H^p[0,1]$$ $H^p[0,1]$,$$1\le p < \infty $$ $1\le p<\infty$. We say that$$D:X(Y)\rightarrow X(Y)$$ $D:X\left(Y\right)\to X\left(Y\right)$is a Haar multiplier if$$D(h_I\otimes h_J) = d_{I,J} h_I\otimes h_J$$ $D(h_I\otimes h_J)=d_{I,J}{h}_I\otimes h_J$, where$$d_{I,J}\in \mathbb {R}$$ $d_{I,J}\in R$, and ask which more elementary operators factor throughD. A decisive role is played by theCapon projection$$\mathcal {C}:\mathcal {V}(\delta ^2)\rightarrow \mathcal {V}(\delta ^2)$$ $C:V\left({\delta }^2\right)\to V\left({\delta }^2\right)$given by$$\mathcal {C} h_I\otimes h_J = h_I\otimes h_J$$ $C{h}_I\otimes h_J=h_I\otimes h_J$if$$|I|\le |J|$$ $\left|I\right|\le \left|J\right|$, and$$\mathcal {C} h_I\otimes h_J = 0$$ $C{h}_I\otimes h_J=0$if$$|I| > |J|$$ $\left|I\right|>\left|J\right|$, as our main result highlights: Given any bounded Haar multiplier$$D:X(Y)\rightarrow X(Y)$$ $D:X\left(Y\right)\to X\left(Y\right)$, there exist$$\lambda ,\mu \in \mathbb {R}$$ $\lambda ,\mu \in R$such that$$\begin{aligned} \lambda \mathcal {C} + \mu ({{\,\textrm{Id}\,}}-\mathcal {C})\text { approximately 1-projectionally factors through }D, \end{aligned}$$ $\begin{array}{c}\lambda C+\mu (\text{Id}-C)\text{approximately 1-projectionally factors through}D,\end{array}$i.e., for all$$\eta > 0$$ $\eta >0$, there exist bounded operatorsA, Bso thatABis the identity operator$${{\,\textrm{Id}\,}}$$ $\text{Id}$,$$\Vert A\Vert \cdot \Vert B\Vert = 1$$ $\Vert A\Vert ·\Vert B\Vert =1$and$$\Vert \lambda \mathcal {C} + \mu ({{\,\textrm{Id}\,}}-\mathcal {C}) - ADB\Vert < \eta $$ $\Vert \lambda C+\mu (\text{Id}-C)-ADB\Vert <\eta$. Additionally, if$$\mathcal {C}$$ $C$is unbounded onX(Y), then$$\lambda = \mu $$ $\lambda =\mu$and then$${{\,\textrm{Id}\,}}$$ $\text{Id}$either factors throughDor$${{\,\textrm{Id}\,}}-D$$ $\text{Id}-D$. 

   more » « less
4. On the stationarity for nonlinear optimization problems with polyhedral constraints

   https://doi.org/10.1007/s10107-023-01979-9  

   di Serafino, Daniela; Hager, William W.; Toraldo, Gerardo; Viola, Marco (May 2023, Mathematical Programming)

   Abstract For polyhedral constrained optimization problems and a feasible point$$\textbf{x}$$ $x$, it is shown that the projection of the negative gradient on the tangent cone, denoted$$\nabla _\varOmega f(\textbf{x})$$ ${\nabla }_{\Omega }f\left(x\right)$, has an orthogonal decomposition of the form$$\varvec{\beta }(\textbf{x}) + \varvec{\varphi }(\textbf{x})$$ $\beta \left(x\right)+\phi \left(x\right)$. At a stationary point,$$\nabla _\varOmega f(\textbf{x}) = \textbf{0}$$ ${\nabla }_{\Omega }f\left(x\right)=0$so$$\Vert \nabla _\varOmega f(\textbf{x})\Vert $$ $\Vert {\nabla }_{\Omega }f\left(x\right)\Vert$reflects the distance to a stationary point. Away from a stationary point,$$\Vert \varvec{\beta }(\textbf{x})\Vert $$ $\Vert \beta \left(x\right)\Vert$and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$ $\Vert \phi \left(x\right)\Vert$measure different aspects of optimality since$$\varvec{\beta }(\textbf{x})$$ $\beta \left(x\right)$only vanishes when the KKT multipliers at$$\textbf{x}$$ $x$have the correct sign, while$$\varvec{\varphi }(\textbf{x})$$ $\phi \left(x\right)$only vanishes when$$\textbf{x}$$ $x$is a stationary point in the active manifold. As an application of the theory, an active set algorithm is developed for convex quadratic programs which adapts the flow of the algorithm based on a comparison between$$\Vert \varvec{\beta }(\textbf{x})\Vert $$ $\Vert \beta \left(x\right)\Vert$and$$\Vert \varvec{\varphi }(\textbf{x})\Vert $$ $\Vert \phi \left(x\right)\Vert$. 

   more » « less
5. Probing light scalars and vector-like quarks at the high-luminosity LHC

   https://doi.org/10.1140/epjc/s10052-025-14085-1  

   Qureshi, U S; Gurrola, A; Flórez, A; Rodriguez, C (April 2025, The European Physical Journal C)

   Abstract A model based on a$$U(1)_{T^3_R}$$ $U{\left(1\right)}_{T_R^3}$extension of the Standard Model can address the mass hierarchy between generations of fermions, explain thermal dark matter abundance, and the muon$$g - 2$$ $g-2$,$$R_{(D)}$$ $R_{\left(D\right)}$, and$$R_{(D^*)}$$ $R_{\left(D^{\ast }\right)}$anomalies. The model contains a light scalar boson$$\phi '$$ ${\varphi }^{\prime }$and a heavy vector-like quark$$\chi _\textrm{u}$$ ${\chi }_{\text{u}}$that can be probed at CERN’s large hadron collider (LHC). We perform a phenomenology study on the production of$$\phi '$$ ${\varphi }^{\prime }$and$${\chi }_u$$ ${\chi }_u$particles from proton–proton$$(\textrm{pp})$$ $\left(\text{pp}\right)$collisions at the LHC at$$\sqrt{s}=13.6$$ $\sqrt{s}=13.6$TeV, primarily through$$g{-g}$$ $g-g$and$$t{-\chi _\textrm{u}}$$ $t-{\chi }_{\text{u}}$fusion. We work under a simplified model approach and directly take the$$\chi _\textrm{u}$$ ${\chi }_{\text{u}}$and$$\phi '$$ ${\varphi }^{\prime }$masses as free parameters. We perform a phenomenological analysis considering$$\chi _\textrm{u}$$ ${\chi }_{\text{u}}$final states to b-quarks, muons, and neutrinos, and$$\phi '$$ ${\varphi }^{\prime }$decays to$$\mu ^+\mu ^-$$ ${\mu }^{+}{\mu }^{-}$. A machine learning algorithm is used to maximize the signal sensitivity, considering an integrated luminosity of 3000$$\text {fb}^{-1}$$ ${\text{fb}}^{-1}$. The proposed methodology can be a key mode for discovery over a large mass range, including low masses, traditionally considered difficult due to experimental constraints. 

   more » « less