Explore the vast range of titles available.

Top‐down entrainment shapes the vertical gradients of sensible heat, latent heat, and CO2 fluxes, influencing the interpretation of eddy covariance (EC) measurements in the unstable atmospheric surface layer (ASL). Using large eddy simulations for convective boundary layer flows, we demonstrate that decreased temperature gradients across the entrainment zone increase entrainment fluxes by enhancing the entrainment velocity, amplifying the asymmetry between top‐down and bottom‐up flux contributions. These changes alter scalar flux profiles, causing flux divergence or convergence and leading to the breakdown of the constant flux layer assumption (CFLA) in the ASL. As a result, EC‐measured fluxes either underestimate or overestimate “true” surface fluxes during divergence or convergence phases, contributing to energy balance non‐closure. The varying degrees of the CFLA breakdown are a fundamental cause for the non‐closure issue. These findings highlight the underappreciated role of entrainment in interpreting EC fluxes, addressing non‐closure, and understanding site‐to‐site variability in flux measurements. Plain Language Summary In the atmosphere over a heated surface, water vapor, carbon dioxide, and heat are transported from both the ground (bottom‐up) and the top of the air column (top‐down). The swirling motion of air within the column helps to even out the distribution of these quantities, known as “scalars.” Scalar fluxes measure how many molecules of these substances cross a unit area over time. At the surface, energy balance and plant processes control heat, water vapor, and carbon dioxide fluxes. However, fluxes at the top of the air column do not follow the same rules and abide by the same constraints as their ground counterpart. This study uses numerical simulations to show that when the temperature difference across a layer at the top of the boundary layer decreases, the boundary layer becomes deeper, increasing the transport of heat from the top. This causes changes in the slopes of flux profiles, disrupting the assumption that fluxes remain constant with height even close to the ground surface. As a result, measurements near the surface often underestimate or overestimate true surface fluxes, contributing to the much‐debated surface energy balance non‐closure problem. Key Points Entrainment‐modulated top‐down transport influences the slopes of the scalar flux profiles in the unstable atmospheric surface layer Variations in scalar flux profiles lead to differing degrees of failure in the constant flux layer assumption (CFLA) for different scalars The failure of the CFLA explains the non‐closure issue in the surface energy balance

It has been hypothesized that submesoscale flows play an important role in the vertical transport of climatically important tracers, due to their strong associated vertical velocities. However, the multi-scale, non-linear, and Lagrangian nature of transport makes it challenging to attribute proportions of the tracer fluxes to certain processes, scales, regions, or features. Here we show that criteria based on the surface vorticity and strain joint probability distribution function (JPDF) effectively decomposes the surface velocity field into distinguishable flow regions, and different flow features, like fronts or eddies, are contained in different flow regions. The JPDF has a distinct shape and approximately parses the flow into different scales, as stronger velocity gradients are usually associated with smaller scales. Conditioning the vertical tracer transport on the vorticity-strain JPDF can therefore help to attribute the transport to different types of flows and scales. Applied to a set of idealized Antarctic Circumpolar Current simulations that vary only in horizontal resolution, this diagnostic approach demonstrates that small-scale strain dominated regions that are generally associated with submesoscale fronts, despite their minuscule spatial footprint, play an outsized role in exchanging tracers across the mixed layer base and are an important contributor to the large-scale tracer budgets. Resolving these flows not only adds extra flux at the small scales, but also enhances the flux due to the larger-scale flows.

Part I of this study examines the relationship among winter cold anomalies over Eurasia, Ural blocking (UB), and the background conditions associated with Arctic warming over the Barents and Kara Seas (BKS) using reanalysis data. It is found that the intensity, persistence, and occurrence region of UB-related Eurasian cold anomalies depend strongly on the strength and vertical shear (VS) of the mean westerly wind (MWW) over mid–high-latitude Eurasia related to BKS warming. Observational analysis reveals that during 1951–2015 UB days are 64% (54%) more frequent during weak MWW (VS) winters, with 26.9 (28.4) days per winter, than during strong MWW (VS) winters. During weak MWW or VS winters, as frequently observed during 2000–15, persistent and large UB-related warming is seen over the BKS together with large and widespread midlatitude Eurasian cold anomalies resulting from increased quasi stationarity and persistence of the UB. By contrast, when the MWW or VS is strong as frequently observed during 1979–99, the cold anomaly is less intense and persistent and confined to a narrow region of Europe because of a rapid westward movement of the strong UB. For this case, the BKS warming is relatively weak and less persistent. The midlatitude cold anomalies are maintained primarily by reduced downward infrared radiation (IR), while the surface heat fluxes, IR, and advection all contribute to the BKS warming. Thus, the large BKS warming since 2000 weakens the meridional temperature gradient, MWW, and VS, which increases quasi stationarity and persistence of the UB (rather than its amplitude) and then leads to more widespread Eurasian cold events and further enhances the BKS warming.

Clouds are assumed to play an important role in the Arctic amplification process. This motivated a detailed investigation of cloud processes, including radiative and turbulent fluxes. Data from the aircraft campaign ACLOUD were analyzed with a focus on the mean and turbulent structure of the cloudy boundary layer over the Fram Strait marginal sea ice zone in late spring and early summer 2017. Vertical profiles of turbulence moments are presented from contrasting atmospheric boundary layers (ABLs) from 4 d. They differ by the magnitude of wind speed, boundary-layer height, stability, the strength of the cloud-top radiative cooling and the number of cloud layers. Turbulence statistics up to third-order moments are presented, which were obtained from horizontal-level flights and from slanted profiles. It is shown that both of these flight patterns complement each other and form a data set that resolves the vertical structure of the ABL turbulence well. The comparison of the 4 d shows that especially during weak wind, even in shallow Arctic ABLs with mixing ratios below 3 g kg−1, cloud-top cooling can serve as a main source of turbulent kinetic energy (TKE). Well-mixed ABLs are generated where TKE is increased and vertical velocity variance shows pronounced maxima in the cloud layer. Negative vertical velocity skewness points then to upside-down convection. Turbulent heat fluxes are directed upward in the cloud layer as a result of cold downdrafts. In two cases with single-layer stratocumulus, turbulent transport of heat flux and of temperature variance are both negative in the cloud layer, suggesting an important role of large eddies. In contrast, in a case with weak cloud-top cooling, these quantities are positive in the ABL due to the heating from the surface. Based on observations and results of a mixed-layer model it is shown that the maxima of turbulent fluxes are, however, smaller than the jump of the net terrestrial radiation flux across the upper part of a cloud due to the (i) shallowness of the mixed layer and (ii) the presence of a downward entrainment heat flux. The mixed-layer model also shows that the buoyancy production of TKE is substantially smaller in stratocumulus over the Arctic sea ice compared to subtropics due to a smaller surface moisture flux and smaller decrease in specific humidity (or even humidity inversions) right above the cloud top. In a case of strong wind, wind shear shapes the ABL turbulent structure, especially over rough sea ice, despite the presence of a strong cloud-top cooling. In the presence of mid-level clouds, cloud-top radiative cooling and thus also TKE in the lowermost cloud layer are strongly reduced, and the ABL turbulent structure becomes governed by stability, i.e., by the surface–air temperature difference and wind speed. A comparison of slightly unstable and weakly stable cases shows a strong reduction of TKE due to increased stability even though the absolute value of wind speed was similar. In summary, the presented study documents vertical profiles of the ABL turbulence with a high resolution in a wide range of conditions. It can serve as a basis for turbulence closure evaluation and process studies in Arctic clouds.

A submesoscale filament of dense water in the oceanic surface layer can undergo frontogenesis with a secondary circulation that has a surface horizontal convergence and downwelling in its center. This occurs either because of the mesoscale straining deformation or because of the surface boundary layer turbulence that causes vertical eddy momentum flux divergence or, more briefly, vertical momentum mixing. In the latter case the circulation approximately has a linear horizontal momentum balance among the baroclinic pressure gradient, Coriolis force, and vertical momentum mixing, that is, a turbulent thermal wind. The frontogenetic evolution induced by the turbulent mixing sharpens the transverse gradient of the longitudinal velocity (i.e., it increases the vertical vorticity) through convergent advection by the secondary circulation. In an approximate model based on the turbulent thermal wind, the central vorticity approaches a finite-time singularity, and in a more general hydrostatic model, the central vorticity and horizontal convergence are amplified by shrinking the transverse scale to near the model’s resolution limit within a short advective period on the order of a day.

Discovery and analysis of submesoscale variability O (0.3–30) km on the continental shelf is made possible by a high-resolution (Δ x = 75 m) Regional Oceanic Modeling System (ROMS) simulation of the Southern California Bight (SCB). This variability is manifest in ubiquitous yet ephemeral coherent structures: fronts, filaments, and vortices. Similar to their open-ocean counterparts, fronts and filaments on the shelf are identified by their strong vertical velocity, surface convergence, cyclonic vorticity, and horizontal density gradient. Life cycles of these features typically last 3–5 days, with the formation dominated by a horizontal advective tendency that increases density and velocity gradients (i.e., frontogenesis). The shape of the coastline and depth of the water column both influence the abundance and spatial orientation of shallow-water fronts and filaments. Closer to shore, fronts and filaments often align themselves parallel to isobaths, and headlands often act as sites of intense vorticity generation through bottom stress. A quasi-steady, approximate momentum balance among rotation, pressure gradient, and vertical mixing—known as turbulent thermal wind (TTW)—often is valid in the strong secondary circulations local to fronts and filaments. However, front and filament circulations subject to strong diurnal variation in surface heating and vertical mixing are inconsistent with steady-state TTW balance. The secondary circulations can induce ephemeral material trapping and substantial vertical heat fluxes on the shelf.

Air‐sea exchanges across oceanic fronts are critical in powering cloud formation, precipitation, and atmospheric storms. Oceanic submesoscale fronts of scales 1–10 km are characterized by strong sea surface temperature (SST) gradients. However, it remains elusive how submesoscale fronts affect the overlying atmosphere due to a lack of high‐resolution observations or models. Based on rare high‐resolution in situ observations in the Kuroshio Extension region, we quantify the air‐sea exchanges across an oceanic submesoscale front. The cross‐front SST and turbulent heat flux gradients reaches 2.4°C/km and 47 W/m2/km, respectively, far stronger than that typically found in mesoscale‐resolving products. The stronger SST gradient drives substantially stronger air‐sea fluxes and vertical mixing than mesoscale fronts, enhancing cloud formations. The intense air‐sea exchanges across submesoscale fronts are confirmed in idealized model simulations, but not resolved in mesoscale‐resolving climate models. Our finding provides essential knowledge for improving simulations of cloud formation, precipitation, and storms in climate models. Plain Language Summary Oceanic fronts, characterized by large sea surface temperature (SST) gradients, are ubiquitous in the global ocean. Through intense heat and moisture release, these oceanic fronts induce large horizontal gradient of sea level pressure or increasing vertical mixing intensity in the lower atmosphere, are critical in powering cloud formation, precipitation, and atmospheric storms, but are sensitive to SST gradients. Oceanic submesoscale fronts of spatial scales 1–10 km are characterized by strong SST gradients. However, our knowledge of how the submesoscale fronts affect the overlying atmosphere is by and large void, due to a lack of high‐resolution observations or models. Here, based on high‐resolution in situ observations and model simulations, we show that submesoscale fronts drive much stronger air‐sea exchanges and vertical mixing as compared to mesoscale fronts, with significant implications for marine atmosphere boundary layer changes and cloud formations. Limited by the coarse resolution, the intense air‐sea exchanges across submesoscale fronts are not resolved in mesoscale‐resolving climate models. These results highlight the importance of submesoscale air‐sea interactions and call for a proper representation of submesoscale air‐sea exchanges in the next generation of climate models. Key Points Observations show strong gradient in sea surface temperature and turbulent heat flux across a submesoscale oceanic front Submesoscale fronts drive substantially stronger air‐sea fluxes and vertical mixing than mesoscale fronts The intense air‐sea exchanges across submesoscale fronts are not resolved in mesoscale‐resolving climate models

Turbulent flow in a weakly convective marine atmospheric boundary layer (MABL) driven by geostrophic winds V g = 10 m s −1 and heterogeneous sea surface temperature (SST) is examined using fine-mesh large-eddy simulation (LES). The imposed SST heterogeneity is a single-sided warm or cold front with jumps Δ θ = (2, −1.5) K varying over a horizontal x distance of 1 km characteristic of an upper-ocean mesoscale or submesoscale front. The geostrophic winds are oriented parallel to the SST isotherms (i.e., the winds are alongfront). Previously, Sullivan et al. examined a similar flow configuration but with geostrophic winds oriented perpendicular to the imposed SST isotherms (i.e., the winds were across-front). Results with alongfront and across-front winds differ in important ways. With alongfront winds, the ageostrophic surface wind is weak, about 5 times smaller than the geostrophic wind, and horizontal pressure gradients couple the SST front and the atmosphere in the momentum budget. With across-front winds, horizontal pressure gradients are weak and mean horizontal advection primarily balances vertical flux divergence. Alongfront winds generate persistent secondary circulations (SC) that modify the surface fluxes as well as turbulent fluxes in the MABL interior depending on the sign of Δ θ . Warm and cold filaments develop opposing pairs of SC with a central upwelling or downwelling region between the cells. Cold filaments reduce the entrainment near the boundary layer top that can potentially impact cloud initiation. The surface-wind–SST-isotherm orientation is an important component of atmosphere–ocean coupling. The results also show frontogenetic tendencies in the MABL.

Mesoscale convective systems (MCSs) contribute significantly to summer precipitation in the tropics and midlatitude. Although soil moisture (SM) effects on convection are globally recognized, its specific role on mature MCSs in East China remains unclear. Using convection‐permitting simulations spanning 22 summers, we find that convective cores within mature MCSs preferentially develop on the drier side of strong SM gradients (∼200 km). This is evidenced by a 2.5‐fold increase in core occurrences downstream of the steepest 10% of SM gradients versus a near‐uniform surface. SM gradients shape sensible heat flux gradients via evapotranspiration, while upstream pre‐storm rain‐producing clouds suppress surface available energy. These processes favor MCSs through enhancing near‐surface temperature gradients which strengthen moisture convergence and zonal wind shear. Our results highlight the critical role of SM gradients in favoring MCS propagation in East China. As climate change intensifies SM heterogeneity, improved land‐surface representation offers potential for advancing rainfall prediction and projection.

The midlatitude poleward atmospheric energy transport increases in radiatively forced simulations of warmed climates across a range of models from comprehensive coupled general circulation models (GCMs) to idealized aquaplanet moist GCMs to diffusive moist energy balance models. These increases have been rationalized from two perspectives. The energetic (or radiative) perspective takes the atmospheric energy budget and decomposes energy flux changes (radiative forcing, feedbacks, or surface fluxes) to determine the energy transport changes required by the budget. The diffusive perspective takes the net effect of atmospheric macroturbulence to be a diffusive energy transport down-gradient, so transport changes can arise from changes in mean energy gradients or turbulent diffusivity. Here, we compare these perspectives in idealized moist, gray-radiation GCM simulations over a wide range of climates. The energetic perspective has a dominant role for radiative forcing in this GCM, with cancellation between the temperature feedback components that account for the GCM’s nonmonotonic energy transport changes in response to warming. Comprehensive CMIP5 simulations have similarities in the Northern Hemisphere to the idealized GCM, although a comprehensive GCM over several CO2 doublings has a distinctly different feedback evolution structure. The diffusive perspective requires a non-constant diffusivity to account for the idealized GCM-simulated changes, with important roles for the eddy velocity, dry static stability, and horizontal energy gradients. Beyond diagnostic analysis, GCM-independent a priori theories for components of the temperature feedback are presented that account for changes without knowledge of a perturbed climate state, suggesting that the energetic perspective is the more parsimonious one.