Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Exploring the impacts of physics and resolution on aqua-planet simulations from a nonhydrostatic global variable-resolution modeling framework

Advances in computing resources are gradually moving regional and global numerical forecasting simulations towards sub-10 km resolution, but global high resolution climate simulations remain a challenge. The non-hydrostatic Model for Prediction Across Scales (MPAS) provides a global framework to achieve very high resolution using regional mesh refinement. Previous studies using the hydrostatic version of MPAS (H-MPAS) with the physics parameterizations of Community Atmosphere Model version 4 (CAM4) found notable resolution dependent behaviors. This study revisits the resolution sensitivity using the non-hydrostatic version of MPAS (NH-MPAS) with both CAM4 and CAM5 physics. A series of aqua-planet simulations at global quasi-uniform resolutions ranging from 240 km to 30 km and global variable resolution simulations with a regional mesh refinement of 30 km resolution over the tropics are analyzed, with a primary focus on the distinct characteristics of NH-MPAS in simulating precipitation, clouds, and large-scale circulation features compared to H-MPAS-CAM4. The resolution sensitivity of total precipitation and column integrated moisture in NH-MPAS is smaller than that in H-MPAS-CAM4. This contributes importantly to the reduced resolution sensitivity of large-scale circulation features such as the inter-tropical convergence zone and Hadley circulation in NH-MPAS compared to H-MPAS. In addition, NH-MPAS shows almost no resolution sensitivity in the simulated westerly jet, in contrast to the obvious poleward shift in H-MPAS with increasing resolution, which is partly explained by differences in the hyperdiffusion coefficients used in the two models that influence wave activity. With the reduced resolution sensitivity, simulations in the refined region of the NH-MPAS global variable resolution configuration exhibit zonally symmetric features that are more comparable to the quasi-uniform high-resolution simulations than those from H-MPAS that displays zonal asymmetry in simulations inside the refined region. Overall, NH-MPAS with CAM5 physics shows less resolution sensitivity compared to CAM4. These results provide a reference for future studies to further explore the use of NH-MPAS for high-resolution climate simulations in idealized and realistic configurations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Sources and pathways of the upscale effects on the Southern Hemisphere jet in MPAS-CAM4 variable-resolution simulations

Impacts of regional grid refinement on large-scale circulations (“upscale effects”) were detected in a previous study that used the Model for Prediction Across Scales-Atmosphere coupled to the physics parameterizations of the Community Atmosphere Model version 4. The strongest upscale effect was identified in the Southern Hemisphere jet during austral winter. This study examines the detailed underlying processes by comparing two simulations at quasi-uniform resolutions of 30 and 120 km to three variable-resolution simulations in which the horizontal grids are regionally refined to 30 km in North America, South America, or Asia from 120 km elsewhere. In all the variable-resolution simulations, precipitation increases in convective areas inside the high-resolution domains, as in the reference quasi-uniform high-resolution simulation. With grid refinement encompassing the tropical Americas, the increased condensational heating expands the local divergent circulations (Hadley cell) meridionally such that their descending branch is shifted poleward, which also pushes the baroclinically unstable regions, momentum flux convergence, and the eddy-driven jet poleward. This teleconnection pathway is not found in the reference high-resolution simulation due to a strong resolution sensitivity of cloud radiative forcing that dominates the aforementioned teleconnection signals. The regional refinement over Asia enhances Rossby wave sources and strengthens the upper level southerly flow, both facilitating the cross-equatorial propagation of stationary waves. Evidence indicates that this teleconnection pathway is also found in the reference high-resolution simulation. Lastly, the result underlines the intricate diagnoses needed to understand the upscale effects in global variable-resolution simulations, with implications for science investigations using the computationally efficient modeling framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Cloud and circulation feedbacks in a near‐global aquaplanet cloud‐resolving model

Abstract A near‐global aquaplanet cloud‐resolving model (NGAqua) with fixed meridionally varying sea‐surface temperature (SST) is used to investigate cloud feedbacks due to three climate perturbations: a uniform 4 K SST increase, a quadrupled‐CO 2 concentration, and both combined. NGAqua has a horizontal resolution of 4 km with no cumulus parameterization. Its domain is a zonally periodic 20,480 km‐long tropical channel, spanning 46°S–N. It produces plausible mean distributions of clouds, rainfall, and winds. After spin‐up, 80 days are analyzed for the control and increased‐SST simulations, and 40 days for those with quadrupled CO 2 . The Intertropical Convergence Zone width and tropical cloud cover are not strongly affected by SST warming or CO 2 increase, except for the expected upward shift in high clouds with warming, but both perturbations weaken the Hadley circulation. Increased SST induces a statistically significant increase in subtropical low cloud fraction and in‐cloud liquid water content but decreases midlatitude cloud, yielding slightly positive domain‐mean shortwave cloud feedbacks. CO 2 quadrupling causes a slight shallowing and a statistically insignificant reduction of subtropical low cloud fraction. Warming‐induced low cloud changes are strongly correlated with changes in estimated inversion strength, which increases modestly in the subtropics but decreases in the midlatitudes. Enhanced clear‐sky boundary layer radiative cooling in the warmer climate accompanies the robust subtropical low cloud increase. The probability distribution of column relative humidity across the tropics and subtropics is compared between the control and increased‐SST simulations. It shows no evidence of bimodality or increased convective aggregation in a warmer climate.

58 GEOSCIENCES↗
DOE OSTI2017

Sensitivity of transpiration to subsurface properties: Exploration with a 1-D model

The amount of moisture transpired by vegetation is critically tied to the moisture supply accessible to the root zone. In a Mediterranean climate, integrated evapotranspiration (ET) is typically greater in the dry summer when there is an uninterrupted period of high insolation. We present a 1-D model to explore the subsurface factors that may sustain ET through the dry season. The model includes a stochastic parameterization of hydraulic conductivity, root water uptake efficiency, and hydraulic redistribution by plant roots. Model experiments vary the precipitation, the magnitude and seasonality of ET demand, as well as rooting profiles and rooting depths of the vegetation. The results show that the amount of subsurface moisture remaining at the end of the wet winter is determined by the competition among abundant precipitation input, fast infiltration, and winter ET demand. The weathered bedrock retains math formula of the winter rain and provides a substantial moisture reservoir that may sustain ET of deep-rooted (>8 m) trees through the dry season. A small negative feedback exists in the root zone, where the depletion of moisture by ET decreases hydraulic conductivity and enhances the retention of moisture. Hence, hydraulic redistribution by plant roots is impactful in a dry season, or with a less conductive subsurface. Suggestions for implementing the model in the CESM are discussed.

58 GEOSCIENCES↗
DOE OSTI2016

The latitudinal variation of geoelectromagnetic disturbances during large ( Dst ≤ ₋100 nT) geomagnetic storms

Geoelectromagnetic disturbances (GMDs) are an important consequence of space weather that can directly impact many types of terrestrial infrastructure. In this paper, we analyze 30 years of SuperMAG magnetometer data from the range of magnetic latitudes 20°≤λ≤75° to derive characteristic latitudinal profiles for median GMD amplitudes. Based on this data, we obtain a parameterization of these latitudinal profiles of different types of GMDs, providing an analytical fit with Dst-dependent parameters. Finally, we also obtain probabilistic estimates for the magnitudes of “100 year” GMDs, finding that $\dot{_B}$ = 6.9 (3.60–12.9) nT/s should be expected at 45°≤λ< 50°, exceeding the 5 nT/s threshold for dangerous inductive heating.

58 GEOSCIENCES↗
DOE OSTI2017

Investigation of representing hysteresis in macroscopic models of two-phase flow in porous media using intermediate scale experimental data

Incorporating hysteresis into models is important to accurately capture the two phase flow behavior when porous media systems undergo cycles of drainage and imbibition such as in the cases of injection and post-injection redistribution of CO 2 during geological CO 2 storage (GCS). In the traditional model of two-phase flow, existing constitutive models that parameterize the hysteresis associated with these processes are generally based on the empirical relationships. This manuscript presents development and testing of mathematical hysteretic capillary pressure—saturation—relative permeability models with the objective of more accurately representing the redistribution of the fluids after injection. The constitutive models are developed by relating macroscopic variables to basic physics of two-phase capillary displacements at pore-scale and void space distribution properties. The modeling approach with the developed constitutive models with and without hysteresis as input is tested against some intermediate-scale flow cell experiments to test the ability of the models to represent movement and capillary trapping of immiscible fluids under macroscopically homogeneous and heterogeneous conditions. The hysteretic two-phase flow model predicted the overall plume migration and distribution during and post injection reasonably well and represented the postinjection behavior of the plume more accurately than the nonhysteretic models. In conclusion, based on the results in this study, neglecting hysteresis in the constitutive models of the traditional two-phase flow theory can seriously overpredict or underpredict the injected fluid distribution during post-injection under both homogeneous and heterogeneous conditions, depending on the selected value of the residual saturation in the nonhysteretic models.

58 GEOSCIENCES↗
DOE OSTI2017

A new process sensitivity index to identify important system processes under process model and parametric uncertainty

Hydrological models are always composed of multiple components that represent processes key to intended model applications. When a process can be simulated by multiple conceptual-mathematical models (process models), model uncertainty in representing the process arises. While global sensitivity analysis methods have been widely used for identifying important processes in hydrologic modeling, the existing methods consider only parametric uncertainty but ignore the model uncertainty for process representation. To address this problem, this study develops a new method to probe multimodel process sensitivity by integrating the model averaging methods into the framework of variance-based global sensitivity analysis, given that the model averaging methods quantify both parametric and model uncertainty. A new process sensitivity index is derived as a metric of relative process importance, and the index includes variance in model outputs caused by uncertainty in both process models and model parameters. For demonstration, the new index is used to evaluate the processes of recharge and geology in a synthetic study of groundwater reactive transport modeling. The recharge process is simulated by two models that converting precipitation to recharge, and the geology process is also simulated by two models of different parameterizations of hydraulic conductivity; each process model has its own random parameters. The new process sensitivity index is mathematically general, and can be applied to a wide range of problems in hydrology and beyond.

Dai, Heng↗
DOE OSTI2017

Experimental evidence for dynamic friction on rock fractures from frequency-dependent nonlinear hysteresis and harmonic generation: Low-Frequency Torsional Dynamic Friction

Frictional properties affect the propagation of high-amplitude seismic waves across rock fractures and faults. Laboratory evidence suggests that these properties can be measured in active seismic surveys, potentially offering a route to characterizing friction in situ. Here, we present experimental results from a subresonance torsional modulus and attenuation apparatus that utilizes micron-scale sinusoidal oscillations to probe the nonlinear stress-strain relation at a range of strain amplitudes and rates. Nonlinear effects are further quantified using harmonic distortion; however, time series data best illuminate underlying physical processes. The low-frequency stress-strain hysteretic loops show stiffening at the sinusoid's static ends, but stiffening is reduced above a threshold frequency. This shape is determined by harmonic generation in the strain; the stress signal has no harmonics, confirming that the fractured sample is the source of the nonlinearity. These qualitative observations suggest the presence of rate-dependent friction and are consistent between fractures in three different rock types. We propose that static friction at the low strain rate part of the cycle, when given sufficient “healing” time at low oscillation frequencies, causes this stiffening cusp shape in the hysteresis loop. While rate-and-state friction is commonly used to represent dynamic friction, it cannot capture static friction or negative slip velocities. So we implement another dynamic friction model, based on the work of Dahl, which describes this process and produces similar results. Since the two models have a similar form, parameterizations of field data could constraint fault model inputs, such as specific location velocity strengthening or weakening properties.

58 GEOSCIENCES↗
DOE OSTI2017

Investigating the dependence of SCM simulated precipitation and clouds on the spatial scale of large-scale forcing at SGP [Investigating the scale dependence of SCM simulated precipitation and cloud by using gridded forcing data at SGP]

Large-scale forcing data, such as vertical velocity and advective tendencies, are required to drive single-column models (SCMs), cloud-resolving models, and large-eddy simulations. Previous studies suggest that some errors of these model simulations could be attributed to the lack of spatial variability in the specified domain-mean large-scale forcing. This study investigates the spatial variability of the forcing and explores its impact on SCM simulated precipitation and clouds. A gridded large-scale forcing data during the March 2000 Cloud Intensive Operational Period at the Atmospheric Radiation Measurement program's Southern Great Plains site is used for analysis and to drive the single-column version of the Community Atmospheric Model Version 5 (SCAM5). When the gridded forcing data show large spatial variability, such as during a frontal passage, SCAM5 with the domain-mean forcing is not able to capture the convective systems that are partly located in the domain or that only occupy part of the domain. This problem has been largely reduced by using the gridded forcing data, which allows running SCAM5 in each subcolumn and then averaging the results within the domain. This is because the subcolumns have a better chance to capture the timing of the frontal propagation and the small-scale systems. As a result, other potential uses of the gridded forcing data, such as understanding and testing scale-aware parameterizations, are also discussed.

58 GEOSCIENCES↗
DOE OSTI2017

Using ARM Observations to Evaluate Climate Model Simulations of Land-Atmosphere Coupling on the U.S. Southern Great Plains

Several independent measurements of warm-season soil moisture and surface atmospheric variables recorded at the ARM Southern Great Plains (SGP) research facility are used to estimate the terrestrial component of land-atmosphere coupling (LAC) strength and its regional uncertainty. The observations reveal substantial variation in coupling strength, as estimated from three soil moisture measurements at a single site, as well as across six other sites having varied soil and land cover types. The observational estimates then serve as references for evaluating SGP terrestrial coupling strength in the Community Atmospheric Model coupled to the Community Land Model. These coupled model components are operated in both a free-running mode and in a controlled configuration, where the atmospheric and land states are reinitialized daily, so that they do not drift very far from observations. Although the controlled simulation deviates less from the observed surface climate than its free-running counterpart, the terrestrial LAC in both configurations is much stronger and displays less spatial variability than the SGP observational estimates. Preliminary investigation of vegetation leaf area index (LAI) substituted for soil moisture suggests that the overly strong coupling between model soil moisture and surface atmospheric variables is associated with too much evaporation from bare ground and too little from the vegetation cover. Lastly, these results imply that model surface characteristics such as LAI, as well as the physical parameterizations involved in the coupling of the land and atmospheric components, are likely to be important sources of the problematical LAC behaviors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Distinct Contributions of Ice Nucleation, Large-Scale Environment, and Shallow Cumulus Detrainment to Cloud Phase Partitioning With NCAR CAM5

Mixed-phase clouds containing both liquid droplets and ice particles occur frequently at high latitudes and in the midlatitude storm track regions. Simulations of the cloud phase partitioning between liquid and ice hydrometeors in state-of-the-art global climate models are still associated with large biases. For this study, the phase partitioning in terms of liquid mass phase ratio (MPR liq , defined as the ratio of liquid mass to total condensed water mass) simulated from the NCAR Community Atmosphere Model version 5 (CAM5) is evaluated against the observational data from A-Train satellite remote sensors. Modeled MPR liq is significantly lower than observations on the global scale, especially in the Southern Hemisphere (e.g., Southern Ocean and the Antarctic). Sensitivity tests with CAM5 are conducted to investigate the distinct contributions of heterogeneous ice nucleation, shallow cumulus detrainment, and large-scale environment (e.g., winds, temperature, and water vapor) to the low MPR liq biases. Our results show that an aerosol-aware ice nucleation parameterization increases the MPR liq especially at temperatures colder than -20°C and significantly improves the model agreements with observations in the Polar regions in summer. The decrease of threshold temperature over which all detrained cloud water is liquid from 268 to 253 K enhances the MPR liq and improves the MPR liq mostly over the Southern Ocean. By constraining water vapor in CAM5 toward reanalysis, modeled low biases in many geographical regions are largely reduced through a significant decrease of cloud ice mass mixing ratio.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Modeling East Asian Dust and Its Radiative Feedbacks in CAM4-BAM

East Asian dust and its radiative feedbacks are analyzed by the use of the fourth version of the Community Atmosphere Model (CAM4) with a bulk aerosol model parameterization (BAM) for the dust size distribution (CAM4-BAM). Two numerical experiments are conducted and intercompared: one with (Active) and one without (Passive) the radiative effects of dust aerosols. This CAM4-BAM captures the main spatial distribution of the dust aerosol optical depth (AOD) and the dust surface concentrations over East Asia, with positive correlations with the local observational data on annual and seasonal means. A comparative analysis of the Active and Passive experiments reveals that consideration of the dust-radiation interaction can significantly reduce dust emissions, loading, transport, and dry and wet depositions over East Asia, which is opposite to the enhanced dust cycle over North Africa. Further analysis of the contrasting dust-radiation feedbacks between North Africa and East Asia shows that over North Africa, the dust radiative forcing significantly increases the surface temperature and 10-m wind speed, whereas it decreases the surface temperature and the surface wind speeds over East Asia. These contrasting radiative effects, in turn, result in distinct dust cycle changes over these two regions. Thus, mechanistic analysis reveals that the radiative contrasts between East Asia and North Africa are mainly due to the differences in their regional surface albedo, dust vertical distribution and size distribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Drivers of Seasonal Variability in Marine Boundary Layer Aerosol Number Concentration Investigated Using a Steady State Approach

Marine boundary layer (MBL) aerosol particles affect the climate through their interaction with MBL clouds. Although both MBL clouds and aerosol particles have pronounced seasonal cycles, the factors controlling seasonal variability of MBL aerosol particle concentration are not well-constrained. In this paper an aerosol budget is constructed representing the effects of wet deposition, free-tropospheric entrainment, primary surface sources, and advection on the MBL accumulation mode aerosol number concentration (N a ). These terms are further parameterized, and by assuming that on seasonal timescales N a is in steady state, the budget equation is rearranged to form a diagnostic equation for N a based on observable variables. Using data primarily collected in the subtropical northeast Pacific during the MAGIC campaign (Marine ARM (Atmospheric Radiation Measurement) GPCI (GCSS Pacific Cross-section Intercomparison) Investigation of Clouds), estimates of both mean summer and winter N a concentrations are made using the simplified steady-state model and seasonal mean observed variables, and are found to match well with the observed N a . To attribute the modeled difference between summer and winter aerosol concentrations to individual observed variables (e.g. precipitation rate, free-tropospheric aerosol number concentration), a local sensitivity analysis is combined with the seasonal difference in observed variables. This analysis shows that despite wintertime precipitation frequency being lower than summer, the higher winter precipitation rate accounted for approximately 60% of the modeled seasonal difference in N a , which emphasizes the importance of marine stratocumulus precipitation in determining MBL aerosol concentrations on longer time scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Comparisons of Ice Water Path in Deep Convective Systems Among Ground-Based, GOES, and CERES-MODIS Retrievals

Retrievals of convective cloud microphysical properties based on passive satellite imagery are difficult. To help quantify their uncertainties, ice water paths (IWPs) retrieved from the NASA Clouds and the Earth's Radiant Energy System project using Geostationary Operational Environmental Satellite (GOES) and Terra/Aqua MODerate-resolution Imaging Spectroradiometer observations are compared with IWPs retrieved from Next-Generation Radar (NEXRAD) observations over a large domain (32°N to 40°N and 105°W to 91°W) during the 2011 Midlatitude Continental Convective Clouds Experiment field campaign. Based on comparisons of pixel-level (4 km × 4 km) daytime IWP retrievals from NEXRAD and GOES, it is found that NEXRAD- and GOES-retrieved mean IWPs are 2.03 and 1.83 kg m -2 , respectively, for ice-phase cloud in thick anvil area. Their mean difference of 0.20 kg m -2 (with 95% confidence interval: 0.14-0.26 kg m -2 ) is within the uncertainty of NEXRAD retrievals. However, the low correlation between pixel-to-pixel comparisons indicates a large variation in GOES-retrieved IWP. For mixed-phase clouds in thick anvil areas, in addition to IWPs, total water paths (TWPs, sum of ice and liquid water path) are estimated with aid of aircraft measurements for NEXRAD retrievals and corrected using a TWP parameterization for GOES retrievals. The mean values of estimated TWPs from NEXRAD (corrected using aircraft in situ measurements) and GOES are similar. GOES and Clouds and the Earth's Radiant Energy System-MODerate-resolution Imaging Spectroradiometer-retrieved IWPs/TWPs generally do not exceed 5 kg m - 2. Large differences and low correlations exist between satellite and NEXRAD retrievals in stratiform rain areas. Possible reasons for the differences between retrievals are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Investigation of Turbulent Entrainment-Mixing Processes With a New Particle-Resolved Direct Numerical Simulation Model

Here, a new particle-resolved three dimensional direct numerical simulation (DNS) model is developed that combines Lagrangian droplet tracking with the Eulerian field representation of turbulence near the Kolmogorov microscale. Six numerical experiments are performed to investigate the processes of entrainment of clear air and subsequent mixing with cloudy air and their interactions with cloud microphysics. The experiments are designed to represent different combinations of three configurations of initial cloudy area and two turbulence modes (decaying and forced turbulence). Five existing measures of microphysical homogeneous mixing degree are examined, modified, and compared in terms of their ability as a unifying measure to represent the effect of various entrainment-mixing mechanisms on cloud microphysics. Also examined and compared are the conventional Damköhler number and transition scale number as a dynamical measure of different mixing mechanisms. Relationships between the various microphysical measures and dynamical measures are investigated in search for a unified parameterization of entrainment-mixing processes. The results show that even with the same cloud water fraction, the thermodynamic and microphysical properties are different, especially for the decaying cases. Further analysis confirms that despite the detailed differences in cloud properties among the six simulation scenarios, the variety of turbulent entrainment-mixing mechanisms can be reasonably represented with power-law relationships between the microphysical homogeneous mixing degrees and the dynamical measures.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

On Which Microphysical Time Scales to Use in Studies of Entrainment-Mixing Mechanisms in Clouds

The commonly used time scales in entrainment-mixing studies are examined in this paper to seek the most appropriate one, based on aircraft observations of cumulus clouds from the RACORO campaign and numerical simulations with the Explicit Mixing Parcel Model. The time scales include: τ evap , the time for droplet complete evaporation; τ phase , the time for saturation ratio deficit (S) to reach 1/e of its initial value; τ satu , the time for S to reach -0.5%; τ react , the time for complete droplet evaporation or S to reach -0.5%. It is found that the proper time scale to use depends on the specific objectives of entrainment-mixing studies. First, if the focus is on the variations of liquid water content (LWC) and S, then τ react for saturation, τ satu and τ phase are almost equivalently appropriate, because they all represent the rate of dry air reaching saturation or of LWC decrease. Second, if one focuses on the variations of droplet size and number concentration, τ react for complete evaporation and τ evap are proper because they characterize how fast droplets evaporate and whether number concentration decreases. Moreover, τ react for complete evaporation and τ evap are always positively correlated with homogeneous mixing degree (ψ), thus the two time scales, especially τ evap , are recommended for developing parameterizations. However, ψ and the other time scales can be negatively, positively, or not correlated, depending on the dominant factors of the entrained air (i.e., relative humidity or aerosols). Third and finally, all time scales are proportional to each other under certain microphysical and thermodynamic conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Observed and Simulated Sensitivities of Spring Greenup to Preseason Climate in Northern Temperate and Boreal Regions: Climate control on spring phenology

Vegetation phenology plays an important role in regulating land-atmosphere energy, water, and trace-gas exchanges. Changes in spring greenup (SG) have been documented in the past half-century in response to ongoing climate change. We use normalized difference vegetation index generated from NOAA's advanced very high resolution radiometer data in the Global Inventory Modeling and Monitoring Study project over the 1982–2005 period, coupled with climate reanalysis (Climate Research Unit-National Centers for Environmental Prediction) to investigate the SG responses to preseason climate change in northern temperate and boreal regions. We compared these observed responses to the simulated SG responses to preseason climate inferred from the Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) over 1982–2005. The observationally inferred SG suggests that there has been an advance of about 1 days per decade between 1982 and 2005 in the northern midlatitude to high latitude, with significant spatial heterogeneity. The spatial heterogeneity of the SG advance results from heterogeneity in the change of the preseason climate as well as varied vegetation responses to the preseason climate across biomes. The SG to preseason temperature sensitivity is highest in forests other than deciduous needleleaf forests, followed by temperate grasslands and woody savannas. The SG in deciduous needleleaf forests, open shrublands, and tundra is relatively insensitive to preseason temperature. Although the extent of regions where the SG is sensitive to preseason precipitation is smaller than the extent of regions where the SG is sensitive to preseason temperature, the biomes that are more sensitive to temperature are also more sensitive to precipitation, suggesting the interactive control of temperature and precipitation. Finally, in the mean, the CMIP5 ESMs reproduced the dominant latitudinal preseason climate trends and SG advances. However, large biases in individual ESMs for the preseason period, climate, and SG sensitivity imply needed model improvements to climate prediction and phenological process parameterizations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Impact of Lateral Mixing in the Ocean on El Nino in a Suite of Fully Coupled Climate Models

Abstract We examine the dependence of the amplitude of the El Nino‐Southern Oscillation (ENSO) on the mixing coefficient parameterizing the lateral mixing of tracers (A Redi ) The value of this coefficient is very uncertain, ranging in Earth System Models between a few hundred and a few thousand m 2 s −1 , with some observational estimates showing even higher values. A suite of simulations is made with two spatially varying distributions of A Redi derived from satellite observations as well as four simulations where a spatially constant A Redi is varied over a factor of 6. Surprisingly, larger values of A Redi result in stronger ENSO variability despite the higher mixing coefficients producing more efficient lateral diffusive damping of anomalies. This is because lateral mixing also warms the cold tongue, increasing vertical temperature gradients and decreasing horizontal temperature gradients. Larger vertical temperature gradients make sea surface temperatures more responsive to atmospheric forcing, while smaller horizontal temperature gradients shift the location of convection and make the atmosphere more responsive to sea surface temperature anomalies. The last effect holds across simulations as well as within individual simulations and thus also helps to explain interdecadal variability in ENSO amplitude. By contrast, a previously proposed anticorrelation between the amplitude interannual and annual variability does not hold across simulations, although it does hold within simulations. The propagation of thermocline depth anomalies is relatively insensitive to the value of A Redi . Properly specifying the lateral mixing along the equator (including distinguishing the impacts of subgridscale turbulence and tropical instability waves) appears essential to simulating ENSO.

54 ENVIRONMENTAL SCIENCES↗