  Earth observation from satellite sensory data poses challenging problems,
where machine learning is currently a key player. In recent years, Gaussian
Process (GP) regression has excelled in biophysical parameter estimation tasks
from airborne and satellite observations. GP regression is based on solid
Bayesian statistics and generally yields efficient and accurate parameter
estimates. However, GPs are typically used for inverse modeling based on
concurrent observations and in situ measurements only. Very often a forward
model encoding the well-understood physical relations between the state vector
and the radiance observations is available though and could be useful to
improve predictions and understanding. In this work, we review three GP models
that respect and learn the physics of the underlying processes in the context
of both forward and inverse modeling. After reviewing the traditional
application of GPs for parameter retrieval, we introduce a Joint GP (JGP) model
that combines in situ measurements and simulated data in a single GP model.
Then, we present a latent force model (LFM) for GP modeling that encodes
ordinary differential equations to blend data-driven modeling and physical
constraints of the system governing equations. The LFM performs multi-output
regression, adapts to the signal characteristics, is able to cope with missing
data in the time series, and provides explicit latent functions that allow
system analysis and evaluation. Finally, we present an Automatic Gaussian
Process Emulator (AGAPE) that approximates the forward physical model using
concepts from Bayesian optimization and at the same time builds an optimally
compact look-up-table for inversion. We give empirical evidence of the
performance of these models through illustrative examples of vegetation
monitoring and atmospheric modeling.
