Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ADVECTION”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

On the Degradation of the Effectiveness of Nonlinear Diffusion Acceleration with Parallel Block Jacobi Splitting

MOOSE's radiation transport module Rattlesnake comprises a first order S N solver that can be accelerated with a flexible Nonlinear Diffusion Acceleration. At each NDA iteration, the streaming and collision operator of the S{sub N} equations are (approximately) inverted using a transport sweep. When executed on multiple processors parallel block Jacobi (PBJ) splitting is utilized lagging the processor sub-domain's incoming angular fluxes making the implemented method asynchronous, i.e. the convergence properties of NDA become dependent on the number of concurrent tasks and details of the partitioning. Within this work, we study the NDA convergence properties for the two-dimensional C5G7 benchmark problem PBJ with p = 1, 2, 4, 8, 12, 16, 24, 32, 48, 60, 72, 96 and 144 concurrent tasks. Rosa studied the convergence properties of PBJ splitting methods with and without transport synthetic acceleration (TSA) in slab geometry and two-dimensional Cartesian geometry, respectively. Rosa's iteration scheme lags both the incoming angular fluxes and the scattering source making it very similar to the method described here. Rosa found by means of Fourier analysis that convergence of unaccelerated PBJ can degrade for optically thin subdomains, i.e. for a large number of processors given a fixed domain size. Depending on the choice of parameters application of TSA can lead to an effective algorithm. In contrast to Rosa's work, this work focuses on PBJ splitting for NDA accelerated solution of the multigroup eigenvalue problem. Later, Zerr and Anistratov used their integral transport matrix method in conjunction with PBJ. This iteration scheme is different as it does not lag the scattering source. Nevertheless, Fourier analysis with Diamond Difference spatial discretization showed that the method degrades with increasing the overall size of the spatial domain while keeping the sub-domain size constant. In the case of an infinite medium, the observed spectral radius becomes unity. These references establish the iterative performance of PBJ for one-group fixed source problems discretized on rectangular cells, but for the Rattlesnake first order NDA method, an emphasis lies on multigroup eigenvalue problems. In this setting, the additional complication arises that the low order equations are not exactly inverted, but only to some iterative criterion. Further, as the low order problem is solved using Jacobian-Free Newton Krylov methods, a suitable preconditioner for the GMRES iterations is required. As the effectiveness of the preconditioners may decrease when more concurrent processes are used, we expect the low order solve to impact the overall effectiveness of the NDA method in parallel execution as well. Therefore, this work focuses on the question of how detrimental the asynchronicity of the high and low order iterative solutions is for the effectiveness of NDA solutions of realistic, unstructured grid, multigroup eigenvalue problems. The iterative convergence properties of Rattlesnake's NDA method for multigroup eigenvalue problems are studied for varying number of concurrent processes. Rattlesnake's NDA method is asynchronous because the high order transport sweep uses PBJ splitting and because the low order solve's preconditioners, algebraic multigrid and additive Schwarz, are asynchronous. We found that up to 144 processes the number of Picard iterations increases by a factor of 2-3 affecting but not destroying NDA's parallel performance for the C5G7 benchmark problem, a representative LWR (light water reactor) configuration. The solution of the NDA low order eigenvalue problem posed greater problems than expected for the solvers currently available in Rattlesnake mainly because of problems to effectively precondition the modified low-order problem. As a matter of fact, the low order problem becomes an advection-diffusion problem. Non-robustness was observed for the AMG preconditioned coarse mesh acceleration for more than 60 concurrent processes, which was fixed by increasing the number of permissible V-cycles per AMG call. For fine mesh acceleration AMG proved to be ineffective and we resorted to ASM. For both AMG and ASM, increasing number of linear iterations per nonlinear are observed when increasing the number of processes, but for ASM this problem appears to be worse. The results from this work also indicate that initial development efforts should be diverted from implementation of a high-order parallel sweeping algorithm to the effective and robust iterative solution of the low order NDA multigroup eigenvalue problem. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Vertical Transport of Water/Contaminant Due to Intermittent Water Influx on the Ground Surface of Vadose Zone - 20088

For movement of moisture and transport of contaminants in the vertical direction, two equations, the Richards equation and the advection-dispersion equation, are solved in a vadose zone for prescribed water and contaminant influx conditions at the ground surface. Two different boundary conditions are considered, one with a single peak of sinusoidal variation and four peaks of finite duration, each of which is sinusoidal in shape. The moisture migration is solved and then the contaminant transport, which is coupled to moisture movement, is solved. In the case of a single influx on the ground surface, at first the hydraulic head and moisture content show a peak. Then above the peak and below the ground surface, the peaks tend to die out, showing a tendency to return to the initial state. The vertical profile of concentration also shows a single peak, which soon becomes stagnant because downward water migration is faster than concentration profile movement. In the case of four intermittent array of peaks, the hydraulic head and moisture content profiles tend to merge and combine together. The concentration profiles show appearance of four peaks clearly. However, although the peak value of concentration remains quite constant, the values may vary with the sorption coefficient because of repeated influx of water and contaminant through the top boundary, the ground surface. The penetration depth decreases with an increasing sorption coefficient. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fundamental Study of Migration Behavior of Cesium under the Unsaturated Condition through One-dimensional Column Simulating Soil around the Disposal Site - 20262

Using a packed column simulating unsaturated soil around the disposal site, flow experiments of cesium ions were conducted in order to examine the difference of retardation effect between saturated and unsaturated conditions. The column was packed with silica sand particles of 350 μm in diameter. In this study, it was confirmed that the X-ray CT images of the column showed that the water saturation was mostly uniform within the column under unsaturated conditions. The values of water saturation in the flow experiments of Cs ions were in a range from 0.85 to 0.93. The retardation coefficients(Rd) were estimated by applying the one-dimensional advection-dispersion equation to the experimental results. The comparison of the Rd-values estimated experimentally and theoretically showed that the experimental value was estimated to be larger than the theoretical value. This means that the theoretical value evaluated through a conventional method would bring the conservative value as an apparent retardation coefficient, at least when the flow direction is from top to bottom as shown in this study. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Choosing the Best Modeling Platform for Radiological Risk Assessment Models - 20468

Radiological risk assessments, in the form of performance or safety assessments, are often required under regulations or guidance for remediation of contaminated land, decommissioning of contaminated buildings or structures, and radioactive waste disposal. These risk assessments are usually supported by fate and transport models that address decay and ingrowth of radionuclides, as well as their movement through engineered systems and the natural environment. These models are often projected thousands, or more, years into the future, largely because the radioactive species change through decay and ingrowth, and hence the magnitude of the radioactive effect changes with time. There are many computer codes that are available to address this type of modeling. They range from addressing specific pathways or processes such as infiltration of water, groundwater, surface water, air, biota, diffusion and advection of water and gases, to those that try to couple all processes together to evaluate the impact of fate and transport through the entire system to places in space and time to which access is assumed. These different types of codes are sometimes separated with the monikers process-level and systems-level codes, although it is often not clear that this separation does justice to the capabilities of the many codes that are available to evaluate fate and transport of radionuclides. The focus of this paper is the latter group of modeling codes. Several systems level modeling codes exist and are used. There are differences between these codes in terms of utility, flexibility, complexity and cost. The purpose of this paper is to compare a few of these codes in the context of work currently being performed by the International Atomic Energy Agency (IAEA) Modeling and Data for Radiological Impact Assessments (MODARIA) II Working Group 1 (WG1). The MODARIA II WG1's main focus is how stakeholder engaged decision analysis can, or should, be applied to radiological contamination problems so that better, longstanding, sustainable, solutions are reached. However, the WG1 also recognizes the potential impact of the modeling tools that are chosen to address radiological risk, which is often a primary objective of decision making for radiological problems. Other objectives might also be important, such as constraining costs, obtaining financing, minimizing impact on ecosystems, saving cultural resources, saving jobs, farmland, environmental justice, etc., in a full decision analysis for a given radiological contamination problem, but none of these other objectives have the same types of complex modeling needs as minimize radiological dose. Consequently, a further focus of WG1 is to evaluate the potential impacts on decision making of the choice of fate and transport, and risk assessment, modeling codes that are used to support decision making. The WG1 will produce a report at the end of 2020 that will focus on an approach to effective decision making and stakeholder engagement. The report will also consider the role that performance assessment modeling should play in the decision-making process, including the impact of the choice of modeling tools or computer codes on risk-informed decision making. Several sites around the World have been made available by Member States for these model comparisons, and several modeling tools have been considered. However, the focus of this paper is on two of the sites, one in Belgium and one in Ukraine, and on three of the tools: NORMALYSA (NORM And Legacy Site Assessment); GoldSim{sup C}, and AMBER{sup C}. The final report from this working group will also cover other modeling tools, including RESRAD, and PC-Cream{sup R}. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

The filamentary structure of mixing fronts and its control on reaction kinetics in porous media flows

We investigate the effective kinetics of a reaction front for mixing limited bimolecular reaction $A+B\rightarrow C$ in a porous medium. While Fickian diffusion predicts a scaling of the cumulative mass produced as $M_C \propto t^{1/2}$, we observe two time regimes in which the total product mass evolves faster then $t^{1/2}$. At early times the invading solute is organized in fingers of high velocity. Reactions take place only at the fingers boundaries whose surface grows linearly in time. We show that this configuration leads to a mass scaling $M_C \propto t^2$. When diffusion mixes reactants and destroy these finger structures, the effective reaction rate slows down and we relate it to the longitudinal advective spreading providing $M_C \propto \sigma_x$. The transition time between these two regimes is characterized by the diffusion time over the transverse fingers cross section.

reactive transport, multiscale methods, smoothed p↗
DOE OSTI2014

Interactions between cumulus convection and its environment as revealed by the MC3E sounding array

This study attempts to understand interactions between midlatitude convective systems and their environments through a heat and moisture budget analysis using the sounding data collected from the Midlatitude Continental Convective Clouds Experiment (MC3E) in central Oklahoma. Distinct large-scale structures and diabatic heating and drying profiles are presented for cases of weaker and elevated thunderstorms as well as intense squall line and supercell thunderstorm events during the campaign. The elevated cell events were nocturnal convective systems occurring in an environment having low convective available potential energy (CAPE) and a very dry boundary layer. In contrast, deeper convective events happened during the morning into early afternoon within an environment associated with large CAPE and a near-saturated boundary layer. As the systems reached maturity, the diagnosed diabatic heating in the latter deep convective cases was much stronger and of greater vertical extent than the former. Both groups showed considerable diabatic cooling in the lower troposphere, associated with the evaporation of precipitation and low-level clouds. The horizontal advection of moisture also played a dominant role in moistening the lower troposphere, particularly for the deeper convective events, wherein the near surface southeasterly flow allows persistent low-level moisture return from the Gulf of Mexico to support convection. The moisture convergence often was present before these systems develop, suggesting a strong correlation between the large-scale moisture convergence and convection. As a result, sensitivity tests indicated that the uncertainty in the surface precipitation and the size of analysis domain mainly affected the magnitude of these analyzed fields rather than their vertical structures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

RACORO continental boundary layer cloud investigations. 2. Large-eddy simulations of cumulus clouds and evaluation with in-situ and ground-based observations

A 60-hour case study of continental boundary layer cumulus clouds is examined using two large-eddy simulation (LES) models. The case is based on observations obtained during the RACORO Campaign (Routine Atmospheric Radiation Measurement [ARM] Aerial Facility [AAF] Clouds with Low Optical Water Depths [CLOWD] Optical Radiative Observations) at the ARM Climate Research Facility's Southern Great Plains site. The LES models are driven by continuous large-scale and surface forcings, and are constrained by multi-modal and temporally varying aerosol number size distribution profiles derived from aircraft observations. We compare simulated cloud macrophysical and microphysical properties with ground-based remote sensing and aircraft observations. The LES simulations capture the observed transitions of the evolving cumulus-topped boundary layers during the three daytime periods, and generally reproduce variations of droplet number concentration with liquid water content (LWC), corresponding to the gradient between the cloud centers and cloud edges at given heights. The observed LWC values fall within the range of simulated values; the observed droplet number concentrations are commonly higher than simulated, but differences remain on par with potential estimation errors in the aircraft measurements. Sensitivity studies examine the influences of bin microphysics versus bulk microphysics, aerosol advection, supersaturation treatment, and aerosol hygroscopicity. Simulated macrophysical cloud properties are found to be insensitive in this non-precipitating case, but microphysical properties are especially sensitive to bulk microphysics supersaturation treatment and aerosol hygroscopicity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Effect of dry large‐scale vertical motions on initial MJO convective onset

Abstract Anomalies of eastward propagating large‐scale vertical motion with ~30 day variability at Addu City, Maldives, move into the Indian Ocean from the west and are implicated in Madden‐Julian Oscillation (MJO) convective onset. Using ground‐based radar and large‐scale forcing data derived from a sounding array, typical profiles of environmental heating, moisture sink, vertical motion, moisture advection, and Eulerian moisture tendency are computed for periods prior to those during which deep convection is prevalent and those during which moderately deep cumulonimbi do not form into deep clouds. Convection with 3–7 km tops is ubiquitous but present in greater numbers when tropospheric moistening occurs below 600 hPa. Vertical eddy convergence of moisture in shallow to moderately deep clouds is likely responsible for moistening during a 3–7 day long transition period between suppressed and active MJO conditions, although moistening via evaporation of cloud condensate detrained into the environment of such clouds may also be important. Reduction in large‐scale subsidence, associated with a vertical velocity structure that travels with a dry eastward propagating zonal wavenumbers 1–1.5 structure in zonal wind, drives a steepening of the lapse rate below 700 hPa, which supports an increase in moderately deep moist convection. As the moderately deep cumulonimbi moisten the lower troposphere, more deep convection develops, which itself moistens the upper troposphere. Reduction in large‐scale subsidence associated with the eastward propagating feature reinforces the upper tropospheric moistening, helping to then rapidly make the environment conducive to formation of large stratiform precipitation regions, whose heating is critical for MJO maintenance.

Powell, Scott W.↗
DOE OSTI2015

Resolution‐dependent behavior of subgrid‐scale vertical transport in the Z hang‐ M c F arlane convection parameterization

Abstract To better understand the behavior of quasi‐equilibrium‐based convection parameterizations at higher resolution, we use a diagnostic framework to examine the resolution‐dependence of subgrid‐scale vertical transport of moist static energy as parameterized by the Zhang‐McFarlane convection parameterization (ZM). Grid‐scale input to ZM is supplied by coarsening output from cloud‐resolving model (CRM) simulations onto subdomains ranging in size from 8 × 8 to 256 × 256 km 2 . Then the ZM‐based parameterization of vertical transport of moist static energy for scales smaller than the subdomain size ( ) are compared to those directly calculated from the CRM simulations ( ) for different subdomain sizes. The ensemble mean decreases by more than half as the subdomain size decreases from 128 to 8 km across while decreases with subdomain size only for strong convection cases and increases for weaker cases. The resolution dependence of is determined by the positive‐definite grid‐scale tendency of convective available potential energy (CAPE) in the convective quasi‐equilibrium (QE) closure. Further analysis shows the actual grid‐scale tendency of CAPE (before taking the positive definite value) and behave very similarly as the subdomain size changes because they are both tied to grid‐scale advective tendencies. We can improve the resolution dependence of significantly by averaging the grid‐scale tendency of CAPE over an appropriately large area surrounding each subdomain before taking its positive definite value. Even though the ensemble mean decreases with increasing resolution, its variability increases dramatically. cannot capture such increase in the variability, suggesting the need for stochastic treatment of convection at relatively high spatial resolution (8 or 16 km).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Anomalous solute transport in saturated porous media: Relating transport model parameters to electrical and nuclear magnetic resonance properties

Abstract The advection‐dispersion equation (ADE) fails to describe commonly observed non‐Fickian solute transport in saturated porous media, necessitating the use of other models such as the dual‐domain mass‐transfer (DDMT) model. DDMT model parameters are commonly calibrated via curve fitting, providing little insight into the relation between effective parameters and physical properties of the medium. There is a clear need for material characterization techniques that can provide insight into the geometry and connectedness of pore spaces related to transport model parameters. Here, we consider proton nuclear magnetic resonance (NMR), direct‐current (DC) resistivity, and complex conductivity (CC) measurements for this purpose, and assess these methods using glass beads as a control and two different samples of the zeolite clinoptilolite, a material that demonstrates non‐Fickian transport due to intragranular porosity. We estimate DDMT parameters via calibration of a transport model to column‐scale solute tracer tests, and compare NMR, DC resistivity, CC results, which reveal that grain size alone does not control transport properties and measured geophysical parameters; rather, volume and arrangement of the pore space play important roles. NMR cannot provide estimates of more‐mobile and less‐mobile pore volumes in the absence of tracer tests because these estimates depend critically on the selection of a material‐dependent and flow‐dependent cutoff time. Increased electrical connectedness from DC resistivity measurements are associated with greater mobile pore space determined from transport model calibration. CC was hypothesized to be related to length scales of mass transfer, but the CC response is unrelated to DDMT.

Swanson, Ryan D.↗
DOE OSTI2015

Multiscale solute transport upscaling for a three‐dimensional hierarchical porous medium

Abstract A laboratory‐generated hierarchical, fully heterogeneous aquifer model (FHM) provides a reference for developing and testing an upscaling approach that integrates large‐scale connectivity mapping with flow and transport modeling. Based on the FHM, three hydrostratigraphic models (HSMs) that capture lithological (static) connectivity at different resolutions are created, each corresponding to a sedimentary hierarchy. Under increasing system ln K variances (0.1, 1.0, 4.5), flow upscaling is first conducted to calculate equivalent hydraulic conductivity for individual connectivity (or unit) of the HSMs. Given the computed flow fields, an instantaneous, conservative tracer test is simulated by all models. For the HSMs, two upscaling formulations are tested based on the advection‐dispersion equation (ADE), implementing space versus time‐dependent macrodispersivity. Comparing flow and transport predictions of the HSMs against those of the reference model, HSMs capturing connectivity at increasing resolutions are more accurate, although upscaling errors increase with system variance. Results suggest: (1) by explicitly modeling connectivity, an enhanced degree of freedom in representing dispersion can improve the ADE‐based upscaled models by capturing non‐Fickian transport of the FHM; (2) when connectivity is sufficiently resolved, the type of data conditioning used to model transport becomes less critical. Data conditioning, however, is influenced by the prediction goal; (3) when aquifer is weakly‐to‐moderately heterogeneous, the upscaled models adequately capture the transport simulation of the FHM, despite the existence of hierarchical heterogeneity at smaller scales. When aquifer is strongly heterogeneous, the upscaled models become less accurate because lithological connectivity cannot adequately capture preferential flows; (4) three‐dimensional transport connectivities of the hierarchical aquifer differ quantitatively from those analyzed for two‐dimensional systems.

Zhang, Mingkan↗
DOE OSTI2015

A hybrid pore‐scale and continuum‐scale model for solute diffusion, reaction, and biofilm development in porous media

Abstract It is a challenge to upscale solute transport in porous media for multispecies bio‐kinetic reactions because of incomplete mixing within the elementary volume and because biofilm growth can change porosity and affect pore‐scale flow and diffusion. To address this challenge, we present a hybrid model that couples pore‐scale subdomains to continuum‐scale subdomains. While the pore‐scale subdomains involving significant biofilm growth and reaction are simulated using pore‐scale equations, the other subdomains are simulated using continuum‐scale equations to save computational time. The pore‐scale and continuum‐scale subdomains are coupled using a mortar method to ensure continuity of solute concentration and flux at the interfaces. We present results for a simplified two‐dimensional system, neglect advection, and use dual Monod kinetics for solute utilization and biofilm growth. The results based on the hybrid model are consistent with the results based on a pore‐scale model for three test cases that cover a wide range of Damköhler (Da = reaction rate/diffusion rate) numbers for both homogeneous (spatially periodic) and heterogeneous pore structures. We compare results from the hybrid method with an upscaled continuum model and show that the latter is valid only for cases of small Damköhler numbers, consistent with other results reported in the literature.

Tang, Youneng↗
DOE OSTI2015

Atmospheric and oceanic conditions and the extremely low Bothnian Bay sea ice extent in 2014/2015

The winter of 2014/2015 was the first for which we can be certain that the northernmost embayment of the Baltic Sea, the Bothnian Bay, remained partially ice-free. We assess atmospheric and oceanic drivers responsible for this extreme event. In terms of atmospheric drivers, mild Bothnian Bay winters are characterized by positive North Atlantic Oscillation (NAO) associated with the advection of air from the North Atlantic. In winter 2014/2015, NAO was at a record high level and resulted in a northerly storm track associated with vigorous southwesterlies that drove ice toward the northeastern end of the bay. Furthermore, we found that the 1985–2015 Baltic Sea ice extent distribution differs from any other preceding 30 winter period since 1720 with a high confidence. Based on climate projections, it is very likely that in the next 50 years partially ice-free Bothnian Bay will go from being an extremely rare event to a new normal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Quantifying hydrate solidification front advancing using method of characteristics

Here, we develop a one-dimensional analytical solution based on the method of characteristics to explore hydrate formation from gas injection into brine-saturated sediments within the hydrate stability zone. Our solution includes fully coupled multiphase and multicomponent flow and the associated advective transport in a homogeneous system. Our solution shows that hydrate saturation is controlled by the initial thermodynamic state of the system and changed by the gas fractional flow. Hydrate saturation in gas-rich systems can be estimated by 1 – cl 0 /cl e when Darcy flow dominates, where cl 0 is the initial mass fraction of salt in brine, and cl e is the mass fraction of salt in brine at three-phase (gas, liquid, and hydrate) equilibrium. Hydrate saturation is constant, gas saturation and gas flux decrease, and liquid saturation and liquid flux increase with the distance from the gas inlet to the hydrate solidification front. The total gas and liquid flux is constant from the gas inlet to the hydrate solidification front and decreases abruptly at the hydrate solidification front due to gas inclusion into the hydrate phase. The advancing velocity of the hydrate solidification front decreases with hydrate saturation at a fixed gas inflow rate. This analytical solution illuminates how hydrate is formed by gas injection (methane, CO 2 , ethane, propane) at both the laboratory and field scales.

58 GEOSCIENCES↗
DOE OSTI2016

Three-dimensional modeling of outcrop-to-outcrop hydrothermal circulation on the eastern flank of the Juan de Fuca Ridge

We present three-dimensional simulations of coupled fluid and heat transport in the ocean crust, to explore patterns and controls on ridge-flank hydrothermal circulation on the eastern flank of the Juan de Fuca Ridge. Field studies have shown that there is large-scale fluid flow in the volcanic ocean crust in this region, including local convection and circulation between two basement outcrops separated by ~50 km. New simulations include an assessment of crustal permeability and aquifer thickness, outcrop permeability, the potential influence of multiple discharging outcrops, and a comparison between two-dimensional (profile) and three-dimensional representations of the natural system. Field observations that help to constrain new simulations include a modest range of flow rates between recharging and discharging outcrops, secondary convection adjacent to the recharging outcrop, crustal permeability determinations made in boreholes, and the lack of a regional seafloor heat flux anomaly as a consequence of advective heat loss from the crust. Three-dimensional simulations are most consistent with field observations when models use a crustal permeability of 3 ×10 -13 to 2 × 10 -12 m 2 , and the crustal aquifer is ≤300 m thick, values consistent with borehole observations. Here, we find fluid flow rates and crustal cooling efficiencies that are an order of magnitude greater in three-dimensional simulations than in two-dimensional simulations using equivalent properties. Lastly, simulations including discharge from an additional outcrop can also replicate field observations but tend to increase the overall rate of recharge and reduce the flow rate at the primary discharge site.

58 GEOSCIENCES↗
DOE OSTI2015

Possible mechanism of abrupt jump in winter surface air temperature in the late 1980s over the Northern Hemisphere

Abstract Possible cause of an abrupt warming in winter mean surface air temperature in the midlatitudes of the Northern Hemisphere in the late 1980s is investigated using observation and reanalysis data. To determine the timing of abrupt warming, we use a regime shift index based on detection of the largest significant differences between the mean values of two contiguous periods. Results show that the abrupt warming occurred in association with a regime shift after the 1980's in which the zonal mean sea level pressure (SLP) is significantly increased (decreased) at the latitude 25–35°N (60–70°N), in the form of north‐south dipole‐like SLP anomaly spanning the subtropics and high latitude. The dipole SLP anomaly can be attributed to a northward expansion of Hadley cell, a poleward broadening and intensification of the Ferrel cell, coupled with a collapse of polar cell. During the abrupt warming, strong anomalous southerly warm advection at the surface was induced by an enhanced and expanded Ferrel circulation, in association with a northward and downward shift of maximum center of northward eddy heat flux over the midlatitudes. An intensification of polar jet subsequent to regime shift may be instrumental in sustaining the warming up to more than 5 years.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗
DOE OSTI2015

The impact of vertical shear on the sensitivity of tropical cyclogenesis to environmental rotation and thermodynamic state

Abstract The impact of vertical wind shear on the sensitivity of tropical cyclogenesis to environmental rotation and thermodynamic state is investigated through idealized cloud‐resolving simulations of the intensification of an incipient vortex. With vertical shear, tropical cyclones intensify faster with a higher Coriolis parameter, f , irrespective of the environmental thermodynamic state. The vertical shear develops a vertically tilted vortex, which undergoes a precession process with the midlevel vortices rotating cyclonically around the surface center. With a higher f , the midlevel vortices are able to rotate continuously against the vertical shear, leading to the realignment of the tilted vortex and rapid intensification. With a lower f , the rotation is too slow such that the midlevel vortices are advected away from the surface center and the intensification is suppressed. The parameter, χ b , measuring the effect from the low‐entropy downdraft air on the boundary layer entropy, is found to be a good indicator of the environmental thermodynamic favorability for tropical cyclogenesis in vertical shear. Without vertical shear, tropical cyclones are found to intensify faster with a lower f by previous studies. We show this dependency on f is sensitive to the environmental thermodynamic state. The thermodynamical favorability for convection can be measured by χ m , which estimates the time it takes for surface fluxes to moisten the midtroposphere. A smaller χ m not only leads to a faster intensification due to a shorter period for moist preconditioning of the inner region but also neutralizes the faster intensification with a lower f due to enhanced peripheral convection.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗