@article{19672,
  abstract     = {Some of the classical models of tropical cyclone intensification predict tropical cyclones to intensify up to a steady intensity, which depends on surface fluxes only, without any relevant role played by convective motions in the troposphere, typically assumed to have a moist adiabatic lapse rate. Simulations performed using the non-hydrostatic, high-resolution model System for Atmosphere Modeling in idealized settings (rotating radiative-convective equilibrium on a doubly periodic domain) show early intensification consistent with these theoretical expectations, but different intensity evolution, with the cyclone undergoing an oscillation in wind speed. This oscillation can be linked to feedbacks between the cyclone intensity and air buoyancy: convective heating, radiative heating, and mixing with warm low stratospheric air warm the mid and upper troposphere of the cyclone stabilizing the air column and thus reducing its intensity. After the intensity decay phase, mid and upper tropospheric cooling, mostly through cold advection from the surroundings, cooled by radiation, rebuilds Convective Available Potential Energy, that peaks just before a new intensification phase. These idealized simulations thus highlight the potentially important interactions between a tropical cyclone, its environment and radiation.},
  author       = {Polesello, Andrea and Charinti, Giousef Alexandros and Meroni, Agostino Niyonkuru and Muller, Caroline J and Pasquero, Claudia},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {4},
  publisher    = {Wiley},
  title        = {{Intensity oscillations of tropical cyclones: Surface versus mid and upper tropospheric processes}},
  doi          = {10.1029/2024MS004613},
  volume       = {17},
  year         = {2025},
}

@article{20705,
  abstract     = {Optical tweezers are widely used as a highly sensitive tool to measure forces on micron-scale particles. One such application is the measurement of the electric charge of a particle, which can be done with high precision in liquids, air, or vacuum. We experimentally investigate how the trapping laser itself can electrically charge such a particle, in our case a ∼1  μ⁢m SiO2 sphere in air. We model the charging mechanism as a two-photon process which reproduces the experimental data with high fidelity.},
  author       = {Stöllner, Andrea and Lenton, Isaac C and Volosniev, Artem and Millen, James and Shibuya, Renjiro and Ishii, Hisao and Rak, Dmytro and Alpichshev, Zhanybek and David, Grégory and Signorell, Ruth and Muller, Caroline J and Waitukaitis, Scott R},
  issn         = {1079-7114},
  journal      = {Physical Review Letters},
  number       = {21},
  publisher    = {American Physical Society},
  title        = {{Using optical tweezers to simultaneously trap, charge, and measure the charge of a microparticle in air}},
  doi          = {10.1103/5xd9-4tjj},
  volume       = {135},
  year         = {2025},
}

@article{20795,
  abstract     = {The tropical climate variability is characterized by various oscillations across a range of timescales. Oscillations that imprint the tropical mean state are generally attributed to slow processes, such as the seasonal cycle or interannual variability. Here, we identify a pronounced tropics-wide intraseasonal oscillation (TWISO) in satellite observations and reanalyses. This oscillation, with a period of 30 to 60 d, is evident across multiple variables and involves interactions between convection, radiation, surface fluxes, and large-scale circulation. It is primarily manifested as convective perturbations in the tropical Indo-Pacific warm pool accompanied by oscillations in the large-scale tropical overturning circulation. Here, we examine the relationship between TWISO, the Madden–Julian Oscillation (MJO), and the instability of radiative-convective equilibrium. Certain phases of TWISO coincide with specific phases of the MJO, suggesting a potential connection between the two. However, although the MJO can amplify the oscillation amplitude of TWISO, it is not essential for TWISO to occur. Finally, due to its broad manifestation across the tropics, TWISO potentially exerts widespread influence on tropical weather and climate at regional scales.},
  author       = {Bao, Jiawei and Bony, Sandrine and Takasuka, Daisuke and Muller, Caroline J},
  issn         = {1091-6490},
  journal      = {Proceedings of the National Academy of Sciences},
  number       = {48},
  publisher    = {National Academy of Sciences},
  title        = {{Tropics-wide intraseasonal oscillations}},
  doi          = {10.1073/pnas.2511549122},
  volume       = {122},
  year         = {2025},
}

@article{20026,
  abstract     = {Deep Convective Systems (DCSs) reaching scales of 100–1000 km play a pivotal role as the primary precipitation source in the tropics. Those systems can have large cloud shields, and thus not only affect severe precipitation patterns but also play a crucial part in modulating the tropical radiation budget. Understanding the complex factors that control how these systems grow and how they will behave in a warming climate remain fundamental challenges. Research efforts have been directed, on one hand, towards understanding the environmental control on these systems, and on the other hand, towards exploring the internal potential of systems to develop and self-aggregate in idealized simulations. However, we still lack understanding on the relative role of the environment and internal feedbacks on DCS mature size and why. The novel high-resolution global SAM simulation from the DYAMOND project, combined with the TOOCAN Lagrangian tracking of DCSs and machine learning tools, offers an unprecedented opportunity to explore this question. We find that a system’s growth rate during the first 2 h of development predicts its final size with a Pearson correlation coefficient of 0.65. Beyond this period, growth rate emerges as the strongest predictor. However, in the early stages, additional factors–such as ice water path heterogeneity, migration distance, interactions with neighboring systems, and deep shear–play a more significant role. Our study quantitatively assesses the relative influence of internal versus external factors on the mature cloud shield size. Our results show that system-intrinsic properties exert a stronger influence than environmental conditions, suggesting that the initial environment does not strictly constrain final system size, particularly for larger systems where internal dynamics dominate.},
  author       = {Abramian, Sophie and Muller, Caroline J and Risi, Camille and Fiolleau, Thomas and Roca, Rémy},
  issn         = {2397-3722},
  journal      = {npj Climate and Atmospheric Science},
  publisher    = {Springer Nature},
  title        = {{How key features of early development shape deep convective systems}},
  doi          = {10.1038/s41612-025-01154-1},
  volume       = {8},
  year         = {2025},
}

@article{18708,
  abstract     = {This study investigates the response of Indian summer monsoon (ISM) precipitation to CO2 removal, with a specific focus on regional and subseasonal variations. Following CO2 removal, monsoon circulation weakens throughout the summer owing to the reduced large-scale meridional temperature gradient around India. Weakened monsoon circulation decreases the local-scale thermodynamic stability within India, following monsoon-onset periods. While the frequency of synoptic-scale ISM low-pressure systems (LPSs) decreases overall, the lower thermodynamic stability causes the LPSs to form and resultantly shift west and south from their typical paths, last longer and move more quickly zonally during August and September. Changes in these rain-producing processes induce distinct regional (Western Ghats, south-central-east India, and Tamil Nadu) and subseasonal precipitation responses. Also, extreme precipitation exhibits similar patterns, but is more strongly affected by changes in LPS. Our results suggest that reliable future projections of regional hydroclimate change require a more accurate understanding of multi-scale precipitation processes.},
  author       = {Paik, Seungmok and Kim, Daehyun and An, Soon Il and Oh, Hyoeun and Shin, Jongsoo and Goswami, Bidyut B and Min, Seung Ki and Mondal, Sanjit Kumar},
  issn         = {2397-3722},
  journal      = {npj Climate and Atmospheric Science},
  publisher    = {Springer Nature},
  title        = {{Exploring causes of distinct regional and subseasonal Indian summer monsoon precipitation responses to CO2 removal}},
  doi          = {10.1038/s41612-024-00858-0},
  volume       = {7},
  year         = {2024},
}

@inproceedings{19005,
  abstract     = {Causal representation learning promises to extend causal models to hidden causal
variables from raw entangled measurements. However, most progress has focused
on proving identifiability results in different settings, and we are not aware of any
successful real-world application. At the same time, the field of dynamical systems
benefited from deep learning and scaled to countless applications but does not allow
parameter identification. In this paper, we draw a clear connection between the two
and their key assumptions, allowing us to apply identifiable methods developed
in causal representation learning to dynamical systems. At the same time, we can
leverage scalable differentiable solvers developed for differential equations to build
models that are both identifiable and practical. Overall, we learn explicitly controllable models that isolate the trajectory-specific parameters for further downstream
tasks such as out-of-distribution classification or treatment effect estimation. We
experiment with a wind simulator with partially known factors of variation. We
also apply the resulting model to real-world climate data and successfully answer
downstream causal questions in line with existing literature on climate change.
Code is available at https://github.com/CausalLearningAI/crl-dynamical-systems.},
  author       = {Yao, Dingling and Muller, Caroline J and Locatello, Francesco},
  booktitle    = {38th Conference on Neural Information Processing Systems},
  location     = {Vancouver, Canada},
  publisher    = {Neural Information Processing Systems Foundation},
  title        = {{Marrying causal representation learning with dynamical systems for science}},
  volume       = {37},
  year         = {2024},
}

@misc{19307,
  abstract     = {This repository contains the data, scripts, SAM codes and files required to reproduce the results of the manuscript "The Unreasonable Efficiency of Total Rain Evaporation Removal in Triggering Convective Self-Aggregation" submitted to the Geophysical Research Letters (GRL).

Brief description of project: This project aims to examine the impact of rain evaporation removal or reduction in the planetary boundary layer (PBL) on convective self aggregation (CSA). Non-rotating radiative-convective equilibrium (RCE) simulations were conducted with the System for Atmospheric Modeling (SAM) cloud resolving model. Rain evaporation in the lowest 1 km was progressively reduced and the effect on CSA was investigated. The physical processes underlying this type of aggregation (referred to in the manuscript as no-evaporation CSA, or NE-CSA) were analyzed and described. 
The default SAM code base (version 6.10.8) can be downloaded from here: http://rossby.msrc.sunysb.edu/~marat/SAM.html},
  author       = {Hwong, Yi-Ling and Muller, Caroline J},
  publisher    = {Zenodo},
  title        = {{Data - The unreasonable efficiency of total rain evaporation removal in triggering convective self-aggregation}},
  doi          = {10.5281/ZENODO.10687169},
  year         = {2024},
}

@article{15047,
  abstract     = {Tropical precipitation extremes and their changes with surface warming are investigated using global storm resolving simulations and high-resolution observations. The simulations demonstrate that the mesoscale organization of convection, a process that cannot be physically represented by conventional global climate models, is important for the variations of tropical daily accumulated precipitation extremes. In both the simulations and observations, daily precipitation extremes increase in a more organized state, in association with larger, but less frequent, storms. Repeating the simulations for a warmer climate results in a robust increase in monthly-mean daily precipitation extremes. Higher precipitation percentiles have a greater sensitivity to convective organization, which is predicted to increase with warming. Without changes in organization, the strongest daily precipitation extremes over the tropical oceans increase at a rate close to Clausius-Clapeyron (CC) scaling. Thus, in a future warmer state with increased organization, the strongest daily precipitation extremes over oceans increase at a faster rate than CC scaling.},
  author       = {Bao, Jiawei and Stevens, Bjorn and Kluft, Lukas and Muller, Caroline J},
  issn         = {2375-2548},
  journal      = {Science Advances},
  number       = {8},
  publisher    = {American Association for the Advancement of Science},
  title        = {{Intensification of daily tropical precipitation extremes from more organized convection}},
  doi          = {10.1126/sciadv.adj6801},
  volume       = {10},
  year         = {2024},
}

@article{15097,
  abstract     = {Global storm-resolving models (GSRMs) use strongly refined horizontal grids compared with the climate models typically used in the Coupled Model Intercomparison Project (CMIP) but employ comparable vertical grid spacings. Here, we study how changes in the vertical grid spacing and adjustments to the integration time step affect the basic climate quantities simulated by the ICON-Sapphire atmospheric GSRM. Simulations are performed over a 45 d period for five different vertical grids with between 55 and 540 vertical layers and maximum tropospheric vertical grid spacings of between 800 and 50 m, respectively. The effects of changes in the vertical grid spacing are compared with the effects of reducing the horizontal grid spacing from 5 to 2.5 km. For most of the quantities considered, halving the vertical grid spacing has a smaller effect than halving the horizontal grid spacing, but it is not negligible. Each halving of the vertical grid spacing, along with the necessary reductions in time step length, increases cloud liquid water by about 7 %, compared with an approximate 16 % decrease for halving the horizontal grid spacing. The effect is due to both the vertical grid refinement and the time step reduction. There is no tendency toward convergence in the range of grid spacings tested here. The cloud ice amount also increases with a refinement in the vertical grid, but it is hardly affected by the time step length and does show a tendency to converge. While the effect on shortwave radiation is globally dominated by the altered reflection due to the change in the cloud liquid water content, the effect on longwave radiation is more difficult to interpret because changes in the cloud ice concentration and cloud fraction are anticorrelated in some regions. The simulations show that using a maximum tropospheric vertical grid spacing larger than 400 m would increase the truncation error strongly. Computing time investments in a further vertical grid refinement can affect the truncation errors of GSRMs similarly to comparable investments in horizontal refinement, because halving the vertical grid spacing is generally cheaper than halving the horizontal grid spacing. However, convergence of boundary layer cloud properties cannot be expected, even for the smallest maximum tropospheric grid spacing of 50 m used in this study.},
  author       = {Schmidt, Hauke and Rast, Sebastian and Bao, Jiawei and Cassim, Amrit and Fang, Shih Wei and Jimenez-De La Cuesta, Diego and Keil, Paul and Kluft, Lukas and Kroll, Clarissa and Lang, Theresa and Niemeier, Ulrike and Schneidereit, Andrea and Williams, Andrew I.L. and Stevens, Bjorn},
  issn         = {1991-9603},
  journal      = {Geoscientific Model Development},
  number       = {4},
  pages        = {1563--1584},
  publisher    = {European Geosciences Union},
  title        = {{Effects of vertical grid spacing on the climate simulated in the ICON-Sapphire global storm-resolving model}},
  doi          = {10.5194/gmd-17-1563-2024},
  volume       = {17},
  year         = {2024},
}

@article{15165,
  abstract     = {Current knowledge suggests a drought Indian monsoon (perhaps a severe one) when the El Nino Southern Oscillation and Pacific Decadal Oscillation each exhibit positive phases (a joint positive phase). For the monsoons, which are exceptions in this regard, we found northeast India often gets excess pre-monsoon rainfall. Further investigation reveals that this excess pre-monsoon rainfall is produced by the interaction of the large-scale circulation associated with the joint phase with the mountains in northeast India. We posit that a warmer troposphere, a consequence of excess rainfall over northeast India, drives a stronger monsoon circulation and enhances monsoon rainfall over central India. Hence, we argue that pre-monsoon rainfall over northeast India can be used for seasonal monsoon rainfall prediction over central India. Most importantly, its predictive value is at its peak when the Pacific Ocean exhibits a joint positive phase and the threat of extreme drought monsoon looms over India.},
  author       = {Goswami, Bidyut B},
  issn         = {1944-8007},
  journal      = {Geophysical Research Letters},
  number       = {5},
  publisher    = {Wiley},
  title        = {{A pre-monsoon signal of false alarms of Indian monsoon droughts}},
  doi          = {10.1029/2023GL106569},
  volume       = {51},
  year         = {2024},
}

@article{15186,
  abstract     = {The elimination of rain evaporation in the planetary boundary layer (PBL) has been found to lead to convective self‐aggregation (CSA) even without radiative feedback, but the precise mechanisms underlying this phenomenon remain unclear. We conducted cloud‐resolving simulations with two domain sizes and progressively reduced rain evaporation in the PBL. Surprisingly, CSA only occurred when rain evaporation was almost completely removed. The additional convective heating resulting from the reduction of evaporative cooling in the moist patch was found to be the trigger, thereafter a dry subsidence intrusion into the PBL in the dry patch takes over and sets CSA in motion. Temperature and moisture anomalies oppose each other in their buoyancy effects, hence explaining the need for almost total rain evaporation removal. We also found radiative cooling and not cold pools to be the leading cause for the comparative ease of CSA to take place in the larger domain.},
  author       = {Hwong, Yi-Ling and Muller, Caroline J},
  issn         = {1944-8007},
  journal      = {Geophysical Research Letters},
  keywords     = {General Earth and Planetary Sciences, Geophysics},
  number       = {6},
  publisher    = {Wiley},
  title        = {{The unreasonable efficiency of total rain evaporation removal in triggering convective self‐aggregation}},
  doi          = {10.1029/2023gl106523},
  volume       = {51},
  year         = {2024},
}

@article{15313,
  abstract     = {Our goal is to investigate fundamental properties of the system of internally cooled convection. The system consists of an upward thermal flux at the lower boundary, a mean temperature lapse-rate and a constant cooling term in the bulk with the bulk cooling in thermal equilibrium with the input heat flux. This simple model represents idealised dry convection in the atmospheric boundary layer, where the cooling mimics the radiative cooling to space notably through longwave radiation. We perform linear stability analysis of the model for different values of the mean stratification to derive the critical forcing above which the fluid is convectively unstable to small perturbations. The dynamic behavior of the fluid system is described and the scaling of various important measured quantities such as the total vertical convective heat flux and the upward mass flux is measured. We introduce a lapse-rate dependent dimensionless Rayleigh-number Ray that determines the behavior of the system, finding that the convective heat-flux and mass-flux scale approximately as Ray0.5 and Ray0.7 respectively. The area-fraction of the domain that is occupied by upward and downward moving fluid and the skewness of the vertical velocity are studied to understand the asymmetry inherent in the system. We conclude with a short discussion on the relevance to atmospheric convection and the scope for further investigations of atmospheric convection using similar simplified approaches.},
  author       = {Agasthya, Lokahith N and Muller, Caroline J},
  issn         = {1007-5704},
  journal      = {Communications in Nonlinear Science and Numerical Simulation},
  publisher    = {Elsevier},
  title        = {{Dynamics and scaling of internally cooled convection}},
  doi          = {10.1016/j.cnsns.2024.108011},
  volume       = {134},
  year         = {2024},
}

@article{17396,
  abstract     = {Wildfire intensity and severity have been increasing in the Iberian Peninsula in recent years, particularly in the Galicia region, due to rising temperatures and accumulating drier combustible vegetation in unmanaged lands. This leads to substantial emissions of air pollutants, notably fine particles (PM2.5), posing a risk to public health. This study aims to assess the impact of local and regional wildfires on PM2.5 levels in Galicia's main cities and their implications for air quality and public health. Over a decade (2013–2022), PM2.5 data during wildfire seasons were analyzed using statistical methods and Lagrangian tracking to monitor smoke plume evolution. The results reveal a notable increase in PM2.5 concentration during the wildfire season (June–November) in Galicia, surpassing health guidelines during extreme events and posing a significant health risk to the population. Regional wildfire analyses indicate that smoke plumes from Northern Portugal contribute to pollution in Galician cities, influencing the seasonality of heightened PM2.5 levels. During extensive wildfires, elevated PM2.5 concentration values persisted for several days, potentially exacerbating health concerns in Galicia. These findings underscore the urgency of implementing air pollution prevention and management measures in the region, including developing effective alerts for large-scale events and improved wildfire management strategies to mitigate their impact on air quality in Galician cities.},
  author       = {Quishpe-Vásquez, César and Oliva, Patricia and López-Barrera, Ellie Anne and Casallas Garcia, Alejandro},
  issn         = {0301-4797},
  journal      = {Journal of Environmental Management},
  publisher    = {Elsevier},
  title        = {{Wildfires impact on PM2.5 concentration in galicia Spain}},
  doi          = {10.1016/j.jenvman.2024.122093},
  volume       = {367},
  year         = {2024},
}

@article{17435,
  abstract     = {The Mediterranean region is experiencing pronounced aridification and in certain areas higher occurrence of intense precipitation. In this work, we analyze the evolution of the precipitation probability distribution in terms of precipitating days (or “wet-days”) and all-days quantile trends, in Europe and the Mediterranean, using the ERA5 reanalysis. Looking at the form of wet-days quantile trends curves, we identify four regimes. Two are predominant: in most of northern Europe the precipitation quantiles all intensify, while in the Mediterranean the low-medium quantiles are mostly decreasing as extremes intensify or decrease. The wet-days distribution is then modeled by a Weibull law with two parameters, whose changes capture the four regimes. Assessing the significance of the parameters' changes over 1950–2020 shows that a signal on wet-days distribution has already emerged in northern Europe (where the distribution shifts to more intense precipitation), but not yet in the Mediterranean, where the natural variability is stronger. We extend the results by describing the all-days distribution change as the wet-days’ change plus a contribution from the dry-days frequency change, and study their relative contribution. In northern Europe, the wet-days distribution change is the dominant driver, and the contribution of dry-days frequency change can be neglected for wet-days percentiles above about 50%. In the Mediterranean, however, the change of precipitation distribution comes from the significant increase of dry-days frequency instead of an intensity change during wet-days. Therefore, in the Mediterranean the increase of dry-days frequency is crucial for all-days trends, even for heavy precipitation.},
  author       = {André, Julie and D'Andrea, Fabio and Drobinski, Philippe and Muller, Caroline J},
  issn         = {2169-8996},
  journal      = {Journal of Geophysical Research: Atmospheres},
  number       = {15},
  publisher    = {Wiley},
  title        = {{Regimes of precipitation change over Europe and the Mediterranean}},
  doi          = {10.1029/2023JD040413},
  volume       = {129},
  year         = {2024},
}

@article{14453,
  abstract     = {Squall lines are substantially influenced by the interaction of low-level shear with cold pools associated with convective downdrafts. Beyond an optimal shear amplitude, squall lines tend to orient themselves at an angle with respect to the low-level shear. While the mechanisms behind squall line orientation seem to be increasingly well understood, uncertainties remain on the implications of this orientation. Roca and Fiolleau (2020, https://doi.org/10.1038/s43247-020-00015-4) show that long lived mesoscale convective systems, including squall lines, are disproportionately involved in rainfall extremes in the tropics. This article investigates the influence of the interaction between low-level shear and squall line outflow on squall line generated precipitation extrema in the tropics. Using a cloud resolving model, simulated squall lines in radiative convective equilibrium amid a shear-dominated regime (super optimal), a balanced regime (optimal), and an outflow dominated regime (suboptimal). Our results show that precipitation extremes in squall lines are 40% more intense in the case of optimal shear and remain 30% superior in the superoptimal regime relative to a disorganized case. With a theoretical scaling of precipitation extremes (C. Muller & Takayabu, 2020, https://doi.org/10.1088/1748-9326/ab7130), we show that the condensation rates control the amplification of precipitation extremes in tropical squall lines, mainly due to its change in vertical mass flux (dynamic component). The reduction of dilution by entrainment explains half of this change, consistent with Mulholland et al. (2021, https://doi.org/10.1175/jas-d-20-0299.1). The other half is explained by increased cloud-base velocity intensity in optimal and superoptimal squall lines.},
  author       = {Abramian, Sophie and Muller, Caroline J and Risi, Camille},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {10},
  publisher    = {Wiley},
  title        = {{Extreme precipitation in tropical squall lines}},
  doi          = {10.1029/2022MS003477},
  volume       = {15},
  year         = {2023},
}

@article{14564,
  abstract     = {Cumulus parameterization (CP) in state‐of‐the‐art global climate models is based on the quasi‐equilibrium assumption (QEA), which views convection as the action of an ensemble of cumulus clouds, in a state of equilibrium with respect to a slowly varying atmospheric state. This view is not compatible with the organization and dynamical interactions across multiple scales of cloud systems in the tropics and progress in this research area was slow over decades despite the widely recognized major shortcomings. Novel ideas on how to represent key physical processes of moist convection‐large‐scale interaction to overcome the QEA have surged recently. The stochastic multicloud model (SMCM) CP in particular mimics the dynamical interactions of multiple cloud types that characterize organized tropical convection. Here, the SMCM is used to modify the Zhang‐McFarlane (ZM) CP by changing the way in which the bulk mass flux and bulk entrainment and detrainment rates are calculated. This is done by introducing a stochastic ensemble of plumes characterized by randomly varying detrainment level distributions based on the cloud area fraction of the SMCM. The SMCM is here extended to include shallow cumulus clouds resulting in a unified shallow‐deep CP. The new stochastic multicloud plume CP is validated against the control ZM scheme in the context of the single column Community Climate Model of the National Center for Atmospheric Research using data from both tropical ocean and midlatitude land convection. Some key features of the SMCM CP such as it capability to represent the tri‐modal nature of organized convection are emphasized.},
  author       = {Khouider, B. and GOSWAMI, BIDYUT B and Phani, R. and Majda, A. J.},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  keywords     = {General Earth and Planetary Sciences, Environmental Chemistry, Global and Planetary Change},
  number       = {11},
  publisher    = {American Geophysical Union},
  title        = {{A shallow‐deep unified stochastic mass flux cumulus parameterization in the single column community climate model}},
  doi          = {10.1029/2022ms003391},
  volume       = {15},
  year         = {2023},
}

@article{14654,
  abstract     = {Two assumptions commonly applied in convection schemes—the diagnostic and quasi-equilibrium assumptions—imply that convective activity (e.g., convective precipitation) is controlled only by the large-scale (macrostate) environment at the time. In contrast, numerical experiments indicate a “memory” or dependence of convection also on its own previous activity whereby subgrid-scale (microstate) structures boost but are also boosted by convection. In this study we investigated this memory by comparing single-column model behavior in two idealized tests previously executed by a cloud-resolving model (CRM). Conventional convection schemes that employ the diagnostic assumption fail to reproduce the CRM behavior. The memory-capable org and Laboratoire de Météorologie Dynamique Zoom cold pool schemes partially capture the behavior, but fail to fully exhibit the strong reinforcing feedbacks implied by the CRM. Analysis of this failure suggests that it is because the CRM supports a linear (or superlinear) dependence of the subgrid structure growth rate on the precipitation rate, while the org scheme assumes a sublinear dependence. Among varying versions of the org scheme, the growth rate of the org variable representing subgrid structure is strongly associated with memory strength. These results demonstrate the importance of parameterizing convective memory, and the ability of idealized tests to reveal shortcomings of convection schemes and constrain model structural assumptions.},
  author       = {Hwong, Yi-Ling and Colin, M. and Aglas, Philipp and Muller, Caroline J and Sherwood, S. C.},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {12},
  publisher    = {Wiley},
  title        = {{Assessing memory in convection schemes using idealized tests}},
  doi          = {10.1029/2023MS003726},
  volume       = {15},
  year         = {2023},
}

@article{12791,
  abstract     = {We investigate the capabilities of Physics-Informed Neural Networks (PINNs) to reconstruct turbulent Rayleigh–Bénard flows using only temperature information. We perform a quantitative analysis of the quality of the reconstructions at various amounts of low-passed-filtered information and turbulent intensities. We compare our results with those obtained via nudging, a classical equation-informed data assimilation technique. At low Rayleigh numbers, PINNs are able to reconstruct with high precision, comparable to the one achieved with nudging. At high Rayleigh numbers, PINNs outperform nudging and are able to achieve satisfactory reconstruction of the velocity fields only when data for temperature is provided with high spatial and temporal density. When data becomes sparse, the PINNs performance worsens, not only in a point-to-point error sense but also, and contrary to nudging, in a statistical sense, as can be seen in the probability density functions and energy spectra.},
  author       = {Clark Di Leoni, Patricio and Agasthya, Lokahith N and Buzzicotti, Michele and Biferale, Luca},
  issn         = {1292-895X},
  journal      = {The European Physical Journal E},
  number       = {3},
  publisher    = {Springer Nature},
  title        = {{Reconstructing Rayleigh–Bénard flows out of temperature-only measurements using Physics-Informed Neural Networks}},
  doi          = {10.1140/epje/s10189-023-00276-9},
  volume       = {46},
  year         = {2023},
}

@article{13256,
  abstract     = {The El Niño-Southern Oscillation (ENSO) and the Indian summer monsoon (ISM, or monsoon) are two giants of tropical climate. Here we assess the future evolution of the ENSO-monsoon teleconnection in climate simulations with idealized forcing of CO2 increment at a rate of 1% year-1 starting from a present-day condition (367 p.p.m.) until quadrupling. We find a monotonous weakening of the ENSO-monsoon teleconnection with the increase in CO2. Increased co-occurrences of El Niño and positive Indian Ocean Dipoles (pIODs) in a warmer climate weaken the teleconnection. Co-occurrences of El Niño and pIOD are attributable to mean sea surface temperature (SST) warming that resembles a pIOD-type warming pattern in the Indian Ocean and an El Niño-type warming in the Pacific. Since ENSO is a critical precursor of the strength of the Indian monsoon, a weakening of this relation may mean a less predictable Indian monsoon in a warmer climate.},
  author       = {Goswami, Bidyut B and An, Soon Il},
  issn         = {2397-3722},
  journal      = {npj Climate and Atmospheric Science},
  publisher    = {Springer Nature},
  title        = {{An assessment of the ENSO-monsoon teleconnection in a warming climate}},
  doi          = {10.1038/s41612-023-00411-5},
  volume       = {6},
  year         = {2023},
}

@article{14752,
  abstract     = {Radiative cooling of the lowest atmospheric levels is of strong importance for modulating atmospheric circulations and organizing convection, but detailed observations and a robust theoretical understanding are lacking. Here we use unprecedented observational constraints from subsidence regimes in the tropical Atlantic to develop a theory for the shape and magnitude of low‐level longwave radiative cooling in clear‐sky, showing peaks larger than 5–10 K/day at the top of the boundary layer. A suite of novel scaling approximations is first developed from simplified spectral theory, in close agreement with the measurements. The radiative cooling peak height is set by the maximum lapse rate in water vapor path, and its magnitude is mainly controlled by the ratio of column relative humidity above and below the peak. We emphasize how elevated intrusions of moist air can reduce low‐level cooling, by sporadically shading the spectral range which effectively cools to space. The efficiency of this spectral shading depends both on water content and altitude of moist intrusions; its height dependence cannot be explained by the temperature difference between the emitting and absorbing layers, but by the decrease of water vapor extinction with altitude. This analytical work can help to narrow the search for low‐level cloud patterns sensitive to radiative‐convective feedbacks: the most organized patterns with largest cloud fractions occur in atmospheres below 10% relative humidity and feel the strongest low‐level cooling. This motivates further assessment of favorable conditions for radiative‐convective feedbacks and a robust quantification of corresponding shallow cloud dynamics in current and warmer climates.},
  author       = {Fildier, B. and Muller, Caroline J and Pincus, R. and Fueglistaler, S.},
  issn         = {2576-604X},
  journal      = {AGU Advances},
  keywords     = {General Earth and Planetary Sciences},
  number       = {3},
  publisher    = {American Geophysical Union},
  title        = {{How moisture shapes low‐level radiative cooling in subsidence regimes}},
  doi          = {10.1029/2023av000880},
  volume       = {4},
  year         = {2023},
}

