How does the mixing time of a geodesic flow on a surface vary with the genus? I have been looking at the numerical behavior of a particular quantity (of no direct importance here, though if you must know the gory details start with figure 17 here) associated to the geodesic flow on a surface of constant negative curvature and genus $g$. The behavior is quantitatively similar for $g = 2,3,4$ and physical intuition based on this quantity suggests that some "intrinsic" timescale--prime candidates are the mixing or relaxation time--should therefore depend on $g$ only weakly or even not at all. 

So: what is known about the behavior
  of the mixing and relaxation (or
  similar) times associated to these flows as $g$
  varies?

 A: Here is another thought that struck me on the way home, that I should have realized earlier.  Suppose that $S$ is a closed hyperbolic surface, of genus $g$. Then the area of $S$ is $-2\pi\chi(S) = 2\pi(2g - 2)$.  Since the area of a disk in the hyperbolic plane is exponential in its radius, it follows that the diameter of $S$ is at least logarithmic in $g$.  The mixing time of a space has to be at least the diameter, right?  So this gives a uniform lower bound on the mixing time. 
Thurston's comment is pointing out that there is no uniform upper bound.  To see this: The injectivity radius is one-half the systole (the length of the shortest closed geodesic).  For a hyperbolic surface, the collar lemma implies that as the injectivity radius goes to zero the diameter goes to infinity (this is the previously mentioned "long thin tube").  Thus the mixing time also has to grow, by the previous paragraph. 
I roughly expect the mixing time can be estimated from the logarithm of the genus and the inverse of the injectivity radius.  One reference for the geometric facts above is Peter Buser's book "Geometry and spectra of compact Riemann surfaces".  
A: OK, so I wanted to elaborate here on Sam's helpful comments. 
Set $N = 8g-4$. An explicit description of the $N$-gon $F$ that I have in mind is 
\begin{equation}
F = D \ \backslash \ \bigcup_{j=1}^{N} \left(\sqrt{a-1} \cdot D + \sqrt{a} e^{2\pi i(j-2g)/N}\right)
\end{equation}
with $a = \sec \frac{2\pi}{N}$. In particular, the nearest point to the origin is at a Euclidean distance $u := \sqrt{a} - \sqrt{a-1}$, so the hyperbolic distance is $d = \int_0^u \frac{dr}{1-r^2} = \frac{1}{2}\log\frac{1+u}{1-u}$, which evidently grows as $\log g$.
A bit more context also: I expect that $t_g f(g) \approx const$, where $t_g$ is whatever timescale and $f(g)$ is the quantity mentioned in the question. If $t_g \sim \log g$ then I'd expect that $f(g) \sim 1/\log g$, which is actually a weak enough dependence on $g$ to not be surprising based on the numerics alluded to in the question.
