@article{21013,
  abstract     = {We have addressed convective self‐aggregation (CSA) in steady and oscillating sea surface temperature (SST) and solar radiation (SOLIN) cloud‐resolving model simulations in a non‐rotating radiative‐convective equilibrium (RCE) framework. Our experiment designs are motivated by land‐ocean heterogeneity of atmospheric convection. The steady and oscillating forcings are idealizations of ocean and land conditions, respectively, based on their differences in heat capacities. In both kinds of simulations, the diurnal mean SST and SOLIN are the same, and both SST and SOLIN are only varied in time (i.e., they are spatially homogeneous at any given time). We find that diurnally oscillating forcing accelerates CSA. Stronger long‐wave cooling in dry regions at night and during the warm SST phase (late afternoon) both allow the long‐wave feedback, known to favor aggregation, to intensify compared to steady forcing simulations. In addition to the long‐wave, reduced short‐wave warming in dry regions (during the day) further enhances radiative cooling there compared to moist regions. Overall, the radiative cooling is enhanced in dry regions compared to neighboring moist convective regions. A dry subsidence is driven by this net radiative (short‐wave plus long‐wave) cooling, consistent with earlier work on CSA. Stronger radiative cooling allows stronger subsidence which allows low‐level circulation to more efficiently transport moisture and energy up‐gradient, driving convection to aggregate faster. We also note a sensitivity of our experimental setup to initial conditions, more so at warmer SST. This stochastic behavior might be critical in reconciling the differences of opinion regarding the response of convection aggregation to oscillating SST forcing.},
  author       = {GOSWAMI, BIDYUT B and Lu, Ziyin and Muller, Caroline J},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {1},
  publisher    = {Wiley},
  title        = {{Convective self‐aggregation in diurnally oscillating sea surface temperature and solar forcing experiments}},
  doi          = {10.1029/2024ms004576},
  volume       = {18},
  year         = {2026},
}

@article{21035,
  abstract     = {According to the scientific consensus, tropical convection must decrease with global warming. This decrease is manifested by a decrease of the mass transported in the upward branch of the atmospheric overturning circulation – the convective mass flux – and a connected decrease of high clouds in the tropics, with implications for climate sensitivity. By using kilometer-scale simulations in radiative-convective equilibrium and a convective tracking algorithm, we show that no such decrease occurs in storms when taken individually and that the mass transport per storm increases instead. Storms can achieve this result by aggregating more surface of the convective cores – the inner part of the storm doing the vertical transport – so that the decrease of tropical convection is actually explained by a decrease in the total number of storms. There is little variation of the mean pressure velocity in the cores of the storms, a robust finding of this study. This remarkable invariance of the mean pressure velocity points to an emerging property of convection that should receive more attention in future studies.},
  author       = {Bolot, Maximilien and Roca, Rémy and Fiolleau, Thomas and Muller, Caroline J},
  issn         = {2397-3722},
  journal      = {npj Climate and Atmospheric Science},
  publisher    = {Springer Nature},
  title        = {{No decrease of tropical convection in individual deep convective systems with global warming}},
  doi          = {10.1038/s41612-025-01285-5},
  volume       = {9},
  year         = {2026},
}

@article{21217,
  abstract     = {This study investigates the mechanisms driving clustered convection and the breakdown of the Intertropical Convergence Zone (ITCZ) over the Western Pacific Warm Pool using high‐resolution cloud‐resolving simulations and machine‐learning sensitivity experiments. Results show that ITCZ breakdown episodes, marked by spatially homogeneous convection and weakened meridional moisture gradients, are triggered primarily by anomalous moisture advection linked to the equatorial Rossby‐wave activity. While large‐scale moisture advection regulates the background convective state strongly, it is the surface and low‐level meridional winds that dominate transitions between clustered and random convection. Simulations demonstrate that moisture alone can sustain convective clustering, but breakdown episodes are more persistent and widespread when coupled with southerly meridional advection. These findings confirm that wave‐driven advection acts as a regulatory mechanism, periodically disrupting convective clustering and reshaping the meridional moisture gradient. This modulation of organization by wave‐induced breakdown events is critical for understanding tropical convection variability and its implications for the climate system.},
  author       = {Casallas Garcia, Alejandro and Mark Tompkins, Adrian and Muller, Caroline J},
  issn         = {1477-870X},
  journal      = {Quarterly Journal of the Royal Meteorological Society},
  publisher    = {Wiley},
  title        = {{Moisture and wind effects of Rossby waves on Western Pacific Intertropical Convergence Zone breakdown events}},
  doi          = {10.1002/qj.70131},
  year         = {2026},
}

@article{18605,
  abstract     = {The response of clouds and moist-convective processes to heat loss to space by long-wave radiative cooling is an important feedback in the Earth's atmosphere. It is known that moist convection increases roughly in equilibrium with radiative cooling, an assumption often made in simplified models of the tropical atmosphere. In this study, we use an idealised two-dimensional model of the atmosphere introduced by Vallis et. al. and incorporate a bulk-cooling term, which is an idealisation of radiative cooling in the atmosphere. We comment briefly on the static stability of the system to dry and moist convection and characteris its moist convective response to changes in the bulk cooling. We find that, while the clear-sky regions of the model respond directly to the change in the cooling term, the regions dominated by moist convective plumes are insensitive to changes in cooling. Similar to previous findings from cloud-resolving models, we too find in our idealised setting that the majority of the increase in convection occurs via an increase in the areal coverage of convection, rather than its intensity. We argue that these small-scale convective processes are an upper bound on how quickly convective intensity can change to stay in equilibrium with radiative cooling.},
  author       = {Agasthya, Lokahith N and Muller, Caroline J and Cheve, Mathis},
  issn         = {1477-870X},
  journal      = {Quarterly Journal of the Royal Meteorological Society},
  number       = {766},
  publisher    = {Wiley},
  title        = {{Moist convective scaling: Insights from an idealised model}},
  doi          = {10.1002/qj.4902},
  volume       = {151},
  year         = {2025},
}

@article{20319,
  abstract     = {The time needed by deep convection to bring the atmosphere back to equilibrium is called convective adjustment timescale or simply adjustment timescale, typically denoted by . In the Community Atmospheric Model|Community Atmosphere Model (CAM),  is the convective available potential energy (CAPE) relaxation timescale and is 1 hr, worldwide. Observational evidence suggests that  is generally longer than 1 hr. Further, continental and oceanic convection are different in terms of the vigor of updrafts and can have different longevities. So using  hour worldwide in CAM has two potential caveats. A longer  improves the simulation of the mean climate. However, it does not address the land‐ocean heterogeneity of atmospheric deep convection. We investigate the prescription of two different CAPE relaxation timescales for land ( hr) and ocean ( to 4 hr). It is arguably an extremely crude parameterization of boundary layer control on atmospheric convection. We contrast a suite of 5‐year‐long simulations with two different  for land and ocean to having one  globally. The choice of longer  over ocean is guided by previous studies and inspired by observational pieces of evidence. Nonetheless, to complement our variable  experiments, we perform a simulation with  hr and  hrs. Most importantly, our key findings are immune to the exact values of prescribed  and . The CAM model, with two  values , improves convective‐stratiform rainfall partitioning and the Madden–Julian oscillation propagation characteristics.},
  author       = {GOSWAMI, BIDYUT B and Polesello, Andrea and Muller, Caroline J},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {9},
  publisher    = {Wiley},
  title        = {{An assessment of representing land‐ocean heterogeneity via CAPE relaxation timescale in the Community Atmospheric Model 6 (CAM6)}},
  doi          = {10.1029/2025ms005035},
  volume       = {17},
  year         = {2025},
}

@article{20590,
  abstract     = {Moist convection is a fundamental process occurring in the Earth's atmosphere. It plays a central role in the weather and climate of the Tropics, where, to first order, the heating of the atmosphere by convection is in balance with the cooling of the atmosphere by the emission of radiation to outer space. In this study, we use a cloud-resolving model in radiative–convective equilibrium with an imposed constant rate of radiative cooling and study the response of moist convection to varying this rate of radiative cooling. In particular, we study two types of simulation: varying air temperature (VAT) simulations, where the air temperature is allowed to adjust to the imposed radiative cooling, and constant air temperature (CAT) simulations, where the surface temperature is tuned to ensure that the atmospheric temperature profile in the domain is constant. We recover the previously known result that, in response to increasing radiative cooling, the area of convection expands rapidly, while the intensity of convection does not change. We find that this response is explained by the increased boundary-layer variability in simulations with greater radiative cooling, which compensates for the decreasing temperature by adding a larger initial velocity close to the cloud base. We also propose a fundamental scaling of the non-dimensional cumulus mass flux in moist convection, which is robust across models of different complexity. We aim to bridge the gap between highly idealised prototypes of moist convection, such as the “Rainy–Bénard convection” introduced by Vallis et al., and comprehensive cloud-resolving models.},
  author       = {Agasthya, Lokahith N and Muller, Caroline J},
  issn         = {1477-870X},
  journal      = {Quarterly Journal of the Royal Meteorological Society},
  publisher    = {Wiley},
  title        = {{Moist convection and radiative cooling: Dynamical response and scaling}},
  doi          = {10.1002/qj.70044},
  year         = {2025},
}

@article{19416,
  abstract     = {Recently, Biagioli and Tompkins (2023, https://doi.org/10.1029/2022ms003231) used a simple stochastic model to derive a dimensionless parameter to predict convective self aggregation (SA) development, which was based on the derivation of the maximum free convective distance ($d_{clr}$) expected in the pre-aggregated, random state. Our goal is to test and further investigate this hypothesis, namely that $d_{clr}$ can predict SA occurrence, using an ensemble of twenty-four distinct combinations of horizontal mixing, planetary boundary layer (PBL), and microphysical parameterizations. We conclude that the key impact of parameterization schemes on SA is through their control of the number of convective cores and their relative spacing, $d_{clr}$, which itself is impacted by cold-pool (CP) properties and mean updraft core size. SA is more likely when the convective core count is small, while CPs modify convective spacing via suppression in their interiors and triggering by gust-front convergence and collisions. Each parameterization scheme emphasizes a different mechanism. Subgrid-scale horizontal turbulent mixing mainly affects SA through the determination of convective core size and thus spacing. The sensitivity to the microphysics is mainly through rain evaporation and the subsequent impact on CPs, while perturbations to the ice cloud microphysics have a limited effect. Non-local PBL mixing schemes promote SA primarily by increasing convective inhibition through inversion entrainment and altering low cloud amounts, leading to fewer convective cores and larger $d_{clr}$. },
  author       = {Casallas Garcia, Alejandro and Tompkins, A.M. and Muller, Caroline J and Thompson, G.},
  issn         = {1942-2466},
  journal      = {Journal of Advances in Modeling Earth Systems},
  number       = {3},
  publisher    = {Wiley},
  title        = {{Sensitivity of self-aggregation and the key role of the free convection distance}},
  doi          = {10.1029/2024MS004791},
  volume       = {17},
  year         = {2025},
}

@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},
}

@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{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{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{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{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},
}

