Harmonic level sets and boundary data This is probably a classic problem, so a good reference book or paper to get me started on this type of question would be great:
Let $\mathbb{D} \subset \mathbb{C}$ be the unit disk with boundary $\mathbb{T}$. Take a function $f : \mathbb{T} \to \mathbb{R}$ and let $u_f$ be the harmonic extension to $\mathbb{D}$ using the Poisson kernel, i.e. the Poisson integral of $f$.

What can be said about the levels sets of $u_f$ in relation to the data $f$?

For example, If the data $f$ has certain behaviour such as, say, oscillatory taking positive and negative values.

 A: (This really should be a comment, but is a bit long.)
Consider elliptic regularity and maximum principle. Interior regularity of harmonic function means that away from the boundary, your function $u_f$ is real analytic, and hence its level sets will also be analytic arcs away from where $\nabla u_f = 0$. Near the boundary the behaviour may degenerate, and the regularity will depend on the regularity of $f$. 
Where $\nabla u_f = 0$, by the maximum principle for harmonic functions, the Hessian $\nabla^2 u_f$, if non-zero, must be indefinite (or you can see that just by noting it is trace free). So the non-degenerate critical points of $u_f$ are saddle points, and have the classical structure with two level-lines intersecting there. 
In general, near a critical point, taking the Taylor expansion of the function $u_f$, you must have 
$$ u_f(x) = u_f(x_0) + \sum_{|\alpha| \geq m} a_\alpha (x-x_0)^\alpha $$
where $\alpha$ are multi-indices, $a_\alpha$ are coefficients, and $m \geq 1$ is the highest number of derivatives to which $u_f$ vanish. The maximum principle states that 
$$ \sum_{|\alpha| = m+1} a_\alpha x^\alpha $$
cannot be a signed function. So the critical point there correspond to an intersection of at least 2 and up to $m+1$ level curves. 
A: Here are some images to complement the nice comments made by Will and Willie:

I should have been more precise. I actually meant: 

What can be said about the geometry of the level regions of $u_f$ in relation to the data $f$?

That is, for $a, b \in \mathbb{R}$ with $a<b$, we define a level region as the set
$$
\{z:a \le u_f(z) \le b\}.
$$
In particular, I am interested what can be said about the level regions for "bad" boundary data. Let say I have some boundary data that causes $u_f$ to blow-up as it approaches the boundary $\mathbb{T}$. For example, here is an image where I "cycle" colors for the level regions, so fast color cycling as one approaches a point $z \in \mathbb{T}$ radially is equivalent to fast blow-up.


It seems like the level regions have very particular shapes. What can be said about the geometry of these regions?

Maybe a good place to start is to assume that we have the estimate
$$
\sup_{z \in \mathbb{D}} | u_f'(z)|(1-|z|^2) < \infty,
$$
so that $u_f$ is a Bloch function?
I am just a novice to these types of questions. References to books, papers, and theorems to get me started would be much appreciated. I plan to start reading Garnett and Marshall Harmonic measure, would that be a good place to start in regards to obtaining these types of results?
A: I know that's been a while now that the question has been asked, but as I'm looking more or less into this topic, I think I should share some of my discoveries in the literature. I'm somewhat amazed at the really small number of occurrences of this theme.


*

*I. De Carli & S. Hudson, Geometric remarks on the level curves of harmonic functions (2010), MR2586969

*L. Flatto, D. Newman & H. Shapiro, The level curves of harmonic functions (1966), MR197755

*W. Boothby, The topology of the level curves of harmonic functions with critical points (1951), MR43456

