Black-Litterman---in what way are expected returns hard to estimate? In the Wikipedia article about the Black-Litterman model, it states that the motivation
behind model is that "it is difficult to come up with reasonable estimates of expected returns." Why is it that expected returns are difficult to estimate, whereas the covariances are not?
Your question relates more broadly to modern portfolio theory, and can be illustrated via mean-variance analysis of a univariate time-series. The extension to the multi variate (normal) setting, is trivial. Below I use the term accuracy, in a non-formal way, relating to the variance of an estimator. In particular, when the variance of an estimator can be reduced, take this to mean the estimator can be made more accurate.

Short story: The variance of the estimator for the variance of returns can be reduced by taking more observations (finer observations - up to a limit), whereas this is not true for the variance of the estimator of the expected returns. In this sense, it is more difficult to accurately estimate expected (mean) returns compared to variance of returns. 

To see why, consider a time-series of returns, and assume the mean (per unit time), $\mu$, and variance (per unit time), $\sigma^{2}$, are constant over non-overlapping time intervals of frequency $h$. Here, $h$ represents daily, or monthly, or quarterly periods, etc. Further assume that our time-series data are observations over time intervals of length $\Delta$, where $\Delta&lt;&lt;h$. Then, define $n=h/\Delta$, as the number of observations of returns during a time interval of $h$.

For the $k^{\text{th}}$ observation interval (of length $\Delta$) during a period of length $h$, without loss of generality, we can assume that the price process for an asset is given by $$S_{k+\Delta}=S_{k}\exp{(\mu\Delta+\sigma \sqrt{\Delta} \epsilon_{k})},$$ where $\epsilon_{k}$ are iid normal. The logarithmic return (over period $\Delta$) can be written, as
$$X_{k}=\mu\Delta+\sigma \sqrt{\Delta} \epsilon_{k}.$$ For mean-variance portfolio theory, we require estimates of the true values of $\mu$ and $\sigma$, since they are of course, unknown. Call these estimators $\hat{\mu}$ and $\hat{\sigma}$. To answer your question, we only need consider the properties of these estimators.

First
$$\hat{\mu}=\sum_{k=1}^{n}X_{k}/h, \quad \mathbb{E}(\hat{\mu})=\mu, \quad \text{Var}(\hat{\mu})=\sigma^{2}/h.$$

Importantly, the variance of $\hat{\mu}$ depends upon the total length of the observation period and not on the number of observations.

Now consider, the (biased) estimator for $\sigma$:
$$\hat{\sigma}=\sum_{k=1}^{n}X_{k}^{2}/h, \quad \mathbb{E}(\hat{\sigma})=\sigma^{2}+\mu^{2}h/n, \quad \text{Var}(\hat{\sigma})=2\sigma^{4}/n+4\mu^{2}h/n^{2}.$$

The bias can be neglected (see reference for details).

Importantly, the accuracy (variance) of the estimator for $\sigma$, $\text{Var}(\hat{\sigma})$, does depend upon the number of observations, $n$ - for fixed $h$. 

The consequence of this is, for fixed $h$, by taking finer (smaller) observations intervals, $\Delta$, the accuracy of the variance estimator can be improved ($\text{Var}(\hat{\sigma})$ can be reduced). This is not true of the variance of the estimator of the expected (mean) return, $\text{Var}(\hat{\mu})$, whose accuracy can be improved only by increasing the observation period, $h$.

The reference for this is


Merton, R. C. On estimating the expected return on the market. J. financ. econ. 8, 323–361 (1980).