Negative Heat Flux Research Articles

Biennial variability of boreal spring surface air temperature over India

In this study, we report significant biennial variability or oscillation (BO) in the boreal spring (March–May) Surface Air Temperature (SAT) over India and unravelled the causative mechanisms. The positive phase of the BO exhibit significant seasonal warming over India, whereas seasonal cooling is observed during the negative phase of BO. Heat wave days are more during the positive phase of BO compared to negative or neutral phase. The positive (negative) phase of BO is generally coherent with the central (eastern) Pacific warming (cooling) years. The anomalous low-level divergence associated with low-level anticyclonic circulation induces less cloudiness and intense surface solar radiation ovr India during the positive phase, favouring surface warming. The evolution of some positive and/or negative phases of BO without any large scale forcing from the equatorial Pacific suggested the possibility of alternate pathways. The strong anomalous upper-level (at 200 hPa) anticyclonic circulation provoked by mid-latitude Rossby waves is found contributing to the positive phase, thereby highlighting the role of dominant mid-latitude pathways in the biennial SAT variability in addition to El Niño forcing. The sinking motion associated with persistent high, and the associated adiabatic compression also supported surface heating during the positive phase of BO. On the other hand, the mid-latitude Rossby wave induced upper-level cyclonic circulation is found contributing to the negative phase. The sinking motion associated with persistent high, and the associated adiabatic compression also supported surface heating during the positive phase of BO. In contrast, negative soil temperature anomalies and high latent heat flux release to the atmosphere supported surface cooling during the negative phase.

Flow patterns and heat transfer of an idealized square city in non-uniform heat flux and different background wind conditions

The city scale buoyancy-driven flow is crucial for the wind and thermal environment in urban areas, exerting a significant impact on the city ventilation, pollutants dispersion and heat removal. The intra-urban variations of hot spots and heat flux are drawing more attention on beating the urban heat. This study aims at analyzing the influences of non-uniform heat flux and background wind on the city-scale heat removal ability and underlying mechanisms. The large eddy simulation model was used and verified with water tank experiments. Results show that positive heat flux and negative heat flux in the rural areas will enhance and depress the strength of the urban heat dome flow respectively compared to the adiabatic condition. The high-speed regions are wider and stronger when the rural area has positive heat flux condition. Under positive rural heat flux condition, the time-average mixed layer height in urban areas across the city center plane is nearly twice that under negative rural heat flux condition. Besides, the city-scale heat transfer coefficients in positive rural heat flux condition are higher than those under negative rural heat flux and adiabatic conditions. For cases of non-uniform heat flux in urban areas, the differences of time-average mixed layer height and heat transfer coefficients are not obvious. The heat transfer coefficients decrease initially and then increase with increasing background wind speed indicating the non-linear interaction between approaching wind and urban heat dome flow. Moreover, the growth rate of the heat transfer coefficient decreases when the non-dimensional wind speed exceeds 1.19.

Environmental conditions affecting air pollution dispersion over Richards Bay, South Africa

Air pollution dispersion over Richards Bay, South Africa is studied using satellite, reanalysis and in-situ measurements. Near-surface concentrations of trace gases and particulates reach high concentrations under dry stable ‘berg’ winds. Temporal records of CO, NO2, O3, PM2.5, SO2 (2000-2023) exhibit no trend and short-term dose exceedances of individual pollutants are rare and confined to winter mornings. Regional point-to-field correlation maps of a Richards Bay winter air pollution index found that anomalous warming of the west Indian Ocean accelerates the jet stream over Madagascar and pulls smoke- and dust-laden north-westerly winds across the South African plateau, implicating long-range transport. At daily and hourly time scales dewpoint temperature and boundary layer height best indicate changes in local air pollution dispersion that impact respiratory health. A case study of 15-17 June 2021 revealed a 10ºC thermal inversion caused by negative sensible heat fluxes of -40 W/m2 through the night, leading to SO2 concentrations above 100 ppb in the morning along the main access road (28.8ºS, 32ºE). Poor dispersion seldom happens in the summer when turbulent mixing and rainfall cleanse the air.

Interpretability of negative latent heat fluxes from eddy covariance measurements in dry conditions

Abstract. It is known from arid and semi-arid ecosystems that atmospheric water vapor can directly be adsorbed by the soil matrix. Soil water vapor adsorption was typically neglected and only recently received attention because of improvements in measurement techniques. One technique rarely explored for the measurement of soil water vapor adsorption is eddy covariance (EC). Soil water vapor adsorption may be detectable as downwardly directed (i.e., negative) EC latent heat (λE) flux measurements under dry conditions, but a systematic assessment of the use of negative λE fluxes from EC flux stations to characterize adsorption is missing. We propose a classification method to characterize soil water vapor adsorption, excluding conditions of dew and fog when λE derived from EC is not trustworthy due to stable atmospheric conditions. We compare downwardly directed λE fluxes from EC with measurements from weighing lysimeters for 4 years in a Mediterranean savanna ecosystem and 3 years in a temperate agricultural site. Our aim is to assess if overnight water inputs from soil water vapor adsorption differ between ecosystems and how well they are detectable by EC. At the Mediterranean site, the lysimeters measured soil water vapor adsorption each summer, whereas at the temperate site, soil water vapor adsorption was much rarer and was measured predominantly under an extreme drought event in 2018. During 30 % of nights in the 4-year measurement period at the Mediterranean site, the EC technique detected downwardly directed λE fluxes of which 88.8 % were confirmed to be soil water vapor adsorption by at least one lysimeter. At the temperate site, downwardly directed λE fluxes were only recorded during 15 % of the nights, with only 36.8 % of half hours matching simultaneous lysimeter measurement of soil water vapor adsorption. This relationship slightly improved to 61 % under bare-soil conditions and extreme droughts. This underlines that soil water vapor adsorption is likely a much more relevant process in arid ecosystems compared to temperate ones and that the EC method was able to capture this difference. The comparisons of the amounts of soil water vapor adsorption between the two methods revealed a substantial underestimation of the EC compared to the lysimeters. This underestimation was, however, comparable with the underestimation in evaporation by the eddy covariance and improved in conditions of higher turbulence. Based on a random-forest-based feature selection, we found the mismatch between the methods being dominantly related to the site's inherent variability in soil conditions, namely soil water status, and soil (surface) temperature. We further demonstrate that although the water flux is very small with mean values of 0.04 or 0.06 mm per night for EC or lysimeter, respectively, it can be a substantial fraction of the diel soil water balance under dry conditions. Although the two instruments substantially differ with regard to the measured ratio of adsorption to evaporation over 24 h with 64 % and 25 % for the lysimeter and EC methods, they are in either case substantial. Given the usefulness of EC for detecting soil water vapor adsorption as demonstrated here, there is potential for investigating adsorption in more climate regions thanks to the greater abundance of EC measurements compared to lysimeter observations.

Methods and metrics to assess the accuracy of heat flux data from a spark ignition IC engine

The measurement of in-cylinder heat transfer can be a valuable diagnostic tool for quantifying and improving IC engine performance. Increased heat transfer out of the combustion chamber negatively impacts overall efficiency, while heat transfer from the metal to the charge can increase knock propensity. The heat flux at various locations in an engine cylinder can be measured using heat flux probes consisting of two thermocouples – a surface thermocouple exposed to the changing in-cylinder temperatures, and a far field thermocouple which has a constant temperature at steady state operating conditions. The buildup of carbon deposits on the surface thermocouple can affect the heat flux measurements and lead to errors in the data. Heat flux measurements were performed on a proprietary single-cylinder research engine with instrumented heat flux probes in the cylinder head and cylinder liner. The engine was run rich with a high PMI fuel to expedite the buildup of carbon deposits. A metric was developed based on the output of the surface thermocouple to determine the cleanliness of the heat flux probes non-intrusively. Heat flux calculations were done using the FFT method and validated using the Cook-Felderman method. An uncertainty analysis was conducted on the heat flux data to verify that the negative heat flux observed was not due to errors in the distance between the two thermocouples or the measurement of temperature. Finally, to bolster confidence in the results, heat flux measurements were made for motoring conditions with varying coolant and intake air temperatures to gauge whether the trends in the output were in the expected direction.

Effects of geometry on the thermal-fluid dynamics of a transition water flow in a square heated asymmetric cavity channel

Cavity channels candidate as a promising passive cooling configuration in thermal management. The superior performance and the limited scientific investigation regarding cavity channels dictate the need for better understanding of the thermal-fluid dynamics. To the best knowledge of the authors, there are few works investigating the effect of the geometry on the thermal-fluid dynamics of cavity channels with no focus on thermal data and quantification of the recirculation effect which can depend on geometrical parameters. In this work, numerical RANS simulations were performed in ANSYS Fluent with the kT-kL-ω turbulence model for a transitional water flow in an asymmetric square cavity channel with different attachment angles in the range of 20°–90°. The optimal heat transfer performance in terms of uniform and high Nusselt numbers is provided by the configuration with an angle in the range of approximately 30°–40°. Recirculation plays an important role in enhancing heat transfer and its strength is associated with the penetration capacity in the boundary layer. Two regions of strong and weak recirculation are identified along the channel; when the recirculation is weak, other mechanisms such as the thinning of the thermal boundary layer play a more important role. By increasing the attachment angle it is shown the local thinning of the thermal boundary layer, an augmentation of the magnitude of the negative turbulent heat flux peak, and a similarity of the conductive sublayer thickness with the recirculation strength.

The Variation Characteristics of Stratospheric Circulation under the Interdecadal Variability of Antarctic Total Column Ozone in Early Austral Spring

Antarctic Total Column Ozone (TCO) gradually began to recover around 2000, and a large number of studies have pointed out that the recovery of the Antarctic TCO is most significant in the austral early spring (September). Based on the Bodeker Scientific Filled Total Column Ozone and ERA5 reanalysis dataset covering 1979–2019, the variation characteristics of the Antarctic TCO and stratospheric circulation for the TCO ‘depletion’ period (1979–1999) and the ‘recovery’ period (2000–2019) are analyzed in September. Results show that: (1) Stratospheric elements significantly related to the TCO have corresponding changes during the two eras. (2) The interannual variability of the TCO and the above-mentioned stratospheric circulation elements in the recovery period are stronger than those in the depletion period. (3) Compared with the depletion period, due to the stronger amplitude of the planetary wave 1, stronger Eliassen–Palm (EP) flux corresponds to EP flux convergence, larger negative eddy heat flux, and positive eddy momentum flux in the stratosphere during the recovery period. The polar temperature rises in the lower and middle stratosphere and the polar vortex weakens in the middle and upper stratosphere, accompanied by the diminished area of PSC. This contributes to the understanding of Antarctic ozone recovery.

Regional Controls on Climate and Weather Variability on the Southwest Coast of Peru

Southwestern Peru has an arid climate typical of subtropical west coasts bordering cold ocean currents. Mountain runoff is barely able to sustain urban needs and motivates this research. Using high-resolution satellite reanalysis products, the meso-scale climate and weather variability are explored via point-to-field regression. A time series spanning 1970–2022 of Tacna area (18 S, 70.2 W, 570 m) rainfall, potential evaporation, wind, and weather parameters were evaluated for thermodynamic and kinematic features. Although sea breezes draw marine air inland, they simultaneously generate low-level divergence and subsidence aloft. Potential evaporation in early summer causes water deficits that are rarely offset by late summer runoff from the Andes Mountains. Winter (May–September) showers from passing cold fronts are more frequent during El Niño. Warming of the tropical east Pacific accelerates subtropical westerly winds that lift over the coastal plains. Quasi-stationary Rossby wave patterns amplify transient troughs at 70 W, but the winter showers rarely exceed 4 mm/day due to low-level stability from negative heat fluxes over cool seawater offshore. Two winter wet spells were studied using satellite and surface data (July 2002, July 2009). Light showers were prominent in elevations from 400 to 900 m. An early summer dry spell was considered (November 2020), wherein southeast winds, coastal upwelling, and low dewpoint temperatures coincided with La Niña conditions. A rain-gauge transect showed that summer convection stays east of the Andes escarpment and seldom benefits the coastal plains. Thus, water resources in Tacna are strained beyond the carrying capacity.

Characterizing passive flow in nuclear prismatic modular reactor core channels: Temperature, velocity, and heat transfer analysis

This study aims to advance understanding and visualization of the passive flow inside the core channels of a dual-channel plenum-to-plenum test facility designed as representative geometry of nuclear prismatic modular reactors core. Simultaneous experimental measurements of local surface and coolant fluid temperatures and flow fields in addition to local heat transfer coefficients were provided by integrating sophisticated flow dynamics and heat transfer techniques at different axial locations (Z/L=0.04,0.28,0.41,0.59,0.77,0.96) along the facility heated channel. Results emphasize the significant dependency of thermal and flow fields, of naturally driven gas, on heating intensity and configuration. Interestingly, although the entire length of the channel is subjected to heating, an inflection in surface and gas temperature in addition to changes in velocity distribution were observed through the top section of the heated channel. This highlights the significance of the complexity of the facility which causes severe interactions, and potential recirculation and back mixing, between rising hot gas with the cold gas in the upper plenum. This is validated by the negative heat flux measurements collected at the top section of the channel indicating heat reversal where heat is transferred from air to the channel. Moreover, densitometric Froude number for different heating profiles were estimated and found to range from 0.66 to 0.75 indicating the existence of a jet flow at the hot channel outlet. The obtained results provide valuable benchmark data for validating further investigations by computational fluid dynamics (CFD) and to perform further investigations of thermal hydraulics in prismatic block reactors.

Geotextile protection of glacier: Observed and simulated impacts on energy and mass balance

The detailed physical processes involved in slowing glacier ablation by material cover remain poorly understood so far. In the present study, using the snow cover model SNOWPACK, the effect of geotextile cover on the energy and mass balance at the tongue of the Urumqi Glacier No. 1 (Chinese Tien Shan) was simulated between July 12, 2022 and August 31, 2022. The mass changes and the energy fluxes with and without material cover were compared. The results indicated that the geotextile covering reduced glacier ablation by approximately 68% compared to the ablation in the uncovered regions. The high solar reflectivity of the geotextile reduced the net short-wave radiation energy available for the melt by 45%. Thermal insulation of the geotextile reduced the sensible heat flux by 15%. In addition, the wet geotextile exerted a cooling effect through long-wave radiation and negative latent heat flux. This cooling effect reduced the energy available for ablation by 20%. Consequently, only 37% of the energy was used for melting compared to that used in the uncovered regions (67%). Sensitivity experiments revealed that the geotextile cover used at a thickness range of 0.045–0.090 m reduced the ice loss by approximately 68%–72%, and a further increase in the thickness of the geotextile cover led to little improvements. A higher temperature and greater wind speed increased glacier ablation, although their effects were small. When the precipitation was set to zero, it led to a significantly increased melt. Overall, the geotextile effectively protected the glacier tongue from rapid melting, and the observed results have provided inspiration for developing an effective and sustainable approach to protect the glaciers using geotextile cover.

Compressible turbulent convection in highly stratified adiabatic background

Buoyancy-driven turbulent convection leads to a fully compressible flow with a prominent top-down asymmetry of first- and second-order statistics when the adiabatic equilibrium profiles of temperature, density and pressure change very strongly across the convection layer. The growth of this asymmetry and the formation of an increasingly thicker stabilized sublayer with a slightly negative mean convective heat flux$J_c(z)$at the top of the convection zone is reported here by a series of highly resolved three-dimensional direct numerical simulations beyond the Oberbeck–Boussinesq and anelastic limits for dimensionless dissipation numbers,$0.1 \le D\le 0.8$, at fixed Rayleigh number$Ra=10^6$and superadiabaticity$\epsilon =0.1$. The highly stratified compressible convection regime appears for$D > D_{crit}\approx 0.65$, when density fluctuations collapse to those of pressure; it is characterized by an up to nearly 50 % reduced global turbulent heat transfer and a sparse network of focused thin and sheet-like thermal plumes falling through the top sublayer deep into the bulk.

Cloud and Precipitation in Low-Latitude Extratropical Cyclones Conditionally Sorted on CYGNSS Surface Latent and Sensible Heat Fluxes

Abstract Extratropical cyclones are the primary cause of midlatitude winds, precipitation, and clouds. Surface latent and sensible heat fluxes from the ocean impact cyclone intensity, but their role in moisture transport leading to cloud and precipitation is still under investigation. While numerical simulations can help establish such links, evaluating and constraining these simulations require observations. To this end, satellite-based observations of cloud and precipitation in low-latitude extratropical cyclones are examined for four distinct classifications based on the Cyclone Global Navigation Satellite System (CYGNSS) latent and sensible heat fluxes, averaged in the post-cold-frontal or warm-sector areas of the cyclones. Using a cyclone compositing approach, contrasts in cloud and precipitation in strong- versus weak-surface-flux conditions are examined. In the post-cold-frontal region, stronger latent or sensible heat fluxes are associated with lower precipitation rates and higher cloud opacity, indicating more vigorous shallow convection. However, larger sensible heat flux cases display larger cloud fraction, while larger latent heat flux cases exhibit lower cloud fraction, which could indicate differing cloud morphologies. In the comma region of the cyclones, clouds and precipitation depend on both cyclone strength and moisture availability. Consistent with this, larger cloud amount and precipitation are found for strong fluxes in the post-cold-frontal region, and weak or negative sensible heat fluxes, indicative of poleward warm advection, in the warm sector. The strong regional differences in the surface heat flux–cloud and precipitation relationships highlight the need for further investigation into moisture supply and transport in cyclones, while providing guidance for future work.

Experimental study on the influence of embankment slope direction on near-surface thermal conditions in permafrost region, Qinghai-Tibet Plateau

Heterogeneous ground thermal conditions caused by slope orientation are significant in Qinghai-Tibet Plateau because of the high altitude and strong solar radiation. Variation in solar radiation may result in contrasting thermal and upper boundary conditions, causing engineering problems at infrastructure with sloped embankments. Quantitative research on conditions on different man-made slopes is insufficient in permafrost regions, and consequently, planning for long-term effects on linear infrastructure is challenging. Net radiation, heat flux, near-surface temperature, soil moisture content, were recorded for one whole year (2021) at a field platform with eight slopes. Slope orientation affected net radiation (Rn) and ground heat flux (G), resulting in different energy balances between slopes. The north facing slope had a lowest Rn (174 W m-2), while the south slope had a maximum (239.8 W m-2). The north slope had a negative annual heat flux (-0.7 W m-2), while a maximum positive flux (1.5 W m-2) was recorded on the southeast slope. The annual mean surface temperature (Ts) was highest on the south slope and lowest on the north slope. The maximum difference of ground temperatures (within 30 cm depth) (Tg), was close to differences in Ts of over 4 °C. The maximum difference in soil moisture content during the thawing season (May–October) was 11% between the W and E slopes. Due to the different ground temperatures and soil moisture conditions, the annual number of freeze-thaw cycles was variable, and the maximum freeze-thaw cycles were 106 times in south slope and the minimum were 18 times in west slope. Those different thermal conditions between different slope orientations were significant, and pertinent to planning and maintenance of infrastructure. The results provide accurate boundary conditions for modeling to support development of new infrastructure that have embankment slopes in any direction, and for ensuring the stability of existing structures on Qinghai-Tibet Plateau.

Inferring the relationship between core-mantle heat flux and seismic tomography from mantle convection simulations

The heat flux pattern at Earth’s core-mantle boundary (CMB) imposes a heterogeneous boundary condition on core dynamics that may profoundly affect the geodynamo. Owing to the expected temperature dependence of seismic velocities, this pattern is classically approximated as proportional to the lowermost layer of seismic tomography models for the global mantle. Two biases however undermine such a simple linear relationship: 1) other contributions than thermal (compositional and mineralogical) influence seismic velocities and 2) the radial average is inherent to tomographic models whereas the local thermal state at the CMB is relevant for the heat flux. We analyze here simulations of thermochemical mantle convection where, owing to their spatial characteristics, specific mantle components are readily identified: hot thermochemical piles (TCPs), “normal” mantle (NM) and, when post-peroskite (pPv) is included, a cold region where this phase is present. Synthetic seismic velocities (i.e. from the mantle simulations) are then computed based on thermal, compositional and mineralogical sensitivities. A formalism to infer the CMB heat flux from these seismic shear velocity anomalies is derived. In this formalism, within each mantle population (i.e. TCPs, NM or pPv) the CMB heat flux vs. seismic anomalies follows a unique fitting function. The transition from one mantle population to another is marked by a jump in the seismic anomaly, i.e. a range of seismic anomalies in between two mantle populations corresponds to a similar CMB heat flux. Applying our formalism to the seismic anomalies from the mantle convection simulations provides far superior fits than the commonly used linear fits. The results highlight reduced negative heat flux anomalies beneath large low shear velocity provinces (LLSVPs), while positive heat flux anomalies are enhanced, both with respect to the classical linear interpretation.

Towards a better understanding of heat transfer and flow mechanisms in a cavity channel with numerical simulations

In this work, both numerical simulations performed with ANSYS Fluent using the kT-kL-ω turbulence model and experimental data are employed to investigate heat transfer in water flow in a uniformly heated (9500–19,350 W/m2) asymmetric channel for transitional Reynolds numbers ranging from about 3640 to 4970. The model predicts the experimental data mostly within about 15% and 16% respectively for the Nusselt number and pressure drops, with larger deviations for the Nusselt number occurring at the highest Reynolds numbers investigated. Simulations reveal the occurrence of a recirculation in the upper part of the channel. A new dimensionless parameter is introduced to characterize the strength of the recirculation with respect to the inertia flow which showed a distribution similar to the Nusselt number, thus suggesting an important role of the recirculation in enhancing the heat transfer. The thickness of both the thermal boundary layer and the conductive sublayer is determined both experimentally and numerically. While the former increases along the channel, the latter undergoes mostly a thinning, which can be attributed to the recirculation region. Peaks of negative turbulent heat flux are attributed to the flow recirculation and explain the minimum observed in the experimental temperature profiles. Among the RANS models, the kT-kL-ω model resulted to be the best in predicting the experimental local Nusselt numbers. This work may aid in achieving a more efficient design of passive configurations for thermal management applications.

The Influences of the Multi‐Scale Sea Surface Temperature Anomalies in the North Pacific on the Jet Stream in Winter

AbstractUsing Climate Forecast System Reanalysis (CFSR) data and numerical simulations, the impacts of the multi‐scale sea surface temperature (SST) anomalies in the North Pacific on the boreal winter atmospheric circulations are investigated. The basin‐scale SST anomaly as the Pacific Decadal Oscillation (PDO) pattern, a narrow meridional band of frontal‐scale smoothed SST anomaly in the subtropical front zone (STFZ) and the spatial dispersed eddy‐scale SST anomalies within the STFZ are the three types of forcings. The results of statistical methods find that all three oceanic forcings may correspond to the winter North Pacific jet changing with the similar pattern. Furthermore, several atmospheric general circulation model simulations are used to reveal the differences and detail processes of the three forcings. The basin‐scale cold PDO‐pattern SST anomaly first causes negative turbulent heat flux anomalies, atmospheric cooling, and wind deceleration in the lower atmosphere. Subsequently, the cooling temperature with an amplified southern lower temperature gradient and baroclinity brings a lagging middle warming because of the enhanced atmospheric eddy heat transport. The poleward and upward development of baroclinic fluctuations eventually causes the acceleration of the upper jet. The smoothed frontal‐ and eddy‐scales SST anomalies in the STFZ cause comparable anomalous jet as the basin‐scale by changing the upward baroclinic energy and Eliassen‐Palm fluxes. The forcing effects of multi‐scales SST anomalies coexist simultaneously in the mid‐latitude North Pacific, which can cause similar anomalous upper atmospheric circulations. This is probably why it is tricky to define the certain oceanic forcing to the specific observed atmospheric circulation variation.

3D Numerical Analysis of a Phase Change Material Solidification Process Applied to a Latent Thermal Energy Storage System

The techniques for releasing thermal energy accumulated in periods of high availability to meet the demand in periods of low energy supply contribute to the continuity of the cycles involved in thermodynamic processes. In this context, phase change materials are capable of absorbing and releasing large amounts of energy in relatively short periods of time and under specific operating conditions. However, phase change materials have low thermal conductivity and need to be coupled with high-thermal-conductivity materials so that the heat flux can be intensified and the energy absorption and release times can be controlled. This work aims to numerically study the solidification process of a phase change material inserted into a triplex tube heat exchanger with finned copper walls to intensify the thermal exchange between the phase change material and the cooling heat transfer fluid, water, that will receive the energy accumulated in the material. This work proposes the 3D numerical modeling of the triplex tube heat exchanger with finned walls and meets the need for numerical models that allow for the analysis of the full geometry of the latent heat thermal energy storage system and the thermal and fluid dynamic phenomena that are influenced by this geometry. Results of the temperature, liquid fractions and velocity fields during phase transformations are presented, analyzed and validated with experimental data, presenting average errors of below 5%. The total material discharge time was approximately 168 min, necessary for the complete solidification of the phase change material, with water injected into the triplex tube heat exchanger at a flow rate of 8.3 L/min and a temperature of 68 °C. The solidification process occurred more slowly in the same direction as the length of the triplex tube heat exchanger, and from 80% of the material in the solid state, the difference between the solidification time for z = 0 and z = 480 mm was 30 min. The fluid dynamic conditions developed in the latent heat thermal energy storage system promoted a maximum negative heat flux of −6423 w/m2 to the annular internal surface and −742 w/m2 to the annular external surface, representing a heat removal process nine times less intense on the external surface. The total energy released to the cooling heat transfer fluid was 239.56 kJ/kg.

Processes Contributing to North American Cold Air Outbreaks Based on Air Parcel Trajectory Analysis

Abstract Wintertime cold air outbreaks are periods of extreme cold, often persisting for several days and spanning hundreds of kilometers or more. They are commonly associated with intrusions of cold polar air into the midlatitudes, but it is unclear whether the air mass’s initial temperature in the Arctic or its cooling as it travels is the determining factor in producing a cold air outbreak. By calculating air parcel trajectories for a preindustrial climate model scenario, we study the role of the origin and evolution of air masses traveling over sea ice and land and resulting in wintertime cold air outbreaks over central North America. We find that not all Arctic air masses result in a cold air outbreak when advected into the midlatitudes. We compare trajectories that originate in the Arctic and result in cold air outbreaks to those that also originate in the Arctic but lead to median temperatures when advected into the midlatitudes. While about one-third of the midlatitude temperature difference can be accounted for by the initial height and temperature in the Arctic, the other two-thirds are a result of differences in diabatic heating and cooling as the air masses travel. Vertical mixing of cold surface air into the air mass while it travels dominates the diabatic cooling and contributes to the cold events. Air masses leading to cold air outbreaks experience more negative sensible heat flux from the underlying surface, suggesting that preconditioning to establish a cold surface is key to producing cold air outbreaks. Significance Statement Wintertime cold air outbreaks can cause temperatures to plummet tens of degrees below freezing over the northern United States, with the potential to damage agriculture, infrastructure, and human health. Accurate predictions under climate change could help mitigate these effects, but there is disagreement over whether cold air outbreaks have declined in line with the already-observed global warming trend or persisted in spite of it. Focusing on cold air outbreaks that originate from the Arctic, we find that there must be additional cooling of the traveling air mass by mixing with very cold surface air as it moves south over North America in order to result in a cold outbreak.

Strong Red Noise Ocean Forcing on Atlantic Multidecadal Variability Assessed from Surface Heat Flux: Theory and Application

Abstract The role of ocean forcing on Atlantic multidecadal variability (AMV) is assessed from the (downward) heat flux–SST relation in the framework of a new stochastic climate theory forced by red noise ocean forcing. Previous studies suggested that atmospheric forcing drives SST variability from monthly to interannual time scales, with a positive heat flux–SST correlation, while heat flux induced by ocean processes can drive SST variability at decadal and longer time scales, with a negative heat flux–SST correlation. Here, first, we develop a theory to show how the sign of heat flux–SST correlation is affected by atmospheric and oceanic forcing with time scale. In particular, a red noise ocean forcing is necessary for the sign reversal of heat flux–SST correlation. Furthermore, this sign reversal can be detected equivalently in three approaches: the low-pass correlation at lag zero, the unfiltered correlation at long (heat flux) lead, and the real part of the heat flux–SST coherence. Second, we develop a new scheme in combination with the theory to assess the magnitude and time scale of the red noise ocean forcing for AMV in the GFDL SPEAR model (Seamless System for Prediction and Earth System Research) and observations. In both the model and observations, the ocean forcing on AMV is in general comparable with the atmospheric forcing, with a 90% probability greater than the atmospheric forcing in observations. In contrast to the white noise atmospheric forcing, the ocean forcing has a persistence time comparable or longer than a year, much longer than the SST persistence of ∼3 months. This slow ocean forcing is associated implicitly with slow subsurface ocean dynamics. Significance Statement A new theoretical framework is developed to estimate the ocean forcing on Atlantic multidecadal variability form heat flux–SST relations in climate models and observation. Our estimation shows the ocean forcing is comparable with the atmospheric forcing and, in particular, has a slow time scale of years.

Vertical mixing and oxygen flux caused by daily sea breezes in a shallow stratified lake

Vertical transport caused by mixing is essential for understanding physical processes in lakes. However, mixing processes in shallow lakes are not well understood because of the lack of turbulence measurements. This study presents observations of vertical mixing and oxygen flux in a shallow lake, Lake Kitaura of the Lake Kasumigaura continuous lake system, which is located along the central eastern coast of the Japanese mainland. Mooring and microstructure surveys were conducted in August 2020. The vertical eddy diffusivity was estimated from the Ellison scale using high-frequency sampled temperature data from the mooring location, and the estimations were consistent with those observed from the microstructure profiler. The estimated eddy diffusivity revealed a mixing structure and oxygen flux in the lake during the study period. The daily cycle of stratification and mixing was caused by daily heating due to solar radiation and by winds from daily sea breezes, respectively. The daily stratification maximum occurred around noon, which suppressed vertical mixing. Vertical mixing was intensified due to sea breezes in the afternoon, which led to movement of oxygen from the surface layer to the bottom layer. The maximum vertical mixing was observed at 18:00. The oxygen concentration did not increase during nighttime when the negative surface heat flux was observed; thus, nighttime cooling may not primarily contribute to the vertical oxygen supply.