Search NASA⌕ Search

SEARCH · Search NASA

Results for “large ensembles”

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 163 records · Page 9

The Tropical Rainfall Measuring Mission (TRMM)

Recognizing the importance of rain in the tropics and the accompanying latent heat release, NASA for the U.S. and NASDA for Japan have partnered in the design, construction and flight of an Earth Probe satellite to measure tropical rainfall and calculate the associated heating. Primary mission goals are: 1) the understanding of crucial links in climate variability by the hydrological cycle, 2) improvement in the large-scale models of weather and climate, and 3) improvement in understanding cloud ensembles and their impacts on larger scale circulations. The linkage with the tropical oceans and landmasses are also emphasized. The Tropical Rainfall Measuring Mission (TRMM) satellite was launched in November 1997 with fuel enough to obtain a four to five year data set of rainfall over the global tropics from 37 deg N to 37 deg S. This paper reports progress from launch date through the spring of 1999. The data system and its products and their access is described, as are the algorithms used to obtain the data. Some exciting early results from TRMM are described. Some important algorithm improvements are shown. These will be used in the first total data reprocessing, scheduled to be complete in early 2000. The reader is given information on how to access and use the data.

Simpson, Joanne↗

Crustal Deformation in Southcentral Alaska: The 1964 Prince William Sound Earthquake Subduction Zone

This article, for Advances in Geophysics, is a summary of crustal deformation studies in southcentral Alaska. In 1964, southcentral Alaska was struck by the largest earthquake (moment magnitude 9.2) occurring in historical times in North America and the second largest earthquake occurring in the world during the past century. Conventional and space-based geodetic measurements have revealed a complex temporal-spatial pattern of crustal movement. Numerical models suggest that ongoing convergence between the North America and Pacific Plates, viscoelastic rebound, aseismic creep along the tectonic plate interface, and variable plate coupling all play important roles in controlling both the surface and subsurface movements. The geodetic data sets include tide-gauge observations that in some cases provide records back to the decades preceding the earthquake, leveling data that span a few decades around the earthquake, VLBI data from the late 1980s, and GPS data since the mid-1990s. Geologic data provide additional estimates of vertical movements and a chronology of large seismic events. Some of the important features that are revealed by the ensemble of studies that are reviewed in this paper include: (1) Crustal uplift in the region that subsided by up 2 m at the time of the earthquake is as much as 1 m since the earthquake. In the Turnagain Arm and Kenai Peninsula regions of southcentral Alaska, uplift rates in the immediate aftermath of the earthquake reached 150 mm/yr , but this rapid uplift decayed rapidly after the first few years following the earthquake. (2) At some other locales, notably those away the middle of the coseismic rupture zone, postseismic uplift rates were initially slower but the rates decay over a longer time interval. At Kodiak Island, for example, the uplift rates have been decreasing at a rate of about 7mm/yr per decade. At yet other locations, the uplift rates have shown little time dependence so far, but are thought not to be sustainable throughout the several hundred year recurrence time for great earthquake. The nearly 10 mm/yr uplift rate at Seldovia on the Kenai Peninsula is an example.

Cohen, Steven C.↗

Evaluation of Arctic Warming in Mid-Pliocene Climate Simulations

Palaeoclimate simulations improve our understanding of the climate, inform us about the performance of climate models in a different climate scenario, and help to identify robust features of the climate system. Here, we analyse Arctic warming in an ensemble of 16 simulations of the mid-Pliocene Warm Period (mPWP), derived from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The PlioMIP2 ensemble simulates Arctic (60–90∘ N) annual mean surface air temperature (SAT) increases of 3.7 to 11.6 ∘C compared to the pre-industrial period, with a multi-model mean (MMM) increase of 7.2 ∘C. The Arctic warming amplification ratio relative to global SAT anomalies in the ensemble ranges from 1.8 to 3.1 (MMM is 2.3). Sea ice extent anomalies range from −3.0 to −10.4×106 km2, with a MMM anomaly of −5.6×106 km2, which constitutes a decrease of 53 % compared to the pre-industrial period. The majority (11 out of 16) of models simulate summer sea-ice-free conditions (≤1×106 km2) in their mPWP simulation. The ensemble tends to underestimate SAT in the Arctic when compared to available reconstructions, although the degree of underestimation varies strongly between the simulations. The simulations with the highest Arctic SAT anomalies tend to match the proxy dataset in its current form better. The ensemble shows some agreement with reconstructions of sea ice, particularly with regard to seasonal sea ice. Large uncertainties limit the confidence that can be placed in the findings and the compatibility of the different proxy datasets. We show that while reducing uncertainties in the reconstructions could decrease the SAT data–model discord substantially, further improvements are likely to be found in enhanced boundary conditions or model physics. Lastly, we compare the Arctic warming in the mPWP to projections of future Arctic warming and find that the PlioMIP2 ensemble simulates greater Arctic amplification than CMIP5 future climate simulations and an increase instead of a decrease in Atlantic Meridional Overturning Circulation (AMOC) strength compared to pre-industrial period. The results highlight the importance of slow feedbacks in equilibrium climate simulations, and that caution must be taken when using simulations of the mPWP as an analogue for future climate change.

Palaeoclimate simulations↗

Atom Interferometry for Detection of Gravitational Waves: NASA Institute for Advanced Concepts (NIAC) - Phase 1 Final Report

This report presents the results of the 2012-2013 NASA Institute for Advanced Concepts (NIAC) Phase 1 "Atom Interferometry for Detection of Gravitational Waves" project. The origin of this GW (Gravitational Wave) detection concept using atoms can be traced to theoretical work that first appeared in 2008 and also to a satellite mission-focused followup study that was done in 2011. The goal of the current project was to explore both theoretical and technical issues surrounding the implementation of this idea, as well as to begin performing proof-of-concept experiments to validate critical aspects of the proposal.The top level trade space for the detector design is driven by the strategy employed to mitigate laser frequency noise, which, if uncontrolled, can mask GW signatures. One of the advantages of the atom interferometric approach is the possibility of single baseline detection (Fig. 1.1), even in the presence of laser noise. This is enabled by the differential measurement between the two ensembles of atoms, which can result in substantial laser noise suppression. The details of this suppression depend on the atomic physics techniques used to implement the atom interferometry. Specifically, we considered the effect on noise suppression that results from using traditional two-photon Raman transitions (with alkali atoms) and also single-photon transitions (with alkaline earth-like atoms).The interferometers shown in Fig 1.1(b) take advantage of single-photon transitions (as opposed to traditional Raman transitions) because using light pulses from one direction at a time allows for near perfect common-mode cancellation of laser phase noise, even for long baselines. This calls for the use of atomic transitions with an (ideally large) optical energy level difference with a long (greater than 1 second) lifetime, such as high-transitions routinely used for optical atomic clocks in species like Sr, Ca and Yb. Notably, large momentum transfer (LMT) atom optics - and the sensitivity enhancement they confer – can still be realized by simply adding additional pairs of alternating pulses to each beam splitter process. Section 3 reports on the theoretical work we performed to justify this GW detection protocol using single-photon transitions. This approach represents a new method for GW detection using atoms that is distinct from the original proposal from 2008. At the system level, we evaluated three architectures, each of which implements a different solution to the laser frequency noise issue. The first two designs are based on two-photon Raman transitions with Rb atoms. One of these is a three-satellite, multiple baseline design while the other is a two-satellite, single baseline design. The third proposal is a two-satellite, single baseline design that uses single-photon transitions with Sr atoms. These three architectures are described in more detail in Section 2. There are a number of known technical issues that we have started to address using ground-based experiments. These issues include atom technology development needs such as, for example, lower ensemble temperature requirements and large momentum transfer (LMT) atom optics. To this end, we have built a 10-meter scale atom drop tower, where we can perform proof-of-principle demonstrations of the proposed AGIS detector in an environment that permits more than 2.5 seconds of free-fall time. This facility allows for demonstration of atom interferometry with long interrogation time (seconds) and large atom wavepacket

Saif, Babak N.↗

Intercomparison of Model Simulations of the Impact of 1997/98 El Nino on South American Summer Monsoon

The simulations of climatology and response of the South American summer monsoon (SASM) to the 1997/98 El Nino are investigated using six atmospheric general circulation models. Results show all models simulate the large-scale features of the SASM reasonably well. However, both stationary and seasonal components of the surface pressure are overestimated, resulting in an excessively strong SASM in the model climatology. The low-level northwesterly jet over eastern foothills of the Andes is not well resolved because of the coarse resolution of the models. Large rainfall simulation biases are found in association with the Andes and the Atlantic ITCZ, indicating model problems in handling steep mountains and parameterization of convective processes. The simulation of the 1997/98 El Nino impact on SASM is examined based on an ensemble of ten two-year (September 1996 - August 1998) integration. Results show that most models can simulate the large-scale tropospheric warming response over the tropical central Pacific, including the dynamic response of Rossby wave propagation of the Pacific-South America (PSA) pattern that influences remote areas. Deficiencies are found in simulating the regional impacts over South America. Model simulation fails to capture the southeastward expansion of anomalously warm tropospheric air. As a result, the upper tropospheric anomalous high over the subtropical Andes is less pronounced, and the enhancement of subtropical westerly jet is displaced 5deg-10deg equatorward compared to the observed. Over the Amazon basin, the shift of Walker cell induced by El Nino is not well represented, showing anomalous easterlies in both upper and lower troposphere.

Zhou, Jiayu↗

Assessing the Impact of Pre-gpm Microwave Precipitation Observations in the Goddard WRF Ensemble Data Assimilation System

The forthcoming Global Precipitation Measurement (GPM) Mission will provide next generation precipitation observations from a constellation of satellites. Since precipitation by nature has large variability and low predictability at cloud-resolving scales, the impact of precipitation data on the skills of mesoscale numerical weather prediction (NWP) is largely affected by the characterization of background and observation errors and the representation of nonlinear cloud/precipitation physics in an NWP data assimilation system. We present a data impact study on the assimilation of precipitation-affected microwave (MW) radiances from a pre-GPM satellite constellation using the Goddard WRF Ensemble Data Assimilation System (Goddard WRF-EDAS). A series of assimilation experiments are carried out in a Weather Research Forecast (WRF) model domain of 9 km resolution in western Europe. Sensitivities to observation error specifications, background error covariance estimated from ensemble forecasts with different ensemble sizes, and MW channel selections are examined through single-observation assimilation experiments. An empirical bias correction for precipitation-affected MW radiances is developed based on the statistics of radiance innovations in rainy areas. The data impact is assessed by full data assimilation cycling experiments for a storm event that occurred in France in September 2010. Results show that the assimilation of MW precipitation observations from a satellite constellation mimicking GPM has a positive impact on the accumulated rain forecasts verified with surface radar rain estimates. The case-study on a convective storm also reveals that the accuracy of ensemble-based background error covariance is limited by sampling errors and model errors such as precipitation displacement and unresolved convective scale instability.

microwave observations↗

A Model-Model and Data-Model Comparison for the Early Eocene Hydrological Cycle

A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Boundary conditions↗

Novae as distance indicators

Nova shells are characteristically prolate with equatorial bands and polar caps. Failure to account for the geometry can lead to large errors in expansion parallaxes for individual novae. When simple prescriptions are used for deriving expansion parallaxes from an ensemble of randomly oriented prolate spheroids, the average distance will be too small by factors of 10 to 15 percent. The absolute magnitudes of the novae will be underestimated and the resulting distance scale will be too small by the same factors. If observations of partially resolved nova shells select for large inclinations, the systematic error in the resulting distance scale could easily be 20 to 30 percent. Extinction by dust in the bulge of M31 may broaden and shift the intrinsic distribution of maximum nova magnitudes versus decay rates. We investigated this possibility by projecting Arp's and Rosino's novae onto a composite B - 6200A color map of M31's bulge. Thirty two of the 86 novae projected onto a smooth background with no underlying structure due to the presence of a dust cloud along the line of sight. The distribution of maximum magnitudes versus fade rates for these unreddened novae is indistinguishable from the distribution for the entire set of novae. It is concluded that novae suffer very little extinction from the filamentary and patchy distribution of dust seen in the bulge of M31. Time average B and H alpha nova luminosity functions are potentially powerful new ways to use novae as standard candles. Modern CCD observations and the photographic light curves of M31 novae found during the last 60 years were analyzed to show that these functions are power laws. Consequently, unless the eruption times for novae are known, the data cannot be used to obtain distances.

Ford, Holland C.↗

The orbit properties of colliding co-orbiting bodies

It is generally assumed that an ensemble of small bodies located in similar Keplarian orbits will, because of collisions, tend to disperse into more and more dissimilar orbits. For example, it is thought that the asteroids may represent the remnants of a few larger bodies that broke up or failed to fully accrete. A proposed experiment using the space station consists of an ensemble of small bodies or particles released gently from a central location in a large chamber, much like the breaking of billiard balls. The particles would them co-orbit and occasionally collide. Their subsequent behavior could be monitored by several video recorders, their linear and angular velocities before and after collision calculated, and their general behavior studied. This experiment could yield results of fundamental importance for theories of the origin of the planets, asteroids, comets, and probably ring systems.

Freeman, J. W.↗

Large-scale motions in a supersonic turbulent boundary layer on a curved surface

This paper presents measurements in a Mach 3 turbulent boundary layer which was subjected to a short region of concave surface curvature. Mean velocity and time-averaged turbulence measurements are presented to indicate the severity of the distortion, although the focus of this study is on how the large-scale organized motions are affected by the distortion. A 12 deg increase in the inclination angle of organized motions is observed, and correlation measurements indicate an increase in the streamwise length scale. VITA conditional sampling is used to identify strong positive streamwise mass-flux gradients as the large-scale structure responsible for the measured structure angle and to determine the ensemble-averaged flowfield associated with the mass-flux gradient structures.

Donovan, J. F.↗

Fast Responses on Pre-Industrial Climate from Present-Day Aerosols in A CMIP6 Multi-Model Study

In this work, we use Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations from 10 Earth system models (ESMs) and general circulation models (GCMs) to study the fast climate responses on pre-industrial climate, due to present-day aerosols. All models carried out two sets of simulations: a control experiment with all forcings set to the year 1850 and a perturbation experiment with all forcings identical to the control, except for aerosols with precursor emissions set to the year 2014. In response to the pattern of all aerosols effective radiative forcing (ERF), the fast temperature responses are characterized by cooling over the continental areas, especially in the Northern Hemisphere, with the largest cooling over East Asia and India, sulfate being the dominant aerosol surface temperature driver for present-day emissions. In the Arctic there is a warming signal for winter in the ensemble mean of fast temperature responses, but the model-to-model variability is large, and it is presumably linked to aerosol-induced circulation changes. The largest fast precipitation responses are seen in the tropical belt regions, generally characterized by a reduction over continental regions and presumably a southward shift of the tropical rain belt. This is a characteristic and robust feature among most models in this study, associated with weakening of the monsoon systems around the globe (Asia, Africa and America) in response to hemispherically asymmetric cooling from a Northern Hemisphere aerosol perturbation, forcing possibly the Intertropical Convergence Zone (ITCZ) and tropical precipitation to shift away from the cooled hemisphere despite that aerosols' effects on temperature and precipitation are only partly realized in these simulations as the sea surface temperatures are kept fixed. An interesting feature in aerosol-induced circulation changes is a characteristic dipole pattern with intensification of the Icelandic Low and an anticyclonic anomaly over southeastern Europe, inducing warm air advection towards the northern polar latitudes in winter.

Pre-industrial climate↗

Second-order closure PBL model with new third-order moments: Comparison with LES data

This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

Canuto, V. M.↗

Regional and Model-Specific Response Types in A Global Gridded Crop Model Ensemble

Crop models are often employed to project crop yields under changing conditions such as global warming and associated management change for adaptation. Multi-model ensembles are promoted to enhance the robustness of projections, but questions remain on what causes often large differences between projections of individual models. Global Gridded Crop Models (GGCMs) are especially exposed to this question when applied for assessing climate change impacts, adaptation, environmental impacts of agricultural production, because their results are used in downstream analyses, such as in integrated assessment or economic modeling for projecting future land-use change. Even though global gridded crop models are often based on detailed field-scale models or have implemented similar modeling principles in other ecosystem models, global-scale models are subject to substantial uncertainties from both model structure and parametrization as well as from calibration and input data quality. AgMIP’s Global Gridded Crop Model Intercomparison (GGCMI) has thus set out to intercompare GGCMs in order to evaluate model performance, describe model uncertainties, identify inconsistencies within the ensemble and underlying reasons, and to ultimately improve models and modeling capacities. In phase 2 of the GGCMI activities, 12 modeling groups followed a modeling protocol that asked for up to 1404 31-year global simulations at 0.5 arc-degree spatial resolution to assess models’ sensitivities to changes in carbon dioxide (C; 4 different levels) temperature (T; 7 different offset levels), water supply (W; 9 levels), and nitrogen (N; 3 levels), the so-called CTWN experiment (Franke et al. 2020; http://dx.doi.org/10.5194/gmd-13-2315-2020). We here present analyses of model response types using impact response surfaces along the C, T, W, and N dimensions, respectively and collectively. Doing so, we can understand differences in simulated responses per driver rather than aggregated changes in yields. We find that models’ sensitivities to the individual driver dimensions are substantially different and often more different across models than across regions. A cluster analysis finds regional and model-specific patterns. There is some agreement across models with respect to the spatial patterns of response types but strong differences in the distribution of response type clusters across models suggests that models need to undergo further scrutiny. We suggest establishing standards in model process evaluation not only against historical dynamics but also against dedicated experiments across the CTWN dimensions.

crop models↗

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗

Multi-Agency Ensemble Forecast of Wildfire Air Quality in the United States: Toward Community Consensus of Early Warning

Wildfires pose increasing risks to human health and properties in North America. Due to large uncertainties in fire emission, transport, and chemical transformation, it remains challenging to accurately predict air quality during wildfire events, hindering our collective capability to issue effective early warnings to protect public health and welfare. Here we present a new real-time Hazardous Air Quality Ensemble System (HAQES) by leveraging various wildfire smoke forecasts from three U.S. federal agencies (NOAA, NASA, and Navy). Compared to individual models, the HAQES ensemble forecast significantly enhances forecast accuracy. To further enhance forecasting performance, a weighted ensemble forecast approach was introduced and tested. Compared to the unweighted ensemble mean, the multilinear regression weighted ensemble reduced fractional bias by 34% in the major fire regions, false alarm rate by 72%, and increased hit rate by 17%. Finally, we improved the weighted ensemble using quantile regression and weighted regression methods to enhance the forecast of extreme air quality events. The advanced weighted ensemble increased the PM2.5 exceedance hit rate by 55% compared to the ensemble mean. Our findings provide insights into the development of advanced ensemble forecast methods for wildfire air quality, offering a practical way to enhance decision-making support to protect public health.

Yunyao Li↗

Influences of Local Sea-Surface Temperatures and Large-scale Dynamics on Monthly Precipitation Inferred from Two 10-year GCM-Simulations

Two parallel sets of 10-year long: January 1, 1982 to December 31, 1991, simulations were made with the finite volume General Circulation Model (fvGCM) in which the model integrations were forced with prescribed sea-surface temperature fields (SSTs) available as two separate SST-datasets. One dataset contained naturally varying monthly SSTs for the chosen period, and the oth& had the 12-monthly mean SSTs for the same period. Plots of evaporation, precipitation, and atmosphere-column moisture convergence, binned by l C SST intervals show that except for the tropics, the precipitation is more strongly constrained by large-scale dynamics as opposed to local SST. Binning data by SST naturally provided an ensemble average of data contributed from disparate locations with same SST; such averages could be expected to mitigate all location related influences. However, the plots revealed: i) evaporation, vertical velocity, and precipitation are very robust and remarkably similar for each of the two simulations and even for the data from 1987-ENSO-year simulation; ii) while the evaporation increased monotonically with SST up to about 27 C, the precipitation did not; iii) precipitation correlated much better with the column vertical velocity as opposed to SST suggesting that the influence of dynamical circulation including non-local SSTs is stronger than local-SSTs. The precipitation fields were doubly binned with respect to SST and boundary-layer mass and/or moisture convergence. The analysis discerned the rate of change of precipitation with local SST as a sum of partial derivative of precipitation with local SST plus partial derivative of precipitation with boundary layer moisture convergence multiplied by the rate of change of boundary-layer moisture convergence with SST (see Eqn. 3 of Section 4.5). This analysis is mathematically rigorous as well as provides a quantitative measure of the influence of local SST on the local precipitation. The results were recast to examine the dependence of local rainfall on local SSTs; it was discernible only in the tropics. Our methodology can be used for computing relationship between any forcing function and its effect(s) on a chosen field.

Sud, Y. C.↗

Mechanisms of Diurnal Precipitation over the United States Great Plains: A Cloud-Resolving Model Simulation

The mechanisms of summertime diurnal precipitation in the US Great Plains were examined with the two-dimensional (2D) Goddard Cumulus Ensemble (GCE) cloud-resolving model (CRM). The model was constrained by the observed large-scale background state and surface flux derived from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Program s Intensive Observing Period (IOP) data at the Southern Great Plains (SGP). The model, when continuously-forced by realistic surface flux and large-scale advection, simulates reasonably well the temporal evolution of the observed rainfall episodes, particularly for the strongly forced precipitation events. However, the model exhibits a deficiency for the weakly forced events driven by diurnal convection. Additional tests were run with the GCE model in order to discriminate between the mechanisms that determine daytime and nighttime convection. In these tests, the model was constrained with the same repeating diurnal variation in the large-scale advection and/or surface flux. The results indicate that it is primarily the surface heat and moisture flux that is responsible for the development of deep convection in the afternoon, whereas the large-scale upward motion and associated moisture advection play an important role in preconditioning nocturnal convection. In the nighttime, high clouds are continuously built up through their interaction and feedback with long-wave radiation, eventually initiating deep convection from the boundary layer. Without these upper-level destabilization processes, the model tends to produce only daytime convection in response to boundary layer heating. This study suggests that the correct simulation of the diurnal variation in precipitation requires that the free-atmospheric destabilization mechanisms resolved in the CRM simulation must be adequately parameterized in current general circulation models (GCMs) many of which are overly sensitive to the parameterized boundary layer heating.

Lee, M.-I.↗

Verification of the SMAP Level-4 Soil Moisture Analysis Using Rainfall Observations in Australia

Global, 3-hourly, 9-km resolution soil moisture estimates are available with a mean latency of ~2.5 days from the NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product. These estimates are based on the assimilation of SMAP radiometer brightness temperature (Tb) observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. Routine monitoring of the L4_SM system's assimilation diagnostics revealed occasionally large observation-minus-forecast Tb differences across eastern central Australia that resulted in large analysis increments (or adjustments) of the model forecast soil moisture. Because this region lacks in situ soil moisture measurements, we developed an alternative approach to assess the veracity of the soil moisture analysis increments in the L4_SM system. Using regional gauge-based precipitation data, we demonstrate that the L4_SM soil moisture increments are correlated with errors in the L4_SM precipitation forcing, suggesting that the SMAP Tb observations contribute valuable information to the L4_SM soil moisture estimates.

Reichle, R.↗