Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUX”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34

PROTOSTELLAR ACCRETION FLOWS DESTABILIZED BY MAGNETIC FLUX REDISTRIBUTION

Magnetic flux redistribution lies at the heart of the problem of star formation in dense cores of molecular clouds that are magnetized to a realistic level. If all of the magnetic flux of a typical core were to be dragged into the central star, the stellar field strength would be orders of magnitude higher than the observed values. This well-known magnetic flux problem can in principle be resolved through non-ideal MHD effects. Two-dimensional (axisymmetric) calculations have shown that ambipolar diffusion, in particular, can transport magnetic flux outward relative to matter, allowing material to enter the central object without dragging the field lines along. We show through simulations that such axisymmetric protostellar accretion flows are unstable in three dimensions to magnetic interchange instability in the azimuthal direction. The instability is driven by the magnetic flux redistributed from the matter that enters the central object. It typically starts to develop during the transition from the prestellar phase of star formation to the protostellar mass accretion phase. In the latter phase, the magnetic flux is transported outward mainly through advection by strongly magnetized low-density regions that expand against the collapsing inflow. The tussle between the gravity-driven infall and magnetically driven expansion leads to a highly filamentary inner accretion flow that is more disordered than previously envisioned. The efficient outward transport of magnetic flux by advection lowers the field strength at small radii, making the magnetic braking less efficient and the formation of rotationally supported disks easier in principle. However, we find no evidence for such disks in any of our rotating collapse simulations. We conclude that the inner protostellar accretion flow is shaped to a large extent by the flux redistribution-driven magnetic interchange instability. How disks form in such an environment is unclear.

79 ASTRONOMY AND ASTROPHYSICS↗

THE EVOLUTION OF SOLAR FLUX FROM 0.1 nm TO 160 {mu}m: QUANTITATIVE ESTIMATES FOR PLANETARY STUDIES

Understanding changes in the solar flux over geologic time is vital for understanding the evolution of planetary atmospheres because it affects atmospheric escape and chemistry, as well as climate. We describe a numerical parameterization for wavelength-dependent changes to the non-attenuated solar flux appropriate for most times and places in the solar system. We combine data from the Sun and solar analogs to estimate enhanced UV and X-ray fluxes for the young Sun and use standard solar models to estimate changing visible and infrared fluxes. The parameterization, a series of multipliers relative to the modern top of the atmosphere flux at Earth, is valid from 0.1 nm through the infrared, and from 0.6 Gyr through 6.7 Gyr, and is extended from the solar zero-age main sequence to 8.0 Gyr subject to additional uncertainties. The parameterization is applied to a representative modern day flux, providing quantitative estimates of the wavelength dependence of solar flux for paleodates relevant to the evolution of atmospheres in the solar system (or around other G-type stars). We validate the code by Monte Carlo analysis of uncertainties in stellar age and flux, and with comparisons to the solar proxies {kappa}{sup 1} Cet and EK Dra. The model is applied to the computation of photolysis rates on the Archean Earth.

79 ASTRONOMY AND ASTROPHYSICS↗

ON THE INJECTION OF HELICITY BY THE SHEARING MOTION OF FLUXES IN RELATION TO FLARES AND CORONAL MASS EJECTIONS

An investigation of helicity injection by photospheric shear motions is carried out for two active regions (ARs), NOAA 11158 and 11166, using line-of-sight magnetic field observations obtained from the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. We derived the horizontal velocities in the ARs from the differential affine velocity estimator (DAVE) technique. Persistent strong shear motions at maximum velocities in the range of 0.6-0.9 km s{sup -1} along the magnetic polarity inversion line and outward flows from the peripheral regions of the sunspots were observed in the two ARs. The helicities injected in NOAA 11158 and 11166 during their six-day evolution period were estimated as 14.16 Multiplication-Sign 10{sup 42} Mx{sup 2} and 9.5 Multiplication-Sign 10{sup 42} Mx{sup 2}, respectively. The estimated injection rates decreased up to 13% by increasing the time interval between the magnetograms from 12 minutes to 36 minutes, and increased up to 9% by decreasing the DAVE window size from 21 Multiplication-Sign 18 to 9 Multiplication-Sign 6 pixel{sup 2}, resulting in 10% variation in the accumulated helicity. In both ARs, the flare-prone regions (R2) had inhomogeneous helicity flux distribution with mixed helicities of both signs and coronal mass ejection (CME) prone regions had almost homogeneous distribution of helicity flux dominated by a single sign. The temporal profiles of helicity injection showed impulsive variations during some flares/CMEs due to negative helicity injection into the dominant region of positive helicity flux. A quantitative analysis reveals a marginally significant association of helicity flux with CMEs but not flares in AR 11158, while for the AR 11166, we find a marginally significant association of helicity flux with flares but not CMEs, providing evidence of the role of helicity injection at localized sites of the events. These short-term variations of helicity flux are further discussed in view of possible flare-related effects. This study suggests that flux motions and spatial distribution of helicity injection are important to understanding the complex nature of the magnetic flux system of the AR, and how it can lead to conditions favorable for eruptive events.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional prominence-hosting magnetic configurations: Creating a helical magnetic flux rope

The magnetic configuration hosting prominences and their surrounding coronal structure is a key research topic in solar physics. Recent theoretical and observational studies strongly suggest that a helical magnetic flux rope is an essential ingredient to fulfill most of the theoretical and observational requirements for hosting prominences. To understand flux rope formation details and obtain magnetic configurations suitable for future prominence formation studies, we here report on three-dimensional isothermal magnetohydrodynamic simulations including finite gas pressure and gravity. Starting from a magnetohydrostatic corona with a linear force-free bipolar magnetic field, we follow its evolution when introducing vortex flows around the main polarities and converging flows toward the polarity inversion line near the bottom of the corona. The converging flows bring the feet of different loops together at the polarity inversion line, where magnetic reconnection and flux cancellation happen. Inflow and outflow signatures of the magnetic reconnection process are identified, and thereby the newly formed helical loops wind around preexisting ones so that a complete flux rope grows and ascends. When a macroscopic flux rope is formed, we switch off the driving flows and find that the system relaxes to a stable state containing a helical magnetic flux rope embedded in an overlying arcade structure. A major part of the formed flux rope is threaded by dipped field lines that can stably support prominence matter, while the total mass of the flux rope is in the order of 4-5× 10{sup 14} g.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurements of EUV coronal holes and open magnetic flux

Coronal holes are regions on the Sun's surface that map the footprints of open magnetic field lines. We have developed an automated routine to detect and track boundaries of long-lived coronal holes using full-disk extreme-ultraviolet (EUV) images obtained by SOHO/EIT, SDO/AIA, and STEREO/EUVI. We measure coronal hole areas and magnetic flux in these holes, and compare the measurements with calculations by the potential field source surface (PFSS) model. It is shown that, from 1996 through 2010, the total area of coronal holes measured with EIT images varies between 5% and 17% of the total solar surface area, and the total unsigned open flux varies between (2-5)× 10{sup 22} Mx. The solar cycle dependence of these measurements is similar to the PFSS results, but the model yields larger hole areas and greater open flux than observed by EIT. The AIA/EUVI measurements from 2010-2013 show coronal hole area coverage of 5%-10% of the total surface area, with significant contribution from low latitudes, which is under-represented by EIT. AIA/EUVI have measured much enhanced open magnetic flux in the range of (2-4)× 10{sup 22} Mx, which is about twice the flux measured by EIT, and matches with the PFSS calculated open flux, with discrepancies in the location and strength of coronal holes. A detailed comparison between the three measurements (by EIT, AIA-EUVI, and PFSS) indicates that coronal holes in low latitudes contribute significantly to the total open magnetic flux. These low-latitude coronal holes are not well measured with either the He I 10830 line in previous studies, or EIT EUV images; neither are they well captured by the static PFSS model. The enhanced observations from AIA/EUVI allow a more accurate measure of these low-latitude coronal holes and their contribution to open magnetic flux.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetohydrostatic equilibrium. II. Three-dimensional multiple open magnetic flux tubes in the stratified solar atmosphere

A system of multiple open magnetic flux tubes spanning the solar photosphere and lower corona is modeled analytically, within a realistic stratified atmosphere subject to solar gravity. This extends results for a single magnetic flux tube in magnetohydrostatic equilibrium, described in Gent et al. Self-similar magnetic flux tubes are combined to form magnetic structures, which are consistent with high-resolution observations. The observational evidence supports the existence of strands of open flux tubes and loops persisting in a relatively steady state. Self-similar magnetic flux tubes, for which an analytic solution to the plasma density and pressure distribution is possible, are combined. We calculate the appropriate balancing forces, applying to the equations of momentum and energy conservation to preserve equilibrium. Multiplex flux tube configurations are observed to remain relatively stable for up to a day or more, and it is our aim to apply our model as the background condition for numerical studies of energy transport mechanisms from the solar surface to the corona. We apply magnetic field strength, plasma density, pressure, and temperature distributions consistent with observational and theoretical estimates for the lower solar atmosphere. Although each flux tube is identical in construction apart from the location of the radial axis, combinations can be applied to generate a non-axisymmetric magnetic field with multiple non-uniform flux tubes. This is a considerable step forward in modeling the realistic magnetized three-dimensional equilibria of the solar atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

INITIATION AND ERUPTION PROCESS OF MAGNETIC FLUX ROPE FROM SOLAR ACTIVE REGION NOAA 11719 TO EARTH-DIRECTED CME

An eruption event launched from the solar active region (AR) NOAA 11719 is investigated based on coronal EUV observations and photospheric magnetic field measurements obtained from the Solar Dynamic Observatory. The AR consists of a filament channel originating from a major sunspot and its south section is associated with an inverse-S sigmoidal system as observed in Atmospheric Imaging Assembly passbands. We regard the sigmoid as the main body of the flux rope (FR). There also exists a twisted flux bundle crossing over this FR. This overlying flux bundle transforms in shape similar to kink-rise evolution, which corresponds with the rise motion of the FR. The emission measure and temperature along the FR exhibits an increasing trend with its rising motion, indicating reconnection in the thinning current sheet underneath the FR. Net magnetic flux of the AR, evaluated at north and south polarities, showed decreasing behavior whereas the net current in these fluxes exhibits an increasing trend. Because the negative (positive) flux has a dominant positive (negative) current, the chirality of AR flux system is likely negative (left handed) in order to be consistent with the chirality of inverse S-sigmoidal FR. This analysis of magnetic fields of the source AR suggests that the cancelling fluxes are prime factors of the monotonous twisting of the FR system, reaching to a critical state to trigger kink instability and rise motion. This rise motion may have led to the onset of the torus instability, resulting in an Earth-directed coronal mass ejection, and the progressive reconnection in the thinning current sheet beneath the rising FR led to the M6.5 flare.

79 ASTRONOMY AND ASTROPHYSICS↗

Role of microstructure on flux expulsion of superconducting radio frequency cavities

The trapped residual magnetic flux during the cool-down due to the incomplete Meissner state is a significant source of radio frequency losses in superconducting radio frequency cavities. Here, in this study, we clearly correlate the niobium microstructure in elliptical cavity geometry and flux expulsion behavior. In particular, a traditionally fabricated Nb cavity half-cell from an annealed poly-crystalline Nb sheet after an 800 °C heat treatment leads to a bi-modal microstructure that ties in with flux trapping and inefficient flux expulsion. This non-uniform microstructure is related to varying strain profiles along the cavity shape. A novel approach to prevent this non-uniform microstructure is presented by fabricating a 1.3 GHz single cell Nb cavity with a cold-worked sheet and subsequent heat treatment leading to better flux expulsion after 800 °C/3 h. Microstructural evolution by electron backscattered diffraction-orientation imaging microscopy on cavity cutouts, and flux pinning behavior by dc-magnetization on coupon samples confirms a reduction in flux pinning centers with increased heat treatment temperature. The heat treatment temperature-dependent mechanical properties and thermal conductivity are reported. The significant impact of cold work in this study demonstrates clear evidence for the importance of the microstructure required for high-performance superconducting cavities with reduced losses caused by magnetic flux trapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Simulations of tokamak boundary plasma turbulence transport in setting the divertor heat flux width

The BOUT++ code has been used to simulate edge plasma electromagnetic (EM) turbulence and transport, and to study the role of EM turbulence in setting the scrape-off layer (SOL) heat flux width λ q . More than a dozen tokamak discharges from C-Mod, DIII-D, EAST, ITER and CFETR have been simulated with encouraging success. The parallel electron heat fluxes onto the target from the BOUT++ simulations of C-Mod, DIII-D, and EAST follow the experimental heat flux width scaling of the inverse dependence on the poloidal magnetic field. Further turbulence statistics analysis shows that the blobs are generated near the pedestal pressure peak gradient region inside the separatrix and contribute to the transport of the particle and heat in the SOL region. Transport simulations indicate two distinct regimes: drift dominant regime and turbulence dominant regime. Goldston's heuristic drift-based (HD) model yields a consistent divertor heat flux width in the drift dominant regime. For C-Mod enhanced D α H-mode discharges, drifts and turbulence are competing in setting the divertor heat flux width, possibly due to its compact machine size and good pedestal confinement. The simulations for ITER and CFETR indicate that divertor heat flux width λ q of the future machines may no longer follows the 1/B pol,OMP HD-based empirical (Eich) scalings and the HD model gives a pessimistic limit of divertor heat flux width. The simulation results show a transition from a drift dominant regime to a turbulence dominant regime from current machines to future machines such as ITER and CFETR for two reasons. (1) The magnetic drift-based radial transport decreases due to large CFETR and ITER machine sizes. (2) The SOL turbulence thermal diffusivity increases due to larger turbulent fluxes ejected from the pedestal into the SOL when operating in a different pedestal structure, from an ELM-free H-mode pedestal regime to a small and grassy ELM regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development and validation of non-axisymmetric heat flux simulations with 3D fields using the HEAT code

A new comprehensive module to simulate heat fluxes from three-dimensional (3D) magnetic fields has been implemented in the HEAT code. Especially compact tokamaks like SPARC require tools to predict and manage large heat fluxes. Existing release versions of HEAT can only simulate axisymmetric heat flux on 3D plasma facing components. The new module uses an M3D-C1 perturbed equilibrium and the MAFOT code to trace field lines of the perturbed 3D magnetic field. Heat flux is then assigned to the resulting footprints via a 3D layer model. The model distinguishes between the scrape-off layer, the magnetic lobes and the private flux region, and employs only 0D parameters like the layer width, diffusive spread and the last closed flux surface position in the perturbed edge to generate a heat flux profile. The magnitude is normalized to the total input power. Resulting heat flux simulations are compared and validated against infrared measurements in the DIII-D tokamak with applied 3D fields; good agreement is found for several cases. The new module can now be applied to the SPARC tokamak; a preliminary result for applied rotating 3D fields is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

OBSERVING FLUX ROPE FORMATION DURING THE IMPULSIVE PHASE OF A SOLAR ERUPTION

Magnetic flux ropes are believed to be an important structural component of coronal mass ejections (CMEs). While there exists much observational evidence of flux ropes after the eruption, e.g., as seen in remote-sensing coronagraph images or in situ solar wind data, the direct observation of flux ropes during CME impulsive phase has been rare. In this Letter, we present an unambiguous observation of a flux rope still in the formation phase in the low corona. The CME of interest occurred above the east limb on 2010 November 3 with footpoints partially blocked. The flux rope was seen as a bright blob of hot plasma in the Atmospheric Imaging Assembly (AIA) 131 A passband (peak temperature {approx}11 MK) rising from the core of the source active region, rapidly moving outward and stretching the surrounding background magnetic field upward. The stretched magnetic field seemed to curve-in behind the core, similar to the classical magnetic reconnection scenario in eruptive flares. On the other hand, the flux rope appeared as a dark cavity in the AIA 211 A passband (2.0 MK) and 171 A passband (0.6 MK); in these relatively cool temperature bands, a bright rim clearly enclosed the dark cavity. The bright rim likely represents the pileup of the surrounding coronal plasma compressed by the expanding flux rope. The composite structure seen in AIA multiple temperature bands is very similar to that in the corresponding coronagraph images, which consists of a bright leading edge and a dark cavity, commonly believed to be a flux rope.

79 ASTRONOMY AND ASTROPHYSICS↗

FINE-SCALE STRUCTURES OF FLUX ROPES TRACKED BY ERUPTING MATERIAL

We present Solar Dynamics Observatory observations of two flux ropes tracked out by material from a surge and a failed filament eruption on 2012 July 29 and August 4, respectively. For the first event, the interaction between the erupting surge and a loop-shaped filament in the east seems to 'peel off' the filament and add bright mass into the flux rope body. The second event is associated with a C-class flare that occurs several minutes before the filament activation. The two flux ropes are, respectively, composed of 85 {+-} 12 and 102 {+-} 15 fine-scale structures, with an average width of about 1.''6. Our observations show that two extreme ends of the flux rope are rooted in opposite polarity fields and each end is composed of multiple footpoints (FPs) of fine-scale structures. The FPs of the fine-scale structures are located at network magnetic fields, with magnetic fluxes from 5.6 Multiplication-Sign 10{sup 18} Mx to 8.6 Multiplication-Sign 10{sup 19} Mx. Moreover, almost half of the FPs show converging motion of smaller magnetic structures over 10 hr before the appearance of the flux rope. By calculating the magnetic fields of the FPs, we deduce that the two flux ropes occupy at least 4.3 Multiplication-Sign 10{sup 20} Mx and 7.6 Multiplication-Sign 10{sup 20} Mx magnetic fluxes, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

NEW VACUUM SOLAR TELESCOPE OBSERVATIONS OF A FLUX ROPE TRACKED BY A FILAMENT ACTIVATION

One main goal of the New Vacuum Solar Telescope (NVST) which is located at the Fuxian Solar Observatory is to image the Sun at high resolution. Based on the high spatial and temporal resolution NVST Hα data and combined with the simultaneous observations from the Solar Dynamics Observatory for the first time, we investigate a flux rope tracked by filament activation. The filament material is initially located at one end of the flux rope and fills in a section of the rope; the filament is then activated by magnetic field cancellation. The activated filament rises and flows along helical threads, tracking the twisted flux rope structure. The length of the flux rope is about 75 Mm, the average width of its individual threads is 1.11 Mm, and the estimated twist is 1π. The flux rope appears as a dark structure in Hα images, a partial dark and partial bright structure in 304 Å, and as a bright structure in 171 Å and 131 Å images. During this process, the overlying coronal loops are quite steady since the filament is confined within the flux rope and does not erupt successfully. It seems that, for the event in this study, the filament is located and confined within the flux rope threads, instead of being suspended in the dips of twisted magnetic flux.

79 ASTRONOMY AND ASTROPHYSICS↗

ELECTRON HEAT FLUX IN THE SOLAR WIND: ARE WE OBSERVING THE COLLISIONAL LIMIT IN THE 1 AU DATA?

Using statistically significant data at 1 AU, it has recently been shown (Bale et al.) that in the solar wind, when the Knudsen number K {sub T} (the ratio between the electron mean free path and the electron temperature scale height) drops below about 0.3, the electron heat flux q intensity rapidly approaches the classical collisional Spitzer-Härm limit. Using a fully kinetic model including the effect of Coulomb collisions and the expansion of the solar wind with heliocentric distance, we observe that the heat flux strength does indeed approach the collisional value for Knudsen numbers smaller than about 0.3 in very good agreement with the observations. However, closer inspection of the heat flux properties, such as its variation with the heliocentric distance and its dependence on the plasma parameters, shows that for Knudsen numbers between 0.02 and 0.3 the heat flux is not conveniently described by the Spitzer-Härm formula. We conclude that even though observations at 1 AU seem to indicate that the electron heat flux intensity approaches the collisional limit when the Knudsen drops below ∼0.3, the collisional limit is not a generally valid closure for a Knudsen larger than 0.01. Moreover, the good agreement between the heat flux from our model and the heat flux from solar wind measurements in the high-Knudsen number regime seems to indicate that the heat flux at 1 AU is not constrained by electromagnetic instabilities as both wave-particle and wave-wave interactions are neglected in our calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

SLIPPING MAGNETIC RECONNECTION TRIGGERING A SOLAR ERUPTION OF A TRIANGLE-SHAPED FLAG FLUX ROPE

We report the first simultaneous activities of the slipping motion of flare loops and a slipping eruption of a flux rope in 131 Å and 94 Å channels on 2014 February 2. The east hook-like flare ribbon propagated with a slipping motion at a speed of about 50 km s{sup –1}, which lasted about 40 minutes and extended by more than 100 Mm, but the west flare ribbon moved in the opposite direction with a speed of 30 km s{sup –1}. At the later phase of flare activity, there was a well developed ''bi-fan'' system of flare loops. The east footpoints of the flux rope showed an apparent slipping motion along the hook of the ribbon. Simultaneously, the fine structures of the flux rope rose up rapidly at a speed of 130 km s{sup –1}, much faster than that of the whole flux rope. We infer that the east footpoints of the flux rope are successively heated by a slipping magnetic reconnection during the flare, which results in the apparent slippage of the flux rope. The slipping motion delineates a ''triangle-shaped flag surface'' of the flux rope, implying that the topology of a flux rope is more complex than anticipated.

79 ASTRONOMY AND ASTROPHYSICS↗

Fermionic condensate in a conical space with a circular boundary and magnetic flux

The fermionic condensate is investigated in a (2+1)-dimensional conical spacetime in the presence of a circular boundary and a magnetic flux. It is assumed that on the boundary the fermionic field obeys the MIT bag boundary condition. For irregular modes, we consider a special case of boundary conditions at the cone apex, when the MIT bag boundary condition is imposed at a finite radius, which is then taken to zero. The fermionic condensate is a periodic function of the magnetic flux with the period equal to the flux quantum. For both exterior and interior regions, the fermionic condensate is decomposed into boundary-free and boundary-induced parts. Two integral representations are given for the boundary-free part for arbitrary values of the opening angle of the cone and magnetic flux. At distances from the boundary larger than the Compton wavelength of the fermion particle, the condensate decays exponentially, with the decay rate depending on the opening angle of the cone. If the ratio of the magnetic flux to the flux quantum is not a half-integer number for a massless field the boundary-free part in the fermionic condensate vanishes, whereas the boundary-induced part is negative. For half-integer values of the ratio of the magnetic flux to the flux quantum, the irregular mode gives a nonzero contribution to the fermionic condensate in the boundary-free conical space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Air–Sea Turbulent Heat Fluxes in Climate Models and Observational Analyses: What Drives Their Variability?

A traditional view is that the ocean outside of the tropics responds passively to atmosphere forcing, which implies that air–sea heat fluxes are mainly driven by atmosphere variability. This paper tests this viewpoint using state-of-the-art air–sea turbulent heat flux observational analyses and a climate model run at different resolutions. It is found herein that in midlatitude ocean frontal zones the variability of air–sea heat fluxes is not predominantly driven by the atmosphere variations but instead is forced by sea surface temperature (SST) variations arising from intrinsic oceanic variability. Meanwhile in most of the tropics and subtropics wind is the dominant driver of heat flux variability, and atmosphere humidity is mainly important in higher latitudes. The predominance of ocean forcing of heat fluxes found in frontal regions occurs on scales of around 700 km or less. Spatially smoothing the data to larger scales results in the traditional atmosphere-driving case, while filtering to retain only small scales of 5° or less leads to ocean forcing of heat fluxes over most of the globe. All observational analyses examined (1° OAFlux; 0.25° J-OFURO3; 0.25° SeaFlux) show this general behavior. A standard resolution (1°) climate model fails to reproduce the midlatitude, small-scale ocean forcing of heat flux: refining the ocean grid to resolve eddies (0.1°) gives a more realistic representation of ocean forcing but the variability of both SST and of heat flux is too high compared to observational analyses.

54 ENVIRONMENTAL SCIENCES↗

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

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.

54 ENVIRONMENTAL SCIENCES↗