-
Using rare event algorithms to understand the statistics and dynamics of extreme heatwave seasons in South Asia
Authors:
Clément Le Priol,
Joy M. Monteiro,
Freddy Bouchet
Abstract:
Computing the return times of extreme events and assessing the impact of climate change on such return times is fundamental to extreme event attribution studies. However, the rarity of such events in the observational record makes this task a challenging one, even more so for "record-shattering" events that have not been previously observed at all. While climate models could be used to simulate su…
▽ More
Computing the return times of extreme events and assessing the impact of climate change on such return times is fundamental to extreme event attribution studies. However, the rarity of such events in the observational record makes this task a challenging one, even more so for "record-shattering" events that have not been previously observed at all. While climate models could be used to simulate such extremely rare events, such an approach entails a huge computational cost: gathering robust statistics for events with return time of centuries would require a few thousand years of simulation.
In this study, we use an innovative tool, rare event algorithm, that allows to sample numerous extremely rare events at a much lower cost than direct simulations. We employ the algorithm to sample extreme heatwave seasons, corresponding to large anomalies of the seasonal average temperature, in a heatwave hotspot of South Asia using the global climate model Plasim. We show that the algorithm estimates the return levels of extremely rare events with much greater precision than traditional statistical fits. It also enables the computation of various composite statistics, whose accuracy is demonstrated through comparison with a very long control run. In particular, our results reveal that extreme heatwave seasons are associated with an anticyclonic anomaly embedded within a large-scale hemispheric quasi-stationary wave-pattern. Additionally, the algorithm accurately represents the intensity-duration-frequency statistics of sub-seasonal heatwaves, offering insights into both seasonal and sub-seasonal aspects of extreme heatwave seasons. This innovative approach could be used in extreme event attribution studies to better constrain the changes in event's probability and intensity with global warming, particularly for events with return times spanning centuries or millennia.
△ Less
Submitted 24 July, 2024; v1 submitted 11 April, 2024;
originally announced April 2024.
-
Synchrotron Microtomography and Neutron Radiography Characterization of the Microstruture and Water Absorption of Concrete from Pompeii
Authors:
Ke Xu,
Anton S. Tremsin,
Jiaqi Li,
Daniela M. Ushizima,
Catherine A. Davy,
Amine Bouterf,
Ying Tsun Su,
Milena Marroccoli,
Anna Maria Mauro,
Massimo Osanna,
Antonio Telesca,
Paulo J. M. Monteiro
Abstract:
There is renewed interest in using advanced techniques to characterize ancient Roman concrete. In the present work, samples were drilled from the "Hospitium" in Pompeii and were analyzed by synchrotron microtomography (uCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were…
▽ More
There is renewed interest in using advanced techniques to characterize ancient Roman concrete. In the present work, samples were drilled from the "Hospitium" in Pompeii and were analyzed by synchrotron microtomography (uCT) and neutron radiography to study how the microstructure, including the presence of induced cracks, affects their water adsorption. The water distribution and absorptivity were quantified by neutron radiography. The 3D crack propagation, pore size distribution and orientation, tortuosity, and connectivity were analyzed from uCT results using advanced imaging methods. The concrete characterization also included classical methods (e.g., differential thermal-thermogravimetric, X-ray diffractometry, and scanning electron microscopy). Ductile fracture patterns were observed once cracks were introduced. When compared to Portland cement mortar/concrete, Pompeii samples had relatively high porosity, low connectivity, and similar coefficient of capillary penetration. In addition, the permeability was predicted from models based on percolation theory and the pore structure data to evaluate the fluid transport properties.
△ Less
Submitted 2 December, 2020; v1 submitted 26 May, 2020;
originally announced May 2020.
-
The Importance of the Ultra-alkaline Volcanic Nature of the Raw Materials to the Ductility of Roman Marine Concrete
Authors:
Jackson MacFarlane,
Tiziana Vanorio,
Paulo J. M. Monteiro
Abstract:
Roman-era concrete is the iconic embodiment of long-term physicochemical resilience. We investigated the basis of this behavior across scales of observations by coupling time-lapse (4-D) tomographic imaging of macroscopic mechanical stressing with structural microscopy and chemical spectroscopy on Roman marine concrete (RMC) from ancient harbors in Italy and Israel. Stress-strain measurements reve…
▽ More
Roman-era concrete is the iconic embodiment of long-term physicochemical resilience. We investigated the basis of this behavior across scales of observations by coupling time-lapse (4-D) tomographic imaging of macroscopic mechanical stressing with structural microscopy and chemical spectroscopy on Roman marine concrete (RMC) from ancient harbors in Italy and Israel. Stress-strain measurements revealed that RMC creeps and exhibits a ductile deformation mode. The permeability of specimens from Italy were found to be low due to increased matrix-aggregate bonding. Structural and chemical imaging shows the presence of well-developed sulfur-rich, fibrous minerals that are intertwined and embedded in a crossbred matrix having the chemical traits of both a calcium-aluminum-silicate-hydrate and a geopolymer. This latter likely reflects the ultra-alkaline volcanic nature of the primary source materials. We hypothesize that the fine interweave of sulfur-rich fibers within this crossbred matrix enhances aggregate bonding, which altogether contributes to the durability of RMC.
△ Less
Submitted 29 January, 2020;
originally announced January 2020.
-
Predicting the Expansion of Concrete Exposed to Sulfate Attack with a Regression Model Based on a Performance Classification
Authors:
Xiangru Jian,
Paulo J. M Monteiro,
Kimberly E. Kurtis
Abstract:
This paper mainly describes the development of a new type of regression model to predict the long-term expansion of concrete subjected to a sulfate-rich environment. The experimental data originated from a long-term (40+ years), nonaccelerated test program performed by the U.S. Bureau of Reclamation (USBR). Expansion data of specimens composed of 54 different mixtures were measured periodically th…
▽ More
This paper mainly describes the development of a new type of regression model to predict the long-term expansion of concrete subjected to a sulfate-rich environment. The experimental data originated from a long-term (40+ years), nonaccelerated test program performed by the U.S. Bureau of Reclamation (USBR). Expansion data of specimens composed of 54 different mixtures were measured periodically throughout the entire test program. In this analysis, the mixtures were first classified into three groups using K-means clustering based on their expansion patterns. Within each group, the expansion rate was predicted as an exclusive regression function of the water-cement ratio (W/C), tricalcium aluminate (C_3 A) content of cement, cement content of cement or time of expansion. Then, a support vector machine (SVM) was employed to determine the classification criteria by relying on the characteristics of the mixture proportions rather than the experimental performance, thereby enabling the model to offer predictions for new mixtures without test data. An analysis of the model indicated that concrete specimens with different mixture proportions, especially with different W/C values or C_3 A contents, are unlikely to share identical expansion patterns and should be considered and predicted separately.separately.
△ Less
Submitted 13 November, 2018; v1 submitted 10 October, 2018;
originally announced October 2018.
-
A Tutorial on Canonical Correlation Methods
Authors:
Viivi Uurtio,
João M. Monteiro,
Jaz Kandola,
John Shawe-Taylor,
Delmiro Fernandez-Reyes,
Juho Rousu
Abstract:
Canonical correlation analysis is a family of multivariate statistical methods for the analysis of paired sets of variables. Since its proposition, canonical correlation analysis has for instance been extended to extract relations between two sets of variables when the sample size is insufficient in relation to the data dimensionality, when the relations have been considered to be non-linear, and…
▽ More
Canonical correlation analysis is a family of multivariate statistical methods for the analysis of paired sets of variables. Since its proposition, canonical correlation analysis has for instance been extended to extract relations between two sets of variables when the sample size is insufficient in relation to the data dimensionality, when the relations have been considered to be non-linear, and when the dimensionality is too large for human interpretation. This tutorial explains the theory of canonical correlation analysis including its regularised, kernel, and sparse variants. Additionally, the deep and Bayesian CCA extensions are briefly reviewed. Together with the numerical examples, this overview provides a coherent compendium on the applicability of the variants of canonical correlation analysis. By bringing together techniques for solving the optimisation problems, evaluating the statistical significance and generalisability of the canonical correlation model, and interpreting the relations, we hope that this article can serve as a hands-on tool for applying canonical correlation methods in data analysis.
△ Less
Submitted 7 November, 2017;
originally announced November 2017.
-
Tropical Vorticity Forcing and Superrotation in the Spherical Shallow Water Equations
Authors:
Suhas DL,
Jai Sukhatme,
Joy M. Monteiro
Abstract:
The response of the nonlinear shallow water equations (SWE) on a sphere to tropical vorticity forcing is examined with an emphasis on momentum fluxes and the emergence of a superrotating (SR) state. Fixing the radiative damping and momentum drag timescales to be of the order of a few days, a state of SR is shown to emerge under steady large-scale and random small-scale vorticity forcing. In the fi…
▽ More
The response of the nonlinear shallow water equations (SWE) on a sphere to tropical vorticity forcing is examined with an emphasis on momentum fluxes and the emergence of a superrotating (SR) state. Fixing the radiative damping and momentum drag timescales to be of the order of a few days, a state of SR is shown to emerge under steady large-scale and random small-scale vorticity forcing. In the first example, the stationary response to a pair of equal and oppositely signed vortices placed on the equator is considered. Here, the equatorial flux budget is dominated by the eddy fluxes and these drive the system into a state of SR. Eventually, the flux associated with the momentum drag increases to balance the eddy fluxes, resulting in a steady state with a SR jet at the equator. The initial value problem with these twin vortices also exhibits SR driven by eddy fluxes. Curiously, this transient solution spontaneously propagates westward and continually circumnavigates the globe. It is worth emphasizing that this SR state does not rely on any particular form of dissipation at large scales, and is in fact observed even in the absence of any large-scale damping mechanism. In the second example, a random small-scale vorticity forcing is applied across the tropics. The statistically steady state obtained is fairly turbulent in nature, but here too, the eddy fluxes dominate, and the system exhibits SR. It is important to note that in both these cases, the direct forcing of the zonal mean zonal flow is zero by construction, and the eddy fluxes at the equator are responsible for its eastward acceleration. Further, in both these examples, the rotational part of the flow dominates the momentum fluxes as well as the stationary zonal mean zonal flow itself. Arguments based on the nature of potential vorticity and enstrophy are put forth to shed some light on these results.
△ Less
Submitted 6 December, 2016;
originally announced December 2016.
-
Quasi-geostrophic dynamics in the presence of moisture gradients
Authors:
Joy M. Monteiro,
Jai Sukhatme
Abstract:
The derivation of a quasi-geostrophic (QG) system from the rotating shallow water equations on a midlatitude beta-plane coupled with moisture is presented. Condensation is prescribed to occur whenever the moisture at a point exceeds a prescribed saturation value. It is seen that a slow condensation time scale is required to obtain a consistent set of equations at leading order. Further, since the…
▽ More
The derivation of a quasi-geostrophic (QG) system from the rotating shallow water equations on a midlatitude beta-plane coupled with moisture is presented. Condensation is prescribed to occur whenever the moisture at a point exceeds a prescribed saturation value. It is seen that a slow condensation time scale is required to obtain a consistent set of equations at leading order. Further, since the advecting wind fields are geostrophic, changes in moisture (and hence, precipitation) occur only via non-divergent mechanisms. Following observations, a saturation profile with gradients in the zonal and meridional directions is prescribed. A purely meridional gradient has the effect of slowing down the dry Rossby waves, through a reduction in the "equivalent gradient" of the background potential vorticity. A large scale unstable moist mode results on the inclusion of a zonal gradient by itself, or in conjunction with a meridional moisture gradient. For gradients that are are representative of the atmosphere, the most unstable moist mode propagates zonally in the direction of increasing moisture, matures over an intraseasonal timescale and has small phase speed.
△ Less
Submitted 13 April, 2016;
originally announced April 2016.