Surface Fluxes Research Articles

The Advective Flux Transport Model: Improving the Far Side with Active Regions Observed by STEREO 304 Å

Observations of the Sun’s photospheric magnetic field are often confined to the Sun–Earth line. Surface flux transport (SFT) models, such as the Advective Flux Transport (AFT) model, simulate the evolution of the photospheric magnetic field to produce magnetic maps over the entire surface of the Sun. While these models are able to evolve active regions that transit the near side of the Sun, new far-side side flux emergence is typically neglected. We demonstrate a new method for creating improved maps of the magnetic field over the Sun’s entire photosphere using data obtained by the Solar TErrestrial RElations Observatory (STEREO) mission. STEREO He ii 304 Å intensity images are used to infer the time, location, and total unsigned magnetic flux of far-side active regions. We have developed an automatic detection algorithm for finding and ingesting new far-side active region emergence into the AFT model. We conduct a series of simulations to investigate the impact of including active region emergence in AFT, both with and without data assimilation of magnetograms. We find that while He ii 304 Å can be used to improve surface flux models, care must taken to mitigate intensity surges from flaring events. We estimate that during Solar Cycle 24's maximum (2011–2015), 4–6 × 1022 Mx of flux is missing from SFT models that do not include far-side data. We find that while He ii 304 Å data alone can be used to create synchronic maps of photospheric magnetic field that resemble the observations, it is insufficient to produce a complete picture without direct magnetic observations from magnetographs.

The Sun’s Non-Potential Corona over Solar Cycle 24

The global magnetic field in the solar corona is known to contain free magnetic energy and magnetic helicity above that of a current-free (potential) state. But the strength of this non-potentiality and its evolution over the solar cycle remain uncertain. Here we model the corona over Solar Cycle 24 using a simplified magneto-frictional model that retains the magnetohydrodynamic induction equation but assumes relaxation towards force-free equilibrium, driven by solar surface motions and flux emergence. The model is relatively conservative compared to some others in the literature, with free energy approximately 20 – 25% of the potential field energy. We find that unsigned helicity is about a factor 10 higher at Maximum than Minimum, while free magnetic energy shows an even greater increase. The cycle averages of these two quantities are linearly correlated, extending a result found previously for active regions. Also, we propose a practical measure of eruptivity for these simulations, and show that this increases concurrently with the sunspot number, in accordance with observed coronal mass ejection rates. Whilst shearing by surface motions generates 50%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$50\\%$\\end{document} or more of the free energy and helicity in the corona, we show that active regions must emerge with their own internal helicity otherwise the eruptivity is substantially reduced and follows the wrong pattern over time.

Sublimation of Snow

Abstract Snow is a vital part of water resources, and sublimation may remove 10%–90% of snowfall from the system. To improve our understanding of the physics that govern sublimation rates, as well as how those rates might change with the climate, we deployed an array of four towers with over 100 instruments from NCAR’s Integrated Surface Flux System from November 2022 to June 2023 in the East River watershed, Colorado, in conjunction with the U.S. Department of Energy’s Surface Atmosphere Integrated Field Laboratory (SAIL) and the National Oceanic and Atmospheric Administration (NOAA)’s Study of Precipitation, the Lower Atmosphere and Surface for Hydrometeorology (SPLASH) campaigns. Mass balance observations, snow pits, particle flux sensors, and terrestrial lidar scans of the evolving snowfield demonstrated how blowing snow influences sublimation rates, which we quantified with latent heat fluxes measured by eddy-covariance systems at heights 1–20 m above the snow surface. Detailed temperature profiles at finer resolutions highlighted the role of the stable boundary layer. Four-stream radiometers indicated the important role of changing albedo in the energy balance and its relationship to water vapor losses. Collectively, these observations span scales from seconds to seasons, from boundary layer turbulence to valley circulation to mesoscale meteorology. We describe the field campaign, highlights in the observations, and outreach and education products we are creating to facilitate cross-disciplinary dialogue and convey relevant findings to those seeking to better understand Colorado River snow and streamflow.

On the Extent of Applicability of Various Non-linear Similarity Functions for Computation of Surface Fluxes under Stable Conditions in Numerical Models

On the Extent of Applicability of Various Non-linear Similarity Functions for Computation of Surface Fluxes under Stable Conditions in Numerical Models

Experiments and modelling of negative triangularity ASDEX Upgrade plasmas in view of DTT scenarios

The paper presents experimental and modelling results of a comparison of negative (NT) and positive (PT) triangularity ASDEX Upgrade (AUG) discharges using the plasma shapes presently foreseen in the DTT tokamak, under construction in Italy. This work is part of a broader effort of investigation to understand whether the good properties observed in NT scenarios in DIII-D and TCV may be extrapolated to the DTT device and more generally to DEMO future operations. The experimental results have shown a practical gain of running these AUG plasmas with only ECRH and mixed NBI+ECRH phases in negative triangularity, even if they access the H-mode. Indeed, the NT electron kinetic profiles recover in all cases the PT electron pressures inside mid-radius due to reduced transport in the region ρtor=0.7−0.9 , while exhibiting lower individual ELM (Edge Localised Mode) energy losses. The ion pressure and expected fusion performance are comparable in the case of similar densities. Integrated modelling has been performed using the transport solver ASTRA and the quasi-linear turbulent model TGLF, investigating the transport properties of these discharges. The modelling reproduces the experiments qualitatively with reasonable accuracy. Nonetheless, the heat transport in NT cases is partially overestimated. This may be because TGLF uses the Miller equilibrium, which approximates the magnetic flux surfaces as up-down symmetric. In the caseof these asymmetric NT shapes, the simulated outer surfaces lose part of the tilt with respect to the z-axis, reducing the upper δ < 0 effect. A numerical test to discern the impact of the geometry by symmetrically flipping the shape has shown a beneficial effect of the negative triangularity on heat transport.

Investigation on the edge cooling threshold of the density limit in the J-TEXT tokamak with limiter and divertor configurations

Cooling of the plasma edge is widely considered to be a key element in the density limit (DL) of tokamaks. This paper investigates the critical edge cooling threshold of the DL by exploring various plasma configurations in the joint-Texas experimental tokamak. Notably, significant differences in the edge electron temperature in the vicinity of the last closed flux surface were observed between the limiter and divertor configurations. However, the electron temperature drops to a similar level in the vicinity of q= 3 surface close to DL, independent of the magnetic field configuration. In addition, to evaluate the reliability of the critical edge cooling threshold, experiments were conducted by increasing the carbon impurity content to enhance the edge cooling rate. These experiments involved two approaches to increase the carbon impurity content: methane injection and penetration of a graphite solid source. Results from these experiments indicate that the temperature threshold of the q = 3 surfaces remains consistent even with a stronger edge cooling rate. The consistency observed in the electron temperature threshold near the q = 3 surface, regardless of magnetic configuration or edge cooling rate, could help refine existing theoretical and simulation works and improve the prediction accuracy of the DL disruption.

Shortwave Radiative Flux Variability Through the Lens of the Pacific Decadal Oscillation

AbstractThe variability of the shortwave radiative fluxes at the surface and top of atmosphere (TOA) is examined in a pre‐industrial modeling setup using the Pacific Decadal Oscillation (PDO) as a possible pacemaker of atmospheric decadal‐scale variability. Within models from the Coupled Model Intercomparison Project—Phase 6, downwelling shortwave radiation at the surface, the net shortwave fluxes at the surface and TOA, as well as cloud radiative effects show remarkably similar patterns associated with the PDO. Through ensemble simulations designed with a pure PDO pattern in the North Pacific only, we show that the PDO relates to about 20%–40% of the unforced year‐to‐year variability of these shortwave fluxes over the Northern Hemispheric continents. The sea surface temperature imprint on shortwave‐flux variability over land is larger for spatially aggregated time series as compared to smaller areas, due to the blurring effect of small‐scale atmospheric noise. The surface and TOA radiative flux anomalies associated with the PDO index range of [−1.64; 1.64] are estimated to reach up to ±6 Wm−2 for North America, ∓3 Wm−2 for India and ±2 Wm−2 for Europe. We hypothesize that the redistribution of clouds in response to a North Pacific PDO anomaly can impact the South Pacific and North Atlantic SSTs.

Flux‐gradient relations and their dependence on turbulence anisotropy

AbstractMonin–Obukhov similarity theory (MOST) is used in virtually all Earth System Models to parametrize the near‐surface turbulent exchanges and mean variable profiles. Despite its widespread use, there is high uncertainty in the literature about the appropriate parametrizations to use. In addition, MOST has limitations in very stable and unstable regimes, over heterogeneous terrain and complex orography, and has been found to represent the surface fluxes incorrectly. A new approach including turbulence anisotropy as a non‐dimensional scaling parameter has recently been developed and has been shown to overcome these limitations and generalize the flux‐variance relations to complex terrain. In this article, we analyze the flux‐gradient relations for five well‐known datasets, ranging from flat and homogeneous to slightly complex terrain. The scaling relations show substantial scatter and highlight the uncertainty in the choice of parametrization even over canonical conditions. We show that, by including information on turbulence anisotropy as an additional scaling parameter, the original scatter becomes well bounded and new formulations can be developed that drastically improve the accuracy of the flux‐gradient relations for wind shear () in unstable conditions and for temperature gradient () in both unstable and stable regimes. This analysis shows that both and are strongly dependent on turbulence anisotropy and allows us finally to settle the extensively discussed free convection regime for , which clearly exhibits a power law when anisotropy is accounted for. Furthermore, we show that the eddy diffusivities for momentum and heat and the turbulent Prandtl number are strongly dependent on anisotropy and that the latter goes to zero in the free convection limit. These results highlight the necessity to include anisotropy in the study of near‐surface atmospheric turbulence and lead the way for theoretically more robust simulations of the boundary layer over complex terrain.

Evaluating nine different air-sea flux algorithms coupled with CAM6

The accuracy of air-sea flux algorithms significantly impacts the simulation performance of climate models. In this study, we evaluated nine widely-used air-sea flux algorithms: those used in the BCC model, the CAM model, the ECMWF model, and the WRF model (with three options), as well as three other commonly-used algorithms: COARE, UA, and VM. An offline comparison revealed significant discrepancies in the flux calculations obtained with different algorithms. Subsequently, we conducted nine sensitivity experiments by coupling these algorithms into the Community Atmosphere Model version 6 (CAM6) to assess the extent and manner in which the selection of flux algorithms influenced the simulated climate. The online simulations introduced complexities due to interactions between the atmosphere and surface fluxes, making the results more intricate. The iterative process of the model created negative feedback between wind speed and momentum flux, as well as between air temperature and heat flux, leading to inconsistencies between online and offline test results. Based on the comprehensive evaluation metric Distance between Indices of Simulation and Observation (DISO), a concept of combined or segmented algorithm was introduced, and a segmented algorithm based on wind speed was found resulting in an overall superior simulation performance compared to the previous nine experiments.

Transport-driven toroidal rotation with general viscosity profile

Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.

An unmodified recycled polyethylene terephthalate (PET) nanofibrous membrane for water-in-oil and oil-in-water emulsions separation

A nanofibrous membrane made of unmodified recycled polyethylene terephthalate (rPET) was fabricated by electrospinning using plastic water bottle as raw materials. Using rPET polymer to fabricate membranes may help limit the discharge of plastic bottle waste into the environment and can treat oily wastewater simultaneously, especially emulsions. The variation in the concentration of rPET polymer and the properties of electrospun nanofibers, such as surface morphology and fiber diameter, had a significant impact on both the efficiency of oil droplet capturing and the surface wettability of the membrane. Water-in-oil emulsions and oil-in-water emulsions were separated using the nanofibrous membranes of different concentrations of rPET polymer, and the separation efficiencies at different concentrations of oil, and oil types were measured. The rPET nanofibrous membrane exhibited a water contact angle of 140.1° and was superoleophilic in air and underwater. The membranes exhibited excellent separation performances: >98% separation efficiency and a flux of more than 6500 L m-2 h-1 for water in hexadecane and octane emulsions; >95 % efficiency and a flux of 1509.70 L m-2 h-1 after three separation cycles of water-in-oil emulsions. For oil-in-water emulsions, the separation efficiencies were higher than 95 %, and flux was more than 3525.22 L m-2 h-1 for hexadecane octane and heptane in water emulsions, with an efficiency of >87.9 % and flux of 4186.20 L m-2 h-1 after three separation cycles. Moreover, the as-prepared membrane exhibited high performance of oil-water separation in harsh conditions, including strong acidic, strong basic, and strong salt solutions, with separation efficiencies of more than 98 % and flux of more than 12,700 L m-2 h-1. However, as water accumulated on the membrane surface, absorption and separation efficiency and flux were lowered significantly. This suggests the need to enhance the efficiency of the membrane and reduce the chance of water accumulation on the membrane surface.

Ocean surface radiation measurement best practices

Ocean surface radiation measurement best practices have been developed as a first step to support the interoperability of radiation measurements across multiple ocean platforms and between land and ocean networks. This document describes the consensus by a working group of radiation measurement experts from land, ocean, and aircraft communities. The scope was limited to broadband shortwave (solar) and longwave (terrestrial infrared) surface irradiance measurements for quantification of the surface radiation budget. Best practices for spectral measurements for biological purposes like photosynthetically active radiation and ocean color are only mentioned briefly to motivate future interactions between the physical surface flux and biological radiation measurement communities. Topics discussed in these best practices include instrument selection, handling of sensors and installation, data quality monitoring, data processing, and calibration. It is recognized that platform and resource limitations may prohibit incorporating all best practices into all measurements and that spatial coverage is also an important motivator for expanding current networks. Thus, one of the key recommendations is to perform interoperability experiments that can help quantify the uncertainty of different practices and lay the groundwork for a multi-tiered global network with a mix of high-accuracy reference stations and lower-cost platforms and practices that can fill in spatial gaps.

The Impact of Extended H2O Cross Sections on Temperate Anoxic Planet Atmospheres: Implications for Spectral Characterization of Habitable Worlds

JWST has created a new era of terrestrial exoplanet atmospheric characterization, and with it, the possibility to detect potential biosignature gases like CH4. Our interpretation of exoplanet atmospheric spectra, and the veracity of these interpretations, will be limited by our understanding of atmospheric processes and the accuracy of input modeling data. Molecular cross sections are essential inputs to these models. The photochemistry of temperate planets depends on photolysis reactions whose rates are governed by the dissociation cross sections of key molecules. H2O is one such molecule; the photolysis of H2O produces OH, a highly reactive and efficient sink for atmospheric trace gases. We investigate the photochemical effects of improved H2O cross sections on anoxic terrestrial planets as a function of host star spectral type and CH4 surface flux. Our results show that updated H2O cross sections, extended to wavelengths >200 nm, substantially impact the predicted abundances of trace gases destroyed by OH. The differences for anoxic terrestrial planets orbiting Sun-like host stars are greatest, showing changes of up to 3 orders of magnitude in surface CO levels, and over an order of magnitude in surface CH4 levels. These differences lead to observable changes in simulated planetary spectra, especially important in the context of future direct-imaging missions. In contrast, the atmospheres of planets orbiting M-dwarf stars are substantially less affected. Our results demonstrate a pressing need for refined dissociation cross-section data for H2O, where uncertainties remain, and other key molecules, especially at mid-UV wavelengths >200 nm.

Massive, long-lived electrostatic potentials in a rotating mirror plasma

Hot plasma is highly conductive in the direction parallel to a magnetic field. This often means that the electrical potential will be nearly constant along any given field line. When this is the case, the cross-field voltage drops in open-field-line magnetic confinement devices are limited by the tolerances of the solid materials wherever the field lines impinge on the plasma-facing components. To circumvent this voltage limitation, it is proposed to arrange large voltage drops in the interior of a device, but coexist with much smaller drops on the boundaries. To avoid prohibitively large dissipation requires both preventing substantial drift-flow shear within flux surfaces and preventing large parallel electric fields from driving large parallel currents. It is demonstrated here that both requirements can be met simultaneously, which opens up the possibility for magnetized plasma tolerating steady-state voltage drops far larger than what might be tolerated in material media.

Just How River‐Like Are Atmospheric Rivers?

AbstractAtmospheric Rivers (ARs) are synoptic‐scale conduits for poleward transport of heat and water, often associated with extreme rainfall. Using NASA surface heat flux observations and climate model simulations, we assess whether ARs are “rivers”, transporting heat and moisture over longer distances, or mostly local convergence. The observations indicate that ARs reduce extratropical surface energy fluxes, even during early development. This damping of surface fluxes during ARs is also simulated in the Goddard Institute for Space Studies ModelE, 2.1 (GISS‐E2.1) nudged to (MERRA2) reanalysis winds. Furthermore, water provenance tracers in GISS‐E2.1 identify the moisture source for ∼7,500 ARs globally during 2018–2022 as farther upstream and equatorward compared to climatology. These results quantitatively show that ARs source relatively less moisture from the surface beneath them and more from a greater distance than during non‐AR times.

Interactions between trade wind clouds and local forcings over the Great Barrier Reef: a case study using convection-permitting simulations

Abstract. Trade wind clouds are ubiquitous across the subtropical oceans, including the Great Barrier Reef (GBR), playing an important role in modulating the regional energy budget. These shallow clouds, however, are by their nature sensitive to perturbations in both their thermodynamic environment and microphysical background. In this study, we employ the Weather Research and Forecasting (WRF) model with a convection-permitting configuration at 1 km resolution to examine the sensitivity of the trade wind clouds to different local forcings over the GBR. A range of local forcings including coastal topography, sea surface temperature (SST), and local aerosol loading is examined. This study shows a strong response of cloud fraction and accumulated precipitation to orographic forcing both over the mountains and upwind over the GBR. Orographic lifting, low-level convergence, and lower troposphere stability are found to be crucial in explaining the cloud and precipitation features over the coastal mountains downwind of the GBR. However, clouds over the upwind ocean are more strongly constrained by the trade wind inversion, whose properties are, in part, regulated by the coastal topography. On the scales considered in this study, the warm-cloud fraction and the ensuant precipitation over the GBR show only a small response to the local SST forcing, with this response being tied to the surface flux and lower troposphere stability. Cloud microphysical properties, including cloud droplet number concentration, liquid water path, and precipitation, are sensitive to the changes in atmospheric aerosol population over the GBR. While cloud fraction shows little responses, a slight deepening of the simulated clouds is evident over the upwind region in correspondence to the increased aerosol number concentration. A downwind effect of aerosol loading on simulated cloud and precipitation properties is further noted.

Direct stellarator coil design using global optimization: application to a comprehensive exploration of quasi-axisymmetric devices

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a ‘global-to-local’ approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighbourhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform and quality of quasi-axisymmetry. The database of stellarators, which comprises approximately 200 000 coil sets, is available online and is called QUASR, for ‘quasi-symmetric stellarator repository’.

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

Abstract This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

Integration of surface-based and space-based atmospheric CO2 measurements for improving carbon flux estimates using a new developed 3-GAS inversion model

The ongoing necessity and developmental trajectory in atmospheric CO2 assimilation system revolve around enhancing the optimization efficacy of assimilation inversions through the synergistic integration of distinct observational datasets. In this study, we developed an new inversion system-Global Greenhouse-Gases Assimilation System (3-GAS), by coupling the global atmospheric transport model GEOS-Chem with the four-dimensional local ensemble transform Kalman filter (4D-LETKF) algorithm. We compiled an integrated dataset utilizing the widely recognized carbon observation data from ObsPack (Observation Package), GOSAT (The Greenhouse Gases Observing Satellite), and OCO-2 (Orbiting Carbon Observatory-2) to estimate global surface fluxes. The integration of multiple observational data significantly amplifies spatial coverage, enhances the variability in flux inversion, and fosters improvements in localized assimilation efficacy. The spatial distribution of the carbon flux from the combined assimilation emerges as a comprehensive outcome, born from the assimilation with individual observational data. This distribution is primarily influenced by the coverage attributes of the original observations and the meticulous selection of observations during the assimilation process. The validation with surface and aircraft observations reveals that the results of combined assimilation possess composite advantages, which are manifested in alleviating certain intrinsic errors inherent to individual assimilation. The combined assimilation improves the accuracy based on root-mean-square error reduced by 9%, 22% and 33% in comparison to surface-only, GOSAT-only and OCO-2-only assimilation. Eventually, predicated upon preliminary inversion results and regional characteristics, we endeavor to propose a more refined filtration scheme and error determination method for the integrated dataset, as the initiative for the efficient assimilation of multi-source data.

Experimental investigation on the boiling heat transfer characteristic of copper surfaces with different hydrophilicity

Surface hydrophilicity has a significant effect on boiling heat transfer. In order to systematically investigate the influence of hydrophilic degree on the boiling heat transfer performance, seven different surfaces with gradually changing contact angles (CAs) were prepared using the chemical etching technology. Meanwhile, Visual studies on pool boiling of four of these different hydrophilic surfaces were conducted using a high-speed camera at atmospheric pressure, with deionized water as the working fluid. The influence laws of the hydrophilicity coefficient and relative wall superheat on the enhancement of boiling heat transfer are systematically summarized. It is shown that the boiling heat transfer coefficient (HTC) and critical heat flux (CHF) of the heating surface increased as the CA decreased. Owing to its vertically and horizontally interlaced nanorods and nanoflower structures, the surface oxidized by NaOH/(NH4)2S2O8 has higher surface roughness and larger phase change heat transfer area, resulting in a super-hydrophilic property of the copper surface and a faster bubble detachment frequency. The CHF and boiling HTC were 1.16 and 1.59 times that of the untreated surface, respectively, and the onset of nucleate boiling was advanced by 3.1 K, greatly enhancing boiling heat transfer and exhibiting better heat transfer performance.