I'm trying to prove that there exists a bijection $f:[0,1]\times [0,1]\longrightarrow [0,1]$ such that $(a,b)\leq_p(c,d)$ implies that $f((a,b))\leq f((c,d))$, where $\leq_p$ is the product order in $[0,1]\times [0,1]$ and $\leq$ is the usual order in $[0,1]$. In fact, I suspect that this function indeed exists, but I can't prove. I was trying something using Zorn's Lemma.
Obs: $(a,b)\leq_p (c,d)$ means $(a\leq c)\wedge (b\leq d)$.
Lemma: There is an order isomorphism from $\mathbb{N}^{\mathbb{N}}$ (under the lexicographic ordering) to $[0, 1)$.
Proof: Define $f : \mathbb{N}^\mathbb{N} \to [0, 1)$ by $$f(n_1, n_2, n_3, \dots) = 1 - 2^{-n_1} - 2^{-n_1-n_2} - 2^{-n_1-n_2-n_3} - \cdots = 1 - \sum_{k \geq 1} 2^{-(n_1 + \cdots + n_k)}.$$ Equivalently, $f(n_1, n_2, n_3, \dots)$ is the number in $[0, 1)$ with binary expansion $0.x_1x_2x_3\dots$, where $x_{n_1}, x_{n_1+n_2}, x_{n_1+n_2+n_3}, \dots = 0$, and all other $x_j = 1$. To show $f$ is surjective, let $x \in [0, 1)$ have binary expansion $0.x_1x_2x_3\dots$, where there are infinitely many $x_j$ with $x_j = 0$ since the expansion does not end in an infinite string of $1$'s. Letting $a_1 < a_2 < a_3 < \cdots$ be the indices of the zeroes in the expansion, so $a_1$ is the smallest number with $x_{a_1} = 0$, $a_2$ is the second smallest, and so on, it immediately follows that $x = f(a_1, a_2 - a_1, a_3 - a_2, \dots)$.
To show $f$ is order-preserving, suppose $(n_1, n_2, \dots) > (n_1', n_2', \dots)$, so there is some unique $j$ with $n_i = n_i'$ for $i < j$ and $n_j > n_j'$. Then we have \begin{align*} f(n_1, n_2, \dots) - f(n_1', n_2', \dots) &= \sum_{k \geq j} 2^{-(n_1' + \cdots + n_k')} - \sum_{k \geq j} 2^{-(n_1 + \cdots + n_k)} \\ &> 2^{-(n_1' + \cdots + n_j')} - \sum_{k \geq j} 2^{-(n_1 + \cdots + n_k)} \\ &\geq 2^{-(n_1' + \cdots + n_j')} - 2^{-(n_1 + \cdots + n_j - 1)} \\ &\geq 0 \end{align*} where the first inequality holds because we are simply removing the rest of the first sum after the first term, and the rest of the sum is positive, the second inequality holds because the $k$-th term of the second sum has $2^{-(n_1 + \cdots + n_k)} \leq 2^{-(n_1 + \cdots + n_j + (k-j))}$, and the final inequality holds because $n_1' + \cdots + n_j' \leq n_1 + \cdots + n_j - 1$. It follows that $f$ is order-preserving, hence also injective, and thus $f$ is an order isomorphism since both $\mathbb{N}^{\mathbb{N}}$ and $[0, 1)$ are linearly ordered. [end of lemma.]
Now we'll use this lemma to find an order-preserving bijection from $[0, 1) \times [0, 1)$ to $[0, 1)$. First, similar to @mathworker21's answer, note that the map $g : \mathbb{N}^{\mathbb{N}} \times \mathbb{N}^{\mathbb{N}} \to \mathbb{N}^{\mathbb{N}}$ given by $$((a_1, a_2, \dots), (b_1, b_2, \dots)) \to (a_1, b_1, a_2, b_2, \dots)$$ is an order-preserving bijection, where again each $\mathbb{N}^\mathbb{N}$ has the lexicographic ordering and $\mathbb{N}^{\mathbb{N}} \times \mathbb{N}^{\mathbb{N}}$ has the product ordering. Also, it follows from the lemma that $f^{-1} \times f^{-1}$ is an order isomorphism from $[0, 1) \times [0, 1)$ to $\mathbb{N}^{\mathbb{N}} \times \mathbb{N}^{\mathbb{N}}$ (both under the product ordering), and thus we conclude that map $h = f \circ g \circ (f^{-1} \times f^{-1})$ is an order-preserving bijection from $[0, 1) \times [0, 1)$ to $[0, 1)$.
Finally, using $h$ we can define $\tilde{h} : [0, 1] \times [0, 1] \to [0, 1]$ by $$\tilde{h}(x, y) = \begin{cases} \frac{1}{3}h(x, y) & x, y < 1 \\ \frac{1 + x}{3} & x < 1, y = 1 \\ \frac{2 + y}{3} & x = 1 \end{cases}$$ which you can check is an order-preserving bijection since $h$ is.
Note: For a more concrete description of $h$: if $x, y \in [0, 1)$ such that $x$ has zeroes in its binary expansion at indices $a_1, a_2, a_3, \dots$, and $y$ has zeroes at indices $b_1, b_2, b_3, \dots$, then $h(x, y)$ is the number in $[0, 1)$ with zeroes at the indices $a_1, a_1 + b_1, a_2 + b_1, a_2 + b_2, a_3 + b_2, a_3 + b_3, \dots$.