Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatiotemporal analysis”

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 91 records · Page 5

Spatiotemporal Methane Emissions from Global Lakes and Reservoirs

Inland aquatic systems, such as lakes and reservoirs, contribute substantially to global methane (CH4) emissions; yet are among the most uncertain components of the total CH4 budget. Lakes and reservoirs have received recent attention as they may generate high CH4 fluxes. Improved quantification of these CH4 fluxes, particularly their spatiotemporal distribution, is key to realistically incorporating them in CH4 modeling and budget studies. Here we report on a new global, gridded (0.25° lat × 0.25° lon) study of lake and reservoir CH4 emissions, accounting for new knowledge regarding lake and reservoir areal extent and distribution, and spatiotemporal emission patterns influenced by diurnal variability, temperature-dependent seasonality, satellite-derived freeze-thaw dynamics, and eco-climatic and physical CH4-centric type classification. The results of this new data set comprise daily CH4 emissions from lake and reservoirs throughout the full annual cycle and are tightly anchored to field observations, in situ measurements, and remote-sensing observations. Results show that reservoirs cover 297 × 103 km2 globally and emit 10.1 Tg CH4 yr-1 from diffusive (1.2 Tg CH4 yr-1) and ebullitive (8.9 Tg CH4 yr-1) emission pathways. On a global scale, CH4 emitting areas of lakes cover 1853 × 103 km2 and emit 37.1 Tg CH4 yr-1 from diffusive (17.1 Tg CH4 yr-1) and ebullitive (23.0 Tg CH4 yr-1) emission pathways and an additional ebullition flux of 5.2 Tg CH4 yr-1 upon ice-melt due to the accumulation of bubbles during the freeze period. This analysis of lakes and reservoir CH4 emission addresses multiple gaps and uncertainties in previous studies and represents an important contribution to studies of the global CH4 budget. The new data sets and methodologies from this study provide a framework to better understand and model the current and future role of lakes and reservoirs in the global CH4 budget and to guide efforts to mitigate inland aquatic system CH4 emissions. This presentation will describe the methodologies applied and the major results of this study focusing on the spatiotemporal distribution of lake and reservoir type classification, processes driving lake and reservoir emission seasonality, and the contribution of individual lake and reservoir types to the annual cycle of inland aquatic CH4 emissions.

Spatiotemporal↗

The Sensitivity of WRF Daily Summertime Simulations over West Africa to Alternative Parameterizations. Part 1: African Wave Circulation

The performance of the NCAR Weather Research and Forecasting Model (WRF) as a West African regional-atmospheric model is evaluated. The study tests the sensitivity of WRF-simulated vorticity maxima associated with African easterly waves to 64 combinations of alternative parameterizations in a series of simulations in September. In all, 104 simulations of 12-day duration during 11 consecutive years are examined. The 64 combinations combine WRF parameterizations of cumulus convection, radiation transfer, surface hydrology, and PBL physics. Simulated daily and mean circulation results are validated against NASA's Modern-Era Retrospective Analysis for Research and Applications (MERRA) and NCEP/Department of Energy Global Reanalysis 2. Precipitation is considered in a second part of this two-part paper. A wide range of 700-hPa vorticity validation scores demonstrates the influence of alternative parameterizations. The best WRF performers achieve correlations against reanalysis of 0.40-0.60 and realistic amplitudes of spatiotemporal variability for the 2006 focus year while a parallel-benchmark simulation by the NASA Regional Model-3 (RM3) achieves higher correlations, but less realistic spatiotemporal variability. The largest favorable impact on WRF-vorticity validation is achieved by selecting the Grell-Devenyi cumulus convection scheme, resulting in higher correlations against reanalysis than simulations using the Kain-Fritch convection. Other parameterizations have less-obvious impact, although WRF configurations incorporating one surface model and PBL scheme consistently performed poorly. A comparison of reanalysis circulation against two NASA radiosonde stations confirms that both reanalyses represent observations well enough to validate the WRF results. Validation statistics for optimized WRF configurations simulating the parallel period during 10 additional years are less favorable than for 2006.

variability↗

An Inversion Analysis of Recent Variability in Natural CO2 Fluxes Using GOSAT and In Situ Observations

About one-half of the global CO2 emissions from fossil fuel combustion and deforestation accumulates in the atmosphere, where it contributes to global warming. The rest is taken up by vegetation and the ocean. The precise contribution of the two sinks, and their location and year-to-year variability are, however, not well understood. We use two different approaches, batch Bayesian synthesis inversion and variational data assimilation, to deduce the global spatiotemporal distributions of CO2 fluxes during 2009-2010. One of our objectives is to assess different sources of uncertainties in inferred fluxes, including uncertainties in prior flux estimates and observations, and differences in inversion techniques. For prior constraints, we utilize fluxes and uncertainties from the CASA-GFED model of the terrestrial biosphere and biomass burning driven by satellite observations and interannually varying meteorology. We also use measurement-based ocean flux estimates and two sets of fixed fossil CO2 emissions. Here, our inversions incorporate column CO2 measurements from the GOSAT satellite (ACOS retrieval, filtered and bias-corrected) and in situ observations (individual flask and afternoon-average continuous observations) to estimate fluxes in 108 regions over 8-day intervals for the batch inversion and at 3 x 3.75 weekly for the variational system. Relationships between fluxes and atmospheric concentrations are derived consistently for the two inversion systems using the PCTM atmospheric transport model driven by meteorology from the MERRA reanalysis. We compare the posterior fluxes and uncertainties derived using different data sets and the two inversion approaches, and evaluate the posterior atmospheric concentrations against independent data including aircraft measurements. The optimized fluxes generally resemble those from other studies. For example, the results indicate that the terrestrial biosphere is a net CO2 sink, and a GOSAT-only inversion suggests a shift in the global sink from the tropics south to the north relative to the prior and to an in-situ-only inversion. We also find a smaller terrestrial sink in higher-latitude northern regions in boreal summer of 2010 relative to 2009.

Carbon dioxide↗

Smooth Pursuit of a Partially Occluded Object

There has long been qualitative evidence that humans can pursue an object defined only by the motion of its parts. We explored this quantitatively using an occluded diamond stimulus. Four subjects (one naive) tracked a line-figure diamond moving along an elliptical path (0.9 Hz) either clockwise (CW) or counterclockwise (CCW) behind either an X-shaped aperture (CROSS) or two vertical rectangular apertures (BARS), which obscured the corners. Although the stimulus consisted of only four line segments (108 cd/square m) moving within a visible aperture (0.2 cd/square m) behind a foreground (38 cd/square m), it is largely perceived as a coherently moving diamond. The inter-saccadic portions of eye-position traces were fit with sinusoids. All subjects tracked object motion with considerable temporal accuracy. The mean phase lag was 5 deg/6 deg (CROSS/BARS) and the mean relative phase between the horizontal and vertical components was +95 deg/+92 deg (CW) and -85 deg/-75 deg (CCW), which is close to perfect. Furthermore, a chi-square analysis showed that 56% of BARS trials were consistent with tracking the correct elliptical shape (p is less than 0.05), although segment motion was purely vertical. These data disprove the main tenet of most models of pursuit: that it is a system that seeks to minimize retinal image motion through negative feedback. Rather, the main drive must be a visual signal which has already integrated spatiotemporal retinal information into an object-motion signal.

Stone, L. S.↗

Intergenerational Inequities in Exposure to Climate Extremes: Young Generations Are Severely Threatened By Climate Change

Under continued global warming, extreme events such as heat waves will continue to rise in frequency, intensity, duration, and spatial extent over the next decades. Younger generations are therefore expected to face more such events across their lifetimes compared with older generations. This raises important issues of solidarity and fairness across generations that have fueled a surge of climate protests led by young people in recent years and that underpin issues of intergenerational equity raised in recent climate litigation. However, the standard scientific paradigm is to assess climate change in discrete time windows or at discrete levels of warming, a “period” approach that inhibits quantification of how much more extreme events a particular generation will experience over its lifetime compared with another. By developing a “cohort” perspective to quantify changes in lifetime exposure to climate extremes and compare across generations (see the first figure), we estimate that children born in 2020 will experience a two- to sevenfold increase in extreme events, particularly heat waves, compared with people born in 1960, under current climate policy pledges. Our results highlight a severe threat to the safety of young generations and call for drastic emission reductions to safeguard their future. Meteorological extremes, hazards, or climate change impacts are mostly studied as they evolve over time under varying emission scenarios and socioeconomic pathways. For example, applying a heat wave indicator (see table S1) to four bias-adjusted global climate models indicates that the land area annually affected by such heat waves will increase from ~15% around 2020 to ~22% by 2100 under a scenario compatible with limiting global warming to 1.5°C, and to ~46% under a scenario in line with current emission reduction pledges (see the first figure). Recent studies extended this approach, studying aspects of climate change as a function of global mean temperature (GMT) increments, highlighting the scenario-independence of several extreme event indicators but remaining, in essence, a comparison of time windows. By contrast, we performed a birth cohort analysis by combining a collection of multimodel extreme event projections with country-scale life expectancy information, gridded population data, and future global temperature trajectories from the Intergovernmental Panel on Climate Change (IPCC) Special Report on Global Warming of 1.5°C (see supplementary materials). By integrating the exposure of an average person in a country or region to extreme events across their lifetime, we encapsulate spatiotemporal changes in climate hazards, population density, cohort size, and life expectancy (see the first figure).

global warming↗

Fractal Characterization of Multitemporal Scaled Remote Sensing Data

Scale is an "innate" concept in geographic information systems. It is recognized as something that is intrinsic to the ingestion, storage, manipulation, analysis, modeling, and output of space and time data within a GIS purview, yet the relative meaning and ramifications of scaling spatial and temporal data from this perspective remain enigmatic. As GISs become more sophisticated as a product of more robust software and more powerful computer systems, there is an urgent need to examine the issue of scale, and its relationship to the whole body of spatiotemporal data, as imparted in GISS. Scale is fundamental to the characterization of geo-spatial data as represented in GISS, but we have relatively little insight on the effects of, or how to measure the effects of, scale in representing multiscaled data; i.e., data that are acquired in different formats (e.g., map, digital) and exist in varying spatial, temporal, and in the case of remote sensing data, radiometric, configurations. This is particularly true in the emerging era of Integrated GISs (IGIS), wherein spatial data in a variety of formats (e.g., raster, vector) are combined with multiscaled remote sensing data, capable of performing highly sophisticated space-time data analyses and modeling. Moreover, the complexities associated with the integration of multiscaled data sets in a multitude of formats are exacerbated by the confusion of what the term "scale" is from a multidisciplinary perspective; i.e., "scale" takes on significantly different meanings depending upon one's disciplinary background and spatial perspective which can lead to substantive confusion in the input, manipulation, analyses, and output of IGISs (Quattrochi, 1993). Hence, we must begin to look at the universality of scale and begin to develop the theory, methods, and techniques necessary to advance knowledge on the "Science of Scale" across a wide number of spatial disciplines that use GISs.

Quattrochi, Dale A.↗

Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology

Changing trends in ecosystem productivity can be quantified using satellite observations of Normalized Difference Vegetation Index (NDVI). However, the estimation of trends from NDVI time series differs substantially depending on analyzed satellite dataset, the corresponding spatiotemporal resolution, and the applied statistical method. Here we compare the performance of a wide range of trend estimation methods and demonstrate that performance decreases with increasing inter-annual variability in the NDVI time series. Trend slope estimates based on annual aggregated time series or based on a seasonal-trend model show better performances than methods that remove the seasonal cycle of the time series. A breakpoint detection analysis reveals that an overestimation of breakpoints in NDVI trends can result in wrong or even opposite trend estimates. Based on our results, we give practical recommendations for the application of trend methods on long-term NDVI time series. Particularly, we apply and compare different methods on NDVI time series in Alaska, where both greening and browning trends have been previously observed. Here, the multi-method uncertainty of NDVI trends is quantified through the application of the different trend estimation methods. Our results indicate that greening NDVI trends in Alaska are more spatially and temporally prevalent than browning trends. We also show that detected breakpoints in NDVI trends tend to coincide with large fires. Overall, our analyses demonstrate that seasonal trend methods need to be improved against inter-annual variability to quantify changing trends in ecosystem productivity with higher accuracy.

Inter-Annual↗

Spatiotemporal Path-Matching for Comparisons Between Ground- Based and Satellite Lidar Measurements

The spatiotemporal sampling differences between ground-based and satellite lidar data can contribute to significant errors for direct measurement comparisons. Improvement in sample correspondence is examined by the use of radiosonde wind velocity to vary the time average in ground-based lidar data to spatially match coincident satellite lidar measurements. Results are shown for the 26 February 2004 GLAS/ICESat overflight of a ground-based lidar stationed at NASA GSFC. Statistical analysis indicates that improvement in signal correlation is expected under certain conditions, even when a ground-based observation is mismatched in directional orientation to the satellite track.

Berkoff, Timothy A.↗

Issues of Long-Term Cryogenic Propellant Storage in Microgravity

Modern multi-layer insulation (MLI) allows to sharply reduce the heat leak into cryogenic propellant storage tanks through the tank surface and, as a consequence, significantly extend the storage duration. In this situation the MLI penetrations, such as support struts, feed lines, etc., become one of the most significant challenges of the tanks heat management. This problem is especially acute for liquid hydrogen (LH2) storage, since currently no efficient cryocoolers exist that operate at very low LH2 temperatures (20K). Even small heat leaks under microgravity conditions and over the period of many months give rise to a complex slowly-developing, large-scale spatiotemporal physical phenomena in a multi-phase liquid-vapor mixture. These phenomena are not well-understood nor can be easily controlled. They can be of a potentially hazardous nature for long-term on-orbital cryogenic torage, propellant loading, tank chilldown, engine restart, and other in-space cryogenic fluid management operations. To support the engineering design solutions that would mitigate these effects a detailed physics-based analysis of heat transfer, vapor bubble formation, growth, motion, coalescence and collapse is required in the presence of stirring jets of different configurations and passive cooling devices such as MLI, thermodynamic vent system, and vapor-cooled shield. To develop physics-based models and correlations reliable for microgravity conditions and long-time scales there is a need for new fundamental data to be collected from on-orbit cryogenic storage experiments. Our report discusses some of these physical phenomena and the design requirements and future studies necessary for their mitigation. Special attention is payed to the phenomena occurring near MLI penetrations.

Muratov, C. B.↗

Leveraging STARE for Co-aligned Data Locality with netCDF and Python MPI

We have leveraged STARE indexing to package partitioned data chunks from diverse datasets into netCDF files, distributed them on a cluster of 16 lightweight nodes with their placements spatiotemporally co-aligned, and demonstrated a few integrative analyses using netCDF parallel I/O and Python MPI, with single-user performance and scalability comparable to, or even better than, that of a parallel array database management system (ADBMS) such as SciDB. However, records of the node location and STARE index ranges for each data chunk, similar to the chunk maps of SciDB, must be maintained and consulted by the I/O and analysis code for coordinating the analytic operations in parallel, in order to achieve the good performance and scalability.

Kwo-Sen Kuo↗

Analyzing Machine Learning Predictions of Passive Microwave Brightness Temperature Spectral Difference Over Snow-Covered Terrain in High Mountain Asia

Snow is an important component of the terrestrial freshwater budget in high mountainAsia (HMA) and contributes to the runoff in Himalayan rivers through snowmelt. Despitethe importance of snow in HMA, considerable spatiotemporal uncertainty exists across the different estimates of snow water equivalent for this region. In order to better estimate snow water equivalent, radiative transfer models are often used in conjunction with microwave brightness temperature measurements. In this study, the efficacy of support vector machines (SVMs), a machine learning technique, to predict passive microwave brightness temperature spectral difference (1Tb) as a function of geophysical variables (snow water equivalent, snow depth, snow temperature, and snow density) is explored through a sensitivity analysis. The use of machine learning (as opposed to radiative transfer models) is a relatively new and novel approach for improving snow water equivalent estimates. The Noah-MP land surface model within the NASALand Information System framework is used to simulate the hydrologic cycle over HMA and model geophysical variables that are then used for SVM training. The SVMsserve as a nonlinear map between the geophysical space (modeled in Noah-MP) andthe observation space (1Tb as measured by the radiometer). Advanced MicrowaveScanning Radiometer-Earth Observing System measured passive microwave brightness temperatures over snow-covered locations in the HMA region are used as training data during the SVM training phase. Sensitivity of well-trained SVMs to each Noah-MP modeled state variable is assessed by computing normalized sensitivity coefficients. Sensitivity analysis results generally conform with the known first-order physics. Input states that increase volume scattering of microwave radiation, such as snow density and snow water equivalent, exhibit a plurality of positive normalized sensitivity coefficients. In general, snow temperature was the most sensitive input to the SVM predictions. The sensitivity of each state is location and time dependent. The signs of normalized sensitivity coefficients that indicate physical irrationality are ascribed to significant cross-correlation between Noah-MP simulated states and decreased SVM prediction capability at specific locations due to insufficient training data. SVM prediction pitfalls do exist that serve to highlight the limitations of this particular machine learning algorithm.

high mountain Asia↗

Leveraging Data Intensive Computing to Support Automated Event Services

A large portion of Earth Science investigations is phenomenon- or event-based, such as the studies of Rossby waves, mesoscale convective systems, and tropical cyclones. However, except for a few high-impact phenomena, e.g. tropical cyclones, comprehensive records are absent for the occurrences or events of these phenomena. Phenomenon-based studies therefore often focus on a few prominent cases while the lesser ones are overlooked. Without an automated means to gather the events, comprehensive investigation of a phenomenon is at least time-consuming if not impossible. An Earth Science event (ES event) is defined here as an episode of an Earth Science phenomenon. A cumulus cloud, a thunderstorm shower, a rogue wave, a tornado, an earthquake, a tsunami, a hurricane, or an EI Nino, is each an episode of a named ES phenomenon," and, from the small and insignificant to the large and potent, all are examples of ES events. An ES event has a finite duration and an associated geolocation as a function of time; its therefore an entity in four-dimensional . (4D) spatiotemporal space. The interests of Earth scientists typically rivet on Earth Science phenomena with potential to cause massive economic disruption or loss of life, but broader scientific curiosity also drives the study of phenomena that pose no immediate danger. We generally gain understanding of a given phenomenon by observing and studying individual events - usually beginning by identifying the occurrences of these events. Once representative events are identified or found, we must locate associated observed or simulated data prior to commencing analysis and concerted studies of the phenomenon. Knowledge concerning the phenomenon can accumulate only after analysis has started. However, except for a few high-impact phenomena. such as tropical cyclones and tornadoes, finding events and locating associated data currently may take a prohibitive amount of time and effort on the part of an individual investigator. And even for these high-impact phenomena, the availability of comprehensive records is still only a recent development. A major reason for the lack of comprehensive ,records for the majority of the ES phenomena is the perception that they do not pose immediate and/or severe threat to life and property and are thus not consistently tracked. monitored, and catalogued. Many phenomena even lack commonly accepted criteria for definitions. However. the lack of comprehensive records is also due to the increasingly prohibitive volume of observations and model data that must be examined. NASA Earth Observing System Data Information System (EOSDIS) alone archives several petabytes (PB) of satellite remote sensing data and steadily increases. All of these factors contribute to the difficulty of methodically identifying events corresponding to a given phenomenon and significantly impede systematic investigations. In the following we present a couple motivating scenarios, demonstrating the issues faced by Earth scientists studying ES phenomena.

Clune, Thomas L.↗

Vertical Structure of a Springtime Smoky and Humid Troposphere Over the Southeast Atlantic From Aircraft and Reanalysis

The springtime atmosphere over the southeast Atlantic Ocean (SEA) is subjected to a consistent layer of biomass burning (BB) smoke from widespread fires on the African continent. An elevated humidity signal is coincident with this layer, consistently proportional to the amount of smoke present. The combined humidity and BB aerosol has potentially significant radiative and dynamic impacts. Here, we use aircraft-based observations from the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) deployments in conjunction with reanalyses to characterize covariations in humidity and BB smoke across the SEA. The observed plume–vapor relationship, and its agreement with the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis version 5 (ERA5) and Copernicus Atmosphere Monitoring Service (CAMS) reanalysis, persists across all observations, although the magnitude of the relationship varies as the season progresses. Water vapor is well represented by the reanalyses, while CAMS tends to underestimate carbon monoxide especially under high BB. While CAMS aerosol optical depth (AOD) is generally overestimated relative to ORACLES AOD, the observations show a consistent relationship between carbon monoxide (CO) and aerosol extinction, demonstrating the utility of the CO tracer to understanding vertical aerosol distribution. We next use k-means clustering of the reanalyses to examine multi-year seasonal patterns and distributions. We identify canonical profile types of humidity and of CO, allowing us to characterize changes in vapor and BB atmospheric structures, and their impacts as they covary. While the humidity profiles show a range in both total water vapor concentration and in vertical structure, the CO profiles primarily vary in terms of maximum concentration, with similar vertical structures in each. The distribution of profile types varies spatiotemporally across the SEA region and through the season, ranging from largely one type in the northeast and southwest to more evenly distributed between multiple types where air masses meet in the middle of the SEA. These distributions follow patterns of transport from the humid, smoky source region (greatest influence in the northeast of the SEA) and the seasonal changes in both humidity and smoke (increasing and decreasing through the season, respectively). With this work, we establish a framework for a more complete analysis of the broader radiative and dynamical effects of humid aerosols over the SEA.

Kristina Pistone↗

An Innovative Infrastructure with a Universal Geo-Spatiotemporal Data Representation Supporting Cost-Effective Integration of Diverse Earth Science Data

The SpatioTemporal Adaptive Resolution Encoding (STARE) is a unifying scheme encoding geospatial and temporal information for organizing data on scalable computing/storage resources, minimizing expensive data transfers. STARE provides a compact representation that turns set-logic functions into integer operations, e.g. conditional sub-setting, taking into account representative spatiotemporal resolutions of the data in the datasets. STARE geo-spatiotemporally aligns data placements of diverse data on massive parallel resources to maximize performance. Automating important scientific functions (e.g. regridding) and computational functions (e.g. data placement) allows scientists to focus on domain-specific questions instead of expending their efforts and expertise on data processing. With STARE-enabled automation, SciDB (Scientific Database) plus STARE provides a database interface, reducing costly data preparation, increasing the volume and variety of interoperable data, and easing result sharing. Using SciDB plus STARE as part of an integrated analysis infrastructure dramatically eases combining diametrically different datasets.

Rilee, Michael Lee↗

Alaska Ecological Conservation: Using NASA Earth Observations to Identify Recent Changes in River Ice Phenology and Its Impacts on Caribou Migration

Each fall, caribou (Rangifer tarandus) in the Western Arctic Herd migrate hundreds of kilometers across northwestern Alaska to acquire seasonal resources and reach wintering grounds. Various rivers intersect migration paths, and caribou can only cross over open water or high-ice extents but are unable to cross rivers in stages of partial freezing. Recent temperature increases in Alaska can alter the timing and duration of ice formation periods, impeding migratory patterns. The Massachusetts NASA DEVELOP team partnered with the National Park Service in Alaska to detect river ice onset dates and formation periods using Landsat 8 Operational Land Imager and Landsat 9 Operational Land Imager-2, Sentinel-1 C-band Synthetic Aperture Radar and Sentinel-2A/B MultiSpectral Instrument imagery in Google Earth Engine. This feasibility analysis measured ice coverage using the Normalized Difference Infrared Index, Relative Difference River Ice, and Vertical-Vertical/Vertical-Horizontal backscatter values based on the spectral and surface characteristics of rivers. The team produced annual freezing timelines, time series plots, and maps at three river stretches to analyze river ice phenology changes. Radar imagery assessed ice coverage more accurately, while optical imagery better identified ice onset dates. Meanwhile indices were unable to robustly establish ice formation thresholds across the study period. These study results can help the National Park Service to better evaluate spatiotemporal migratory shifts and contextualize recent regional caribou declines.

river phenology↗

The second-moment climatology of the GATE rain rate data

The first part of this paper presents the description of the GARP (Global Atmospheric Research Program) Atlantic Tropical Experiment (GATE) 1 rain-rate data and its two-dimensional spectral and correlation characteristics, which has made it possible to accomplish the following: to show the concentration of a significant power along the frequency axis in the spatiotemporal spectra; to detect a diurnal cycle (which has a range of variation of about 3.4-5.4 mm/n) as one of the sources of bias in the rain statistics of satellite data; to study the distinction between the north-south and east-west transport of spatial rain-rate field and character of its anisotropy; to evaluate the scales of the distinction between second-moment estimates associated with ground and satellite samples; and to determine the appropriate spatial and temporal scales of simple linear stochastic models fitted to averaged rain-rate fields. The second part of this paper is devoted to an analysis of the diffusion of the rain rate by establishing a relationship between the parameters of the multivariate autoregressive model and the coefficients of a diffusion equation. This analysis led to the use of rain data to estimate the rain advection velocity as well as other coefficients of the diffusion equation of the corresponding field. The results obtained can be used for comparison with corresponding estimates of other sources of data (satellite, Tropical Oceans Global Atmosphere Coupled Ocean - Atmosphere Response Experiment (TOGA, COARE) or simulated by physical models), for generating multiple samples of any size, for solving the inverse problems of some of the hydrodynamic equations, and in some other areas of rain data analysis and modeling.

Polyak, Ilya↗

TRMM-TMI Satellite Observed Soil Moisture and Vegetation Density (1998-2005) Show Strong Connection with El Nino in Eastern Australia

Spatiotemporal patterns in soil moisture and vegetation water content across mainland Australia were investigated from 1998 through 2005, using TRMMITMI passive microwave observations. The Empirical Orthogonal Function technique was used to extract dominant spatial and temporal patterns in retrieved estimates of moisture content for the top 1-cm of soil (theta) and vegetation moisture content (via optical depth tau). The dominant temporal theta and tau patterns were strongly correlated to El Nino/Southern Oscillation (ENSO) in spring (3 = 0.90), and to a progressively lesser extent autumn, summer and winter. The Indian Ocean Dipole (IOD) index also explained part of the variation in spring 8 and z. Cluster analysis suggested that the regions most affected by ENS0 are mainly located in eastern Australia. The results suggest that the drought conditions experienced in eastern Australia since 2000 an clearly expressed in these satellite observations have a strong connection with ENSO patterns.

Liu, Yi↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗